Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,540,647 G→C 6.1% R648R (CGC→CGG 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,6470GC 6.1% 245.6 / 4.3 78R648R (CGC→CGGnarZnitrate reductase 2 (NRZ), alpha subunit
Reads supporting (aligned to +/- strand):  ref base G (33/40);  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 = 1.00e+00

GCAGGCCATACGTAAATCGCCCGATTTCAGCGCCTGGGCGGTAAATTCTGTGGGGGATAATCCGGCCTTGTCGGCTTCAGCTTTAATCCCGAGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAGATGACCGCTGTACTTACTGGCATCGGCGAGCGGTGAAAGTAACTCCTGAGCAGAGACTTTTTCATAGCGCCATTGGCTGGAATGATTGTAGAAAAACGAGGTGCTGTTCATTTGGCGCGGTG  >  NC_000913/1540503‑1540789
                                                                                                                                                |                                                                                                                                              
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:363537/150‑1 (MQ=255)
         aCGTAAATCGCCCGATTTCAGCGCCTGGGCGGTAAATTCTGTGGGGGATAATCCGGCCTTGTCGGCTTCAGCTTTAATCCCGAGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAaatc                                                                                                                                  <  1:57698/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                                                                                                                >  1:1360775/1‑150 (MQ=255)
                           cAGCGCCTGGGCGGTAAATTCTGTGGGGGATAATCCGGCCTTGTCGGCTTCAGCTTTAATCCCGAGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAGATGACCGCTGta                                                                                                                >  2:244621/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                                                                                     <  2:776295/150‑1 (MQ=255)
                                                      ggATAATCCGGCCTTGTCGGCTTCAGCTTTAATCCCGAGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAGATGACCGCTGTACTTACTGGCATCGGCGAGCGGTGAAAg                                                                                     >  1:1279855/1‑150 (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:163198/150‑1 (MQ=255)
                                                                                              gggTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAGATGACCGCTGTACTTACTGGCATCGGCGAGCGGTGAAAGTAACTCCTGAGCAGAGACTTTTTCATAGCGCCATTGGCTg                                             <  1:602273/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                >  2:1361735/1‑150 (MQ=255)
                                                                                                                           tAGCCAGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAGATGACCGCTGTACTTACTGGCATCGGCGAGCGGTGAAAGTAACTCCTGAGCAGAGACTTTTTCATAGCGCCATTGGCTGGAATGATTGTAGAAAAACGAGGTGCTGtt                <  1:430246/150‑1 (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)
                                                                                                                                                |                                                                                                                                              
GCAGGCCATACGTAAATCGCCCGATTTCAGCGCCTGGGCGGTAAATTCTGTGGGGGATAATCCGGCCTTGTCGGCTTCAGCTTTAATCCCGAGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAGATGACCGCTGTACTTACTGGCATCGGCGAGCGGTGAAAGTAACTCCTGAGCAGAGACTTTTTCATAGCGCCATTGGCTGGAATGATTGTAGAAAAACGAGGTGCTGTTCATTTGGCGCGGTG  >  NC_000913/1540503‑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: