Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,540,645 G→C 6.5% A649G (GCC→GGC) ‡ narZ ← nitrate reductase 2 (NRZ), alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,540,6450GC 6.5% 231.8 / 6.5 78A649G (GCC→GGC) ‡narZnitrate reductase 2 (NRZ), alpha subunit
Reads supporting (aligned to +/- strand):  ref base G (33/39);  new base C (2/3);  total (35/43)
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.71e-01

GTTGTTCGCAGGCCATACGTAAATCGCCCGATTTCAGCGCCTGGGCGGTAAATTCTGTGGGGGATAATCCGGCCTTGTCGGCTTCAGCTTTAATCCCGAGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAGATGACCGCTGTACTTACTGGCATCGGCGAGCGGTGAAAGTAACTCCTGAGCAGAGACTTTTTCATAGCGCCATTGGCTGGAATGATTGTAGAAAAACGAGGTGCTGTTCATTTGGCGCGGTG  >  NC_000913/1540496‑1540789
                                                                                                                                                     |                                                                                                                                                
gttgttCGCAGGCCATACGTAAATCGCCCGATTTCAGCGCCTGGGCGGTAAATTCTGTGGGGGATAATCCGGCCTTGTCGGCTTCAGCTTTAATCCCGAGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCgg                                                                                                                                                  >  1:999613/1‑150 (MQ=255)
gttggtCGCAGGCCATACGTAACTCGCCCGATTTCAGCGCCTGGGCGGTAAATTCTGTGGGGGATAACCCGGCCTTGTCGGCTTCAGCTTTAATCCCGAGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCgg                                                                                                                                                  <  2:1590391/150‑1 (MQ=255)
       gCAGTCCATACGTAAATCTCACGATTTCAGCGCCTGGGCGGTAAATTCTGTGGGGGATAATCCGGCCTTGTCGGCTTCAGCTTTAATCCCGAGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGCGAAc                                                                                                                                           <  2:1041285/150‑1 (MQ=255)
           gCCTTACGTACATCGCCCGATTTCAGCGCCTGGGCGGTAAATTCTGTGGGGGATAATCCGGCCTTGTCGGCTTCAGCTTTAATCCCGAGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCAACGTTCGCGCCGAACGtatg                                                                                                                                      <  1:708311/150‑4 (MQ=255)
           gCCATACGTAAATCGCCCGATTTCAGCGCCTGGGCGGTAAATTCTGTGGGGGATAATCCGGCCTTGTCGGCTTCAGCTTTAATCCCGAGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTg                                                                                                                                       >  2:1468219/1‑150 (MQ=255)
            ccATACGTAAATCGCCCGATTTCAGCGCCTGGGCGGTAAATTCTGTGGGGGATAATCCGGCCTTGTCGGCTTCAGCTTTAATCCCGAGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGCGCCGAACGtatg                                                                                                                                      <  2:680124/150‑4 (MQ=255)
             cATACGTAAATCGCCCGATTTCAGCGCCTGGGCGGTAAATTCTGTGGGGGATAATCCGGCCTTGTCGGCTTCAGCTTTAATCCCGAGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGaa                                                                                                                                     >  1:264109/1‑150 (MQ=255)
              aTACGTAAATCGCCCGATTTCAGCGCCTGGGCGGTAAATTCTGTGGGGGATAATCCGGCCTTGTCGGCTTCAGCTTTAATCCCGAGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGaaa                                                                                                                                    >  1:1624949/1‑150 (MQ=255)
              aTACGTAAATCGCCCGATTTCAGCGCCTGGGCGGTAAATTCTGTGGGGGATAATCCGGCCTTGTCGGCTTCAGCTTTAATCCCGAGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGCGCCGAACGtatggg                                                                                                                                    <  1:184122/150‑6 (MQ=255)
                aCGTAAATCGCCCGATTTCAGCGCCTGGGCGGTAAATTCTGTGGGGGATAATCCGGCCTTGTCGGCTTCAGCTTTAATCCCGAGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAaatc                                                                                                                                  <  1:57698/150‑1 (MQ=255)
                aCGTAAATCGCCCGATTTCAGCGCCTGGGCGGTAAATTCTGTGGGGGATAATCCGGCCTTGTCGGCTTCAGCTTTAATCCCGAGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAaatc                                                                                                                                  <  1:363537/150‑1 (MQ=255)
                 cGTAAATCGCCCGATTTCAGCGCCTGGGCGGTAAATTCTGTGGGGGATAATCCGGCCTTGTCGGCTTCAGCTTTAATCCCGAGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAattca                                                                                                                                 >  1:289060/1‑150 (MQ=255)
                  gTAAATCGCCCGATTTCAGCGCCTGGGCGGTAAATTCTGTGGGGGATAATCCGGCCTTGTCGGCTTCAGCTTTAATCCCGAGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAaatcaa                                                                                                                                >  1:1012042/1‑150 (MQ=255)
                   tAAATCGCCCGATTTCAGCGCCTGGGCGGTAAATTCTGTGGGGGATAATCCGGCCTTGTCTGCTTCAGCTTTAATCCCGAGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAaatcaat                                                                                                                               <  2:1163725/150‑1 (MQ=255)
                   tAAATCGCCCGATTTCAGCGCCTGGGCGGTAAATTCTGTGGGGGATAATCCGGCCTTGTCGGCTTCAGCTTTAATCCCGAGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAaatcaat                                                                                                                               <  1:1348388/150‑1 (MQ=255)
                           ccGATTTCAGCGCCTGGGCGGTAAATTCTGTGGGGGATAATCCGGCCTTGTCGGCTTCAGCTTTAATCCCGAGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAGATGAc                                                                                                                       <  1:1526360/150‑1 (MQ=255)
                           ccGATTTCAGCGCCTGGGCGGTAAATTCTGTGGGGGATAATCCGGCCTTGTCGGCTTCAGCTTTAATCCCGAGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAGATGAc                                                                                                                       <  1:621507/150‑1 (MQ=255)
                             gcgtgctgcggCTGGGCGGTATATTTTGTGGGGGATACTCCGGCCTGGTCGGCTTCAGCTTTACTCGCGAGGGGGTTACGCCCCAGCTGCGGCGCAGAAGGTAGGCAGCGCATACGTTCGGCGCGATCGTTGAAATCATTCAGATGACCg                                                                                                                     <  1:111019/143‑1 (MQ=255)
                               tttCAGCGCCTGGGCGGTAAATTCTGTGGGGGATAATCCGGCCTTGTCGGCTTCAGCTTTAATCCCGAGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAGATGACCGCt                                                                                                                   >  2:1092078/1‑150 (MQ=255)
                                 tCAGCGCCTGGGCGGTAAATTGTGTGGGGGATAAACCGGCCTTGTCGGCTTCAGCTTTAATCCCGAGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAGATGACCGCTGt                                                                                                                 <  2:1340868/150‑1 (MQ=255)
                                 tCAGCGCCTGGGCGGTAAATTCTGTGGGGGATAATCCGGCCTTGTCGGCTTCAGCTTTAATCCCGAGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAGATGACCGCTGt                                                                                                                 >  1:388729/1‑150 (MQ=255)
                                  cAGCGCCTGGGCGGTAAATTCTGTGGGGGATAATCCGGCCTTGTCGGCTTCAGCTTTAATCCCGAGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAGATGACCGCTGta                                                                                                                >  2:244621/1‑150 (MQ=255)
                                  cAGCGCCTGGGCGGTAAATTCTGTGGGGGATAATCCGGCCTTGTCGGCTTCAGCTTTAATCCCGAGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAGATGACCGCTGta                                                                                                                >  1:1360775/1‑150 (MQ=255)
                                   aGCGCCTGGGCGGTAAATTCTGTGGGGGATAATCCGGCCTTGTCGGCTTCAGCTTTAATCCCGAGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAGATGACCGCTGtac                                                                                                               >  1:889222/1‑150 (MQ=255)
                                   aGCGCCTGGGCGGTAAATTCTGTGGGGGATAATCCGGCCTTGTCGGCTTCAGCTTTAATCCCGAGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAGATGACCGCTGtac                                                                                                               <  1:448998/150‑1 (MQ=255)
                                   aGCGCCAGGGCGGTAAATTCTGTGGGGGATAATCCGGCCTTGTCGGCTTCAGCTTTAATCCCGAGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGAGCGAACGTTGAAATCAATCAGATGACCGCTGtaa                                                                                                               >  2:1436125/1‑149 (MQ=255)
                                         tGGGCGGTAAATTCTGTGGGGGATAATCCGGCCTTGTCGGCTTCAGCTTTAATCCCGAGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAAATCAATCATATGACCGCTGTACTTACTg                                                                                                         >  2:1436327/1‑150 (MQ=255)
                                          gggCGGTAAATTCTGTGGGGGATAATCCGGCCTTGTCGGCTTCAGCTTTAATCCCGAGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAGATGACCGCTGTACTTACTgg                                                                                                        <  2:911014/150‑1 (MQ=255)
                                          gggCGGTAAATTCTGTGGGGGATAATCCGGCCTTGTCGGCTTCAGCTTTAATCCCGAGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAGATGACCGCTGTACTTACTgg                                                                                                        <  1:570313/150‑1 (MQ=255)
                                              ggTAAATTCTGTGGGGGATAATCCGGCCTTGTCGGCTTCAGCTTTAATCCCGAGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAGATGACCGCTGTACTTACTGGCATc                                                                                                    <  2:971845/150‑1 (MQ=255)
                                               gTAAATTCTGTGGGGGATAATCCGGCCTTGTCGGCTTCAGCTTTAATCCCGAGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAGATGACCGCTGTCCTTACTGGCATCg                                                                                                   >  1:620127/1‑150 (MQ=255)
                                               gTAAATTCTGTGGGGGATAATCCGGCCTTGTCGGCTTCAGCTTTAATCCCGAGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAGATGACCGCTGTACTTACTGGCATCg                                                                                                   >  2:449863/1‑150 (MQ=255)
                                                           ggggATAATCCGGCCTTGTCGGCTTCAGCTTTAATCCCGAGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAGATGACCGCTGTACTTACTGGCATCGGCGAGCGGTGaa                                                                                       <  1:244621/150‑1 (MQ=255)
                                                             ggATAATCCGGCCTTGTCTGCTTCAGCTTTAATCCCGAGCGTGTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAGATGACCGCTGTACTTACTGGGATCGGCGAGCGGTGAAAg                                                                                     <  2:782907/150‑1 (MQ=255)
                                                             ggATAATCCGGCCTTGTCGGCTTCAGCTTTAATCCCGAGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAGATGACCGCTGTACTTACTGGCATCGGCGAGCGGTGAAAg                                                                                     >  1:1279855/1‑150 (MQ=255)
                                                             ggATAATCCGGCCTTGTCGGCTTCAGCTTTAATCCCGAGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAGATGACCGCTGTACTTACTGGCATCGGCGAGCGGTGAAAg                                                                                     <  2:776295/150‑1 (MQ=255)
                                                             ggATAATCCGGCCTTGTCGGCTTCAGCTTTAATCCCGAGCGGGTTACGCCCCAGCTGAGGCGCAGAAGGTAGGCAGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAGATGACCGCTGTACTTACTGGCATCGGCGAGCGGTGAAAg                                                                                     <  2:1507521/150‑1 (MQ=255)
                                                               aTAATCCGGCCTTGTCGGCTTCAGCTTTAATCCCGAGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAGATGACCGCTGTACTTACTGGCATCGGCg                                                                                                >  1:295210/1‑137 (MQ=255)
                                                               aTAATCCGGCCTTGTCGGCTTCAGCTTTAATCCCGAGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAGATGACCGCTGTACTTACTGGCATCGGCg                                                                                                <  2:295210/137‑1 (MQ=255)
                                                                   tCCGGCCTTGTCGGCTTCAGCTTTAATCCCGAGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGAGAACGTTGAAATCAATCAGATGACCGCTGTACTTACTGGCATCGGCGAGCGGTGAAAGTAACTc                                                                               >  1:1163605/1‑150 (MQ=255)
                                                                         cTTGTAGGCTTCAGCTTTAATCACGAGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAGATGACCGCTGTACTTACTGGCATCGGCGAGCGGTGAAAGTAACTCCTGAGc                                                                         <  2:1095532/150‑1 (MQ=255)
                                                                                    cAGCTTTAATCCCGAGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAGATGACCGCTGTACTTACTGGCATCGGCGAGCGGTGAAAGTAACTCCTGAGCAGAGACttttt                                                              >  2:1615122/1‑150 (MQ=255)
                                                                                     aGCTTTAATCCCGAGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAGATGACCGCTGTACTTACTGGCATCGGCGAGCGGTGAAAGTAACTCCTGAGCAGAGACTTTTTc                                                             >  2:10952/1‑150 (MQ=255)
                                                                                       cTTTATTGCCTAGCGGGTTACGCCCCAGCTGCGGCGCAGATGGTCGCCAGCCCATACGTTCGTCGCGAACGTTGAATTCAATCAGATGACCGATTTAGTTACTTGCATCTGCGAGCGGTGAAAGTATCTCCTGTGCAGAGACTTTTTCAt                                                           <  1:1058435/150‑1 (MQ=255)
                                                                                        tttAATCCCGAGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAGATGACCGCTGTACTTACTGGCATCGGCGAGCGGTGAAAGTAACTCCTGAGCAGAGACTTTTTCATa                                                          >  1:1651363/1‑150 (MQ=255)
                                                                                           aaTCCCGAGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAGATGACCGCTGTACTTACTGGCATCGGCGAGCGGTGAAAGTAACTCCTGAGCAGAGACTTTTTCATAgcg                                                       >  2:367446/1‑150 (MQ=255)
                                                                                            aTCCCGAGCGGGTTATGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGCGATCGTTGATATCTATCAGATGACCGCTGTACTTACTGGCATCGGCGAGCGGTGCAAGTAACTCCTGCGCTGAGACTTTTTCATAgcgc                                                      >  2:845646/1‑150 (MQ=255)
                                                                                               ccGAGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCCGCCCATACGTTCGGCGCGAACGTTGATATCAATCAGATGACCGCTGTACTTACTGGCATCGGCGAGCGGTGAAAGTACCTCCTGAGCAGAGACTTTTGCATAGCGCCag                                                   >  1:1575663/1‑149 (MQ=255)
                                                                                                cGAGCGGGTTACGCCCCAGCTGCGGCGCAAAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAGATGACCGCGGTACTTACTGGCATCGGCGCGCGGTGAAAGTAACTCCTGAGCAGAGACTTTTTCATAGCGCCAtt                                                  >  1:1489531/1‑150 (MQ=255)
                                                                                                    cGGGTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAGATGACCGCTGTACTTACTGGCATCGGCGAGCGGTGAAAGTAACTCCTGAGCAGAGACTTTTTCATAGCGCCATTGGCt                                              <  1:1326001/150‑1 (MQ=255)
                                                                                                     ggggTACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAGATGACCGCTGTACTTACTGGCATCGGCGAGCGGTGAAAGTAACTCCTGAGCAGAGACTTTTTCATAGCGCCATTGGCTg                                             <  1:388227/150‑1 (MQ=255)
                                                                                                     gggTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAGATGACCGCTGTACTTACTGGCATCGGCGAGCGGTGAAAGTAACTCCTGAGCAGAGACTTTTTCATAGCGCCATTGGCTg                                             <  1:602273/150‑1 (MQ=255)
                                                                                                     gggTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAGATGACCGCTGTACTTACTGGCATCGGCGAGCGGTGAAAGTAACTCCTGAGCAGAGACTTTTTCATAGCGCCATTGGCTg                                             <  1:163198/150‑1 (MQ=255)
                                                                                                     cttctgcgccgcAGCTGGGGCGTAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAGATGACCGCTGTACTTACTGGCATCGGCGAGCGGTGAAAGTAACTCCTGAGCAGAGACTTTTTCATAGCGCCATTGGCTg                                             >  1:639816/7‑150 (MQ=255)
                                                                                                      ggTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGCGAGCGTTGAAATCAATCAGATGACCGCTGTCCTTACTGGCATCGGCGCGCGGTGAAAGTAACTCCTGAGCAGAGACTTTTTCATAGCGCCATTGGCTgg                                            >  1:40236/1‑150 (MQ=255)
                                                                                                      ggTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAGATGACCGCTGTACTTACTGGCATCGGCGAGCGGTGAAAGTAACTCCTGAGCAGAGACTTTTTCATAGCGCCATTGGCTgg                                            <  2:790663/150‑1 (MQ=255)
                                                                                                      ggTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCAATACGTTCGGCGCGAACGTTGAAATCAATCAGATGACCGCTGTACTTACTGGCATCGGCGAGCGGTGAAAGTAACTCCTGAGCAGAGACTTTTTCATAGCGCCATTGGCTgg                                            <  1:1685935/150‑1 (MQ=255)
                                                                                                          aCGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAGATGACCGCTGTACTTACTGGCATCGGCGAGCGGTGAAAGTAACTCCTGAGCAGAGACTTTTTCATAGCGCCATTGGCTGGAATg                                        <  2:520811/150‑1 (MQ=255)
                                                                                                           cGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAGATGACCGCTGTACTTACTGGCATCGGCGAGCGGTGAAAGTAACTCCTGAGCAGCGACTTTTTCATAGCGCCATTGGCTGGAATGa                                       >  2:1524142/1‑150 (MQ=255)
                                                                                                             ccccAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAGATGACCGCTGTACTTACTGGCATCGGCGAGCGGTGAAAGTAACTCCTGAGCAGAGACTTTTTCATAGCGCCATTGGCTGGAATGAtt                                     <  1:564496/150‑1 (MQ=255)
                                                                                                              cccAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAGATGACCGCTGTACTTACTGGCATCGGCGAGCGGTGAAAGTAACTCCTGAGCAGAGACTTTTTCATAGCGCCATTTGCTGGAATGATTg                                    >  2:1729984/1‑150 (MQ=255)
                                                                                                              cccAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAGATGACCGCTGTACTTACTGGCATCGGCGAGCGGTGAAAGTAACTCCTGAGCAGAGACTTTTTCATAGCGCCATTGGCTGGAATGATTg                                    <  1:521971/150‑1 (MQ=255)
                                                                                                               ccAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAGATGACCGCTGTACTTACTGGCATCGGCGAGCGGTGAAAGTAACTCCTGAGCAGAGACTTATTCATAGCGCCATTGGCTGGAATGATTGt                                   <  2:133212/150‑1 (MQ=255)
                                                                                                                cAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAGATGACCGCTGTACTTACTGGCATCGGCGAGCGGTGAAAGTAACTCCTGAGCAGAGACTTTTTCATAGCGCCATTGGCTGGAATGATTGTa                                  <  2:1640302/150‑1 (MQ=255)
                                                                                                                 aGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAGATGACCGCTGTACTTACTGGCATCGGCGAGCGGTGAAAGTAACTCCTGAGCAGAGACTTTTTCATAGCGCCATTGGCTGGAATGATTGTAg                                 <  1:449863/150‑1 (MQ=255)
                                                                                                                 aGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAGATGACCGCTGTACTTACTGGCATCGGCGAGCGGTGAAAGTAACTCCTGAGCAGAGACTTTTTCATAGCGCCATTGGCTGGAATGATTGTAg                                 <  1:592997/150‑1 (MQ=255)
                                                                                                                  gCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAGATGACCGCTGTACTTACTGGCATCGGCGAGCGGTGAAAGTAACTCCTGAGCAGAGACTTTTTCATAGCGCCATTGGCTGGAATGATTGTAGa                                >  1:843336/1‑150 (MQ=255)
                                                                                                                    tGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAGATGACCGCTGTACTTACTGGCATCGGCGAGCGGTGAAAGTAACTCCTGAGCAGAGACTTTTTCATAGCGCCATTGGCGGGAATGATTGTAGaaa                              >  1:139663/1‑150 (MQ=255)
                                                                                                                           cataaGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAGATGACCGCTGTACTTACTGGCATCGGCGAGCGGTGAAAGTAACTCCTGAGCAGAGACTTTTTCATAGCGCCATTGGCTGGAATGATTGTAGAAAAACGAgg                       <  1:1313158/147‑1 (MQ=255)
                                                                                                                           cagccggTAGCCAGACCCTCCGTTCGGCGCGAACGTTGAACTCCCGCAGCTGACCGCGGTACTTACTGGCATCGGCGAGCGGGGACAGTACCTCCTGCGCCGCGACTTTTTCATCGCGCCATTGGCGGGAAGGATTGTAGAAACACGAgg                       <  2:689830/145‑1 (MQ=255)
                                                                                                                              aaGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAGATGACCGCTGTACTTACTGGCATCGGCGAGCGGTGAAAGTAACTCCTGAGCAGAGACTTTTTCATAGCGCCATTGGCTGGAATGATTGTAGAAAAACGAGGTGc                    <  1:429576/150‑1 (MQ=255)
                                                                                                                               aGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAGATGACCGCTGTACTTACTGGCATCGGCGAGCGGTGAAAGTAACTCCTGAGCAGAGACTTTTTCATAGCGCCATTGGCTGGAATGATTGTAGAAAAACGAGGTGCt                   >  1:1624449/1‑150 (MQ=255)
                                                                                                                                 gcagccagccCATACGTTCGGCGCGAACGTTGAAATCAATCAGATGACCGCTGTACTTACTGGCAGCGGCGAGCGGTGAAAGTAACTCCTGAGCCGAGACTTTTTCATAGCGCCATTGGCTGGAATGATTGTAGAAAAACGAGGTGCTGt                 <  2:1039566/148‑1 (MQ=255)
                                                                                                                                  tAGCCAGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAGATGACCGCTGTACTTACTGGCATCGGCGAGCGGTGAAAGTAACTCCTGAGCAGAGACTTTTTCATAGCGCCATTGGCTGGAATGATTGTAGAAAAACGAGGTGCTGtt                <  1:430246/150‑1 (MQ=255)
                                                                                                                                  tAGCCAGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAGATGACCGCTGTACTTACTGGCATCGGCGAGCGGTGAAAGTAACTCCTGAGCAGAGACTTTTTCATAGCGCCATTGGCTGGAATGATTGTAGAAAAACGAGGTGCTGtt                >  2:1361735/1‑150 (MQ=255)
                                                                                                                                      cagccCATACGTTCGGCGCGAACGTGTAAATCAATAAGATGACCGCTGTACTTACTGGCATCTGCGAGCGGTGAAAGTAACTCCTGAGCAGAGACTTTTTCATAGCGCCATTGGCTGGAATGATTGTAGAAAAACGAGGTGCTGTTCAtt            <  2:347590/150‑1 (MQ=255)
                                                                                                                                       agccCATACGTTCGGCGCGAACGTTGAAATCAATCAGATGACCGCTGTACTTACTGGCATCGGCGAGCGGTGAAAGTAACTCCTGAGCAGAGACTTTTTCATAGCGCCATTGGCTGGAATGATTGTAGAAAAACGAGGTGCTGTTCAttt           >  1:195511/1‑150 (MQ=255)
                                                                                                                                         tttcaaCGTTCGCGCCGAACGTTGAAATCAATCAGATGACCGCTGTACTTACTGGCATCGGCGAGCGGTGAAAGTAACTCCTGAGCAGAGACTTTTTCATAGCGCCATTGGCTGGAATGATTGTAGAAAAACGAGGTGCTGTTCATTTgg         >  2:1724685/6‑150 (MQ=255)
                                                                                                                                          ttcaaCGTTCGCGCCGAACGTTGAAATCAATCAGATGACCGCTGTACTTACTGGCATCGGCGAGCGGTGAAAGTAACTCCTGAGCAGAGACTTTTTCATAGCGCCATTGGCTGGAATGATTGTAGAAAAACGAGGTGCTGTTCATTTGgc        >  2:360612/5‑150 (MQ=255)
                                                                                                                                                gtgcGGCGCGTACGTTGAAATCAATCAGATGACCGCTGTACTTACTGGCATCGGCGAGCGGTGTAAGTAACTCCTGAGCAGAGAGTTTTTCATAGCGCCATTGGCTGGAATGATTGTAGAAAAACGAGGTGCTGTTCATTTGGCGCggtg  <  2:293184/147‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                
GTTGTTCGCAGGCCATACGTAAATCGCCCGATTTCAGCGCCTGGGCGGTAAATTCTGTGGGGGATAATCCGGCCTTGTCGGCTTCAGCTTTAATCCCGAGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAGATGACCGCTGTACTTACTGGCATCGGCGAGCGGTGAAAGTAACTCCTGAGCAGAGACTTTTTCATAGCGCCATTGGCTGGAATGATTGTAGAAAAACGAGGTGCTGTTCATTTGGCGCGGTG  >  NC_000913/1540496‑1540789

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