Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,942,711 A→G intergenic (‑144/+207) gcvA ← / ← ygdI glycine cleavage system transcriptional activator; autorepressor/DUF903 family verified lipoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,942,7110AG98.7% 248.1 / ‑4.4 74intergenic (‑144/+207)gcvA/ygdIglycine cleavage system transcriptional activator; autorepressor/DUF903 family verified lipoprotein
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (36/37);  total (36/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.79e-01

GACATGGCTATTAAACTTTGATGTTAAATGAATTTAACAATTAGATCACACTATGTAACCTATTAGTTTTTTTAATCTGAGCCATTATAAATTGTCCGTTGAGCTTCTACCAGCAAATACCTATAGTGGCGGCACTTCCTGAGCCGGAACGAAAAGTTTTATCGGAATGCGTGTTCTGGTGAACTTTTGGCTTACGGTTGTGATGTTGTGTTGTTGTGTTTGCAATTGGTCTGCGATTCAGACCATGGTAGCAAAGCTACCTTTTTTCACTTCCTGTACATTTACCCTGTCTGTCC  >  NC_000913/2942563‑2942858
                                                                                                                                                    |                                                                                                                                                   
gACATGGCTATTAAACTTTGATGTTAAATGAATTTAACAATTAGATCACACTATGTAACCTATTAGTTTTTTTAATCTGAGCCATTATAAATTGTCCGTTGAGCTTCTACCAGCAAATACCTATAGTGGCGGCACTTCCTGAGCCGGAGCg                                                                                                                                                   >  2:791230/1‑151 (MQ=255)
  cATGGCTATTAAACTTTGATGTTAAATGAATTTAACAATTAGATCACACTATGTAACCTATTAGTTTTTTTAATCTGAGCCATTATAAATTGTCCGTTGAGCTTCTACCAGCAAATACCTATAGTGGCGGCACTTCCTGAGCCGGAGCGaa                                                                                                                                                 <  2:133013/151‑1 (MQ=255)
  cATGGCTATTAAACTTTGATGTTAAATGAATTTAACAATTAGATCACACTATGTAACCTATTAGTTTTTTTAATCTGAGCCATTATAAATTGTCCGTTGAGCTTCTACCAGCAAATACCTATAGTGGCGGCACTTCCTGAGCCGGAGCGa                                                                                                                                                  >  2:755526/1‑150 (MQ=255)
   aTGGCTATTAAACTTTGATGTTAAATGAATTTAACAATTAGATCACACTATGTAACCTATTAGTTTTTTTAATCTGAGCCATTATAAATTGTCCGTTGAGCTTCTACCAGCAAATACCTATAGTGGCGGCACTTCCTGAGCCGGAGCGaaa                                                                                                                                                <  2:285323/151‑1 (MQ=255)
   aTGGCTATTAAACTTTGATGTTAAATGAATTTAACAATTAGATCACACTATGTAACCTATTAGTTTTTTTAATCTGAGCCATTATAAATTGTCCGTTGAGCTTCTACCAGCAAATACCTATAGTGGCGGCACTTCCTGAGCCGGAGCGaa                                                                                                                                                 >  1:115989/1‑150 (MQ=255)
    tGGCTATTAAACTTTGATGTTAAATGAATTTAACAATTAGATCACACTATGTAACCTATTAGTTTTTTTAATCTGAGCCATTATAAATTGTCCGTTGAGCTTCTACCAGCAAATACCTATAGTGGCGGCACTTCCTGAGCCGGAGCGaaa                                                                                                                                                >  1:446770/1‑150 (MQ=255)
     ggCTATTAAACTTTGATGTTAAATGAATTTAACAATTAGATCACACTATGTAACCTATTAGTTTTTTTAATCTGAGCCATTATAAATTGTCCGTTGAGCTTCTACCAGCAAATACCTATAGTGGCGGCACTTCCTGAGCCGGAGCGaaaa                                                                                                                                               <  2:105166/150‑1 (MQ=255)
        tATTAAACTTTGATGTTAAATGAATTTAACAATTAGATCACACTATGTAACCTATTAGTTTTTTTAATCTGAGCCATTATAAATTGTCCGTTGAGCTTCTACCAGCAAATACCTATAGTGGCGGCACTTCCTGAGCCGGAGCGAAAAGttt                                                                                                                                           >  1:611600/1‑151 (MQ=255)
         aTTAAACTTTTATGTTAAATGAATTTAACAATTAAGTCACTCTATGTAACCTATTAGTTTTTTTAATCTGAGCCATTATAAATTGTCCGTTGAGCTTCTACCAGCAAATACCTATAGTGGCGGCACTTCCTGAGCCGGAACGAAAAGttt                                                                                                                                           <  2:197532/150‑1 (MQ=255)
         aTTAAACTTTGATGTTAAATGAATTTAACAATTAGATCACACTATGTAACCTATTAGTTTTTTTAATCTGAGCCATTATAAATTGTCCGTTGAGCTTCTACCAGCAAATACCTATAGTGGCGGCACTTCCTGAGCCGGAGCGAAAAGtttt                                                                                                                                          >  1:498450/1‑151 (MQ=255)
             aaCTTTGATGTTAAATGAATTTAACAATTAGATCACACTATGTAACCTATTAGTTTTTTTAATCTGAGCCATTATAAATTGTCCGTTGAGCTTCTACCAGCAAATACCTATAGTGGCGGCACTTCCTGAGCCGGAGCGAAAAGTTTTATCg                                                                                                                                      <  2:323597/151‑1 (MQ=255)
              aCTTTGATGTTAAATGAATTTAACAATTAGATCACACTATGTAACCTATTAGTTTTTTTAATCTGAGCCATTATAAATTGTCCGTTGAGCTTCTACCAGCAAATACCTATAGTGGCGGCACTTCCTGAGCCGGAGCGAAAAGTTTTATCgg                                                                                                                                     >  1:231940/1‑151 (MQ=255)
               cTTTGATGTTAAATGAATTTAACAATTAGATCACACTATGTAACCTATTAGTTTTTTTAATCTGAGCCATTATAAATTGTCCGTTGAGCTTCTACCAGCAAATACCTATAGTGGCGGCACTTCCTGAGCCGGAGCGAAAAGTTTTATCgg                                                                                                                                     >  2:634762/1‑150 (MQ=255)
                   gATGTTAAATGAATTTAACAATTAGATCACACTATGTAACCTATTAGTTTTTTTAATCTGAGCCATTATAAATTGTCCGTTGAGCTTCTACCAGCAAATACCTATAGTGGCGGCACTTCCTGAGCCGGAGCGAAAAGTTTTATCGGAATGc                                                                                                                                <  1:658246/151‑1 (MQ=255)
                     tGTTAAATGAATTTAACAATTAGATCACACTATGTAACCTATTAGTTTTTTTAATCTGAGCCATTATAAATTGTCCGTTGAGCTTCTACCAGCAAATACCTATAGTGGCGGCACTTCCTGAGCCGGAGCGAAAAGTTTTATCGGAATGCgt                                                                                                                              >  2:385237/1‑151 (MQ=255)
                       ttAAATGAATTTAACAATTAGATCACACTATGTAACCTATTAGTTTTTTTAATCTGAGCCATTATAAATTGTCCGTTGAGCTTCTACCAGCAAATACCTATAGTGGCGGCACTTCCTGAGCCGGAGCGAAAAGATTTATCGGAATGCgtg                                                                                                                             >  1:213885/1‑150 (MQ=255)
                        tAAATGAATTTAACAATTAGATCACACTATGTAACCTATTAGTTTTTTTAATCTGAGCCATTATAAATTGTCCGTTGAGCTTCTACCAGCAAATACCTATAGTGGCGGCACTTCCTGAGCCGGAGCGAAAAGTTTTATCGGAATGCGTGtt                                                                                                                           >  1:231574/1‑151 (MQ=255)
                        tAAATGAATTTAACAATTAGATCACACTATGTAACCTATTAGTTTTTTTAATCTGAGCCATTATAAATTGTCCGTTGAGCTTCTACCAGCAAATACCTATAGTGGCGGCACTTCCTGAGCCGGAGCGAAAAGTTTTATCGGAATGCGTGtt                                                                                                                           >  1:258127/1‑151 (MQ=255)
                             gAATTTAACAATTAGATCACACTATGTAACCTATTAGTTTTTTTAATCTGAGCCATTATAAATTGTCCGTTGAGCTTCTACCAGCAAATACCTATAGTGGCGGCACTTCCTGAGCCGGAGCGAAAAGTTTTATCGGAATGCGTGTTCTgg                                                                                                                       <  1:488423/150‑1 (MQ=255)
                             gAATTTAACAATTAGATCACACTATGTAACCTATTAGTTTTTTTAATCTGAGCCATTATAAATTGTCCGTTGAGCTTCTACCAGCAAATACCTATAGTGGCGGCACTTCCTGAGCCGGAGCGAAAAGTTTTATCGGAATGCGTGTTCTgg                                                                                                                       >  1:430690/1‑150 (MQ=255)
                               aTTTAACAATTAGATCACACTATGTAACCTATTAGTTTTTTTAATCTGAGCCATTATAAATTGTCCGTTGAGCTTCTACCAGCAAATACCTATAGTGGCGGCACTTCCTGAGCCGGAGCGAAAAGTTTTATCGGAATGCGTGTTCTGGTGa                                                                                                                    <  1:296597/151‑1 (MQ=255)
                                    aCAATTAGATCACACTATGTAACCTATTAGTTTTTTTAATCTGAGCCATTATAAATTGTCCGTTGAGCTTCTACCAGCAAATACCTATAGTGGCGGCACTTCCTGAGCCGGAGCGAAAAGTTTTATCGGAATGCGTGTTCTGGTGAACttt                                                                                                               >  2:245280/1‑151 (MQ=255)
                                    aCAATTAGATCACACTATGTAACCTATTAGTTTTTTTAATCTGAGCCATTATAAATTGTCCGTTGAGCTTCTACCAGCAAATACCTATAGTGGCGGCACTTCCTGAGCCGGAGCGAAAAGTTTTATCGGAATGCGTGTTCTGGTGAACttt                                                                                                               >  2:494824/1‑151 (MQ=255)
                                         tAGATCACACTATGTAACCTATTAGTTTTTTTAATCTGAGCCATTATAAATTGTCCGTTGAGCTTCTACCAGCAAATACCTATAGTGGCGGCACTTCCTGAGCCGGAGCGAAAAGTTTTATCGGAATGCGTGTTCTGGTGAACTTTTGGCt                                                                                                          >  2:57273/1‑151 (MQ=255)
                                                 acTATGTAACCTATTAGTTTTTTTAATCTGAGCCATTATAAATTGTCCGTTGAGCTTCTACCAGCAAATACCTATAGTGGCGGCACTTCCTGAGCCGGAGCGAAAAGTTTTATCGGAATGCGTGTTCTGGTGAACTTTTGGCTTACGGTtg                                                                                                  >  2:798564/1‑151 (MQ=255)
                                                    aTGTAACCTATTAGTTTTTTTAATCTGAGCCATTATAAATTGTCCGTTGAGCTTCTACCAGCAAATACCTATAGTGGCGGCACTTCCTGAGCCGGAGCGAAAAGTTTTATCGGAATGCGTGTTCTGGTGAACTTTTGGCTTACGGTtgtg                                                                                                <  2:152505/150‑1 (MQ=255)
                                                           cTATTAGTTTTTTTAATCTGAGCCATTATAAATTGTCCGTTGAGCTTCTACCAGCAAATACCTATAGTGGCGGCACTTCCTGAGCCGGAGCGAGAAGTTTTATCGGAATGCGTGTTCTGGTGAACTTTTGGCTTACGGTTGTGAtgttgt                                                                                         >  2:708950/1‑150 (MQ=255)
                                                             aTTAGTTTTTTTAATCTGAGCCATTATAAATTGTCCGTTGAGCTTCTACCAGCAAATACCTATAGTGGCGGCACTTCCTGAGCCGGAGCGAGAAGTTTTATCGGAATGCGTGTTCTGGTGAACTTTTGGCTTACGGTTGTGAtgttgt                                                                                         <  1:708906/148‑1 (MQ=255)
                                                             aTTAGTTTTTTTAATCTGAGCCATTATAAATTGTCCGTTGAGCTTCTACCAGCAAATACCTATAGTGGCGGCACTTCCTGAGCCGGAGCGAAAAGTTTTATCGGAATGCGTGTTCTGGTGAACTTTTGGCTTACGGTTGTGAtgttgtgtt                                                                                      <  1:45071/151‑1 (MQ=255)
                                                             aTTAGTTTTTTTAATCTGAGCCATTATAAATTGTCCGTTGAGCTTCTACCAGCAAATACCTATAGTGGCGGCACTTCCTGAGCCGGAGCGAAAAGTTTTATCGGAATGCGTGTTCTGGTGAACTTTTGGCTTACGGTTGTGAtgttgtgt                                                                                       <  2:79377/150‑1 (MQ=255)
                                                             aTTAGTTTTTTTAATCTGAGCCATTATAAATTGTCCGTTGAGCTTCTACCAGCAAATACCTATAGTGGCGGCACTTCCTGAGCCGGAGCGAAAAGTTTTATCGGAATGCGTGTTCTGGTGAACTTTTGGCTTACGGTTGTGAtgttgtg                                                                                        <  2:654317/149‑1 (MQ=255)
                                                               tAGTTTTTTTAATCTGAGCCATTATAAATTGTCCGTTGAGCTTCTACCAGCAAATACCTATAGTGGCGGCACTTCCTGAGCCGGAGCGAAAAGTTTTATCGGAATGCGTGTTCTGGTGAACTTTTGGCTTACGGTTGTGAtgttgtgttgt                                                                                    >  2:48635/1‑151 (MQ=255)
                                                               tAGTTTTTTTAATCTGAGCCATTATAAATTGTCCGTTGAGCTTCTACCAGCAAATACCTATAGTGGCGGCACTTCCTGAGCCGGAGCGAAAAGTTTTATCGGAATGCGTGTTCTGGTGAACTTTTGGCTTACGGTTGTGAtgttgtgttg                                                                                     <  2:540727/150‑1 (MQ=255)
                                                                 gTTTTTTTAATCTGAGCCATTATAAATTGTACGTTGAGCTTCTACCAGCGAATACCTATAGTGGCGGCACGTCCTGAGACGGGGGGAAAGGCTTTATCGGAATGCGTGTTCTGGTGGACTTTTGGCTTATTGTTGTGAtgttgtgttgt                                                                                    >  1:148720/1‑149 (MQ=255)
                                                                    tttttAATCTGAGCCATTATAAATTGTCCGTTGAGCTTCTACCAGCAAATACCTATAGTGGCGGCACTTCCTGAGCCGGAGCGAAAAGTTTTATCGGAATGCGTGTTCTGGTGAACTTTTGGCTTACGGTTGTGAtgttgtgttgttgtg                                                                                <  2:38308/150‑1 (MQ=255)
                                                                      tttAATCTGAGCCATTATAAATTGTCCGTTGAGCTTCTACCAGCAAATACCTATAGTGGCGGCACTTCCTGAGCCGGAGCGAAAAGTTTTATCGGAATGCGTGTTCTGGTGAACTTTTGGCTTACGGTTGTGATGTTGTGTTGTTGTGttt                                                                             >  1:387955/1‑151 (MQ=255)
                                                                        tAATCTGAGCCATTATAAATTGTCCGTTGAGCTTCTACCAGCAAATACCTATAGTGGCGGCACTTCCTGAGCCGGAGCGAAAAGTTTTATCGGAATGCGTGTTCTGGTGAACTTTTGGCTTACGGTTGTGATGTTGTGTTGTTGTGTTTGc                                                                           >  2:224417/1‑151 (MQ=255)
                                                                         aaTCTGAGCCATTATAAATTGTCCGTTGAGCTTCTACCAGCAAATACCTATAGTGGCGGCACTTCCTGAGCCGGAGCGAAAAGTTTTATCGGAATGCGTGTTCTGGTGAACTTTTGGCTTACGGTTGTGATGTTGTGTTGTTGTGTTTGc                                                                           >  1:735897/1‑150 (MQ=255)
                                                                              gAGCCATTATAAATTGTCCGTTGAGCTTCTACTAGCAAATACCTATAGTGGCGGCACTTCCTGAGCCGGAGCGAAAAGTTTTATCGGAATGCGTGTTCTGGTGAACTTTTGGCTTACGGTTGTGATGTTGTGTTGTTGTGTTTGCAATTg                                                                      <  1:585192/150‑1 (MQ=255)
                                                                              gAGCCATTATAAATTGTCCGTTGAGCTTCTACCAGCAAATACCTATAGTGGCGGCTCTTCCTGAGCCGGAGCGAAAAGTTTTATCGGAATGCGTGTTCTGGTGAACTTTTGGCTTACGGTTGTGATGTTGTGTTGTTGTGTTTGCAAt                                                                        <  1:406476/148‑1 (MQ=255)
                                                                              gAGCCATTATAAATTGTCCGTTGAGCTTCTACCAGCAAATACCTATAGTGGCGGCACTTCCTGAGCCGGAGCGAAAAGTTTTATCGGAATGCGTGTTCTGGTGAACTTTTGGCTTACGGTTGTGATGTTGTGTTGTTGTGTTTGCAATTgg                                                                     <  2:855228/151‑1 (MQ=255)
                                                                                 ccATTATAAATTGTCCGTTGAGCTTCTACCAGCAAATACCTATAGTGGCGGCACTTCCTGAGCCGGAGCGAAAAGTTTTATCGGAATGCGTGTTCTGGTGAACTTTTGGCTTACGGTTGTGATGTTGTGTTGTTGTGTTTGCAATTGGTCt                                                                  >  1:442327/1‑151 (MQ=255)
                                                                                 ccATTATAAATTGTCCGTTGAGCTTCTACCAGCAAATACCTATAGTGGCGGCACTTCCTGAGCCGGAGCGAAAAGTTTTATCGGAATGCGTGTTCTGGTGAACTTTTGGCTTACGGTTGTGATGTTGTGTTGTTGTGTTTGCAATTGGTCt                                                                  <  1:45826/151‑1 (MQ=255)
                                                                                    ttATAAATTGTCCGTTGAGCTTCTACCAGCAAATACCTATAGTGGCGGCACTTCCTGAGCCGGAGCGAAAAGTTTTATCGGAATGCGTGTTCTGGTGAACTTTTGGCTTACGGTTGTGATGTTGTGTTGTTGTGTTTGCAATTGGTCTGCg                                                               >  2:714549/1‑151 (MQ=255)
                                                                                      ataAATTGTCCGTTGAGCTTCTACCAGCAAATACCTATAGTGGCGGCACTTCCTGAGCCGGAGCGAAAAGTTTTATCCGAATGCGTGTTCTGGTGAACTTTTGGCTTACGGTTGTGATGTTGTGTTGTTGTGTTTGCAATTGGTCTGCGAt                                                             >  1:657709/1‑151 (MQ=255)
                                                                                          aTTGTCCGTTGAGCTTCTACCAGCAAATACCTATAGTGGCGGCACTTCCTGAGCCGGAGCGAAAAGTTTTATCGGAATGCGTGTTCTGGTGAACTTTTGGCTTACGGTTGTGATGTTGTGTTGTTGTGTTTGCAATTGGTCTGCGATTCAg                                                         <  1:351758/151‑1 (MQ=255)
                                                                                             gTCCGTTGAGCTTCTACCAGCAAATACCTATAGTGGCGGCACTTCCTGAGCCGGAGCGAAAAGTTTTATCGGAATGCGTGTTCTGGTGAACTTTTGGCTTACGGTTGTGATGTTGTGTTGTTGTGTTTGCAATTGGTCTGCGATTCAGAcc                                                      >  2:521313/1‑151 (MQ=255)
                                                                                             gTCCGTTGAGCTTCTACCAGCAAATACCTATAGTGGCGGCACTTCCTGAGCCGGAGCGAAAAGTTTTATCGGAATGCGTGTTCTGGTGAACTTTTGGCTTACGGTTGTGATGTTGTGTTGTTGTGTTTGCAATTGGTCTGCGATTCAGAc                                                       <  1:195779/150‑1 (MQ=255)
                                                                                                cGTTGAGCTTCTACCAGCAAATACCTATAGTGGCGGCACTTCCTGAGCCGGAGCGAAAAGTTTTATCGGAATGCGTGTTCTGGTGAACTTTTGGCTTACGGTTGTGATGTTGTGTTGTTGTGTTTGCAATTGGTCTGCGATTCAGACCATg                                                   >  2:473956/1‑151 (MQ=255)
                                                                                                cGTTGAGCTTCTACCAGCAAATACCTATAGTGGCGGCACTTCCTGAGCCGGAGCGAAAAGTTTTATCGGAATGCGTGTTCTGGTGAACTTTTGGCTTACGGTTGTGATGTTGTGTTGTTGTGTTTGCAATTGGTCTGCGATTCAGACCATg                                                   >  2:32690/1‑151 (MQ=255)
                                                                                                     aGCTTCTACCAGCAAATACCTATAGTGGCGGCACTTCCTGAGCCGGAGCGAAAAGTTTTATCGGAATGCGTGTTCTGGTGAACTTTTGGCTTACGGTTGTGATGTTGTGTTGTTGTGTTTGCAATTGGTCTGCGATTCAGACCATGGTAGc                                              >  1:27497/1‑151 (MQ=255)
                                                                                                      gCTTCTACCAGCAAATACCTATAGTGGCGGCACTTCCTGAGCCGGAGCGAAAAGTTTTATCGGAATGCGTGTTCTGGTGAACTTTTGGCTTACGGTTGTGATGTTGTGTTGTTGTGTTTGCAATTGGTCTGCGATTCAGACCATGGTAGCa                                             <  2:490097/151‑1 (MQ=255)
                                                                                                         tCTACCAGCAAATACCTATAGTGGCGGCACTTCCTGAGCCGGAGCGAAAAGTTTTATCGGAATGCGTGTTCTGGTGAACTTTTGGCTTACGGTTGTGATGTTGTGTTGTTGTGTTTGCAATTGGTCTGCGATTCAGACCATGGTAGCaa                                            <  2:95391/149‑1 (MQ=255)
                                                                                                          cTACCAGCAAATACCTATAAGGGCGGCACTTCCTCAGCCAGAGGTAAAAGTTTTATCGGAATGAGTGTTCTAGTGAACTTTTGACTTACTGTTGTGATGTTGTGTTGTTGGATGTGCAATTAGAATGCGATTCAGTCCAAGGTAGCaaa                                           <  1:554155/149‑1 (MQ=255)
                                                                                                               aGCAAATACCTATAGTGGCGGCACTTCCTGAGCCGGAGCGAAAAGTTTTATCGGAATGCGTGTTCTGGTGAACTTTTGGCTTACGGTTGTGATGTTGTGTTGTTGTGTTTGCAATTGGTCTGCGATTCAGACCATGGTAGCAAAGCTACCt                                    >  2:217857/1‑151 (MQ=255)
                                                                                                                gCAAATACCTATAGTGGCGGCACTTCCTGAGCCGGAGCGAAAAGTTTTATCGGAATGCGTGTTCTGGTGAACTTTTGGCTTACGGTTGTGATGTTGTGTTGTTGTGTTTGCAATTGGTCTGCGATTCAGACCATGGTAGCAAAGCTGCCtt                                   <  1:199695/151‑1 (MQ=255)
                                                                                                                      aCCTATAGTGGCGGCACTTCCTGAGCCGGAGCGAAAAGTTTTATCGGAATGCGTGTTCTGGTGAACTTTTGGCTTACGGTTGTGATGTTGTGTTGTTGTGTTTGCAATTGGTCTGCGATTCAGACCATGGTAGCAAAGCTACCTTTTTTCa                             <  1:156404/151‑1 (MQ=255)
                                                                                                                       ccTATAGTGGCGGCACTTCCTGAGCCGGAGCGAAAAGTTTTATCGGAATGCGTGTTCTGGTGAACTTTTGGCTTACGGTTGTGATGTTGTGTTGTTGTGTTTGCAATTGGTCTGCGATTCAGACCATGGTAGCAAAGCTACCTTTTTTc                              <  1:245256/149‑1 (MQ=255)
                                                                                                                       ccTATAGTGGCGGCACTTCCTGAGCCGGAGCGAAAAGTTTTATCGGAATGCGTGTTCTGGTGAACTTTTGGCTTACGGTTGTGATGTTGTGTTGTTGTGTTTGCAATTGGTCTGCGATTCAGACCATGGTAGCAAAGCTACCTTTTTTCa                             <  1:394057/150‑1 (MQ=255)
                                                                                                                       ccTATAGTGGCGGCACTTCCTGAGCCGGAGCGAAAAGTTTTATCGGAATGCGTGTTCTGGTGAACTTTTGGCTTACGGTTGTGATGTTGTGTTGTTGTGTTTGCAATTGGTCTGCGATTCAGACCATGGTAGCAAAGCTACCTTTTTTCAc                            >  1:718919/1‑151 (MQ=255)
                                                                                                                          ataGTGGCGGCACTTCCTGAGCCGGAGCGAAAAGTTTTATCGGAATGCGTGTTCTGGTGAACTTTTGGCTTACGGTTGTGATGTTGTGTTGTTGTGTTTGCAATTGGTCTGCGATTCAGACCATGGTAGCAAAGCTACCTTTTTTCACtt                          >  1:109366/1‑150 (MQ=255)
                                                                                                                          ataGTGGCGGCACTTCCTGAGCCGGAGCGAAAAGTTTTATCGGAATGCGTGTTCTGGTGAACTTTTGGCTTACGGTTGTGATGTTGTGTTGTTGTGTTTGCAATTGGTCTGCGATTCAGACCATGGTAGCAAAGCTACCTTTTTTCACtt                          <  1:340157/150‑1 (MQ=255)
                                                                                                                          ataGTGGCGGCACTTCCTGAGCCGGAGCGAAAAGTTTTATCGGAATGCGTGTTCTGGTGAACTTTTGGCTTACGGTTGTGATGTTGTGTTGTTGTGTTTGCAATTGGTCTGCGATTCAGACCATGGTAGCAAAGCTACCTTTTTTCACTTc                         >  1:508073/1‑151 (MQ=255)
                                                                                                                          ataGTGGCGGCACTTCCTGAGCCGGAGCGAAAAGTTTTATCGGAATGCGTGTTCTGGTGAACTTTTGGCTTACGGTTGTGATGTTGTGTTGTTGTGTTTGCAATTGGTCTGCGATTCAGACCATGGTAGCAAAGCTACCTTTTTTCACTTc                         >  1:622884/1‑151 (MQ=255)
                                                                                                                            aGTGGCGGCACTTCCTGAGCCGGAGCGAAAAGTTTTATCGGAATGCGTGTTCTGGTGAACTTTTGGCTTACGGTTGTGATGTTGTGTTGTTGTGTTTGCAATTGGTCTGCGATTCAGACCATGGTAGCAAAGCTACCTTTTTTCACTTCCt                       <  2:309017/151‑1 (MQ=255)
                                                                                                                               ggcggcACTTCCTGAGCCGGAGCGAAAAGTTTTATCGGAATGCGTGTTCTGGTGAACTTTTGGCTTACGGTTGTGATGTTGTGTTGTTGTGTTTGCAATTGGTCTGCGATTCAGACCATGGTAGCAAAGCTACCTTTTTTCACTTCCTGt                     <  1:798508/150‑1 (MQ=255)
                                                                                                                                 cggcACTTCCTGAGCCGGAGCGAAAAGTTTTATCGGAATGCGTGTTCTGGTGAACTTTTGGCTTACGGTTGTGATGTTGTGTTGTTGTGTTTGCAATTGGTCTGCGATTCAGACCATGGTAGCAAAGCTACCTTTTTTCACTTCCTGt                     <  2:661329/148‑1 (MQ=255)
                                                                                                                                    cACTTCCTGAGCCGGAGCGAAAAGTTTTATCGGAATGCGTGTTCTGGTGAACTTTTGGCTTACGGTTGTGATGTTGTGTTGTTGTGTTTGCAATTGGTCTGCGATTCAGACCATGGTAGCAAAGCTACCTTTTTTCACTTCCTGTACAttt               <  1:380585/151‑1 (MQ=255)
                                                                                                                                     aCTTCCTGAGCCGGAGCGAAAAGTTTTATCGGAATGCGTGTTCTGGTGAACTTTTGGCTTACGGTTGTGATGTTGTGTTGTTGTGTTTGCAATTGGTCTGCGATTCAGACCATGGTAGCAAAGCTACCTTTTTTCACTTCCTGTACATTTa              <  2:389975/151‑1 (MQ=255)
                                                                                                                                        tcatgggCCTGAGCGATAAGTTTTATCGGAATGCGTCTTCTGGTGAACTTTTGGCTTACGGTTGTGATGTGGTGTTGTTGTGTTTGCAAATGGTCTCCGATTCAGACCATGGTAGCAAAGCTACCTTTTTTCACTTCCCATCCATTTAc             <  1:194909/143‑1 (MQ=255)
                                                                                                                                        tCCTGAGCCGGAGCGAAAAGTTTTATCGGAATGCGTGTTCTGGTGAACTTTTGGCTTACGGTTGTGATGTTGTGTTGTTGTGTTTGCAATTGGTCTGCGATTCAGACCATGGTAGCAAAGCTACCTTTTTTCACTTCCTGTACATTTAcc            <  1:74375/150‑1 (MQ=255)
                                                                                                                                               ccGGAGCGAAAAGTTTTATCGGAATGCGTGTTCTGGTGAACTTTTGGCTTACGGTTGTGATGTTGTGTTGTTGTGTTTGCAATTGGTCTGCGATTCAGACCATGGTAGCAAAGCTACCTTTTTTCACTTCCTGTACATTTACCctgtctgt    <  2:231598/151‑1 (MQ=255)
                                                                                                                                                cGGAGCGAAAAGTTTTATCGGAATGCGTGTTCTGGTGAACTTTTGGCTTACGGTTGTGATGTTGTGTTGTTGTGTTTGCAATTGGTCTGCGATTCAGACCATGGTAGCAAAGCTACCTTTTTTCACTTCCTGTACATTTACCCtgtctgtc   <  2:738186/151‑1 (MQ=255)
                                                                                                                                                 ggAGCGAAAAGTTTTATCGGAATGCGTGTTCTGGTGAACTTTTGGCTTACGGTTGTGATGTTGTGTTGTTGTGTTTGCAATTGGTCTGCGATTCAGACCATGGTAGCAAAGCTACCTTTTTTCACTTCCTGTACATTTACCCTGTCTGTcc  >  1:632682/1‑151 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
GACATGGCTATTAAACTTTGATGTTAAATGAATTTAACAATTAGATCACACTATGTAACCTATTAGTTTTTTTAATCTGAGCCATTATAAATTGTCCGTTGAGCTTCTACCAGCAAATACCTATAGTGGCGGCACTTCCTGAGCCGGAACGAAAAGTTTTATCGGAATGCGTGTTCTGGTGAACTTTTGGCTTACGGTTGTGATGTTGTGTTGTTGTGTTTGCAATTGGTCTGCGATTCAGACCATGGTAGCAAAGCTACCTTTTTTCACTTCCTGTACATTTACCCTGTCTGTCC  >  NC_000913/2942563‑2942858

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