Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,762,985 A→T 5.1% N62I (AAC→ATC)  PP_3326 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,762,9850AT 5.1% 239.5 / 3.8 79N62I (AAC→ATC) PP_3326hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (40/34);  new base T (2/2);  total (42/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.06e-01

GATGGTCAGCGTCAACGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGC‑GGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATGGGTCAATCAGGATCGCACCCA  >  NC_002947/3762839‑3763129
                                                                                                                                                  |                                                                                                                                                 
gATGGTCAGCGTCAACGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACAt                                                                                                                                                <  2:1066306/150‑1 (MQ=255)
  tGGTCAGCGTCAACGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGa                                                                                                                                              <  1:1898035/150‑1 (MQ=255)
    ggcAGCGTCAACGCAAGCGCCCACGGCCTGTGGACCGAGGCGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCGGAGGATACCGCCCAGAAGGCGGAAAAGATCAAGCGCGCCGGGGAGTTCGACGTGACCGGGAACATGTTc                                                                                                                                            <  2:113821/148‑1 (MQ=255)
      cAGCGTCAACGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATccc                                                                                                                                          >  1:1016965/1‑150 (MQ=255)
       aGCGTCAACGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCg                                                                                                                                         >  2:733853/1‑150 (MQ=255)
          gTCAACGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGt                                                                                                                                       >  1:1912254/1‑149 (MQ=255)
          gTCAACGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACAGGATCCCGGt                                                                                                                                       <  2:1912254/149‑1 (MQ=255)
           tCAACGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCa                                                                                                                                     >  1:743616/1‑150 (MQ=255)
              aCGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGACGATTACGCCTACAAGGCGGAAAAGATCGAGCGCGCCTGGGCGCTCGACGTGACCGGGAACATGATCCCGGTCACgg                                                                                                                                  >  1:1867715/1‑150 (MQ=255)
               cGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCAGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGt                                                                                                                                 <  2:1098446/150‑1 (MQ=255)
               cGCAGGCGCCCACGGCCTGTGCACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGt                                                                                                                                 <  1:138889/150‑1 (MQ=255)
               cGCAGGCGCCCACGGCCTGCGGACCGAGGAGCGACGTGGCAATATCGAGGCAGTGTGCGACCACGGGCGCTATGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGt                                                                                                                                 <  1:224894/150‑1 (MQ=255)
                   ggCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAg                                                                                                                             <  1:459950/150‑1 (MQ=255)
                        ccACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCt                                                                                                                        >  1:1788015/1‑150 (MQ=255)
                          aCGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTgg                                                                                                                      >  2:894977/1‑150 (MQ=255)
                            ggCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGGGCCAGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGC‑g                                                                                                                   >  2:1159436/1‑150 (MQ=255)
                                tgtgGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGC‑GGATc                                                                                                               >  1:1551042/1‑150 (MQ=255)
                                tgtgGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGC‑GGATc                                                                                                               >  2:721026/1‑150 (MQ=255)
                                    gACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGC‑GGATCATGc                                                                                                           >  2:1827976/1‑150 (MQ=255)
                                        gaggagCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGC‑GGATCATGCCGtg                                                                                                       >  2:568086/1‑150 (MQ=255)
                                           gagCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGATCATGTTCCCGGTc                                                                                                                                      <  1:458002/117‑1 (MQ=255)
                                           gagCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGATCATGTTCCCGGTc                                                                                                                                      >  2:458002/1‑117 (MQ=255)
                                             gCGACGTGGCAATATCGAGGCCGTGTACGTCCACGGCGCCGCGGGAAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTCCGACGTGACCCGGAACATGATCTCGGTCACGGACTCTCGCCTGGC‑GGCTCATGCCGGGCTGCa                                                                                                  >  2:403211/1‑150 (MQ=255)
                                             gCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGC‑GGATCATGCCGTGCTGCa                                                                                                  >  2:84226/1‑150 (MQ=255)
                                             gCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGC‑GGATCATGCCGTGCTGCa                                                                                                  >  1:372907/1‑150 (MQ=255)
                                                aCGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGC‑GGATCATGCCGTGCTGCAAcc                                                                                               <  1:1807226/150‑1 (MQ=255)
                                                 cGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGC‑GGATCATGCCGTGCTGCAACCa                                                                                              <  1:1071589/150‑1 (MQ=255)
                                                     gCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGCAAACCGAATGCCGGTCCTGGTCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGC‑GGATCATGCCGTGCTGCAACCACTTa                                                                                          <  1:431508/150‑1 (MQ=255)
                                                       aaTATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGTC‑GGATCATGCCGTGGTGCAACCACTTAcc                                                                                        >  1:1219960/1‑150 (MQ=255)
                                                           tCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGt                                                                                                                                       <  1:455927/100‑1 (MQ=255)
                                                           tCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGt                                                                                                                                       >  2:455927/1‑100 (MQ=255)
                                                              aGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGC‑GGATCATGCCGTGCTGCAACCACTTACCCCGCCTg                                                                                 >  1:947821/1‑150 (MQ=255)
                                                              aGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGCTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGC‑GGATCATGCCGTGCTGCAACCACTTACCCCGCCTg                                                                                 >  2:1408472/1‑150 (MQ=255)
                                                                 ccGTGTACGGCCACGGCGCCGAGGCTAACGCCTACACGGCGGAAAAGATCAAACGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGC‑GGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCgg                                                                              >  1:1184737/1‑150 (MQ=255)
                                                                         ggCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGC‑GGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGcc                                                                      <  2:1867715/150‑1 (MQ=255)
                                                                           ccACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGC‑GGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGt                                                                    <  2:1016965/150‑1 (MQ=255)
                                                                                        gATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGC‑GGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAAc                                                       >  1:1747646/1‑150 (MQ=255)
                                                                                        gATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGC‑GGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAAc                                                       >  2:1104121/1‑150 (MQ=255)
                                                                                          tAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGC‑GGATCATGCCGTGCTGCGACCACTTACCCCGCCTGCGGCa                                                                            <  2:1877501/127‑1 (MQ=255)
                                                                                          tAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGC‑GGATCATGCCGTGCTGCGACCACTTACCCCGCCTGCGGCa                                                                            >  1:1877501/1‑127 (MQ=255)
                                                                                           aaCGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGC‑GGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACggg                                                    <  1:721026/150‑1 (MQ=255)
                                                                                            aCGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGC‑GGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGc                                                   >  1:608109/1‑150 (MQ=255)
                                                                                             cGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGC‑GGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGcc                                                  >  1:75231/1‑150 (MQ=255)
                                                                                              gCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGC‑GGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCt                                                 >  2:1254666/1‑150 (MQ=255)
                                                                                               ccTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGC‑GGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTgctgct                                                              >  2:265011/1‑136 (MQ=255)
                                                                                               ccTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGC‑GGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTgctgct                                                              <  1:265011/136‑1 (MQ=255)
                                                                                                cTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGC‑GGATCATGCCGTGCTGCAACCAc                                                                                             >  1:847753/1‑104 (MQ=255)
                                                                                                cTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGC‑GGATCATGCCGTGCTGCAACCAc                                                                                             <  2:847753/104‑1 (MQ=255)
                                                                                                cTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGC‑GGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTg                                               <  2:840958/150‑1 (MQ=255)
                                                                                                   cAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCGAGC‑GGCTCATGCCGTGCTGCAAGCATTGTCGCTGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGTTc                                            >  1:1377546/1‑150 (MQ=255)
                                                                                                    aaGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGC‑GGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTca                                           >  2:603076/1‑150 (MQ=255)
                                                                                                      ggCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGC‑GGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTcaca                                         <  2:372907/150‑1 (MQ=255)
                                                                                                        cGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGC‑GGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAgc                                       <  2:1702225/150‑1 (MQ=255)
                                                                                                            aaaGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGC‑GGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGccc                                   >  2:1393589/1‑150 (MQ=255)
                                                                                                                aTCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAgc                                                                                                                            <  1:417546/58‑1 (MQ=255)
                                                                                                                aTCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAgc                                                                                                                            >  2:417546/1‑58 (MQ=255)
                                                                                                                aTCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTgg                                                                                                                      <  2:1804126/64‑1 (MQ=255)
                                                                                                                aTCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTgg                                                                                                                      >  1:1804126/1‑64 (MQ=255)
                                                                                                                aTCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGC‑GGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGAc                               >  1:1784843/1‑150 (MQ=255)
                                                                                                                  cAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGC‑GGATCATGCCGTGCTGCAACCACTTACCCAGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACaa                             <  2:264280/150‑1 (MQ=255)
                                                                                                                        cgcCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGC‑GGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATg                       <  2:1551042/150‑1 (MQ=255)
                                                                                                                        cgcCTGGGCGTTCGACGTCACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGC‑GGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATg                       <  2:576491/150‑1 (MQ=255)
                                                                                                                         gccgcggCGTTCGACGTGACCGGGACCATGATCCCGGTCACGGTCGAGCGCCTGGC‑GGATCATGCCGTGCTGCAACTACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGACACAGCGCCCGGACAAGCAATgg                      <  2:947821/145‑1 (MQ=255)
                                                                                                                         gcCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGC‑gg                                                                                                                  >  2:826975/1‑58 (MQ=255)
                                                                                                                         gcCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGC‑gg                                                                                                                  <  1:826975/58‑1 (MQ=255)
                                                                                                                         gcCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGC‑GGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATgg                      >  1:732149/1‑150 (MQ=255)
                                                                                                                           cTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGC‑GGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATGGGt                    <  1:1393589/150‑1 (MQ=255)
                                                                                                                           cTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGC‑GGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATGGGt                    <  2:732149/150‑1 (MQ=255)
                                                                                                                           cTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAG‑AACTGTCTGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATGGGt                    <  1:593791/150‑1 (MQ=255)
                                                                                                                            tGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGC‑GGATCATGCCGTGCTGCAACCACTTATCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATGGGTc                   >  2:1183976/1‑150 (MQ=255)
                                                                                                                             aggCGCTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCt                                                                                                                        >  1:1380978/2‑49 (MQ=255)
                                                                                                                             aggCGCTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCt                                                                                                                        <  2:1380978/48‑1 (MQ=255)
                                                                                                                                cGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGC‑GGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATGGGTCAATc               >  2:1862255/1‑150 (MQ=255)
                                                                                                                                  ttCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGC‑GGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATGGGTCAATCAg             <  1:1506401/150‑1 (MQ=255)
                                                                                                                                       cGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGC‑GGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATGGGTCAATCAGGATCg        <  1:1159436/150‑1 (MQ=255)
                                                                                                                                          gACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGC‑GGATCATGCCGTGCTGCAACCCCTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATGGGTCAATCAGGATCGCAc     <  1:259543/150‑1 (MQ=255)
                                                                                                                                           aCCGGGATCATGTTCCCGGTCAC‑GTCGAACGCCCGGC‑GGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGc                                                                       <  2:1141565/82‑1 (MQ=21)
                                                                                                                                           aCCGGGATCATGTTCCCGGTCAC‑GTCGAACGCCCGGC‑GGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGc                                                                       >  1:1141565/1‑82 (MQ=21)
                                                                                                                                             cGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGC‑GGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATGGGTCAATCAGGATCGCACCCa  >  2:1314780/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
GATGGTCAGCGTCAACGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGC‑GGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATGGGTCAATCAGGATCGCACCCA  >  NC_002947/3762839‑3763129

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