Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB143 465,344 C→T 63.1% G214D (GGC→GAC)  rpoD ← RNA polymerase sigma factor RpoD

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB143465,3440CT63.1% 54.9 / 70.3 65G214D (GGC→GAC) rpoDRNA polymerase sigma factor RpoD
Reads supporting (aligned to +/- strand):  ref base C (15/9);  new base T (27/14);  total (42/23)
Fisher's exact test for biased strand distribution p-value = 7.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.23e-01

AGCCAGGGCCTGCAGGGCCTCGATGCTTTCCTTGTGGTTGCGACCATTTTTCTTCAGGACCTTGGAGGTTGCCTGAAGCTGATCGGATACCGCACCGAAGCGTTGGGCTGCGACGACCGGATCCGGGCCGCTCTCGGCCTCTTCCTCGTCGTCACCGCTTTCGGATTCTTCTTCGTCGTCGTCGGACTCTTCCTTCGCCGCAGCGGCCTTGGCACCTGGGATCGGCACTTCTTCGGTTGGCGCGGCAATGTTGTCGTCAG  >  ACB143/465217‑465476
                                                                                                                               |                                                                                                                                    
aGCCAGGGCCTGCAGGGCCTCGATGCTTTCCTTGTGGTTGCGTCCATTTTTCTTCAGGACCTTGGAGGTTGCCTGAAGCTGATCGGATACCGCACCGAAGCGTTGGGCTGCGACGACCGGATCCGGGTCGCTcgcg                                                                                                                              >  5:360090/1‑136 (MQ=255)
aGCCAGGGCCTGCAGGGCCTCGATGCTTTCCTTGTGGTTGCGACCATTTTTCTTCAGGACCTTGGAGGTTGCCTGATGCTGATCGGATACCGCACCGCAGCGTTGGTCTGCGACGACCGGATCCGGGTCGcggaag                                                                                                                              >  5:250453/1‑131 (MQ=255)
aGCCAGGGCCTGCAGGGCCTCGATGCTTTCCTTGTGGTTGCGACCATTTTTCTTCAGGACCTTGGAGGTTGCCTGAAGCTGATCGGATACCGCACCGAAGCGTTGGGCTGCGACGACCGGATCCGGGTCGCTCTCg                                                                                                                              >  3:362973/1‑136 (MQ=255)
aGCCAGGGCCTGCAGGGCCTCGATGCTTTCCTTGTGGTTGCGACCATTTTTCTTCAGGACCTTGGAGGTTGCCTGAAGCTGATCGGATACCGCACCGAAGCGTTGGGCTGCGACGACCGGATCCGGGTCGCTCTCg                                                                                                                              >  4:157817/1‑136 (MQ=255)
aGCCAGGGCCTGCAGGGCCTCGATGCTTTCCTTGTGGTTGCGACCATTTTTCTTCAGGACCTTGGAGGTTGCCTGAAGCTGATCGGATACCGCACCGAAGCGTTGGGCTGCGACGACCGGATCCGGGTCGCTCTCg                                                                                                                              >  3:143849/1‑136 (MQ=255)
aGCCAGGGCCTGCAGGGCCTCGATGCTTTCCTTGTGGTTGCGACCATTTTTCTTCAGGACCTTGGAGGTTGCCTGAAGCTGATCGGATACCGCACCGAAGCGTTGGGCTGCGACGACCGGATCCGGGTCGATCTCg                                                                                                                              >  4:88740/1‑136 (MQ=255)
aGCCAGGGCCTGCAGGGCCTCGATGCTTTCCTTGTGGTTGCGACCATTTTTCTTCAGGACCTTGGAGGTTGCCTGAAGCTGATCGGATACCGCACCGAAGCGTTGGGCTGCGACGACCGGATCCGGGCCGCTCTCg                                                                                                                              >  3:211624/1‑136 (MQ=255)
aGCCAGGGCCTGCAGGGCCTCGATGCTTTCCTTGTGGTTGCGACCATTTTTCTTCAGGACCTTGGAGGTTGCCTGAAGCTGATCGGATACCGCACCGAAGCGTTGGGCTGCGACGACCGGATCCGGGCCGCTCTCg                                                                                                                              >  4:160523/1‑136 (MQ=255)
 gCCAGGGCCTGCAGGGCCTCGATGCTTTCCTTGTGGTTGCGACCATTTTTCTTCAGGACCTTGGAGGTTGCCTGAAGCTGATCGGATACCGCACCGAAGCGTTGGGCTGCGACGACCGGATCCGGGTCGCTCTCgg                                                                                                                             <  8:220087/136‑1 (MQ=255)
     gggCCTGCAGGGCCTCGATGCTTTCCTTGTGGTTGCGACCATTTTTCTTCAGGACCTTGGAGGTTGCCTGAAGCTGATCGGATACCGCACCGAAGCGTTGGGCTGCGACGACCGGATCCGGGTCGCTCTCGGCctc                                                                                                                         >  3:279515/1‑136 (MQ=255)
     gggCCTGCAGGGCCTCGATGCTTTCCTTGTGGTTGCGACCATTTTTCTTCAGGACCTTGGAGGTTGCCTGAAGCTGATCGGATACCGCACCGAAGCGTTGGGCTGCGACGACCGGATCCGGGTCGCTCTCGGCctc                                                                                                                         >  4:231704/1‑136 (MQ=255)
          tGCAGGGCCTCGATGCTTTCCTTGTGGTTGCGACCATTTTTCTTCAGGACCTTGGAGGTTGCCTGAAGCTGATCGGATACCGCACCGAAGCGTTGGGCTGCGACGACCGGATCCGGGTCGCTCTGGGCCTGGTCCt                                                                                                                    >  1:198420/1‑136 (MQ=255)
          tGCAGGGCCTCGATGCTTTCCTTGTGGTTGCGACCATTTTTCTTCAGGACCTTGGAGGTTGCCTGAAGCTGATCGGATACCGCACCGAAGCGTTGGGCTGCGACGACCGGATCCGGGTCGCTCTCGGCCTCTTccg                                                                                                                    >  3:403092/1‑135 (MQ=255)
          tGCAGGGCCTCGATGCTTTCCTTGTGGTTGCGACCATTTTTCTTCAGGACCTTGGAGGTTGCCTGAAGCTGATCGGATACCGCACCGAAGCGTTGGGCTGCGACGACCGGATCCGGGTCGCTCTCGGCCTCTTCCt                                                                                                                    >  7:359328/1‑136 (MQ=255)
          tGCAGGGCCTCGATGCTTTCCTTGTGGTTGCGACCATTTTTCTTCAGGACCTTGGAGGTTGCCTGAAGCTGATCGGATACCGCACCGAAGCGTTGGGCTGCGACGACCGGATCCGGGCCGCTCTCGGCCTCTTcc                                                                                                                     >  4:10647/1‑135 (MQ=255)
          tGCAGGGCCTCGATGCTTTCCTTGTGGTTGCGACCATTTTTCTTCAGGACCTTGGAGGTTGCCTGAAGCTGATCGGATACCGCACCGAAGCGTTGGGCTGCGACGACCGGATCCGGGCCGCTCTCGGCCTCTTCCt                                                                                                                    >  3:303508/1‑136 (MQ=255)
          tGCAGGGCCTCGATGCTTTCCTTGTGGTTGCGACCATTTTTCTTCAGGACCTTGGAGGTTGCCTGAAGCTGATCGGATACCGCACCGAAGCGTTGGGCTGCGACGACCGGATCCGGGCCGCTCTCGGCCTCTTCCt                                                                                                                    >  2:216566/1‑136 (MQ=255)
          tGCAGGGCCTCGATGCTTTCCTTGTGGTTGCGACCATTTTTCTTCAGGACCTTGGAGGTTGCCTGAAGCTGATCGGATACCGCACCGAAGCGTTGGGCTGCGACGACCGGATCCGGGCCGCTCTCGGCCTCTTCCt                                                                                                                    >  7:266561/1‑136 (MQ=255)
              gggCCTCGATGCTTTCCTTGTGGTTGCGACCATTTTTCTTCAGGACCTTGGAGGTTGCCTGAAGCTGATCGGATACCGCACCGAAGCGTTGGGCTGCGACGACCGGATCCGGGTCGCTCTCGGCCTCTTCCtcgtc                                                                                                                >  2:14009/1‑136 (MQ=255)
               ggCCTCGATGCTTTCCTTGTGGTTGCGACCATTTTTCTTCAGGACCTTGGAGGTTGCCTGAAGCTGATCGGATACCGCACCGAAGCGTTGGGCTGCGACGACCGGATCCGGGTCGCTCTCGGCCTCTTCCtcgtcg                                                                                                               >  5:26939/1‑136 (MQ=255)
               ggCCTCGATGCTTTCCTTGTGGTTGCGACCATTTTTCTTCAGGACCTTGGAGGTTGCCTGAAGCTGATCGGATACCGCACCGAAGCGTTGGGCTGCGACGACCGGATCCGGGTCGCTCTCGGCCTCTTCCtcgtcg                                                                                                               >  5:375668/1‑136 (MQ=255)
                gCCTCGATGCTTTCCTTGTGGTTGCGACCATTTTTCTTCAGGACCTTGGAGGTTGCCTGAAGCTGATCGGATACCGCACCGAAGCGTTGGGCTGCGACGACCGGATCCGGGTCGCTCTCGGCCTCTTCCTcgtcgt                                                                                                              <  2:404428/136‑1 (MQ=255)
                   tCGATGCTTTCCTTGTGGTTGCGACCATTTTTCTTCAGGACCTTGGAGGTTGCCTGAAGCTGATCGGATACCGCACCGAAGCGTTGGGCTGCGACGACCGGATCCGGGTCGCTCTCGGCCACTTCCTCGTCCTCAc                                                                                                           >  4:97968/1‑136 (MQ=255)
                    cGATGCTTTCCTTGTGGTTGCGACCATTTTTCTTCAGGACCTTGGAGGTTGCCTGAAGCTGATCGGATACCGCACCGAAGCGTTGGGCTGCGACGACCGGATCCGGGTCGCTCTCGGCCTCTTCCTCGTCGTCAcc                                                                                                          >  5:372743/1‑136 (MQ=255)
                    cGATGCTTTCCTTGTGGTTGCGACCATTTTTCTTCAGGACCTTGGAGGTTGCCTGAAGCTGATCGGATACCGCACCGAAGCGTTGGGCTGCGACGACCGGATCCGGGCCGCTCTCGGCCTCTTCCTCGTCGTCAcc                                                                                                          >  7:371631/1‑136 (MQ=255)
                    cGATGCTTTCCTTGTGGTTGCGACCATTTTTCTTCAGGACCTTGGAGGTTGCCTGAAGCTGATCGGATACCGCACCGAAGCGTTGGGCTGCGACGACCGGATCCGGGCCGCTCTCGGCCTCTTCCTCGTCGTCAcc                                                                                                          >  7:272975/1‑136 (MQ=255)
                               ttGTGGTTGCGACCATTTTTCTTCAGGACCTTGGAGGTTGCCTGAAGCTGATCGGATACCGCACCGAAGCGTTGGGCTGCGACGACCGGATCCGGGTCGCTCTCGGCCTCTTCCTCGTCGTCACCGCTTTCGGAt                                                                                                >  6:390653/1‑135 (MQ=255)
                                    gTTGCGACCATTTTTCTTCAGGACCTTGGAGGTTGCCTGAAGCTGATCGGATACCGCACCGAAGCGTTGGGCTGCGACGACCGGATCCGGGCCGCTCTCGGCCTCTTCCTCGTCGTCACCGCTTTCGGAttcttct                                                                                          >  3:342757/1‑136 (MQ=255)
                                     ttGCGACCATTTTTCTTCAGGACCTTGGAGGTTGCCTGAAGCTGATCGGATACCGCACCGAAGCGTTGGGCTGCGACGACCGGATCCGGGCCGCTCTCGGCCTCTTCCTCGTCGTCACCGCTTTCGGAttcttctt                                                                                         <  1:271194/136‑1 (MQ=255)
                                       gCGACCATTTTTCTTCAGGACCTTGGAGGTTGCCTGAAGCTGATCGGATACCGCACCGAAGCGTTGGGCTGCGACGACCGGATCCGGGCCGCTCTCGGCCTCTTCCTCGTCGTCACCGCTTTCGGAttcttcttc                                                                                        <  8:138984/135‑1 (MQ=255)
                                             aTTTTTCTTCAGGACCTTGGAGGTTGCCTGAAGCTGATCGGATACCGCACCGAAGCGTTGGGCTGCGACGACCGGATCCGGGTCGCTCTCGGCCTCTTCCTCGTCGTCACCGCTTTCGGATTCTTCTtcgtcgtcg                                                                                 <  5:375544/136‑1 (MQ=255)
                                                  tcttcAGGACCTTGGAGGTTGCCTGAAGCTGATCGGATACCGCACCGAAGCGTTGGGCTGCGACGACCGGATCCGGGCCGCTCTCGGCCTCTTCCTCGTCGTCACCGCTTTCGGATTCTTCTTcgtcggcgtcgg                                                                             <  3:83702/135‑1 (MQ=255)
                                                      cAGGACCTTGGAGGTTGCCTGAAGCTGATCGGATACCGCACCGAAGCGTTGGGCTGCGACGACCGGATCCGGGCCGCTCTCGGCCTCTTCCTCGTCGTCACCGCTTTCGGATTCTTCTTCGTCGTCGTCGGAc                                                                           <  3:10647/133‑1 (MQ=255)
                                                          aCCTTGGAGGTTGCCTGAAGCTGATCGGATACCGCACCGAAGCGTTGGGCTGCGACGACCGGATCCGGGTCGCTCTCGGCCTCTTCCTCGTCGTCACCGCTTTCGGATTCTTCTTCGTCGTCGTCGGACTcttcct                                                                    <  1:347663/136‑1 (MQ=255)
                                                            cTTGGAGGTTGCCTGAAGCTGATCGGATACCGCACCGAAGCGTTGGGCTGCGACGACCGGATCCGGGTCGCTCTCGGCCTCTTCCTCGTCGTCACCGCTTTCGGATTCTTCTTCGTCGTCGTCGGACTcttccttc                                                                  >  1:21900/1‑136 (MQ=255)
                                                            cTTGGAGGTTGCCTGAAGCTGATCGGATACCGCACCGAAGCGTTGGGCTGCGACGACCGGATCCGGGTCGCTCTCGGCCTCAACCTCGTAGTTCCCGCTTTAGGATTATTCCTAGTCGTGGTGCGACTgttacttc                                                                  >  5:335279/1‑136 (MQ=255)
                                                                 aGGTTGCCTGAAGCTGATCGGATACCGCACCGAAGCGTTGGGCTGCGACGACCGGATCCGGGTCGCTCTCGGCCTCTTCCTCGTCGTCACCGCGTTAGGATTCTTCTTCGTCGTCGTAGGACTCTTACTTcgccgc                                                             >  7:135039/1‑136 (MQ=255)
                                                                         tGAAGCTGATCGGATACCGCACCGAAGCGTTGGGCTGCGACGACCGGATCCGGGTCGCTCTCGGCCGCTTCCGCGTCGTCACCGCTTTCGGATTGTTATTAGTCCTAGTAGGGCGCTTGCTTCGCCGCACCgggcg                                                     >  3:305561/1‑133 (MQ=255)
                                                                         tGAAGCTGATCGGATACCGCACCGAAGCGTTGGGCTGCGACGACCGGATCCGGGCCGCTCTCGGCCTCTTCCTCGTCGTCACCGCTTTCGGATTCTTCTTCGTCGTCGTCGGAGTCTTCCTTCGCCGCAGCGGCCt                                                     >  3:102841/1‑136 (MQ=255)
                                                                         tGAAGCTGATCGGATACCGCACCGAAGCGTTGGGCTGCGACGACCGGATCCGGGCCGCTCTCGGCCTCTTCCTCGTCGTCACCGCTTTCGGATTCTTCTTCGTCGTCGTCGGACTCTTCCTTCGCCGCAGCGGCCt                                                     >  4:404076/1‑136 (MQ=255)
                                                                          gAAGCTGATCGGATACCGCACCGAAGCGTTGGGCTGCGACGACCGGATCCGGGCCGCTCTCGGCCTCTTCCTCGTCGTCACCGCTTTCGGATTCTTCTTCGTCGTCGTCGGACTCTTCCTTCGCCGCACCGGCCtt                                                    >  5:44292/1‑136 (MQ=255)
                                                                                 aTCGGATACCGCACCGAAGCGTTGGGCTGCGACGACCGGATCCGGGTCGCTCTCGGCCTCTTCCTCGTCGTCACCGCTTTCGGATTCTTCTTCGTCGTCGTCGGACTCTTCCTTCGCCGCAGCGGCCTTGGCACCt                                             <  7:359460/136‑1 (MQ=255)
                                                                                    ggATACCGCACCGAAGCGTTGGGCTGCGACGCCCGGATCCGGGCCGCTCTCGGCCTCTTCCTCGTCGTCACCGCTTTCGGATTCTTCTTCGTCGTCGTCGGACTCTTCCTTCGCCGCAGCGGCCTTGGCACCTggg                                          >  1:236880/1‑136 (MQ=255)
                                                                                    ggATACCGCACCGAAGCGTTGGGCTGCGACGACCGGATCCGGGTCGCTCTCGGCCTCTTGCTCGTCGTCACCGCTTTCGGATTCTTCTTCGTCGTCGTAGGCCTCATACTTAGCCGCAGCGGCCTTGGCACCTggg                                          >  2:2627/1‑136 (MQ=255)
                                                                                    ggATACCGCACCGAAGCGTTGGGCTGCGACGACCGGATCCGGGTCGCTCTCGGCCTCTTCCTCGTCGTCCCCGCTTTCGGATTCTTCTTCGTCGTCGTCGGACTCTTCCTTCGCCGCAGCGGCCTTGGCACCTggg                                          >  7:293114/1‑136 (MQ=255)
                                                                                    ggATACCGCACCGAAGCGTTGGGCTGCGACGACCGGATCCGGGTCGCTCTCGGCCTCTTCCTCGTCGTCACCGCTTTCGGATTCTTCTTCGTCGTCGTCGGACTCTTCCTTCGCCGCCGCGGCCTTGGCACCTgg                                           >  1:32005/1‑135 (MQ=255)
                                                                                      aTACCGCACCGAAGCGTTGGGCTGCGACGACCGGATCCGGGCCGCTCTCGGCCTCTTCCTCGTCGTCACCGCTTTCGGATTCTTCTTCGTCGTCGTCGGACTCTTCCTTCGCCGCAGCGGCCTTGGCACCTGGGAt                                        <  1:39897/136‑1 (MQ=255)
                                                                                               cGAAGCGTTGGGCTGCGACGACCGGATCCGGGTCGCTCTAGGCCTCTTCATAGTCGTCACCGCTTTAGGATTCTTCTTCGTCGTCGTCGGACTAGTCCTTCGCCGCAGCGGCCGTGGCACCGGGGATCGGCAcgtc                               >  8:180302/1‑136 (MQ=255)
                                                                                                   gCGTTGGGCTGCGACGACCGGATCCGGGTCGCTCTCGGCCTCTTCCTCGGCGTCCCCGCGTTCGGATTCTTCTTCGTCGTCGTCGGCCTCTTCCTTCGCCTCAGCCGCCGTGGGAACAGGGATGGGCCCGTCTTCg                           >  1:131500/1‑136 (MQ=255)
                                                                                                   gCGTTGGGCTGCGACGACCGGATCCGGGCCGCTCTCGGCCTCTTCCTCGTCGTCACCGCTTTCGGATTCTTCATCGTCGTCGTCGGACTCTTGCTTCGCCGCAGCGGCCTTGGCCCCTGGGGACGGGACTTCTTCg                           >  1:92064/1‑136 (MQ=255)
                                                                                                   gCGTTGGGCTGCGACGACCGGATCCGGGCCGCTCTCGGCCTCTTCCTCGGCGTCACCGCTTTCGGATTCTTGTTCGTCGTCGTGGGACGCTTCCTTCGCCGCAGCGGCCTTGGGACCTGGGATCGGCACGTCTTCg                           >  2:210511/1‑136 (MQ=255)
                                                                                                     gTTGGGCTGCGACGACCGGATCCGGGTCGCTCTCGGCCTCTTCCTCGTCGTCACCGCTTTCGGATTCTTCTTCGTCGTCGTCGGACTCTTCCTTCGCCGCAGCGGCCTTGGCACCTTGGATCGGCACTTCTTCGGt                         <  8:153984/136‑1 (MQ=255)
                                                                                                         ggCTGCGACGACCGGATCCGGGTCGCTGTCGGCCTCTTCCTCGTCGTCACCGCTTTCGGATTCTTCTTCGTCGTCGTCGGACGCTTCCTTCGCCGCAGCGGCCTTGGCACCTGGGATCGGCACTTCTTCGGTTGgc                     >  2:355428/1‑136 (MQ=255)
                                                                                                             gCGACGACCGGATCCGGGTCGCTCTCGGCCTCTTCCTCGTCGTCACCGCTTTCGGATTCTTCTTCGTCGTCGTCGGACTCTTCCTTCGCCGCAGCGGCCTTGGCACCTGGGATCGGCACTTCTTCGGTTGGCGCgg                 <  3:357550/136‑1 (MQ=255)
                                                                                                              cgacgaCCGGATCCGGGCCGCTCTCGGCCTCTTCCTCGTCGTCACCGCTTTCGGATTCTTCTTCGTCGTCGTCGGACTCTTCCTTCGCCGCAGCGGCCTTGGCACCTGGGATCGGCACTTCTTCGGTTGGcgcg                  <  8:366595/134‑1 (MQ=255)
                                                                                                              cgacgaCCGGATCCGGGCCGCTCTCGGCCTCTTCCTCGTCGTCACCGCTTTCGGATTCTTCTTCGTCGTCGTCGGACTCTTCCTTCGCCGCAGCGGCCTTGGCACCTGGGATCGGCACTTCTTCGGTTGGcgcg                  <  6:304264/134‑1 (MQ=255)
                                                                                                                 cgacCGGATCCGGGTCGCTCTCGGCCTCTTCCTCGTCGTCACCGCTTTCGGATTCTTCTTCGTCGTCGTCGGACTCTTCCTTCGCCGCAGCGGCCTTGGCACCTGGGATCGGCACTTCTTCGGTTGGCGCGGCAAt             <  4:279515/136‑1 (MQ=255)
                                                                                                                 cgacCGGATCCGGGTCGCTCTCGGCCTCTTCCTCGTCGTCACCGCTTTCGGATTCTTCTTCGTCGTCGTCGGACTCTTCCTTCGCCGCAGCGGCCTTGGCACCTGGGATCGGCACTTCTTCGGTTGGCGCGGCAAt             <  4:143849/136‑1 (MQ=255)
                                                                                                                 cgacCGGATCCGGGTCGCTCTCGGCCTCTTCCTCGTCGTCACCGCTTTCGGATTCTTCTTCGTCGTCGTCGGACTCTTCCTTCGCCGCAGCGGCCTTGGCACCTGGGATCGGCACTTCTTCGGTTGGCGCGGCAAt             <  1:101205/136‑1 (MQ=255)
                                                                                                                  gacCGGATCCGGGTCGCTCTCGGCCTCTTCCTCGTCGTCACCGCTTTCGGATTCTTCTTCGTCGTCGTCGGACTCTTCCTTCGCCGCAGCGGCCTTGGCACCTGGGATCGGCACTTCTTCGGTTGGCGCGGCAAtg            <  6:26939/136‑1 (MQ=255)
                                                                                                                      ggATCCGGGTCGCTCTCGGCCTCTTCCTCGTCGTCACCGCTTTCGGATTCTTCTTCGTCGTCGTCGGACTCTTCCTTCGCCGCAGCGGCCTTGGCACCTGGGATCGGCACTTCTTCGGTTGGCGCGGCAAtgttgt        <  3:81165/136‑1 (MQ=255)
                                                                                                                       gATGCGGGTCGGTGTCGGCCTCTTCCTCGTCGTCACCGCTTAGGGTTTCTTCTTCGTCGTCGTCGGACTCTTCCTTCGCCGCAGCGGCCTTGGCACCTGGGATCGGCACTTCTTCGGTTGGCGCGGCATTGTTgtc       <  7:354809/136‑1 (MQ=255)
                                                                                                                        aTCCGGGCCGCTCTCGGCCTCTTCCTCGTCGTCACCGCTTTCGGATTCTTCTTCGTCGTCGTCGGACTCTTCCTTCGCCGCAGCGGCCTTGGCACCTGGGATCGGCACTTCTTCGGTTGGCGCGGCAATGTTgtcg      <  4:102688/136‑1 (MQ=255)
                                                                                                                            gggTCGCTCTCGGCCTCTTCCTCGTCGTCACCGCTTTCGGATTCTTCTTCGTCGTCGTCGGACTCTTCCTTCGCCGCAGCGGCCTTGGCACCTGGGATCGGCACTTCTTCGGTTGGCGCGGCAATGTTGTCGTCAg  <  1:47930/136‑1 (MQ=255)
                                                                                                                            gggCCGCTCTCGGCCTCTTCCTCGTCGTCACCGCTTTCGGATTCTTCTTCGTCGTCGTCGGACTCTTCCTTCGCCGCAGCGGCCTTGGCACCTGGGATCGGCACTTCTTCGGTTGGCGCGGCAATGTTGTCGTCAg  <  4:59620/136‑1 (MQ=255)
                                                                                                                               |                                                                                                                                    
AGCCAGGGCCTGCAGGGCCTCGATGCTTTCCTTGTGGTTGCGACCATTTTTCTTCAGGACCTTGGAGGTTGCCTGAAGCTGATCGGATACCGCACCGAAGCGTTGGGCTGCGACGACCGGATCCGGGCCGCTCTCGGCCTCTTCCTCGTCGTCACCGCTTTCGGATTCTTCTTCGTCGTCGTCGGACTCTTCCTTCGCCGCAGCGGCCTTGGCACCTGGGATCGGCACTTCTTCGGTTGGCGCGGCAATGTTGTCGTCAG  >  ACB143/465217‑465476

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: