Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,449,607:1 +G 5.7% coding (443/822 nt) yfcO ← DUF2544 domain‑containing protein YfcO

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,449,6071.G 5.7% 436.3 / 21.0 106coding (443/822 nt)yfcODUF2544 domain‑containing protein YfcO
Reads supporting (aligned to +/- strand):  ref base . (54/46);  new base G (4/2);  total (58/48)
Fisher's exact test for biased strand distribution p-value = 6.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.12e-01

Only 100 of 109 total aligned reads displayed.

TCCCATCTTTCATCAAATACGTTGAAAGCGTTTTCTCAACCACACCGCCCGTCGACTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAACCTCAAATGTATAATACCCGCAGTTCCACCCAGACGGTAATCCA  >  NC_000913/2449458‑2449756
                                                                                                                                                      |                                                                                                                                                     
tCCCATCTTTCATCAAATACGTTGAAAGCGTTTTCTCAACCACACCGCCCGTCGACTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAc                                                                                                                                                        <  2:62304/150‑1 (MQ=255)
   cATCTTTCATCAAATACGTTGAAAGCGTTTTCTCAACCACACCGCCCGTCGACTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTAt                                                                                                                                                    >  1:1430935/1‑150 (MQ=255)
   cATCTTTCATCAAATACGTTGAAAGCGTTTTCTCAACCACACCGCCCGTCGACTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTAt                                                                                                                                                    >  2:1493103/1‑150 (MQ=255)
       tttCATCAAATACGTTGAAAGCGTTTTCTCAACCACACCGCCCGTCGACTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAgg                                                                                                                                                >  2:582925/1‑150 (MQ=255)
        ttCATCAAATACGTTGAAAGCGTTTTCTCAACCACACCGCCCGTCGACTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACGTATCAgg                                                                                                                                                >  1:86297/1‑150 (MQ=255)
         tcatcaAATACGTTGAAAGCGTTTTCTCAACCACACCGCCCGTCGACTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGcc                                                                                                                                              <  2:430909/150‑1 (MQ=255)
          catcaAATACGTTGAAAGCGTTTTCTCAACCACACCGCCCGTCGACTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGccc                                                                                                                                             >  2:1047269/1‑150 (MQ=255)
          catcaAATACGTTGAAAGCGTTTTCTCAACCACACCGCCCGTCGACTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGccc                                                                                                                                             >  1:1205214/1‑150 (MQ=255)
          catcaAATACGTTGAAAGCGTTTTCTCAACCACACCGCCCGTCGACTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGccc                                                                                                                                             >  1:1009976/1‑150 (MQ=255)
          aatcaAATACGTTGAAAGCGTTTTCTCAACCACACCGCCCGTCGACTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGccc                                                                                                                                             <  1:414161/149‑1 (MQ=255)
              aaaTACGTTGAAAGCGTTTTCTCAACCACACCGCCCGTCGACTGCAGGGTGAGCAGTTTTGACTTTCTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGa                                                                                                                                         >  2:58887/1‑150 (MQ=255)
               aaTACGTTGAAAGCGTTTTCTCAACCACACCGCCCGTCGACTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTTTCGTGACTATCAGGCCCCTGAt                                                                                                                                        <  2:392268/150‑1 (MQ=255)
               aaTACGTTGAAAGCGTTTTCTCAACCACACCGCCCGTCGACTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGAt                                                                                                                                        <  2:961722/150‑1 (MQ=255)
                  aCGTTGAAAGCGTTTTCTCAACCACACCGCCCGTCGCCTGCAGGGTGAGCAGTTTTGACTAACTGCCATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACg                                                                                                                                     >  2:1585798/1‑150 (MQ=255)
                      tGAAAGCGTTTTCTCAACCACACCGCCCGTCGACTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCAc                                                                                                                                 <  2:1461894/150‑1 (MQ=255)
                       gAAAGCGTTTTCTCAACCACACCGCCCGTCGACTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACg                                                                                                                                >  1:109962/1‑150 (MQ=255)
                           gCGTTTTCTCAACCACACCGCCCGTCGACTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGcccc                                                                                                                            >  1:601097/1‑150 (MQ=255)
                            cGTTTTCTCAACCACACCGCCCGTCGACTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGccccc                                                                                                                           <  2:1120371/150‑1 (MQ=255)
                               tttCTCAACCACACCGCCCGTCGACTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTgg                                                                                                                        <  2:1292013/150‑1 (MQ=255)
                                 tctcAACCACACCGCCCGTCGACTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGGCAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGgag                                                                                                                      >  1:1546781/1‑150 (MQ=255)
                                  ctcAACCACACCGCCCGTCGACTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGgaga                                                                                                                     >  1:1032610/1‑150 (MQ=255)
                                   tcAACCACACCGCCCGTCGACTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGagag                                                                                                                    >  2:1014595/1‑150 (MQ=255)
                                     aaCCACACCGCCCGTCGACTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGc                                                                                                                  <  1:975040/150‑1 (MQ=255)
                                      aCCACACCGCCCGTCGACTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGGTACGTCACGCCCCCTGGAGAGGCt                                                                                                                 >  1:547553/1‑150 (MQ=255)
                                       ccACACCGCCCGTCGACTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCtt                                                                                                                >  1:1049890/1‑150 (MQ=255)
                                       ccACACCGCCCGTCGACTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCtt                                                                                                                >  2:623003/1‑150 (MQ=255)
                                             cGCCCGTCGACTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGt                                                                                                          <  1:672188/150‑1 (MQ=255)
                                              gCCCGTCGACTGCAGGGTGAGGAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGtc                                                                                                         <  2:1029937/150‑1 (MQ=255)
                                               cccGTCGACTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCag                                                                                                             <  1:756860/145‑1 (MQ=255)
                                               cccGTCGACTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCag                                                                                                             >  2:756860/1‑145 (MQ=255)
                                                 cGTCGACTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGtcatt                                                                                                      <  2:658101/150‑1 (MQ=255)
                                                  gTCGACTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGtcattc                                                                                                     >  2:1058229/1‑150 (MQ=255)
                                                       cTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCataat                                                                                                <  2:1292877/150‑1 (MQ=255)
                                                       cTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCataat                                                                                                <  1:1218300/150‑1 (MQ=255)
                                                         gCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCAtaataa                                                                                              >  1:1308341/1‑150 (MQ=255)
                                                         gCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCAtaataa                                                                                              >  1:1017217/1‑150 (MQ=255)
                                                         gCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCAtaataa                                                                                              >  2:667890/1‑150 (MQ=255)
                                                           aGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATaacc                                                                                            >  2:552142/1‑150 (MQ=255)
                                                              gTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATaaccaac                                                                                         >  1:336564/1‑150 (MQ=255)
                                                              gTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATaaccaac                                                                                         >  1:1549660/1‑150 (MQ=255)
                                                               tGAGCAGTTTTGACTTACTGACATAATTCTCATTCTATGAAACGTCATAAGGTTGAGCAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATaaccaacc                                                                                        <  2:115401/150‑1 (MQ=255)
                                                                gAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAaccaacca                                                                                       >  1:65381/1‑150 (MQ=255)
                                                                  gCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATg                                                                                     <  2:1269182/150‑1 (MQ=255)
                                                                          gACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGaa                                                                             >  2:365768/1‑150 (MQ=255)
                                                                           aCTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGaaa                                                                            <  1:582925/150‑1 (MQ=255)
                                                                            cTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGaaaa                                                                           >  1:1345569/1‑150 (MQ=255)
                                                                            cTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGaaaa                                                                           >  1:246505/1‑150 (MQ=255)
                                                                              tACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGt                                                                         >  2:676335/1‑150 (MQ=255)
                                                                               aCTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGGCAAATCGTGTCGTGACTATCAGGCCTCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTa                                                                        >  1:938552/1‑150 (MQ=255)
                                                                                  gACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAg                                                                               <  1:663513/140‑1 (MQ=255)
                                                                                  gACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGa                                                                     >  1:329008/1‑150 (MQ=255)
                                                                                  gACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGa                                                                     >  1:543340/1‑150 (MQ=255)
                                                                                  gACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGa                                                                     >  1:1304636/1‑150 (MQ=255)
                                                                                    cATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGAtt                                                                   <  1:1040633/150‑1 (MQ=255)
                                                                                       aaTTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCaa                                                                >  2:920739/1‑150 (MQ=255)
                                                                                         ttCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAAtc                                                              >  1:379984/1‑150 (MQ=255)
                                                                                         ttCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTCACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAAtc                                                              >  1:590624/1‑150 (MQ=255)
                                                                                            tcATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATctct                                                           <  2:1549660/150‑1 (MQ=255)
                                                                                              aTTCCATTAAACGTCATAAGGTTGAACACAAACGGATGTACAAATCGTGTCGTGACTATCAGGCCCCTTATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAAAAACCAACCATGGAAAAGAAAAGTAAGATTCAATCTCTGc                                                         <  2:1009976/150‑1 (MQ=255)
                                                                                              aTTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGc                                                         <  2:65381/150‑1 (MQ=255)
                                                                                              aTTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATCGACAAGAAAAGTAAGATTCAATCTATGc                                                         >  1:727696/1‑150 (MQ=255)
                                                                                               ttCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCa                                                        <  1:1014595/150‑1 (MQ=255)
                                                                                               ttCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCa                                                        <  1:408693/150‑1 (MQ=255)
                                                                                                tCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCag                                                       <  2:71390/150‑1 (MQ=255)
                                                                                                tCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCag                                                       <  2:246505/150‑1 (MQ=255)
                                                                                                  cATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAAtct                                                             >  2:868190/1‑142 (MQ=255)
                                                                                                  cATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAAtct                                                             <  1:868190/142‑1 (MQ=255)
                                                                                                  cATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCagag                                                     >  1:143298/1‑150 (MQ=255)
                                                                                                  cATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCagag                                                     <  2:601097/150‑1 (MQ=255)
                                                                                                   aTGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAgaga                                                    >  2:316816/1‑150 (MQ=255)
                                                                                                    tGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGAc                                                   <  1:1027789/150‑1 (MQ=255)
                                                                                                      aaaCGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACaa                                                 >  2:676909/1‑150 (MQ=255)
                                                                                                          gTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAAc                                             <  1:676909/150‑1 (MQ=255)
                                                                                                              tAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAACCTCa                                         >  2:1322838/1‑150 (MQ=255)
                                                                                                              tAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAACCTCa                                         >  1:1184321/1‑150 (MQ=255)
                                                                                                               aaGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAACCTCaa                                        >  2:730835/1‑150 (MQ=255)
                                                                                                                aGGTTGAACAGCAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCAGGAGAGGCTTCAGAGTCATTCATAATCACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAACCTCaaa                                       >  2:742367/1‑150 (MQ=255)
                                                                                                                 ggTTGAACAGAAACGGATGGGCAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAACCTCAAAt                                      <  2:22624/150‑1 (MQ=255)
                                                                                                                  gTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAACCTCAAATg                                     <  1:809945/150‑1 (MQ=255)
                                                                                                                  gTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAACCTCAAATg                                     <  1:1346969/150‑1 (MQ=255)
                                                                                                                  gTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAACCTCAAATg                                     >  1:1584891/1‑150 (MQ=255)
                                                                                                                        cAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAACCTCAAATGTataat                               >  1:768700/1‑150 (MQ=255)
                                                                                                                               ggATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAACCTCAAATGTATAATACCCGCa                        <  2:125157/150‑1 (MQ=255)
                                                                                                                                gATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAACCTCAAATGTATAATACCCGCAg                       <  2:1345569/150‑1 (MQ=255)
                                                                                                                                 aTGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAACCTCAAATGTATAATACCCGCAGt                      >  2:1190439/1‑150 (MQ=255)
                                                                                                                                 aTGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAACCTCAAATGTATAATACCCGCAGt                      <  2:1017217/150‑1 (MQ=255)
                                                                                                                                    gACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAACCTCAAATGTATAATACCCGCAGTTcc                   <  2:1032610/150‑1 (MQ=255)
                                                                                                                                         aTCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGGCAATAACCTCAAATGTATAATACCCGCAGTTCCACCCa              <  1:1058229/150‑1 (MQ=255)
                                                                                                                                             gggcgtgACGTATCAGGGGCCTGATAGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAACCTCAAATGTATAATACCCGCAGTTCCACCCAGACg          >  1:918922/4‑150 (MQ=255)
                                                                                                                                              gtcgtgACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAACCTCAAATGTATAATACCCGCAGTTCCACCCAGACgg         <  2:143298/150‑1 (MQ=255)
                                                                                                                                               tcgtgACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAACCTCAAATGTATAATACCCGCAGTTCCACCCAGACGGt        <  2:1502553/150‑1 (MQ=255)
                                                                                                                                                cgtgACGTATCAGGGGCCTGATAGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAACCTCAAATGt                                    >  1:754160/1‑121 (MQ=255)
                                                                                                                                                cgtgACGTATCAGGGGCCTGATAGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAACCTCAAATGt                                    <  2:754160/121‑1 (MQ=255)
                                                                                                                                                cgtgACGTATCAGGGGCCTGATAGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAACCTCAAATGTATAATACCCGCAGTTCCAACCAGACGGTa       >  2:445383/1‑150 (MQ=255)
                                                                                                                                                cgtgACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACGAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATACCCTCAAATGTATAATACCCGCAGTTCCACCCAGACGGTa       >  2:1440718/1‑150 (MQ=255)
                                                                                                                                                cgtgACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCGGAGACAAt                                                <  2:1051006/109‑1 (MQ=255)
                                                                                                                                                cgtgACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCGGAGACAAt                                                >  1:1051006/1‑109 (MQ=255)
                                                                                                                                                 gtgACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAACCTCAAATGTATAATACCCGCAGTTCCACCCAGACGGTaa      >  1:1090753/1‑150 (MQ=255)
                                                                                                                                                     cTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAACCTCAAATGTATAATACCCGCAGTTCCACCCAGACGGTAATCCa  <  2:1476909/150‑1 (MQ=255)
                                                                                                                                                     cTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAACCTCAAATGTATAATACCCGCAGTTCCACCCAGACGGTAATCCa  >  2:244593/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                                     
TCCCATCTTTCATCAAATACGTTGAAAGCGTTTTCTCAACCACACCGCCCGTCGACTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAACCTCAAATGTATAATACCCGCAGTTCCACCCAGACGGTAATCCA  >  NC_000913/2449458‑2449756

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: