Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 321,440 A→G 6.6% T278A (ACT→GCT)  aguA → agmatine deiminase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947321,4400AG 6.6% 264.9 / 5.4 90T278A (ACT→GCT) aguAagmatine deiminase
Reads supporting (aligned to +/- strand):  ref base A (42/42);  new base G (3/3);  total (45/45)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.25e-01

AAGTGTTACTGGCCTGGACCGATGATTCCAACGACCCCAACTATGCGCGCTGCCACGCGGCCATGGATGTGCTGAAGAACACCCGCGACGCCAAGGGCCGCGAGTTCATCGTGCACAAGATGCCGATCCCGGGCCCGCTGTTTGCCACTGCCGAAGAGTGCGCCGGTGTCGACCAGGTCGCCGGCAGCCAGGAGCGTGACCCGTCGGTGCGCCTGGCGGGTTCGTACGTGAACTTCCTGATTGTTAACGGTGGCATCATTGCCCCAAGCTTCGACGACCCGGCAGATGCCGAA  >  NC_002947/321294‑321586
                                                                                                                                                  |                                                                                                                                                  
aaGTGTTACTGGCCTGGACCGATGATTCCAACGACCCCAACTATGCGCGCTGCCACGCGGCCATGGATGTGCTGAAGAACACCCGCGACGCCAAGGGCCGCGAGTTCATCGTGCACAAGATGCCGATCCCGGGCCCGCTGTTTGCCACTg                                                                                                                                                 >  2:32087/1‑150 (MQ=255)
  gtgtTACTGGCCTGGACCGACGATTCCAACGACCCCAACTATGCGCGCTGCCACGCGGCCATGGATGTGCTGAAGAACACCCGCGACGCCAAGGGCCGCGAGTTCATCGTGCACAAGATGCCGATCCCGGGCCCGCTGTTTGCCACTGcc                                                                                                                                               >  2:55196/1‑150 (MQ=255)
   tgtTACTGGCCTGGACCGATGATTCCAACGACCCCAACTATGCGCGCTGCCACGCGGCCATGGATGTGCTGAAGAACACCCGCGACGCCAAGGGCCGCGAGTTCATCGTGCACAAGATGCCGATCCCGGGCCCGCTGTTTGCCACTGCCg                                                                                                                                              >  1:486737/1‑150 (MQ=255)
   tgtTACTGGCCTGGACCGATGATTCCAACGACCCCAACTATGCGCGCTGCCACGCGGCCATGGATGTGCTGAAGAACACCCGCGACGCCAAGGGCCGCGAGTTCATCGTGCACAAGATGCCGATCCCGGGCCCGCTGTTTGCCACTGCCg                                                                                                                                              >  1:1444990/1‑150 (MQ=255)
     ttACTGGCCTGGACCGATGATTCCAACGACCCCAACTATGCGCGCTGCCACGCGGCCATGGATGTGCTGAAGAACACCCGCGACGCCAAGGGCCGCGAGTTCATCGTGCACAAGATGCCGATCCCGGGCCCGCTGTTTGCCACTGCCGaa                                                                                                                                            >  1:1260433/1‑150 (MQ=255)
       aCTGGCCTGGACCGATGATTCCAACGACCCCAACTATGCGCGCTGCCACGCGGCCATGGATGTGCTGAAGAACACCCGCGACGCCAAGGGCCGCGAGTTCATCGTGCACAAGATGCCGATCCCGGGCCCGCTGTTTGCCACTGCCGAaga                                                                                                                                          <  2:1410467/150‑1 (MQ=255)
       aCTGGCCTGGACCGATGATTCCAACGACCCCAACTATGCGCGCTGCCACGCGGCCATGGATGTGCTGAAGAACACCCGCGACGCCAAGGGCCGCGAGTTCATCGTGCACAAGATGCCGATCCCGGGCCCGCTGTTTGCCACTGCCGAaga                                                                                                                                          <  2:638700/150‑1 (MQ=255)
               ggACCGATGATTCCAACGACCCCAACTATGCGCGCTGCCACGCGGCCATGGATGTGCTGAAGAACACCCGCGACGCCAAGGGCCGCGAGTTCATCGTGCACAAGATGCCGATCCCGGGCCCGCTGTTTGCCACTGCCGAAGAGTGCGCCg                                                                                                                                  >  2:1707390/1‑150 (MQ=255)
                 cccGATGATTCCAACGACCCCAACTATGCGCGCTGCCACGCGGCCATGGATGTGCTGAAGAACACCCGCGACGCCAAGGGCCGCGAGTTCATCGTGCACAAGATGCCGATCCCGGGCCCGCTGTTTGCCACTGCCGAAGAGTGCGCCGgt                                                                                                                                <  2:1328541/149‑1 (MQ=255)
                 aCCGATGATTCCAACGACCCCAACTATGCGCGCTGCCACGCGGCCATGGATGTGCTGAAGAACACCCGCGCCACCTCTTTCCGCGAGTTTGCGGGTCACAAGATGCCGATCCCGGGCCCGCTGTTTGCCACTGCCGAAGAGTGCGCCGgt                                                                                                                                <  1:1178211/150‑1 (MQ=255)
                 aCCGATGATTCCAACGACCCCAACTATGCGCGCTGCCACGCGGCCATGGATGTGCTGAAGAACACCCGCGACGCCAAGGGCCGCGAGTTCATCGTGCACAAGATGCCGATCCCGGGCCCGCTGTTTGCCACTGCCGAAGAGTGCGCCGgt                                                                                                                                <  1:1730809/150‑1 (MQ=255)
                    gatgatTCCAACGACCCCAACTATGCGCGCTGCCACGCGGCCATGGATGTGCTGAAGAACACCCGCGACGCCAAGGGCCGCGAGTTCATCGTGCACAAGATGCCGATCCCGGGCCCGCTGTTTGCCACt                                                                                                                                                  >  2:232448/1‑129 (MQ=255)
                    gatgatTCCAACGACCCCAACTATGCGCGCTGCCACGCGGCCATGGATGTGCTGAAGAACACCCGCGACGCCAAGGGCCGCGAGTTCATCGTGCACAAGATGCCGATCCCGGGCCCGCTGTTTGCCACt                                                                                                                                                  <  1:232448/129‑1 (MQ=255)
                         ttCCAACGACCCCAACTATGCGCGCTGCCACGCGGCCATGGATGTGCTGAAGAACACCCGCGACGCCAAGGGCCGCGAGTTCATCGTGCACAAGATGCCGATCCCGGGCCCGCTGTTTGCCACTGCCGAAGAGTGCGCCGGTGTCGACCa                                                                                                                        <  1:671908/150‑1 (MQ=255)
                         ttCCAACGACCCCAACTATGCGCGCTGCCACGCGGCCATGGATGTGCTGAAGAACACCCGCGACGCCAAGGGCCGCGAGTTCATCGTGCACAAGATGCCGATCCCGGGCCCGCTGTTTGCCACTGCCGAAGAGTGCGCCGGTGTCGACCa                                                                                                                        >  1:419347/1‑150 (MQ=255)
                          tCCAACGACCCCAACTATGCGCGCTGCCACGCGGCCCTGGCTGTGCCTAAGTACACCCGCGACGCCAAGGGCCGCGAGTTCATCGTGCACAAGATGCCGATCCCGGGCCCGCTGTTTGCCACTGCCGAAGAGTGCGCCGGTGTCGACCAg                                                                                                                       <  1:1564442/150‑1 (MQ=255)
                            cAACGACCCCAACTATGCGCGCTGCCACGCGGCCATGGATGTGCTGAAGAACACCCGCGACGCCAAGGGCCGCGAGTTCATCGTGCACAAGATGCCGATCCCGGGCCCGCTGTTTGCCACTGCCGAAGAGTGCGCCGGTGTCGACCAgg                                                                                                                      <  2:462392/149‑1 (MQ=255)
                            cAACGACCCCAACTATGCGCGCTGCCACGCGGCCATGGATGTGCTGAAGAACACCCGCGACGCCAAGGGCCGCGAGTTCATCGTGCACAAGATGCCGATCCCGGGCCCGCTGTTTGCCACTGCCGAAGAGTGCGCCGGTGTCGACCAgg                                                                                                                      >  1:462392/1‑149 (MQ=255)
                            cAACGACCCCAACTATGCGCGCTGCCACGCGGCCATGGATGTGCTGAAGAACACCCGCGACGCCAAGGGCCGCGAGTTCATCGTGCACAAGATGCCGATCCCGGGCCCGCTGTTTGCCACTGCCGAAGAGTGCGCCGGTGTCGACCAGGt                                                                                                                     >  1:898640/1‑150 (MQ=255)
                             aaCGACCCCAACTATGCGCGCTGCCACGCGGCCATGGATGTGCTGAAGAACACCCGCGACGCCAAGGGCCGCGAGTTCATCGTGCACAAGATGCCGATCCCGGGCCCGCTGTTTGCCACTGCCGAAGAGt                                                                                                                                        >  1:288558/1‑130 (MQ=255)
                             aaCGACCCCAACTATGCGCGCTGCCACGCGGCCATGGATGTGCTGAAGAACACCCGCGACGCCAAGGGCCGCGAGTTCATCGTGCACAAGATGCCGATCCCGGGCCCGCTGTTTGCCACTGCCGAAGAGt                                                                                                                                        <  2:288558/130‑1 (MQ=255)
                             aaCGACCCCAACTATGCGCGCTGCCACGCGGCCATGGATGTGCTGAAGAACACCCGCGACGCCAAGGGCCGCGAGTTCATCGTGCACAAGATGCCGATCCCGGGCCCGCTGTTTGCCACTGCCGAAGAGTGCGCCGGTGTCGACCAGGTc                                                                                                                    <  1:1363546/150‑1 (MQ=255)
                             aaCGACCCCAACTATGCGCGCTGCCACGCGGCCATGGATGTGCTGAAGAACACCCGCGACGCCAAGGGCCGCGAGTTCATCGTGCACAAGATGCCGATCCCGGGCCCGCTGTTTGCCACTGCCGAAGAGTGCGCCGGGGTCGACCAGGTc                                                                                                                    >  1:391966/1‑150 (MQ=255)
                             aaCGACCCCAACTATGCGCGCTGCCACGCGGCCATGGATGTGCTGAAGAACACCCGCGACGCCAAGGGCCGCGAGTTCATCGTGCACAAGATGCCGATCCCGGGCCCGCTGTTTGCCACTACCGAAGAGTGCGCCGGTGTCGACCAGGTc                                                                                                                    >  1:1838598/1‑150 (MQ=255)
                                gACCCCAACTATGCGCGCTGCCACGCGGCCATGGATGTGCTGAAGAACACCCGCGACGCCAAGGGCCGCGAGTTCATCGTGCACAAGATGCCGATCCCGGGCCCGCTGTTTGCCACTGCCGAAGAGTGCGCCGGTGTCGACCAGGTCGcc                                                                                                                 >  1:167441/1‑150 (MQ=255)
                                gACCCCAACTATGCGCGCTGCCACGCGGCCATGGATGTGCTGAAGAACACCCGCGACGCCAAGGGCCGCGAGTTCATCGTGCACAAGATGCCGATCCCGGGCCCGCTGTTTGCCACTGCCGAAGAGTGCGCCGGTGTCGACCAGGTCGcc                                                                                                                 >  1:1423183/1‑150 (MQ=255)
                                    ccAACTATGCGCGCTGCCACGCGGCCATGGATGTGCTGAAGAACACCCGCGACGCCAAGGGCCGCGAGTTCATCGTGCACAAGATGCCGATCCCGGGCCCGCTGTTTGCCACTGCCGAaga                                                                                                                                          <  2:1439982/121‑1 (MQ=255)
                                    ccAACTATGCGCGCTGCCACGCGGCCATGGATGTGCTGAAGAACACCCGCGACGCCAAGGGCCGCGAGTTCATCGTGCACAAGATGCCGATCCCGGGCCCGCTGTTTGCCACTGCCGAaga                                                                                                                                          >  1:1439982/1‑121 (MQ=255)
                                     cAACTATGCGCGCTGCCACGCGGCCATGGATGTGCTGAAGAACACCCGCGACGCCAAGGGCCGCGAGTTCATCGTGCACAAGATGCCGATCCCGGGCCCGCTGTTTGCCACTGCCGAAGAGTGCGCCGGTGTCGAc                                                                                                                          >  2:1826500/1‑136 (MQ=255)
                                     cAACTATGCGCGCTGCCACGCGGCCATGGATGTGCTGAAGAACACCCGCGACGCCAAGGGCCGCGAGTTCATCGTGCACAAGATGCCGATCCCGGGCCCGCTGTTTGCCACTGCCGAAGAGTGCGCCGGTGTCGAc                                                                                                                          <  1:1826500/136‑1 (MQ=255)
                                              gcgcTGCCACGCGGCCATGGATGTGCTGAAGAACACCCGCGACGCCAAGGGCCGCGAGTTCATCGTGCACAAGATGCCGATCCCGGGCCCGCTGTTTGCCACTGCCGAAGAGTGCGCCGGTGTCGACCAGGTCGCCGGCAGCCAGGAGCg                                                                                                   >  2:1912575/1‑150 (MQ=255)
                                               cgcTGCCACGCGGCCATGGATGTGCTGAAGAACACCCGCGACGCCAAGGGCCGCGAGTTCATCGTGCACAAGATGCCGATCCCGGGCCCGCTGTTTGCCACTGCCGAAGAGTGCGCCGGTGTCGACCAGGTCGCCGGCAGCCAGGAGCGt                                                                                                  <  1:1342053/150‑1 (MQ=255)
                                                gcTGCCACGCGGCCATGGATGTGCTGAAGAACACCCGCGACGCCAAGGGCCGCGAGTTCATCGTGCACAAGATGCCGATCCCGGGCCCGCTGTTTGCCACTGCCGAAGAGTGCGCCGGTGTCGACCAGGTCGCCGGCAGCCAGGAGCGTg                                                                                                 <  2:419347/150‑1 (MQ=255)
                                                  tGCCACGCGGCCATGGATGTGCTGAAGGACACCCGCGACGCCAAGGGCCGCGAGTTCATCGTGCACAAGATGCCGATCCCGGGCCCGCTGTTTGCCACTGCCGAAGAGTGCGCCGGTGTCGACCAGGTCGCCGGCAGCCAGGAGCGTGAc                                                                                               <  2:1880993/150‑1 (MQ=255)
                                                              aTGGATGTGCTGAAGAACACCCGCGACGCCAAGGGCCGCGAGTTCATCGTGCACAAGATGCCGATCCCGGGCCCGCTGTTTGCCACTGCCGAAGAGTGCGCCGGTGTCGACCAGGTCGCCGGCAGCCAGGAGCGTGACCCGTCGGTgcgc                                                                                   >  2:77148/1‑150 (MQ=255)
                                                              aTGGATGTGCTGAAGAACACCCGCGACGCCAAGGGCCGCGAGTTCATCGTGCACAAGATGCCGATCCCGGGCCCGCTGTTTGCCACTGCCGAAGAGAGCGCCGGTGTCGACCAGGTCGCCGGCAGCCAGGAGCGTGACCCGTCGGTgcgc                                                                                   >  2:547576/1‑150 (MQ=255)
                                                                  aTGTGCTGAAGAACACCCGCGACGCCAAGGGCCGCGAGTTCATCGTGCACAAGATGCCGATCCCGGGCCCGCTGTTTGCCACTGCCGAAGAGTGCGCCGGTGTCGACCAGGTCGcc                                                                                                                 <  2:666937/116‑1 (MQ=255)
                                                                  aTGTGCTGAAGAACACCCGCGACGCCAAGGGCCGCGAGTTCATCGTGCACAAGATGCCGATCCCGGGCCCGCTGTTTGCCACTGCCGAAGAGTGCGCCGGTGTCGACCAGGTCGcc                                                                                                                 >  1:666937/1‑116 (MQ=255)
                                                                     tgCTGAAGAACACCCGCGACGCCAAGGGCCGCGAGTTCATCGTGCACAAGATGCCGATCCCGGGCCCGCTGTTTGCCACTGCCGAAGAGTGCGCCGGTGTCGACCAGGTCGCCGGCAGCCAGGAGCGTGACCCGTCGGTGCGCCTGGCgg                                                                            <  2:167441/150‑1 (MQ=255)
                                                                       cTGAAGAACACCCGCGACGCCAAGGGCCGCGAGTTCATCGTGCACAAGATGCCGATCCCGGGCCCGCTGTTTGCCACTGCCGAAGAGTGCGCCGGTGTCGACCAGGTCGCCGGCAGCCAGGAGCGTGACCCGTCGGTGCGCCTGGCGGGt                                                                          >  2:1247304/1‑150 (MQ=255)
                                                                       cTGAAGAACACCCGCGACGCCAAGGGCCGCGAGTTCATCGTGCACAAGATGCCGATCCCGGGCCCGCTGTTTGCCACTGCCGAAGAGTGCGCCGGTGTCGACCAGGTCGCCGGCAGCCAGGAGCGTGACCCGTCGGTGCGCCTGGCGGGt                                                                          >  2:869315/1‑150 (MQ=255)
                                                                         gaagaaCACCCGCGACGCCAAGGGCCGCGAGTTCATCGTGCACAAGATGCCGATCCCGGGCCCGCTGTTTGCCACTGCCGAAGAGTGCGCCGGTGTCGACCAGGTCGCCGGCAGCCAGGAGCGTGACCCGTCGGTGCGCCTGGCGGGTTc                                                                        <  1:32087/150‑1 (MQ=255)
                                                                            gaaCACCCGCGACGCCAAGGGCCGCGAGTTCATCGTGCACAAGATGCCGATCCCGGGCCCGCTGTTTGCCACTGCCGAAGAGTGCGCCGGTGTCGACCAGGTCGCCGGCAGCCAGGAGCGTGACCCGt                                                                                           >  1:433779/1‑128 (MQ=255)
                                                                            gaaCACCCGCGACGCCAAGGGCCGCGAGTTCATCGTGCACAAGATGCCGATCCCGGGCCCGCTGTTTGCCACTGCCGAAGAGTGCGCCGGTGTCGACCAGGTCGCCGGCAGCCAGGAGCGTGACCCGt                                                                                           <  2:433779/128‑1 (MQ=255)
                                                                                 cgcgctACGCCAAGGGCCGCGCGTTCCTCGTGCACACGCTGGCGCTCCCGTGCCCGCTGTTTGCCACTGCCGAAGAGTGCGCCGGTGTCGACCAGGTCGCCGGCAGCCAGGAGCGTGACCCGTCGGTGCGCCTGGCGGGTTCGTACGTGa                                                                <  2:1260433/148‑1 (MQ=255)
                                                                                 cccGCGACGCCAAGGGCCGCGAGTTCATCGTGCACAAGATGCCGATCCCGGGCCCGCTGTTTGCCACTGCCGAAGAGTGCGCCGGTGTCGACCAGGTCGCCGGCAGCCAGGAGCGTGACCCGTCGGTGCGCCTGGCGGGTTCGTACGTGa                                                                >  2:433284/1‑150 (MQ=255)
                                                                                  ccGCGACGCCAAGGGCCGCGAGTTCATCGTGCACAAGATGCCGATCCCGGGCCCGCTGTTTGCCGCTGCCGAAGAGTGCGCCGGTGTCGACCAGGTCGCCGGCAGCCAGGAg                                                                                                     >  1:1544479/1‑112 (MQ=255)
                                                                                  ccGCGACGCCAAGGGCCGCGAGTTCATCGTGCACAAGATGCCGATCCCGGGCCCGCTGTTTGCCGCTGCCGAAGAGTGCGCCGGTGTCGACCAGGTCGCCGGCAGCCAGGAg                                                                                                     >  1:1544970/1‑112 (MQ=255)
                                                                                  ccGCGACGCCAAGGGCCGCGAGTTCATCGTGCACAAGATGCCGATCCCGGGCCCGCTGTTTGCCGCTGCCGAAGAGTGCGCCGGTGTCGACCAGGTCGCCGGCAGCCAGGAg                                                                                                     <  2:1546389/112‑1 (MQ=255)
                                                                                  ccGCGACGCCAAGGGCCGCGAGTTCATCGTGCACAAGATGCCGATCCCGGGCCCGCTGTTTGCCGCTGCCGAAGAGTGCGCCGGTGTCGACCAGGTCGCCGGCAGCCAGGAg                                                                                                     <  2:1544970/112‑1 (MQ=255)
                                                                                  ccGCGACGCCAAGGGCCGCGAGTTCATCGTGCACAAGATGCCGATCCCGGGCCCGCTGTTTGCCGCTGCCGAAGAGTGCGCCGGTGTCGACCAGGTCGCCGGCAGCCAGGAg                                                                                                     >  1:1546389/1‑112 (MQ=255)
                                                                                  ccGCGACGCCAAGGGCCGCGAGTTCATCGTGCACAAGATGCCGATCCCGGGCCCGCTGTTTGCCGCTGCCGAAGAGTGCGCCGGTGTCGACCAGGTCGCCGGCAGCCAGGAg                                                                                                     <  2:1544479/112‑1 (MQ=255)
                                                                                  ccGCGACGCCAAGGGCCGCGAGTTCATCGTGCACAAGATGCCGATCCCGGGCCCGCTGTTTGCCACTGCCGAAGAGTGCGCCGGTGTCGACCAGGTCGCCGGCAGCCAGGAGCGTGACCCGTCGGTGCGCCTGGCGGGTTCGTACGt                                                                  <  2:257725/147‑1 (MQ=255)
                                                                                  ccGCGACGCCAAGGGCCGCGAGTTCATCGTGCACAAGATGCCGATCCCGGGCCCGCTGTTTGCCACTGCCGAAGAGTGCGCCGGTGTCGACCAGGTCGCCGGCAGCCAGGAGCGTGACCCGTCGGTGCGCCTGGCGGGTTCGTACGt                                                                  >  1:257725/1‑147 (MQ=255)
                                                                                  ccGCGACGCCAAGGGCCGCGAGTTCATCGTGCACAAGATGCCGATCCCGGGCCCGCTGTTTGCCACTGCCGAAGAGTGCGCCGGTGTCGACCAGGTCGCCGGCAGCCAGGAGCGTGACCCGTCGGTGCGCCTGGCGGGTTCGTACGTGaa                                                               <  2:1583579/150‑1 (MQ=255)
                                                                                        cGCCAAGGGCCGCGAGTTCATCGTGCACAAGATGCCGATCCCGGGCCCGCTGTTTGCCACTGCCGAAGAGTGCGCCGGTGTCGACCAGGTCGCCGGCAGCCAGGAGCGTGACCCGTCGGTGCGCCTGGCGGGTTCGTACGTGAACTTCCt                                                         <  2:262136/150‑1 (MQ=255)
                                                                                        cGCCAAGGGCCACGAGTTCATCGTGCACAAGATGCCGATCCCGGGCCCGCTGTTTGCCACTGCCGAAGAGTGCGCCGGTGTCGACCAGGTCGCCGGCAGCCAGGAGCGTGACCCGTCGGTGCGCCTGGCGGGTTCGTACGTGAACTTCCt                                                         <  1:1238938/150‑1 (MQ=255)
                                                                                         gCCAAGGGCCGCGAGTTCATCGTGCACAAGATGCCGATCCCGGGCCCGCTGTTTGCCACTGCCGAAGAGTGAGAAGGTTTCGACCAGGTCGCCGGCAGCCAGGAGCGTGACCCGTCGGTGCGCCTGGCGGGTTCGTACGTGAACTTCCTg                                                        >  1:1669403/1‑150 (MQ=255)
                                                                                            aaGGGCCGCGAGTTCATCGTGCACAAGATGCCGATCCCGGGCCCGCTGTTTGCCACTGCCGAAGAGTGCGCCGGTGTCGACCAGGTCGCCGGCAGCCAGGAGCGTGACCCGTCGGTGCGCCTGGCGGGTTCGTACGt                                                                  >  2:1216114/1‑137 (MQ=255)
                                                                                            aaGGGCCGCGAGTTCATCGTGCACAAGATGCCGATCCCGGGCCCGCTGTTTGCCACTGCCGAAGAGTGCGCCGGTGTCGACCAGGTCGCCGGCAGCCAGGAGCGTGACCCGTCGGTGCGCCTGGCGGGTTCGTACGt                                                                  <  1:1216114/137‑1 (MQ=255)
                                                                                            aaGGGCCGCGAGTTCATCGTGCACAAGATGCCGATCCCGGGCCCGCTGTTTGCCACTGCCGAAGAGTGCGCCGGTGTCGACCAGGTCGCCGGCAGCCAGGAGCGTGACCCGTCGGTGCGCCTGGCGGGTTCGTACGTGAACTTCCTGAtt                                                     >  2:389412/1‑150 (MQ=255)
                                                                                              gggCCGCGAGTTCATCGTGCACAAGATGCCGATCCCGGGCCCGCTGTTTGCCACTGCCGAAGAGTGCGCCGGTGTCGACCAGGTCGCCGGCAGCCAGGAGCGTGACCCGTCGGTGCGCCTGGCGGGTTCGTACGTGAACTTCCTGATTGt                                                   >  1:1835465/1‑150 (MQ=255)
                                                                                                     gAGTTCATCGTGCACAAGATGCCGATCCCGGGCCCGCTGTTTGCCACTGACGAAGAGTGCGCCGGTGTCGACCAGGTCGCCGGCAGCCAGGAGCGTGACCCGTCGGTGCGCCTGGCGGGTTCGTACGTGAACTTCCTGATTGTTAACGGt                                            <  2:982354/150‑1 (MQ=255)
                                                                                                     gAGTTCATCGTGCACAAGATGCCGAGCCCGGGCCCGCTGTTTGCCACTGCCGAAGAGTGCGCCGGTGTCGACCAGGTCGCCGGCAGCCAGGAGCGTGACCCGTCGGTGCGCCTGGCGGGTTCGTACGTGAACTTCCTGATTGTTAACGGt                                            <  2:562125/150‑1 (MQ=255)
                                                                                                      aGTTCATCGTGCACAAGATGCCGATCCCGGGCCCGCTGTTTGCCACTGCCGAAGAGTGCGCCGGTGTCGACCAGGTCGCCGGCAGCCAGGAGCGTGACCCGTCGGTGCGCCTGGCGGGTTCGTACGTGAACTTCCTGATTGTTAACGGTg                                           <  1:869315/150‑1 (MQ=255)
                                                                                                          cATCGTGCACAAGATGCCGATCCCGGGCCCGCTGTTTGCCACTGCCGAAGAGTGCGCCGGTGTCGACCAGGTCGCCGGCAGCCAGGAGCGTGACCCGTCGGTGCGCCTGGCGGGTTCGTACGTGAACTTCCTGATTGTTAACGGTGGcat                                       >  2:162055/1‑150 (MQ=255)
                                                                                                            tCGTGCACAAGATGCCGATCCCGGGCCCGCTGTTTGCCACTGCCGAAGAGTGCGCCGGTGTCGACCAGGTCGCCGGCAGCCAGGAGCGTGACCCGTCGGTGCGCCTGGc                                                                              >  2:440309/1‑109 (MQ=255)
                                                                                                            tCGTGCACAAGATGCCGATCCCGGGCCCGCTGTTTGCCACTGCCGAAGAGTGCGCCGGTGTCGACCAGGTCGCCGGCAGCCAGGAGCGTGACCCGTCGGTGCGCCTGGc                                                                              <  1:440309/109‑1 (MQ=255)
                                                                                                              gTGCACAAGATGCCGATCCCGGGCCCGCTGTTTGCCACTGCCGAAGAGTGCGCCGGTGTCGACCAGGTCGCCGGCAGCCAGGAGCGTGACCCGTCGGTGCGCCTGGCGGGTTCGTACGTGAACTTCCTGATTGTTAACGGTgg                                          >  1:1133362/1‑143 (MQ=255)
                                                                                                              gTGCACAAGATGCCGATCCCGGGCCCGCTGTTTGCCACTGCCGAAGAGTGCGCCGGTGTCGACCAGGTCGCCGGCAGCCAGGAGCGTGACCCGTCGGTGCGCCTGGCGGGTTCGTACGTGAACTTCCTGATTGTTAACGGTgg                                          <  2:1133362/143‑1 (MQ=255)
                                                                                                                 cacaAGATGCCGATCCCGGGCCCGCTGTTTGCCACTGCCGAAGAGTGCGCCGGTGTCGACCAGGTCGCCGGCAGCCAGGAGCGTGACCCGTCGGTGCGCCTGGCGGGTTCGTACGTGAACTTCCTGATTGTTAACGGTGGCATCATTGcc                                <  1:744040/150‑1 (MQ=255)
                                                                                                                 cacaAGATGCCGATCCCGGGCCCGCTGTTTGCCACTGCCGAAGAGTGCGCCGGTGTCGACCAGGTCGCCGGCAGCCAGGAGCGTGACCCGTCGGTGCGCCTGGCGGGTTCGTACGTGAACTTCCTGATTGTTAACGGTGGCATCATTGcc                                <  1:1912575/150‑1 (MQ=255)
                                                                                                                 cacaAGATGCCGATCCCGGGCCCGCTGTTTGCCACTGCCGAAGAGTGCGCCGGTGTCGACCAGGTCGCCGGCAGCCAGGAGCGTGACCCGTCGGTGCGCCTGGCGGGTTCGTACGTGAACTTCCTGATTGTTAACGGTGGCATCATTGcc                                >  2:744040/1‑150 (MQ=255)
                                                                                                                   caAGATGCCGATCCCGGGCCCGCTGTTTGCCACTGCCGAAGAGTGCGCCGGTGTCGACCAGGTCGCCGGCAGCCAGGAGCGTGACCCGTCGGTGCGCCTGGCGGGTTCGTACGTGAACTTCCTGATTGTTAACGGTGGCATCATTGcccc                              <  2:76679/150‑1 (MQ=255)
                                                                                                                       aTGCCGATCCCGGGCCCGCTGTTTGCCACTGCCGAAGAGTGCGCCGGTGTCGACCAGGTCGCCGGCAGCCAGGAGCGTGACCCGTCGGTGCGCCTGGCGGGTTCGTACGt                                                                  >  1:745207/1‑110 (MQ=255)
                                                                                                                       aTGCCGATCCCGGGCCCGCTGTTTGCCACTGCCGAAGAGTGCGCCGGTGTCGACCAGGTCGCCGGCAGCCAGGAGCGTGACCCGTCGGTGCGCCTGGCGGGTTCGTACGt                                                                  <  2:745207/110‑1 (MQ=255)
                                                                                                                          ccGATCCCGGGCCCGCTGTTTGCCACTGGCGAAGAGTGCGCCGGTGTCGACCAGGTCGCCGGCAGCCAGGGGCGTGACCCGTCGGTGCGCCGTGCGGGTACGTACGTGAACTTCCTGATTGTTAACGGTGGCATCATTGCCCCAAGCTTc                       >  2:998387/1‑150 (MQ=255)
                                                                                                                          ccGATCCCGGGCCCGCTGTTTGCCACTGCCGAAGAGTGCGCCGGTGTCGACCAGGTCGCCGGCAGCCAGGAGCGTGACCCGTCGGTGCGCCTGGCGGGTTCGTACGTGAACTTCCTGATTGTTAACGGTGGCATCATTGCCCCAAGCTTc                       <  1:162055/150‑1 (MQ=255)
                                                                                                                           cGATCCCGGGCCCGCTGTTTGCCACTGCCGAAGAGTGCGCCGGTGTCGACCAGGTCGCCGGCAGCCAGGAGCGTGACCCGTCGGTGCGCCTGGCGGGTTCGTACGTGAACTTCCTGATTGTTAACGGTGGCATCATTGCCCCAAGCTTcg                      >  2:1823940/1‑150 (MQ=255)
                                                                                                                            gATCCCGGGCCCGCTGTTTGCCACTGCCGAAGAGTGCGCCGGTGTCGACCAGGTCGCCGGCAGCCAGGAGCGTGACCCGTCGGTGCGCCTGGCGGGTTCGTACGTGAACTTCCTGATTGTTAACGGTGGCATCATTGCCCCAAGCTTcga                     <  1:1157835/150‑1 (MQ=255)
                                                                                                                            cggcccGGGCCCGCTGTTTGCCACTGGCGAAGAGTGCGCCGGTGTCGACCAGGTCGCCGGCAGACAGGAGCGTGACCCGTCGGTGCGCCTGGCGGGTTCGTACGTGAACTTCCTGATTGTTAACGGTGGCATCATTGCCCCAAGCTTcga                     <  1:28032/147‑1 (MQ=255)
                                                                                                                                 cGGGCCCGCTGTTTGCCACTGCCGAAGCGCGACACGGTGTTGATCAGGTCGCCGGCAGCCAGGAGCGTGACCCGTCGGTGCGCCTGGCGGGTTCGTACGTGAACTTCCTGATTGTTAACGGTGGCATCATTGCCCCAAGCTTCGACGAcc                <  2:1590703/150‑1 (MQ=255)
                                                                                                                                        gCTGTTTGCCACTGCCGAAGAGTGCGCCGGTGTCGACCAGGTCGCCGGCAGCCAGGAGCGTGACCCGTCGGTGCGCCTGGCGGGTTCGTACGTGAACTTCCTGATTGTTAACGGTGGCATCATTGcc                                >  2:183795/1‑127 (MQ=255)
                                                                                                                                        gCTGTTTGCCACTGCCGAAGAGTGCGCCGGTGTCGACCAGGTCGCCGGCAGCCAGGAGCGTGACCCGTCGGTGCGCCTGGCGGGTTCGTACGTGAACTTCCTGATTGTTAACGGTGGCATCATTGcc                                <  1:183795/127‑1 (MQ=255)
                                                                                                                                        gCTGTTTGCCACTGCCGAAGAGTGCGCCGGTGTCGACCAGGTCGCCGGCAGCCAGGAGCGTGACCCGTCGGTGCGCCTGGCGGGTTCGTACGTGAACTTCCTGATTGTTAACGGTGGCATCATTGCCCCAAGCTTCGACGACCCGGCAGa         >  1:1693515/1‑150 (MQ=255)
                                                                                                                                          tGTTTGCCACTGCCGAAGAGTGCGCCGGTGTCGACCAGGTCGCCGGCAGCCAGGAGCGTGACCCGTCGGTGCGCCTGGCGGGTTCGTACGTGAACTTCCTGATTGTTAACGGTGGCATCATTGCCCCAAGCTTCGACGACCCGGCAGATg       >  2:1144386/1‑150 (MQ=255)
                                                                                                                                            tttGCCACTGCCGAAGAGTGCGCCGGTGTCGACCAGGTCGCCGGCAGCCAGGAGCGTGACCCGTCGGTGCGCCTGGCGGGTTCGTACGTGAACTTCCTGCGCCATTGCGGTGGCATCATTGCCCCAAGCTTAAACGACCCGGCAGATGcc     >  1:1152118/1‑150 (MQ=255)
                                                                                                                                               gCCACTGCCGAAGAGTGCGCCGGTGTCGACCAGGTCGCCGGCAGCCAGGAGCGTGACCCGTCGGTGCGCCTGGCGGGTTCGTACGTGAACTTCCTGATTGTTAACGGTGGCATCATTGCCCCAAGCTTCGACGACCCGGCAGATGCCGaa  <  2:1835465/150‑1 (MQ=255)
                                                                                                                                                 cACTGCCGAAGAGTGCGCCGGTGTCGACCAGGTCGCCGGCAGCCAGGAGCGTGACCCGTCGGTGCGCCTGGCGGGTTCGTACGTGAACTTc                                                           >  2:503260/1‑91 (MQ=255)
                                                                                                                                                 cACTGCCGAAGAGTGCGCCGGTGTCGACCAGGTCGCCGGCAGCCAGGAGCGTGACCCGTCGGTGCGCCTGGCGGGTTCGTACGTGAACTTc                                                           <  1:503260/91‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
AAGTGTTACTGGCCTGGACCGATGATTCCAACGACCCCAACTATGCGCGCTGCCACGCGGCCATGGATGTGCTGAAGAACACCCGCGACGCCAAGGGCCGCGAGTTCATCGTGCACAAGATGCCGATCCCGGGCCCGCTGTTTGCCACTGCCGAAGAGTGCGCCGGTGTCGACCAGGTCGCCGGCAGCCAGGAGCGTGACCCGTCGGTGCGCCTGGCGGGTTCGTACGTGAACTTCCTGATTGTTAACGGTGGCATCATTGCCCCAAGCTTCGACGACCCGGCAGATGCCGAA  >  NC_002947/321294‑321586

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: