Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 808,438 T→A 6.4% K274* (AAA→TAA)  bioA ← adenosylmethionine‑8‑amino‑7‑oxononanoate aminotransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913808,4380TA 6.4% 352.4 / 10.1 109K274* (AAA→TAA) bioAadenosylmethionine‑8‑amino‑7‑oxononanoate aminotransferase
Reads supporting (aligned to +/- strand):  ref base T (53/49);  new base A (4/3);  total (57/52)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.82e-01

Only 100 of 111 total aligned reads displayed.

GCTGGCGTTTGCTGCCGCGCAGGCCAGCGGATTGCCCATAAAAGTTGGCCCATGCATAAAGCAACCGGCTTCACCGTTACTGATGGTTTCTGCAACCTCGCGCGTGGTGAGTGTGGCGGAAAGGGTCATTGTGCCGCCGGTTAAGGCTTTACCGAGGCACAAAATGTCCGGCGCGATTTCTGCATGTTCACAGGCAAACAGTTTCCCGGTACGACCAAATCCAGTGGCGATCTCGTCGGCAATCAGCAAGATACCTTCGCGATCGCATATTTTGCGGATTCGTTTTAACCATTCC  >  NC_000913/808289‑808583
                                                                                                                                                     |                                                                                                                                                  
gCTGGCGTTTGCTGCCGCGCAGGCCAGCGGATTGCCGATAAAAGTTGGCCCATGCATAAAGCAACCGGCTTCACCGTTACTGATGGTTTCTGCAACCTCGCGCGTGGTGAGTGTGGCGGAAAGGGTCATTGTGCCGCCGGTTAAGGCttt                                                                                                                                                    <  2:865556/150‑1 (MQ=255)
   ggCGTTTGCTGCCGCGCAGGCCAGCGGATTGCCCATAAAAGTTGGCCCATGCATAAAGCAACCGGCTTCACCGTTACTGATGGTTTCTGCAACCTCGCGCGTGGTGAGTGTGGCGGAAAGGGTCATTGTGCCGCCGGTTAAGGCTTTAcc                                                                                                                                                 >  1:1100024/1‑150 (MQ=255)
   ggCGTTTGCTGCCGCGCAGGCCAGCGGATTGCCCATAAAAGTTGGCCCATGCATAAAGCAACCGGCTTCACCGTTACTGATGGTTTCTGCAACCTCGCGCGTGGTGAGTGTGGCGGAAAGGGTCATTGTGCCGCCGGTTAAGGCTTTAcc                                                                                                                                                 >  1:1498176/1‑150 (MQ=255)
      gTTTGCTGCCGCGCAGGCCAGCGGATTGCCCATAAAAGTTGGCCCATGCATAAAGCAACCGGCTTCACCGTTACTGATGGTTTCTGCAACCTCGCGCGTGGTGAGTGTGGCGGAAAGGGTCATTGTGCCGCCGGTTAAGGCTTTACCGAg                                                                                                                                             >  1:1621269/1‑150 (MQ=255)
       tttGCTGCCGCGCAGGCCAGCGGATTGCCCATAAAAGTTGGCCCATGCATAAAGCAACCGGCTTCACCGTTACTGATGGTTTCTGCAACCTCGCGCGTGGTGAGTGTGGCGGAAAGGGTCATTGTGCCGCCGGTTAAGGCTTTACCGAgg                                                                                                                                            <  2:602764/150‑1 (MQ=255)
        ttGCTGCCGCGCAGGCCAGCGGATTGCCCATAAAAGTTGGCCCATGCATAAAGCAACCGGCTTCACCGTTACTGATGGTTTCTGCAACCTCGCGCGTGGTGAGTGTGGCGGAAAGGGTCATTGTGCCGCCGGTTAAGGCTTTACCGAGGc                                                                                                                                           <  2:645152/150‑1 (MQ=255)
         tgctgcCGCGCAGGCCAGCGGATTGCCCATAAAAGTTGGCCCATGCATAAAGCAACCGGCTTCACCGTTACTGATGGTTTCTGCAACCTCGCGCGTGGTGAGTGTGGCGGAAAGGGTCATTGTGCCGCCGGTTAAGGCTTTACCGAGGca                                                                                                                                          >  2:930299/1‑150 (MQ=255)
          gctgcCGCGCAGGCCAGCGGATTGCCCATAAAAGTTGGCCCATGCATAAAGCAACCGGCTTCACCGTTACTGATGGTTTCTGCAACCTCGCGCGTGGTGAGTGTGGCGGAAAGGGTCATTGTGCCTCGGTAAAGCCTTAACCggcggcac                                                                                                                                        <  1:1034007/150‑8 (MQ=255)
          gctgcCGCGCAGGCCAGCGGATTGCCCATAAAAGTTGGCCCATGCATAAAGCAACCGGCTTCACCGTTACTGATGGTTTCTGCAACCTCGCGCGTGGTGAGTGTGGCGGAAAGGGTCATTGTGCCTCGGTAAAGCCTTAACCggcggcac                                                                                                                                        <  1:139925/150‑8 (MQ=255)
           ctgcCGCGCAGGCCAGCGGATTGCCCATAAAAGTTGGCCCATGCATAAAGCAACCGGCTTCACCGTTACTGATGGTTTCTGCAACCTCGCGCGTGGTGATTGTGGCGGAAAGGGTCATTGTGCCGCCGGTTAAGGCTTTACCGAGGcaca                                                                                                                                        >  2:1323606/1‑150 (MQ=255)
           ctgcCGCGCAGGCCAGCGGATTGCCCATAAAAGTTGGCCCATGCATAAAGCAACCGGCTTCACCGTTACTGATGGTTTCTGCAACCTCGCGCGTGGTGAGTGTGGCGGAAAGGGTCATTGTGCCGCCGGTTAAGGCTTTACCGAGGcaca                                                                                                                                        >  1:787715/1‑150 (MQ=255)
            tgcCGCGCAGGCCAGCGGATTGCCCATAAAAGTTGGCCCATGCATAAAGCAACCGGCTTCACCGTTACTGATGGTTTCTGCAACCTCGCGCGTGGTGAGTGTGGCGGAAAGGGTCATTGTGCCGCCGGTTAAGGCTTTACCGAGGCACaa                                                                                                                                       >  1:878003/1‑150 (MQ=255)
             gcCGCGCAGGCCAGCGGATTGCCCATAAAAGTTGGCCCATGCATAAAGCAACCGGCTTCACCGTTACTGATGGTTTCTGCAACCTCGCGCGTGGTGAGTGTGGCGGAAAGGGTCATTGTGCCGCCGGTTAAGGCTTTACCGAGGCACaaa                                                                                                                                      >  1:1104580/1‑150 (MQ=255)
                  gcAGGCCAGCGGATTGCCCATAAAAGTTGGCCCATGCATAAAGCAACCGGCTTCACCGTTACTGATGGTTTCTGCAACCTCGCGCGTGGTGAGTGTGGCGGAAAGGGTCATTGTGCCGCCGGTTAAGGCTTTACCGAGGCACAAAATGTc                                                                                                                                 <  1:1178781/150‑1 (MQ=255)
                        cAGCGGATTGCCCATAAAAGTTGGCCCATGCATAAAGCAACCGGCTTCACCGTTACTGATGGTTTCTGCAACCTCGCGCGTGGTGAGTGTGGCGGAAAGGGTCATTGTGCCGCCGGTTAAGGCTTTACCGAGGCACAAAATGTCCGgcgc                                                                                                                           >  2:824660/1‑150 (MQ=255)
                         aGCAGATTGCCCATAAAAGTTGGCCCATGCATAAAGCAACCGGCTTCACCGTTACTGATGGTTTCTGCAACCTCGCGCGTGGTGAGTGTGGCGGAAAGGGTCATTGTGCCGCCGGTTAAGGCTTTACCGAGGCACAAAATGTCCGGcgcg                                                                                                                          <  2:1380152/150‑1 (MQ=255)
                            ggATTGCCCATAAAAGTTGGCCCATGCATAAAGCAACCGGCTTCACCGTTACTGATGGTTTCTGCAACCTCGCGCGTGGTGAGTGTGGCGGAAAGGGTCATTGTGCCGCCGGTTAAGGCTTTACCGAGGCACAAAATGTCCGGCGCGAtt                                                                                                                       >  1:1648966/1‑150 (MQ=255)
                                 gaacATAAAAGATGCCCCATACATAAACCAAAAGGCTTCACCGATACTGATTGTTTCTGCAACCTCGCGCGTGGTGAGTGTGGCGGAAAGGGTCATTGTGCCGCCGGTTAAGGCTTTACCGAGGCACAAAATGTCCGGCGCGATTTCTGc                                                                                                                  >  2:230772/4‑150 (MQ=255)
                                 gCCCATAAAAGTTGGCCCATGCATAAAGCAACCGGCTTCACCGTTACTGATGGTTTCTGCAACCTCGCGCGTGGTGAGTGTGGCGGAAAGGGTCATTGTGCCGCCGGTTAAGGCTTTACCGAGGCACAAAATGTCCGGCGCGATTTCTGc                                                                                                                  >  1:1545705/1‑150 (MQ=255)
                                 gCCCATAAAAGTTGGCCCATGCATAAAGCAACCGGCTTCACCGTTACTGATGGTTTCTGCAACCTCGCGCGTGGTGAGTGTGGCGGAAAGGGTCATTGTGCCGCCGGTTAAGGCTTTACCGAGGCACAAAATGTCCGGCGCGATTTCTGc                                                                                                                  >  2:1647364/1‑150 (MQ=255)
                                 gCCCATAAAAGTTGGCCCATGCATAAAGCAACCGGCTTCACCGTTACTGATGGTTTCTGCAACCTCGCGCGTGGTGAGTGTGGCGGAAAGGGTCATTGTGCCGCCGGTTAAGGCTTTACCGAGGCACAAAATGTCCGGCGCGATTTCTGc                                                                                                                  >  2:1559335/1‑150 (MQ=255)
                                  cccATAAAAGTTGGCCCATGCATAAAGCAACCGGCTTCACCGTTACTGATGGTTTCTGCAACCTCGCGCGTGGTGAGTGTGGCGGAAAGGGTCATTGTGCCGCCGGTTAAGGCTTTACCGAGGCACAAAATGTCCGGCGCGATTTCTGCa                                                                                                                 >  2:795879/1‑150 (MQ=255)
                                  cccATAAAAGTTGGCCCATGCATAAAGCAACCGGCTTCACCGTTACTGATGGTTTCTGCAACCTCGCGCGTGGTGAGTGTGGCGGAAAGGGTCATTGTGCCGCCGGTTAAGGCTTTACCGAGGCACAAAATGTCCGGCGCGATTTCTGCa                                                                                                                 <  1:1229893/150‑1 (MQ=255)
                                  cccATAAAAGTCGGCCCATGCATAAAGCAACCGGCTTCACCGTTACTGATGGTTTCTGCAACCTCGCGCGTGGTGAGTGTGGCGGAAAGGGTCATTGTGCCGCCGGTTAAGGCTTTACCGAGGCACAAAATGTCCGGCGCGATTTCTGCa                                                                                                                 <  1:903254/150‑1 (MQ=255)
                                   ccATAAAAGTTGGCCCATGCATAAAGCAACCGGCTTCACCGTTACTGATGGTTTCTGCAACCTCGCGCGTGGTGAGTGTGGCGGAAAGGGTCATTGTGCCGCCGGTTAAGGCTTTACCGAGGCACAAAATGTCCGGCGCGATTTCTGCAt                                                                                                                >  2:465558/1‑150 (MQ=255)
                                    cATAAAAGTTGGCCCATGCATAAAGCAACCGGCTTCACCGTTACTGATGGTTTCTGCAACCTCGCGCGTGGTGAGTGTGGCGGAAAGGGTCATTGTGCCGCCGGTTAAGGCTTTACCGAGGCACAAAATGTCCGGCGCGATTTCTGCATg                                                                                                               >  1:1373275/1‑150 (MQ=255)
                                    cATAAAAGTTGGCCCATGCATAAAGCAACCGGCTTCACCGTTACTGATGGTTTCTGCAACCTCGCGCGTGGTGAGTGTGGCGGAAAGGGTCATTGTGCCGCCGGTTAAGGCTTTACCGAGGCACAAAATGTCCGGCGCGATTTCTGCATg                                                                                                               >  1:1513297/1‑150 (MQ=255)
                                     aTAAAAGTTGGCCCATGCATAAAGCAACCGGCTTCACCGTTACTGATGGTTTCTGCAACCTCGCGCGTGGTGAGTGTGGCGGAAAGGGTCATTGTGCCGCCGGTTAAGGCTTTACCGAGGCACAAAATGTCCGGCGCGATTTCTGCATGt                                                                                                              >  2:1244546/1‑150 (MQ=255)
                                      tAAAAGTTGGCCCATGCATAAAGCAACCGGCTTCACCGTTACTGATGGTTTCTGCAACCTCGCGCGTGGTGAGTGTGGCGGAAAGGGTCATTGTGCCGCCGGTTAAGGCTTTACCGAGGCACAAAATGTCCGGCGCGATTTCTGCATGtt                                                                                                             >  2:927400/1‑150 (MQ=255)
                                       aaaaGTTGGCCCATGCATAAAGCAACCGGCTTCACCGTTACTGATGGTTTCTGCAACCTCGCGCGTGGTGAGTGTGGCGGAAAGGGTCATTGTGCCGCCGGTTAAGGCTTTACCGAGGCACAAAATGTCCGGCGCGATTTCTGCATGTTc                                                                                                            <  1:732205/150‑1 (MQ=255)
                                        aaaGTTGGCCCATGCATAAAGCAACCGGCTTCACCGTTACTGATGGTTTCTGCAACCTCGCGCGTGGTGAGTGTGGCGGAAAGGGTCATTGTGCCGCCGGTTAAGGCTTTACCGAGGCACAAAATGTCCGGCGCGATTTCTGCATGTTca                                                                                                           >  1:346715/1‑150 (MQ=255)
                                            ttGGCCCATGCATAAAGCAACCGGCTTCACCGTTACTGATGGTTTCTGCAACCTCGCGCGTGGTGAGTGTGGCGGAAAGGGTCATTGTGCCGCCGGTTAAGGCTTTACCGAGGCACAAAATGTCCGGCGCGATTTCTGCATGTTCACAgg                                                                                                       <  1:1244546/150‑1 (MQ=255)
                                             tGGCCCATGCATAAAGCAACCGGCTTCACCGTTACTGATGGTTTCTGCAACCTCGCGCGTGGTGAGTGTGGCGGAAAGGGTCATTGTGCCGCCGGTTAAGGCTTTACCGAGGCACAAAATGTCCGGCGCGATTTCTGCATGTTCACAGGc                                                                                                      <  1:639212/150‑1 (MQ=255)
                                               gCCCATGCATAAAGCAACCGGCTTCACCGTTACTGATGGTTTCTGCAACCTCGCGCGTGGTGAGTGTGGCGGAAAGGGTCATTGTGCCGCCGGTTAAGGCTTTACCGAGGCACAAAATGTCCGGCGCGATTTCTGCATGTTCACAGGCaa                                                                                                    <  1:1236385/150‑1 (MQ=255)
                                                 ccATGCATAAAGCAACCGGCTTCACCGTTACTGATGGTTTCTGCAACCTCGCGCGTGGTGAGTGTGGCGGAAAGGGTCATTGTGCCGCCGGTTAAGGCTTTACCGAGGCACAAAATGTCCGGCGCGATTTCTGCATGTTCACAGGCAAAc                                                                                                  <  2:297283/150‑1 (MQ=255)
                                                 ccATGCATAAAGCAACCGGCTTCACCGTTACTGATGGTTTCTGCAACCTCGCGCGTGGTGAGTGTGGCGGAAAGGGTCATTGTGCCGCCGGTTAAGGCTTTACCGAGGCACAAAATGTCCGGCGCGATTTCTGCATGTTCACAGGCAAAc                                                                                                  <  1:1617904/150‑1 (MQ=255)
                                                    tGCATAAAGCAACCGGCTTCACCGTTACTGATGGTTTCTGCAACCTCGCGCGTGGTGAGTGTGGCGGAAAGGGTCATTGTGCCGCCGGTTAAGGCTTTACCGAGGCACAAAATGTCCGGCGCGATTTCTGCATGTTCACAGGCAAACAGt                                                                                               >  2:101834/1‑150 (MQ=255)
                                                        tAAAGCAACCGGCTTCACCGTTACTGATGGTTTCTGCAACCTCGCGCGTGGTGAGTGTGGCGGAAAGGGTCATTGTGCCGCCGGTTAAGGCTTTACCGAGGCACAAAATGTCCGGCGCGATTTCTGCATGTTCACAGGCAAACAGTTTcc                                                                                           <  2:344810/150‑1 (MQ=255)
                                                         aaaGCAACCGGCTTCACCGTTACTGATGGTTTCTGCAACCTCGCGCGTGGTGAGTGTGGCGGAAAGGGTCATTGTGCCGCCGGTTAAGGCTTTACCGAGGCACAAAATGTCCGGCGCGATTTCTGCATGTTCACAGGCAAACAGTTTccc                                                                                          <  2:1648966/150‑1 (MQ=255)
                                                         aaaGCAACCGGCTTCACCGTTACTGATGGTTTCTGCAACCTCGCGCGTGGTGAGTGTGGCGGAAAGGGTCATTGTGCCGCCGGTTAAGGCTTTACCGAGGCACAAAATGTCCGGCGCGATTTCTGCATGTTCACAGGCAAACAGTTTccc                                                                                          <  1:1270275/150‑1 (MQ=255)
                                                         aaaGCAACCGGCTTCACCGTTACTGATGGTTTCTGCAACCTCGAGCGTGGTGAGTGTGGCGGAAAGGGTCATTGTGCCGCCGGTTAAGGCTTTACCGACGCACAAAATGTCCGGCGCGATTTCTGCATGTTCACAGGCAAACAGTTTccc                                                                                          >  2:1428645/1‑150 (MQ=255)
                                                          aaGCAACCGGCTTCACCGTTACTGATGGTTTCTGCAACCTCGCGCGTGGTGAGTGTGGCGGAAAGGGTCATTGTGCCGCCGGTTAAGGCTTTACCGAGGCACAAAATGTCCGGCGCGATTTCTGCATGTTCACAGGCAAACAGTTTCCCg                                                                                         >  2:611680/1‑150 (MQ=255)
                                                           aGCAACCGGCTTCACCGTTACTGATGGTTTCTGCAACCTCGCGCGTGGTGAGTGTGGCGGAAAGGGTCATTGTGCCGCCGGTTAAGGCTTTACCGAGGCACAAAATGTCCGGCGCGATTTCTGCATGTTCACAGGCAAACAGTTTCCCgg                                                                                        <  2:1498176/150‑1 (MQ=255)
                                                            gCAACCGGCTTCACCGTTACTGATGGTTTCTGCAACCTCGCGCGTGGTGAGTGTGGCGGAAAGGGTCATTGTGCCGCCGGTTAAGGCTTTACCGAGGCACAAAATGTCCGGCGCGATTTCTGCATGTTCACAGGCAAACAGTTTCCCGGt                                                                                       <  2:376370/150‑1 (MQ=255)
                                                            gCAACCGGCTTCACCGTTACTGATGGTTTCTGCAACCTCGCGCGTGGTGAGTGTGGCGGAAAGGGTCATTGTGCCGCCGGTTAAGGCTTTACCGAGGCACAAAATGTCCGGCGCGATTTCTGCATGTTCACAGGCAAACAGTTTCCCGGt                                                                                       <  2:1246503/150‑1 (MQ=255)
                                                             cAACCGGCTTCACCGTTACTGATGGTTTCTGCAACCTCGCGCGTGGTGAGTGTGGCGGAAAGGGTCATTGTGCCGCCGGTTAAGGCTTTACCGAGGCACAAAATGTCCGGCGCGATTTCTGCATGTTCACAGGCAAACAGTTTCCCGGTa                                                                                      <  2:1621269/150‑1 (MQ=255)
                                                             cAACCGGCTTCACCGTTACTGATGGTTTCTGCAACCTCGCGCGTGGTGAGTGTGGCGGAAAGGGTCATTGTGCCGCCGGTTAAGGCTTTACCGAGGCACAAAATGTCCGGCGCGATTTCTGCATGTTCACAGGCAAACAGTTTCCCGGTa                                                                                      >  2:133382/1‑150 (MQ=255)
                                                               aCCGGCTTCACCGTTACTGATGGTTTCTGCAACCTCGCGCGTGGTGAGTGTGGCGGAAAGGGTCATTGTGCCGCCGGTTAAGGCTTTACCGAGGCACAAAATGTCCGGCGCGATTTCTGCATGTTCACAGGCAAACAGTTTCCCGGTACg                                                                                    <  2:908871/150‑1 (MQ=255)
                                                               aCCGGCTTCACCGTTACTGATGGTTTCTGCAACCTCGCGCGTGGTGAGTGTGGCGGAAAGGGTCATTGTGCCGCCGGTTAAGGCTTTACCGAGGCACAAAATGTCCGGCGCGATTTCTGCATGTTCACAGGCAAACAGTTTCCCGGTACg                                                                                    <  1:1123082/150‑1 (MQ=255)
                                                               aCCGGCTTCACCGTTACTGATGGTTTCTGCAACCTCGCGCGTGGTGAGTGTGGCGGAAAGGGTCATTGTGCCGCCGGTTAAGGCTTTACCGAGGCACAAAATGTCCGGCGCGATTTCTGCATGTTCACAGGCAAACAGTTTCCCGGTACg                                                                                    >  2:308359/1‑150 (MQ=255)
                                                                  ggCTTCACCGTTACTGATGGTTTCTGCAACCTCGCGCGTGGTGAGTGTGGCGGAAAGGGTCATTGTGCCGCCGGTTAAGGCTTTACCGAGGCACAAAATGTCCGGCGCGATTTCTGCATGTTCACAGGCAAACAGTTTCCCGGTACGAcc                                                                                 <  1:328996/150‑1 (MQ=255)
                                                                    cTTCACCGTTACTGATGGTTTCTGCAACCTCGCGCGTGGTGAGTGTGGCGGAAAGGGTCATTGTGCCGCCGGTTAAGGCTTTACCGAGGCACAAAATGTCCGGCGCGATTTCTGCATGTTCACAGGCAAACAGTTTCCCGGTACGACCaa                                                                               >  1:686186/1‑150 (MQ=255)
                                                                      tCACCGTTACTGATGGTTTCTGCAACCTCGCGCGTGGTGAGTGTGGCGGAAAGGGTCATTGTGCCGCCGGTTAAGGCTTTACCGAGGCACAAAATGTCCGGCGCGATTTCTGCATGTTCACAGGCAAACAGTTTCCCGGTACGACCAAAt                                                                             >  2:1068797/1‑150 (MQ=255)
                                                                      tCACCGTTACTGATGGTTTCTGCAACCTCGCGCGTGGTGAGTGTGGCGGAAAGGGTCATTGTGCCGCCGGTTAAGGCTTTACCGAGGCACAAAATGTCCGGCGCGATTTCTGCATGTTCACAGGCAAACAGTTTCCCGGTACGACCAAAt                                                                             <  1:546870/150‑1 (MQ=255)
                                                                         ccGTTACTGATGGTTTCTGCAACCTCGCGCGTGGTGAGTGTGGCGGAAAGGGTCATTGTGCCGCCGGTTAAGGCTTTACCGAGGCACAAAATGTCCGGCGCGATTTCTGCATGTTCACAGGCAAACAGTTTCCCGGTACGACCAAATCCa                                                                          >  2:1278838/1‑150 (MQ=255)
                                                                          cGTTACTGATGGTTTCTGCAACCTCGCGCGTGGTGAGTGTGGCGGAAAGGGTCATTGTGCCGCCGGTTAAGGCTTTACCGAGGCACAAAATGTCCGGCGCGATTTCTGCATGTTCACAGGCAAACAGTTTCCCGGTACGACCAAATCCAg                                                                         >  1:48463/1‑150 (MQ=255)
                                                                            ttACTGATTGTTTCTGCAACCTCGCGCGTGGTGAGTGTGGCGGAAAGGGTCATTGTGCCGCCGGTTAAGGCTTTACCGAGGCACAAAATGTCCGGCGCGATTTCTGCATGTTCACAGGCAAACAGTTTCCCGGTACGACCAAATCCAGTg                                                                       <  1:996291/150‑1 (MQ=255)
                                                                                    ggTTTCTGCAACCTCGCGCGTGGTGAGTGTGGCGGAAAGGGTCATTGTGCCGCCGGTTAAGGCTTTACCGAGGCACAAAATGTCCGGCGCGATTTCTGCATGTTCACAGGCAAACAGTTTCCCGGTACGACCAAATCCAGTGGCGAtctc                                                               >  1:777677/1‑150 (MQ=255)
                                                                                        tCTGCAACCTCGCGCGTGGTGAGTGTGGCGGAAAGGGTCATTGTGCCGCCGGTTAAGGCTTTACCGAGGCACAAAATGTCCGGCGCGATTTCTGCATGTTCACAGGCAAACAGTTTCCCGGTACGACCAAATCCAGTGGCGATCtcgtcg                                                           <  1:1647364/150‑1 (MQ=255)
                                                                                         cTGCAACCTCGCGCGTGGTGAGTGTGGCGGAAAGGGTCATTGTGCCGCCGGTTAAGGCTTTACCGAGGCACAAAATGTCCGGCGCGATTTCTGCATGTTCACAGGCAAACAGTTTCCCGGTACGACCAAATCCAGTGGCGATCTCGTCgg                                                          >  2:1107791/1‑150 (MQ=255)
                                                                                         cTGCAACCTCGCGCGTGGTGAGTGTGGCGGAAAGGGTCATTGTGCCGCCGGTTAAGGCTTTACCGAGGCACAAAATGTCCGGCGCGATTTCTGCATGTTCACAGGCAAACAGTTTCCCGGTACGACCAAATCCAGTGGCGATCTCGTCgg                                                          >  2:803884/1‑150 (MQ=255)
                                                                                         cTGCAACCTCGCGCGTGGTGAGTGTGGCGGAAAGGGTCATTGTGCCGCCGGTTAAGGCTTTACCGAGGCACAAAATGTCCGGCGCGATTTCTGCATGTTCACAGGCAAACAGTTTCCCGGTACGACCAAATCCAGTGGCGATCTCGTCgg                                                          <  1:1313690/150‑1 (MQ=255)
                                                                                          tGCAACCTCGCGCGTGGTGAGTGTGGCGGAAAGGGTCATTGTGCCGCCGGTTAAGGCTTTACCGAGGCACAAAATGTCCGGCGCGATTTCTGCATGTTCACAGGCAAACAGTTTCCCGGTACGACCAAATCCAGTGGCGATCTCGTCGGc                                                         <  1:308359/150‑1 (MQ=255)
                                                                                          tGCAACCTCGCGCGTGGTGAGTGTGGCGGAAAGGGTCATTGTGCCGCCGGTTAAGGCTTTACCGAGGCACAAAATGTCCGGCGCGATTTCTGCATGTTCACAGGCAAACAGTTTCCCGGTACGACCAAATCCAGTGGCGATCTCGTCGGc                                                         <  1:410770/150‑1 (MQ=255)
                                                                                           gCAACCTCGCGCGTGGTGAGTGTGGCGGAAAGGGTCATTGTGCCGCCGGTTAAGGCTTTACCGAGGCACAAAATGTCCGGCGCGATTTCTGCATGTTCACAGGCAAACAGTTTCCCGGTACGACCAAATCCAGTGGCGATCTCGTCGGCa                                                        >  2:1121946/1‑150 (MQ=255)
                                                                                           gCAACCTCGCGCGTGGTGAGTGTGGCGGAAAGGGTCATTGTGCCGCCGGTTAAGGCTTTACCGAGGCACAAAATGTCCGGCGCGATTTCTGCATGTTCACAGGCAAACAGTTTCCCGGTACGACCAAATCCAGTGGCGATCTCGTCGGCa                                                        >  2:327749/1‑150 (MQ=255)
                                                                                            cAACCTCGCGCGTGGTGAGTGTGGCGGAAAGGGTCATTGTGCCGCCGGTTAAGGCTTTACCGAGGCACAAAATGTCCGGCGCGATTTCTGCATGTTCACAGGCAAACAGTTTCCCGGTACGACCAAATCCAGTGGCGATCTCGTCGGCaa                                                       >  1:919522/1‑150 (MQ=255)
                                                                                               ccTCGCGCGTGGTGAGTGTGGCGGAAAGGGTCATTGTGCCGCCGGTTAAGGCTTTACCGAGGCACAAAATGTCCGGCGCGATTTCTGCATGTTCACAGGCAAACAGTTTCCCGGTACGACCAAATCCAGTGGCGATCTCGTCGGCAATCa                                                    >  2:262656/1‑150 (MQ=255)
                                                                                               ccTCGCGCGTGGTGAGTGTGGCGGAAAGGGTCATTGTGCCGCCGGTTAAGGCTTTACCGAGGCACAAAATGTCCGGCGCGATTTCTGCATGTTCACAGGCAAACAGTTTCCCGGTACGACCAAATCCAGTGGCGATCTCGTCGGCAATCa                                                    <  1:1123171/150‑1 (MQ=255)
                                                                                                  cgcgcgTGGTGAGTGTGGCGGAAAGGGTCATTGTGCCGCCGGTTAAGGCTTTACCGAGGCACAAAATGTCCGGCGCGATTTCTGCATGTTCACAGGCAAACAGTTTCCCGGTACGACCAAATCCAGTGGCGATCTCGTCGGCAATCAGCa                                                 >  2:1049929/1‑150 (MQ=255)
                                                                                                        tggtgAGTGTGGCGGAAAGGGTCATTGTGCCGCCGGTTAAGGCTTTACCGAGGCACAAAATGTCCGGCGCGATTTCTGCATGTTCACAGGCAAACAGTTTCCCGGTACGACCAAATCCAGTGGCGATCTCGTCGGCAATCAGCAAGATAc                                           >  1:661090/1‑150 (MQ=255)
                                                                                                             aGTGTGGCGGAAAGGGTCATTGTGCCGCCGGTTAAGGCTTTACCGAGGCACAAAATGTCCGGCGCGATTTCTGCATGTTCACAGGCAAACAGTTTCCCGGTACGACCAAATCCAGTGGCGATCTCGTCGGCAATCAGCAAGATACCTTcg                                      >  1:651542/1‑150 (MQ=255)
                                                                                                              gtgtGGCGGAAAGGGTCATTGTGCCGCCGGTTAAGGCTTTACCGAGGCACAAAATGTCCGGCGCGATTTCTGCATGTTCACAGGCAAACAGTTTCCCGGTACGACCAAATCCAGTGGCGATCTCGTCGGCAATCAGCAAGATACCTTcgc                                     <  1:927400/150‑1 (MQ=255)
                                                                                                                gtgGCGGAAAGGGTCATTGTGCCGCCGGTTAAGGCTTTACCGAGGCACAAAATGTCCGGCGCGATTTCTGCATGTTCACAGGCAAACAGTTTCCCGGTACGACCAAATCCAGTGGCGATCTCGTCGGCAATCAGCAAGATACCTTCGCGa                                   >  2:970843/1‑150 (MQ=255)
                                                                                                                  ggCGGAAAGGGTCATTGTGCCGCCGGTTAAGGCTTTACCGAGGCACAAAATGTCCGGCGCGATTTCTGCATGTTCACAGGCAAACAGTTTACCGGTACGACCAAATCCAGTGGAGATCTCGTCGGCAATCAGCAAGATACCTGCGCGATc                                 >  2:570693/1‑150 (MQ=255)
                                                                                                                   gCGGAAAGGGTCATTGTGCCGCCGGTTAAGGCTTTACCGAGGCACAAAATGTCCGGCGCGATTTCTGCATGTTCACAGGCAAACAGTTTCCCGGTACGACCAAATCCAGTGGCGATCTCGTCGGCAATCAGCAAGATACCTTCGCGATCg                                <  1:176529/150‑1 (MQ=255)
                                                                                                                   gCGGAAAGGGTCATTGTGCCGCCGGTTAAGGCTTTACCGAGGCACAAAATGTCCGGCGCGATTTCTGCATGTTCACAGGCAAACAGTTTCCCGGTACGACCAAATCCAGTGGCGATCTCGTCGGCAATCAGCAAGATACCTTCGCGATCg                                >  1:178561/1‑150 (MQ=255)
                                                                                                                     ggAAAGGGTCATTGTGCCGCCGGTTAAGGCTTTACCGAGGCACAAAATGTCCGGCGCGATTTCTGCATGTTCACAGGCAAACAGTTTCCCGGTACGACCAAATCCAGTGGCGATCtcgtcg                                                           >  2:1668808/1‑121 (MQ=255)
                                                                                                                     ggAAAGGGTCATTGTGCCGCCGGTTAAGGCTTTACCGAGGCACAAAATGTCCGGCGCGATTTCTGCATGTTCACAGGCAAACAGTTTCCCGGTACGACCAAATCCAGTGGCGATCtcgtcg                                                           <  1:1668808/121‑1 (MQ=255)
                                                                                                                      gAAAGGGTCATTGTGCCGCCGGTTAAGGCTTTACCGAGGCACAAAATGTCCGGCGCGATTTCTGCATGTTCACAGGCAAACAGTTTCCCGGTACGACCAAATCCAGTGGCGATCTCGTCGGCAATCAGCAAGATACCTTCGCGATCGCat                             <  1:327749/150‑1 (MQ=255)
                                                                                                                          gggTCATTGTGCCGCCGGTTAAGGCTTTACCGAGGCACAAAATGTCCGGCGCGATTTCTGCATGTTCACAGGCAAACAGTTTCCCGGTACGACCAAATCCAGTGGCGATCTCGTCGGCAATCAGCAAGATACCTTCGCGATCGCATAttt                         >  1:856873/1‑150 (MQ=255)
                                                                                                                            cattttGTGCCTCGGTAAAGCCTTAACCGGCGGCACAAAATGTCCGGCGCGATTTCTGCATGTTCACAGGCAAACAGTTTCCCGGTACGACCAAATCCAGTGGCGATCTCGTCGGCAATCAGCAAGATACCTTCGCGATCGCATATTTTg                       >  2:282545/5‑150 (MQ=38)
                                                                                                                            cattttGTGCCTCGGTAAAGCCTTAACCGGCGGCACAAAATGTCCGGCGCGATTTCTGCATGTTCACAGGCAAACAGTTTCCCGGTACGACCAAATCCAGTGGCGATCTCGTCGGCAATCAGCAAGATACCTTCGCGATCGCATATTTTg                       >  2:209283/5‑150 (MQ=38)
                                                                                                                              ttttGTGCCTCGGTAAAGCCTTAACCGGCGGCACAAAATGTCCGGCGCGATTTCTGCATGTTCACAGGCAAACAGTTTCCCGGTACGACCAAATCCAGTGGCGATCTCGTCGGCAATCAGCAAGATACCTTCGCGATCGCATATTTTGCg                     >  1:1451388/3‑150 (MQ=38)
                                                                                                                              cATTGTGCCGCCGGTTAAGGCTTTACCGAGGCACAAAATGTCCGGCGCGATTTCTGCATGTTCACAGGCAAACAGTTTCCCGGTACGACCAAATCCAGTGGCGATCTCGTCGGCAATCAGCAAGATACCTTCGCGATCGCATATTTTGCg                     <  1:1121946/150‑1 (MQ=255)
                                                                                                                              cATTGTGCCGCCGGTTAAGGCTTTACCGAGGCACAAAATGTCCGGCGCGATTTCTGCATGTTCACAGGCAAACAGTTTCCCGGTACGACCAAATCCAGTGGCGATCTCGTCGGCAATCAGCAAGATACCTTCGCGATCGCATATTTTGCg                     <  1:879351/150‑1 (MQ=255)
                                                                                                                              cATTGTGCCGCCGGTTAAGGCTTTACCGAGGCACAAAATGTCCGGCGCGATTTCTGCATGTTCACAGGCAAACAGTTTCCCGGTACGACCAAATCCAGTGGCGATCTCGTCGGCAATCAGCAAGATACCTTCGCGATCGCATATTTTGCg                     <  1:824660/150‑1 (MQ=255)
                                                                                                                              cATTGTGCCGCCGGTTAAGGCTTTACCGAGGCACAAAATGTCCGGCGCGATTTCTGCATGTTCACAGGCAAACAGTTTCCCGGTACGACCAAATCCAGTGGCGATCTCGTCGGCAATCAGCAAGATACCTTCGCGATCGCATATTTTGCg                     <  2:434958/150‑1 (MQ=255)
                                                                                                                              cATTGTGCCGCCGGTTAAGGCTTTACCGAGGCACAAAATGTCCGGCGCGATTTCTGCATGTTCACAGGCAAACAGTTTCCCGGTACGACCAAATCCAGTGGCGATCTCGTCGGCAATCAGCAAGATACCTTCGCGATCGCATATTTTGCg                     <  2:48463/150‑1 (MQ=255)
                                                                                                                                   gtgcctcgGTAAAGCCTTAACCGGCGGCACAAAATGTCCGGCGCGATTTCTGCATGTTCACAGGCAAACAGTTTCCCGGTACGACCAAATCCAGTGGCGATCTCGTCGGCAATCAGCAAGATACCTTCGCGATCGCATATTTTGCGGAtt                 >  2:293061/7‑150 (MQ=255)
                                                                                                                                    gccgccGGTTAAGGCTTTACCGAGGCACAAAATGTCCGGCGCGATTTCTGCATGTTCACAGGCAAACAGTTTCCCGGTACGACCAAATCCAGTGGCGATCTCGTCGGCAATCAGCAAGATACCTTCGCGATCGCATATTTTCCGGATTCg               >  1:1665246/1‑150 (MQ=255)
                                                                                                                                     ccgccgGTTAAGGCTTTACCGAGGCACAAAATGTCCGGCGCGATTTCTGCATGTTCACAGGCAAACAGTTTCCCGGTACGACCAAATCCAGTGGCGATCTCGTCGGCAATCAGCAAGATACCTTCGCGATCGCATATTTTGCGGATTCGt              >  2:543046/1‑150 (MQ=255)
                                                                                                                                     ccgccgGTTAAGGCTTTACCGAGGCACAAAATGTCCGGCGCGATTTCTGCATGTTCACAGGCAAACAGTTTCCCGGTACGACCAAATCCAGTGGCGATCTCGTCGGCAATCAGCAAGATACCTTCGCGATCGCATATTTTGCGGATTCGt              <  1:130868/150‑1 (MQ=255)
                                                                                                                                     ccgccgGTTAAGGCTTTACCGAGGCACAAAATGTCCGGCGCGATTTCTGCATGTTCACAGGCAAACAGTTTCCCGGTACGACCAAATCCAGTGGCGATCTCGTCGGCAATCAGCAAGATACCTTCGCGATCGCATATTTTGCGGATTCGt              >  2:112585/1‑150 (MQ=255)
                                                                                                                                        ccgGTTAAGGCTTTACCGAGGCACAAAATGTCCGGCGCGATTTCTGCATGTTCACAGGCAAACAGTTTCCCGGTACGACCAAATCCAGTGGCGATCTCGTCGGCAATCAGCAAGATACCTTCGCGATCGCATATTTTGCGGATTCGtttt           <  2:878003/150‑1 (MQ=255)
                                                                                                                                           gTTAAGGCTTTACCGAGGCACAAAATGTCCGGCGCGATTTCTGCATGTTCACAGGCAAACAGTTTCCCGGTACGACCAAATCCAGTGGCGATCTCGTCGGCAATCAGCAAGATACCTTCGCGATCGCATATTTTGCGGATTCGTTTTAAc        >  1:441871/1‑150 (MQ=255)
                                                                                                                                            ttAAGGCTTTACCGAGGCACAAAATGTCCGGCGCGATTTCTGAATGTTCACAGGCAAACAGTTTCCCGGTACGACCAAATCCAGTGGCGATCTCGTCGGCAATCAGCAAGATACCTTCGCGATCGCATATTTTGCGGATTCGTTTTAAcc       >  2:1445144/1‑150 (MQ=255)
                                                                                                                                                 gCTTTACCGAGGCACAAAATGTCCGGCGCGATTTCTGCATGTTCACAGGCAAACAGTTTCCCGGTACGACCAAATCCAGTGGCGATCTCGTCGGCAATCAGCAAGATACCTTCGCGATCGCATATTTTGCGGATTCGTTTTAACCATTcc  <  2:625286/150‑1 (MQ=255)
                                                                                                                                                 gCTTTACCGAGGCACAAAATGTCCGGCGCGATTTCTGCATGTTCACAGGCAAACAGTTTCCCGGTACGACCAAATCCAGTGGCGATCTCGTCGGCAATCAGCAAGATACCTTCGCGATCGCATATTTTGCGGATTCGTTTTAACCATTcc  <  1:348764/150‑1 (MQ=255)
                                                                                                                                                 gCTTTACCGAGGCACAAAATGTCCGGCGCGATTTCTGCATGTTCACAGGCAAACAGTTTCCCGGTACGACCAAATCCAGTGGCGATCTCGTCGGCAATCAGCAAGATACCTTCGCGATCGCATATTTTGCGGATTCGTTTTAACCATTcc  <  1:133382/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
GCTGGCGTTTGCTGCCGCGCAGGCCAGCGGATTGCCCATAAAAGTTGGCCCATGCATAAAGCAACCGGCTTCACCGTTACTGATGGTTTCTGCAACCTCGCGCGTGGTGAGTGTGGCGGAAAGGGTCATTGTGCCGCCGGTTAAGGCTTTACCGAGGCACAAAATGTCCGGCGCGATTTCTGCATGTTCACAGGCAAACAGTTTCCCGGTACGACCAAATCCAGTGGCGATCTCGTCGGCAATCAGCAAGATACCTTCGCGATCGCATATTTTGCGGATTCGTTTTAACCATTCC  >  NC_000913/808289‑808583

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