Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,008,400 +GAAA coding (483/516 nt) b4555 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,008,3981.A98.3% 425.8 / ‑0.6 118coding (481/516 nt)b4555hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (0/2);  new base A (50/66);  total (50/68)
Fisher's exact test for biased strand distribution p-value = 5.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.10e-01
*NC_0009134,008,3982.A98.3% 427.6 / ‑0.6 118coding (481/516 nt)b4555hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (0/2);  new base A (50/66);  total (50/68)
Fisher's exact test for biased strand distribution p-value = 5.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.36e-01
*NC_0009134,008,3983.G98.3% 371.7 / ‑0.9 118coding (481/516 nt)b4555hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (0/2);  new base G (50/66);  total (50/68)
Fisher's exact test for biased strand distribution p-value = 5.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.90e-01
*NC_0009134,008,3984.A98.3% 423.3 / ‑0.6 118coding (481/516 nt)b4555hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (0/2);  new base A (50/66);  total (50/68)
Fisher's exact test for biased strand distribution p-value = 5.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.52e-01

Only 100 of 132 total aligned reads displayed.

GCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAGGCCATGCGCTGGAAGGCTGGCG‑‑‑‑AA‑‑‑‑CTGGCTACAGCTGAAGGATTGATTGCGTT‑‑GTGACGGGATGAACGACGCCGAAATGCTG  >  NC_000913/4008309‑4008458
                                                                                            ||||                                                                      
gCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCG‑‑‑‑a                                                                       <  1:712134/93‑1 (MQ=255)
gCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAAGCCATGCGCTGGAAGCGGTGGCG‑‑‑‑a                                                                       >  1:7950130/1‑93 (MQ=255)
 cGATCGCGTCAATGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGGGTCAAAAGGCCATGCGCTGGAAGCGGTGGCG‑‑‑‑aa                                                                      >  2:8947093/1‑93 (MQ=255)
 cGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGTGCTGGAAGCGGTGGCG‑‑‑‑aa                                                                      <  2:6364701/93‑1 (MQ=255)
 cGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCG‑‑‑‑aa                                                                      <  1:2607125/93‑1 (MQ=255)
 cGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCG‑‑‑‑aa                                                                      <  2:8640424/93‑1 (MQ=255)
  gATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGCTTCAAAAGGCCATGCGCTGGAAGCGGTGGCG‑‑‑‑aag                                                                     >  1:6295637/1‑92 (MQ=255)
   aTCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCCTGCGCTGGAAGCGGTGGCG‑‑‑‑aaga                                                                    >  1:9955995/1‑91 (MQ=255)
    tCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGGAGCGGTGGCG‑‑‑‑aagaa                                                                   >  1:6062766/1‑90 (MQ=255)
     cgcgTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCG‑‑‑‑aagaaa                                                                  <  2:8301252/93‑5 (MQ=255)
     cgcgTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCG‑‑‑‑aagaaa                                                                  <  2:2434993/93‑5 (MQ=255)
                                     tCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCGtt                                 <  2:6362410/93‑1 (MQ=255)
                                     tCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCGtt                                 <  1:7177307/93‑1 (MQ=255)
                                      cTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCGttt                                >  2:3112731/1‑92 (MQ=255)
                                      cTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCGttt                                >  2:5960747/1‑92 (MQ=255)
                                      cTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCGttt                                <  1:1148820/93‑2 (MQ=255)
                                      cTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCGttt                                >  2:2016160/1‑92 (MQ=255)
                                       tGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCGtttg                               >  1:9334969/1‑91 (MQ=255)
                                       tGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCGtttg                               >  1:6088940/1‑91 (MQ=255)
                                       tGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCGtttg                               <  2:3652544/93‑3 (MQ=255)
                                        ggAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCGtttgg                              <  2:5823044/93‑4 (MQ=255)
                                        ggAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCGtttgg                              <  2:2132598/93‑4 (MQ=255)
                                        ggAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCGtttgg                              >  1:1532792/1‑90 (MQ=255)
                                         gAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCGtttggt                             <  1:6573968/93‑5 (MQ=255)
                                         gAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCGtttggt                             <  1:6889041/93‑5 (MQ=255)
                                         gAAGTGATGGCGGGCGGCGTGTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCGtttggt                             <  2:861771/93‑5 (MQ=255)
                                               aTGGCGGGCGGCGTTTCGAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACggg                       <  1:8977973/93‑1 (MQ=255)
                                               aTGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACggg                       >  1:750883/1‑93 (MQ=255)
                                               aTGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACggg                       >  2:3448880/1‑93 (MQ=255)
                                                tGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGa                      <  2:2881758/93‑1 (MQ=255)
                                                tGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGa                      <  2:8409837/93‑1 (MQ=255)
                                                 ggcgggcgGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGAt                     <  2:1818924/93‑1 (MQ=255)
                                                 ggcgggcgGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGAt                     <  1:3963700/93‑1 (MQ=255)
                                                 ggcgggcgGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGAt                     >  1:3153169/1‑93 (MQ=255)
                                                 ggcgggcgGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGAt                     >  2:8863266/1‑93 (MQ=255)
                                                 ggcgggcgGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGAt                     <  2:4646743/93‑1 (MQ=255)
                                                  gcgggcggCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATg                    <  2:3031378/93‑1 (MQ=255)
                                                  gcgggcggCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATg                    <  2:8626408/93‑1 (MQ=255)
                                                    gggcggCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGaa                  <  1:9500316/93‑1 (MQ=255)
                                                    gggcggCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGaa                  <  1:3076046/93‑1 (MQ=255)
                                                    gggcggCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGaa                  <  2:7651843/93‑1 (MQ=255)
                                                    gggcggCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGaa                  >  1:1242716/1‑93 (MQ=255)
                                                     ggcggcGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAac                 <  1:10002085/93‑1 (MQ=255)
                                                     ggcggcGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAac                 <  2:4349308/93‑1 (MQ=255)
                                                     ggcggcGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAac                 >  1:8626919/1‑93 (MQ=255)
                                                     ggcggcGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAac                 >  1:8626918/1‑93 (MQ=255)
                                                     ggcggcGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAac                 <  2:7374816/93‑1 (MQ=255)
                                                     ggcggcGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAac                 >  1:5145751/1‑93 (MQ=255)
                                                      gcggcgTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAacg                >  1:2907039/1‑93 (MQ=255)
                                                      gcggcgTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAacg                <  2:5280233/93‑1 (MQ=255)
                                                      gcggcgTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCCTGAAGGAGTGTATTGCGGTTGGTGACGGGATGAacg                >  1:6401998/1‑93 (MQ=255)
                                                       cggcgTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAacga               <  2:641982/93‑1 (MQ=255)
                                                       cggcgTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAacga               >  1:4095124/1‑93 (MQ=255)
                                                       cggcgTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAacga               <  1:10033062/93‑1 (MQ=255)
                                                        ggcgTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAacgac              <  2:1182648/93‑1 (MQ=255)
                                                        ggcgTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAacgac              >  1:2739865/1‑93 (MQ=255)
                                                        ggcgTTTCAAAAGGCCATGCGCTGGAAGCGGGGGCGAAGAAA‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAacgac              >  1:2868802/1‑93 (MQ=255)
                                                         gcgTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAacgacg             >  2:3783878/1‑93 (MQ=255)
                                                         gcgTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAacgacg             >  2:5876638/1‑93 (MQ=255)
                                                         gcgTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAacgacg             >  1:9759049/1‑93 (MQ=255)
                                                          cgTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGc            <  1:648139/93‑1 (MQ=255)
                                                          cgTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGc            <  1:4552484/93‑1 (MQ=255)
                                                           gTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGcc           <  1:6408410/93‑1 (MQ=255)
                                                           gTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGcc           >  1:7202564/1‑93 (MQ=255)
                                                           gTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGcc           >  1:2422851/1‑93 (MQ=255)
                                                           gTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGcc           >  1:8630491/1‑93 (MQ=255)
                                                           gTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCGTTNGGTGACGGGATGAACGACGcg           >  1:8760513/1‑92 (MQ=255)
                                                           gTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCG‑‑‑‑AAGAATCTGGGCTACAGCCTGAAGGATTGTATTGCGTTGGGTGACGGGATGAACGACGcc           <  2:9151104/93‑1 (MQ=255)
                                                             ttcaaaaGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGCCGCCGa         >  2:3275039/5‑93 (MQ=255)
                                                             ttcaaaaGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGa         >  1:2638946/5‑93 (MQ=255)
                                                             ttcaaaaGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGa         >  1:4004406/5‑93 (MQ=255)
                                                             ttcaaaaGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGa         >  1:4615904/5‑93 (MQ=255)
                                                              tcaaaaGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGaa        >  2:5931017/4‑93 (MQ=255)
                                                              tcaaaaGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGaa        <  2:776121/90‑1 (MQ=255)
                                                              tcaaaaGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGaa        >  2:6112688/4‑93 (MQ=255)
                                                              tcaaaaGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGaa        >  2:7328160/4‑93 (MQ=255)
                                                              tcaaaaGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGaa        >  2:7064425/4‑93 (MQ=255)
                                                              tcaaaaGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGaa        <  1:1680756/90‑1 (MQ=255)
                                                               taaaaGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGaaa       <  1:8741687/91‑1 (MQ=255)
                                                               caaaaGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGaaa       <  1:2353815/91‑1 (MQ=255)
                                                               caaaaGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGaaa       <  2:6411942/91‑1 (MQ=255)
                                                               caaaaGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGaaa       <  2:6578561/91‑1 (MQ=255)
                                                               caaaaGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGaaa       <  2:318095/91‑1 (MQ=255)
                                                               caaaaGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGaaa       <  2:6578560/91‑1 (MQ=255)
                                                                aaaaGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAAt      >  2:6351003/2‑93 (MQ=255)
                                                                aaaaGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAAt      >  1:865711/2‑93 (MQ=255)
                                                                aaaaGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAAt      >  1:798070/2‑93 (MQ=255)
                                                                 aaaGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATg     >  2:2286238/1‑93 (MQ=255)
                                                                 aaaGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATg     >  2:6404630/1‑93 (MQ=255)
                                                                 aaaGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATg     >  1:160295/1‑93 (MQ=255)
                                                                 aaaGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATg     <  2:5419476/93‑1 (MQ=255)
                                                                 aaaGGCCATGCGCTGGAAGCGGTGGCG‑‑‑‑AAGAATCTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATg     <  1:181652/93‑1 (MQ=255)
                                                                  aaGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGc    <  2:7954981/93‑1 (MQ=255)
                                                                  aaGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGc    <  1:2939580/93‑1 (MQ=255)
                                                                   aGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGc  >  2:7669372/1‑92 (MQ=255)
                                                                   aGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCt   >  1:5588573/1‑93 (MQ=255)
                                                                   aGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCt   <  1:2133690/93‑1 (MQ=255)
                                                                   aGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCt   >  1:5588572/1‑93 (MQ=255)
                                                                   aGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCt   <  1:9889411/93‑1 (MQ=255)
                                                                    ggCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTg  <  2:7244493/93‑1 (MQ=255)
                                                                                            ||||                                                                      
GCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAGGCCATGCGCTGGAAGGCTGGCG‑‑‑‑AA‑‑‑‑CTGGCTACAGCTGAAGGATTGATTGCGTT‑‑GTGACGGGATGAACGACGCCGAAATGCTG  >  NC_000913/4008309‑4008458

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 20 ≤ ATCG/ATCG < 28 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: