Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,060,265 A→C 7.0% L217R (CTG→CGG)  cbl ← ssuEADCB/tauABCD operon transcriptional activator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,060,2650AC 7.0% 273.9 / 9.6 97L217R (CTG→CGG) cblssuEADCB/tauABCD operon transcriptional activator
Reads supporting (aligned to +/- strand):  ref base A (51/39);  new base C (2/5);  total (53/44)
Fisher's exact test for biased strand distribution p-value = 2.39e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.22e-01

ATCAAAAAGATGCCGCGTATCCAGGCGGATTAAATTCTCTTCCTCTTGTTCGCCACTGGATTGCTCGGCAACTAATCCGATCCCAAGCCCAAGAGCAACATAGGTTTTAATGACATCAGAATCCTGCGCACTTAATACAATATCTGCCAGCAAACCTTTGCGGGCAAATGCGTCATCAATACGTGAGCGCCCCGTAATCCCCTGTCGGTAAGTGATTAACGGCCACTTCGCTATTGATTCCAGCGTCAATGGTGAAATTTGCGTCAAGGGATGATCGTGTGGAACAAGCAAACTATGG  >  NC_000913/2060117‑2060414
                                                                                                                                                    |                                                                                                                                                     
aTCAAAAAGATGCCGCGTATCCAGGCGGATTAAATTCTCTTCCTCTTGTTCGCCACTGGATTGCTCGGCAACTAATCCGATCCCAAGCCCAAGAGCAACATAGGTTTTAATGACATCCGAATCCTGCGCACTTAATACAATATCTGCCAg                                                                                                                                                      <  1:1663864/150‑1 (MQ=255)
aTCAAAAAGATGCCGCGTATCCAGGCGGATTAAATTCTCTTCCTCTTGTTCGCCACTGGATTGCTCGGCAACTAATCCGATCCCAAGCCCAAGAGCAACATAGGTTTTAATGACATCAGAATCCTGCGCACTTAATACAATATCTGCCAg                                                                                                                                                      >  2:597905/1‑150 (MQ=255)
aTCAAAAAGATGCCGCGTATCAAGGCGGATTAAATTCTCTTCCTCTTGTTCGCCACTGGATTGCTCGGCAACTAATCCGATCGCAAGCCCAAGCGCAACCTAGGTTTTAATGACATCAGAATCCTGCGCACTTAAGTCAATATCTGGCAg                                                                                                                                                      >  1:75068/1‑150 (MQ=255)
 tCAAAAAGATGCCGCGTATCCAGGCGGATTAAATTCTCTTCCTCTTGTTCGCCACTGGATTGCTCGGCAACTAATCCGATCCCAAGCCCAAGAGCAACATAGGTTTTAATGACATCAGAATCCTGCGCACTTAATACAATATCTGCCAGc                                                                                                                                                     <  1:746265/150‑1 (MQ=255)
  cAAAAAGATGCCGCGTATCCAGGCGGATTAAATTCTCTTCCTCTTGTTCGCCACTGGATTGCTCGGCAACTACTCCGATCCCACGCCCAAGAGCAACATAGCTTTTAATGACATCAGAATCCTGCGCACGTACTACAATATCTGCCAGCa                                                                                                                                                    >  1:263988/1‑150 (MQ=255)
   aaaaaGATGCCGCGTATCCAGGCGGATTAAATTCTCTTCCTCTTGTTCGCCACTGGATTGCTCGGCAACTAATCCGATCCCAAGCCCAAGAGCAACATAGGTTTTAATGACATCAGAATCCTGCGCACTTAATACAATATCTGCCAGCaa                                                                                                                                                   >  1:630422/1‑150 (MQ=255)
      aaGATGCCGCGTATCCAGGCGGATTAAATTCTCTTCCTCTTGTTCGCCACTGGATTGCTCGGCAACTAATCCGATCCCAAGCCCAAGAGCAACATAGGTTTTAATGCCATCAGAATCCTGCGCACTTAATACAATATCTGCCAGCAAAcc                                                                                                                                                >  2:593240/1‑150 (MQ=255)
        gATGCCGCGTATCCAGGCGGATTAAATTCTCTTCCTCTTGTTCGCCACTGGATTGCTCGGCAACTAATCCGATCCCAAGCCCAAGAGCAACATAGGTTTTAATGACATCAGAATCCTGCGCACTTAATACAATATCTGCCAGCAAACCtt                                                                                                                                              >  1:1504634/1‑150 (MQ=255)
             cgcgTATCCAGGCGGATTAAATTCTCTTCCTCTTGTTCGCCACTGGATTGCTCGGCAACTAATCCGATCCCAAGCCCAAGAGCAACATAGGTTTTAATGACATCAGAATCCTGCGCACTTAATACAATATCTGCCAGCAAACCTTTGCgg                                                                                                                                         <  2:1142429/150‑1 (MQ=255)
             cgcgTATCCAGGCGGATTAAATTCTCTTCCTCTTGTTCGCCACTGGATTGCTAGGCAACTAATCCGATCGCAATCCCAATAGCAACATAGGTTTTAATGACATCAGAATCCTGCGCACTTAATACAATATCTGCCAGCAAACCTTTGCgg                                                                                                                                         <  1:59465/150‑1 (MQ=255)
              gcgTATCCAGGCGGATTAAATTCTCTTCCTCTTGTTCGCCACTGGATTGCTCGGCAACTAATCCGATCCCAAGCCCAAGAGCAACATAGGTTTTAATGACATCAGAATCCTGCGCACTTAATACAATATCTGCCAGCAAACCTTTGCggg                                                                                                                                        <  1:542810/150‑1 (MQ=255)
                gTATCCAGGCGGATTAAATTCTCTTCCTCTTGTTCGCCACTGGATTGCTCGGCAACTAATCCGATCCCAAGCCCAAGAGCAACATAGGTTTTAATGACATCAGAATCCTGCGCACTTAATACAATATCTGCCCGCAAAGGTTTGCTGGCa                                                                                                                                      <  1:1459568/150‑1 (MQ=255)
                 tATCCAGGCGGATTAAATTCTCTTCCTCTTGTTCGCCACTGGATTGCTCGGCAACTAATCCGATCCCAAGCCCAAGAGCAACATAGGTTTTAATGACATCAGAATCCTGCGCACTTAATACAATATCTGCCCGCAAAGGTTTGCTGGCag                                                                                                                                     <  1:578768/150‑2 (MQ=255)
                 tATCCAGGCGGATTAAATTCTCTTCCTCTTGTACGCCACTGGATTGCTCGGAAACTAATCCGATCCCAAGCCCAAGAGCAACATAGGTTTTAATGACATCAGAATCCTCCGCACTTAATACAATATCTGCCCGCAAAGGTTTGCTGGCag                                                                                                                                     <  2:1696699/150‑2 (MQ=255)
                  aTCCAGGCGGATTAAATTCTCTTCCTCTTGTTCGACACTGGATTGCTCGGCAACTAATCCGATCCCAAGCCCAAGAGCAACATAGGTTTTAATGACATCAGAATCCAGCGCACTTAATACATTATCTGCCAGCAAACCTTTGCGGGCaaa                                                                                                                                    >  2:793618/1‑150 (MQ=255)
                     cAGGCGGATTAAATTCTCTTCCTCTTGTTCGCCACTGGATTGCTCGGCAACTAATCCGATCCCAAGCCCAAGAGCAACATAGGTTTTAATGACATCAGAATCCTGCGCACTTAATACAATATCTGCCAGCAAACCTTTGCGGGCAAATGc                                                                                                                                 >  2:1093555/1‑150 (MQ=255)
                      aGGCGGATTAAATTCTCTTCCTCTTGTTCGCCACTGGATTGCTCGGCAACTAATCCGATCCCAAGCCCAAGAGCAACATAGGTTTTAATGACATCAGAATCCTGCGCACTTAATACAATATCTGCCAGCAAACCTTTGCGGGCAAATGCg                                                                                                                                <  2:1663763/150‑1 (MQ=255)
                      aGGCGGATTAAATTCTCTTCCTCTTGTTCGCCACTGGAGTGCTCGGCAACTAATCCGATCCCAAGCCCAAGAGCAACATAGGTTTTAATGACATCAGAATCCTGCGCACTTAATACAATATCTGCCAGCAAACCTTTGCGGGCAAATGCg                                                                                                                                <  1:1443014/150‑1 (MQ=255)
                       ggCGGATTAAATTCTCTTCCTCTTGTTCGCCACTGGATTGCTCGGCAACTAATCCGATCCCAAGCCCAAGAGCAACATAGGTTTTAATGACATCAGAATCCTGCGCACTTAATACAATATCTGCCAGCAAACCTTTGCGGGCAAATGCGt                                                                                                                               <  1:685106/150‑1 (MQ=255)
                       ggCGGATTAAATTCTCTTCCTCTTGTTCGCCACTGGATTGCTCGGCAACTAATCCGATCCCAAGCCCAAGAGCAACATAGGTTTTAATGACATCAGAATCCTGCGCACTTAATACAATATCTGCCAGCAAACCTTTGCGGGCAAATGCGt                                                                                                                               >  1:1066649/1‑150 (MQ=255)
                       ggCGGATTAAATTCTCTTCCTCTTGTTCGCCACTGGATTGCTCGGCAACTAATCCGATCCCAAGCCCAAGAGCAACATAGGTTTTAATGACATCAGAATCCTGCGCACTTAATACAATATCTGCCAGCAAACCTTTGCGGGCAAATGCGt                                                                                                                               <  2:741484/150‑1 (MQ=255)
                        gCGGATTAAATTCTCTTCCTCTTGTTCGCCACTGGATTGCTCGGCAACTAATCCGATCCCAAGCCCAAGAGCAACATAGGTTTTAATGACATCAGAATCCTGCGCACTTAATACAATATCTGCCCGCAAAGGTTTGCTGGCAGAtattgt                                                                                                                              <  1:1555351/150‑6 (MQ=255)
                         cGGATTAAATTCTCTTCCTCTTGTTCGCCACTGGATTGCTCGGCAACTAATCCGATCCCAAGCCCAAGAGCAACATAGGTTTTAATGACATCAGAATCCTGCGCACTTAATACAATATCTGCCAGCAAACCTTTGCGGGCAAATGCGtca                                                                                                                             >  1:742948/1‑150 (MQ=255)
                           gATTAAATTCTCTTCCTCTTGTTCGCCACTGGATTGCTCGGCAACTAATCCGATCCCAAGCCCATGAGCAACATAGGTTTTAATGACATCAGAATCCTGCGCACTTAATACAATATCTGCCAGCAAACCTTTGCGGGCAAATGCGtcatc                                                                                                                           >  2:661969/1‑150 (MQ=255)
                                 aTTCTCTTCCTCTTGTTCGCCACAGGATTGCTCGGCAACTAATCCGATCCCAAGCCCAAGAGCAACATAGGTTTTAATGACATCAGAATCCTGCGCACTTAATACAAAATCTGCCAGCAAACCTTTGCGGGCAAATGCGTAATCAATACg                                                                                                                     >  1:1695903/1‑150 (MQ=255)
                                   tctcTTCCTCTTGTTCGCCACTGGATTGCTCGGCAACTAATCCGATCCCAAGCCCAAGAGCAACATAGGTTTTAATGACATCAGAATCCTGCGCACTTAATACAATATCTGCCAGCAAACCTTTGCGGGCAAATGCGTCATCAATACGTg                                                                                                                   <  2:12548/150‑1 (MQ=255)
                                     tctTCCTCTTGTTCGCCACTGGATTGCTCGGCAACTAATCCGATCCCAAGCCCAAGAGCAACATAGGTTTTAATGACATCAGAATCCTGCGCACTTAATACAATATCTGCCAGCAAACCTTTGCGGGCAAATGCGTCATCAATACGTGAg                                                                                                                 <  1:1639756/150‑1 (MQ=255)
                                         ccTCTTGTTCGCCACTGGATTGCTCGGCAACTAATCCGATCCCAAGCCCAAGAGCAACATAGGTTTTAATGACAGCAGAATCCTGCGCACTTACTACAATATCTGCCAGCAAACCTTTGCGGGCAAATGCGTCATCAATACGTGAGCGcc                                                                                                             >  1:1049779/1‑150 (MQ=255)
                                          ctctTGTTCGCCACTGGATTGCTCGGCAACTAATCCGATCCCAAGCCCAAGAGCAACATAGGTTTTAATGACATCAGAATCCTGCGCACTTAATACAATATCTGCCAGCAAACCTTTGCGGGCAAATGCGTCATCAATACGTGAGCGccc                                                                                                            >  2:824863/1‑150 (MQ=255)
                                            ctTGTTCGCCACTGGATTGCTCGGCAACTAATCCGATCCCAAGGCCAAGAGCAACATATGTTTTAATGACATCAGAATCCTGCGCACTTAATACAATATCTGCCAGCAAACCTTTGCGGGCAAATGCGTCATCAATACGTGAGCGCCCCg                                                                                                          <  2:931242/150‑1 (MQ=255)
                                            ctTGTTCGCCACTGGATTGCTCGGAAACTAATCCGATCCCAAGCCCAAGAGCAACATAGGTTTTAATGTCATCAGAATCCTGCGCACTTAATACAATATCTGCCAGCAAACCTTTGCGGGCAAATGCGTCATCAATACGTGAGCGcccct                                                                                                          <  2:1398326/150‑2 (MQ=255)
                                               gTTCGCCACTGGATTGCTCGGCAACTAATCCGATCCCAAGCCCAAGAGCAACATAGGTTTTAATGACATCAGAATCCTGCGCACTTAATACAATATCTGCCAGCAAACCTTTGCGGGCAAATGCGTCATCAATACGTGAGCGCCCCGTaa                                                                                                       >  2:1370892/1‑150 (MQ=255)
                                               gTTCGCCACTGGATTGCTCGGCAACTAATCCGATCCCAAGCCCAAGAGCAACATAGGTTTTAATGACATCAGAATCCTGCGCACTTAATACAATATCTGCCAGCAAACCTTTGCGGGCAAATGCGTCATCAATACGTGAGCGCCCCGTaa                                                                                                       >  2:525608/1‑150 (MQ=255)
                                               gTTCGCCACTGGATTGCTCGGCAACTAATCCGATCCCAAGCCCAAGAGCAACATAGGTTTTAATGACATCAGAATCCTGCGCACTTAATACAATATCTGCCAGCAAACCTTTGCGGGCAAATGCGTCATCAATACGTGAGCGCCCCGTaa                                                                                                       >  1:879155/1‑150 (MQ=255)
                                                  cGCCACTGGATTGCTCGGCAACTAATCCGATCCCAAGCCCAAGAGCAACATAGGTTTTAATGACATCAGAATCCTGCGCACTTAATACAATATCTGCCAGCAAACCTTTGCGGGCAAATGCGTCATCAATACGTGAGCGCCCCGTAATcc                                                                                                    >  2:1324311/1‑150 (MQ=255)
                                                     cACTGGATTGCTCGGCAACTAATCCGATCCCAAGCCCAAGAGCAACATAGGTTTTAATGACATCAGAATCCTGCGCACTTAATACAATATCTGCCAGCAAACCTTTGCGGGCTAATGCGTCATCAATACGTGAGCGCCCCGTAATCCCCt                                                                                                 >  1:251589/1‑150 (MQ=255)
                                                        tGGATTGCTCGGCAACTAATCCGATCCCAAGCCCACGAGCAACATAGGTTTTAATGACATCAGAATCCTGCGCACTTAATACAATATCTGCCAGCAAACCTTTGCGGGCAAATGCGTCTCAATACGTGAGCGCCCCGTAATCCCCTGTCg                                                                                             >  2:1070021/1‑150 (MQ=255)
                                                          gATTGCTCGGCAACTAATCCGCTCCCAAGCCCAAGAGCAACATAGGTTTTAATGACATCAGAATCCTGCGCAGTTAATACAATATCTGCCAGCAAACCTTTGCGGGCAAATGCGTCATCAATCCGTGAGCGCCCCGTAATCCCCTGTCgg                                                                                            >  2:1076413/1‑150 (MQ=255)
                                                          gATTGCTCGGCAACTAATCCGATCCCAAGCCCAAGAGCAACATAGGTTTTAATGACATCAGAATCCTGCGCACTTAATACAATATCTGCCAGCAAACCTTTGCGGGCAAATGCGTCATCAATACGTGAGCGCCCCGTAATCCCCTGTCgg                                                                                            <  1:1500445/150‑1 (MQ=255)
                                                           aTTGCTCGGCAACTAATCCGATCCCAAGCCCAAGAGCAACATAGGTTTTAATGACATCAGAATCCTGCGCACTTAATACAATATCTGCCAGCAAACCTTTGCGGGCAAATGCGTCATCAATACGTGAGCGCCCCGTAATCCCCTGTCGGt                                                                                           >  2:814871/1‑150 (MQ=255)
                                                            ttGCTCGGCAACTAATCCGATCCCAAGCCCAAGAGCAACATAGGTTTTAATGACATCAGAATCCTGCGCACTTAATACAATATCTGCCAGCAAACCTTTGCGGGCAAATGCGTCATCAATACGTGAGCGCCCCGTAATCCCCTGTCGGTa                                                                                          >  1:1583998/1‑150 (MQ=255)
                                                             tGCTCGGCAACTAATCCGATCCCAAGCCCAAGAGCAACATAGGTTTTAATGACATCAGAATCCTGCGCACTTAATACAATATCTGCCAGCAAACCTTTGCGGGCAAATGCGTCATCAATACGTGAGCGCCCCGTAATCCCCTGTCGGTaa                                                                                         <  2:742948/150‑1 (MQ=255)
                                                             tGCTCGGCAACTAATCCGATCCCAAGCCCAAGAGCAACATAGGTTTTAATGACATCAGAATCCTCCGCACTTAATACAATAACTGCCAGCAAACCTTTGCGGGCAAATGCGTCATCAATACGTGAGCGCCCCGTAATCCCCTGACGGTaa                                                                                         >  1:315467/1‑150 (MQ=255)
                                                               cTCGGCAACTAATCCGATCCCAAGCCCAAGAGCAACATAGGTTTTAATGACATCAGAATCCTGCGCACTTAATACAATATCTGCCAGCAAACCTTTGCGGGCAAATGCGTCATCAATACGTGAGCGCCCCGTAATCCCCTGTCGGTAAGt                                                                                       >  1:1391636/1‑150 (MQ=255)
                                                                   gCAACTAATCCGATCCCAAGCCCAAGAGCAACATAGGTTTTAATGACATCAGAATCCTGCGCACTTAATACAATATCTGCCAGCAAACCTTTGCGGGCAAATGCGTCATCAATACGTGAGCGCCCCGTAATCCCCTGTCGGTAAGTGAtt                                                                                   >  2:1677625/1‑150 (MQ=255)
                                                                   gCAACTAATACGATCCCAAGCCCAAGAGCAACATAGGTTTTAATGACATCAGAATCCTGCGCACTTAATACAATATCTGCCAGCAAACCTTTGCGGGCAAATGCGTCATCAATACGTGAGCGCCCCGTAATCCCCTGTCGGTAAGTGAtt                                                                                   <  1:1073260/150‑1 (MQ=255)
                                                                    cAACTAATCCGATCCCAAGCCCAAGAGCAACATAGGTTTTAATGACATCAGAATCCTGCGCACTTAATACAATATCTGCCAGCAAACCTTTGCGGGCAAATGCGTCATCAATACGTGAGCGCCCCGTAATCCCCTGTCGGTAAGTGATTa                                                                                  >  2:19335/1‑150 (MQ=255)
                                                                      aCTAATCCGATCCCAAGCCCAAGAGCAACATAGGTTTTAATGACATCAGAATCCTGCGCACTTAATACAATATCTGCCAGCAAACCTTTGCGGGCAAATGCGTCATCAATACGTGAGCGCCCCGTAATCCCCTGTCGGTAAGTGATTAAc                                                                                >  2:410344/1‑150 (MQ=255)
                                                                       cTAATCCGATCCCAAGCCCAAGAGCAACATAGGTTTTAATGACATCAGAATCCTGCGCACTTAATACAATATCTGCCAGCAAACCTTTGCGGGCAAATGCGTCATCAATACGTGAGCGCCCCGTAATCCCCTGTCGGTAAGTGATTAACg                                                                               <  1:1093555/150‑1 (MQ=255)
                                                                       cTAATCCGATCCCAAGCCCAAGAGCAACATAGGTTTTAATGACATCAGAATCCTGCGCACTTAATACAATATCTGCCAGCAAACCTTTGCGGGCAAATGCGTCATCAATACGTGAGCGCCCCGTAATCCCCTGTCGGTAAGTGATTAACg                                                                               >  2:267470/1‑150 (MQ=255)
                                                                       cTAATCCGATCCCAAGCCCAAGAGCAACATAGGTTTTAATGACATCAGAATCCTGCGCACTTAATACAATATCTGCCAGCAAACCTTTGCGGGCAAATGCGTCATCAATAAGTGAGCGCCCCGGAATCCCCTGTAGGTAAGTGATTaaag                                                                               >  1:1453836/1‑148 (MQ=255)
                                                                       cTAATCCGATCCCAAGCCCAAGAGCAACATAGGTTTTAATGACATCAGAATCCTGCACACTTAATACAATATCTGCCAGCAAACCTTTGCGGGCAAATGCGTCATCAATACGTGAGCGCCCCGTAATCCCCTGTCGGTAAGTGATTAACg                                                                               >  1:876728/1‑150 (MQ=255)
                                                                             cGATCCCCAGCCCAAGAGCAACATAGGTTTTAATGACATCAGAATCCTGCGCACTTAATACAATATCTGCCAGCAAACCTTTGCGGGCAAATGCGTCATCAATACGTGAGCGCCCCGTAATCCCCTGTCGGTAAGATAATAACGGCCACt                                                                         >  2:455167/1‑150 (MQ=255)
                                                                             cGATCCCAAGCCCAAGAGCAACATAGGTTTTAATGACATCAGAATCCTGCGCACTTAATACAATATCTGCCAGCAAACCTTTGCGGGCAAATGCGTCATCAATACGTGAGCGCCCCGTAATCCCCTGTCGGTAAGTGATTAACGGCCACt                                                                         >  1:361754/1‑150 (MQ=255)
                                                                                    aagcccaagAGCAACATAGGTTTTAATGACATCAGAATCCTGCGCACTTAATACAATATCTGCCAGCAAACCTTTGCGGGCAAATGCGTCATCAATACGTGAGCGCCCCGTAATCCCCTGTCGGTAAGTGATTAACGGCCACTTCGCTAt                                                                  <  2:920876/150‑1 (MQ=255)
                                                                                       cccaagAGCAACATAGGTTTTAATGACATCAGAATCCTGCGCACTTAATACAATATCTGCCAGCAAACCTTTGCGGGCAAATGCGTCATCAATACGTGAGCGCCCCGTAATCCCCTGTCGGTAAGTGATTAACGGCCACTTCGCTATTGa                                                               <  1:525608/150‑1 (MQ=255)
                                                                                        ccaagAGCAACATAGGTTTTAATGACATCAGAATCCTGCGCACTTAATACAATATCTGCCAGCAAACCTTTGCGGGCAAATGCGTCATCAATACGTGAGCGCCCCGTAATCCCCTGTCGGTAAGTGATTAACGGCCACTTCGCTATTGAt                                                              >  2:717081/1‑150 (MQ=255)
                                                                                        ccaagAGCAACATAGGTTTTAATGACATCAGAATCCTGCGCACTTAATACAATATCTGCCAGCAAACCTTTGCGGGCAAATGCGTCATCAATACGTGAGCGCCCCGTAATCCCCTGTCGGTAAGTGATTAACGGCCACTTCGCTATTGAt                                                              <  2:1671033/150‑1 (MQ=255)
                                                                                           agagCAAGATAGGTTTTATTGACATCAGAATCCTGCTCCCTTAATACAATATCAGCTAGCAAACCTTTGCGGGCAAATGCGTCATCAATAGGTGAGCGCCCCGTAATCCCCTGTCGGTAAGTGATTAACGGCCACTTCGCTATTGATTcc                                                           <  2:860956/150‑1 (MQ=255)
                                                                                           agagCAACATAGGTTTTAATGACATCAGAATCCTGCGCACTTAATACAATATCTGCCAGCAAACCTTTGCGGGCAAATGCGTCATCAATACGTGAGCGCCCCGTAATCCCCTGTCGGTAAGTGATTAACGGCCACTTCGCTATTGATTcc                                                           >  1:657526/1‑150 (MQ=255)
                                                                                           agagCAACATAGGTTTTAATGACATCAGAATCCTGCGCACTTAATACAATATCTGCCAGCAAACCTTTGCGGGCAAATGCGTCATCAATACGTGAGCGCCCCGTAATCCCCTGTCGGTAAGTGATTAACGGCCACTTCGCTATTGATTcc                                                           <  1:717081/150‑1 (MQ=255)
                                                                                           agagCAACATAGGTTTTAATGACATCAGAATCCTGCGCACTTAATACAATATCTGCCAGCAAACCTTTGCGGGCAAATGCGTCATCAATACGTGAGCGCCCCGTAATCCCCTGTCGGTAAGTGATTAACGGCCACTTCGCTATTGATTcc                                                           <  1:1647713/150‑1 (MQ=255)
                                                                                             agCAACATAGGTTTTAATGACATCAGAATCCTGCGCACTTAATACAATATCTGCCAGCAAACCTTTGCGGGCAAATGCGTCATCAATACGTGAGAGCCCCGTAATCCCCTGTCGGTAAGTGATTAACGGCCACTTCGCTATTGATTCCAg                                                         >  1:161781/1‑150 (MQ=255)
                                                                                               cAACATAGGTTTTAATGACATCAGAATCCTGCGCACTTAATACAATATCTGCCAGCAAACCTTTGCGGGCAAATGCGTCATCAATACGTGAGCGCCCCGTAATCCCCTGTCGGTAAGTGATTAACGGCCACTTCGCTATTGATTCCAGCg                                                       <  1:1373955/150‑1 (MQ=255)
                                                                                               cAACATAGGTTTTAATGACATCAGAATCCTGCGCACTTAATACAATATCTGCCAGCAAACCTTTGCGGGCAAATGCGTCATCAATACGTGAGCGCCCCGTAATCCCCTGTCGGTAAGTGATTAACGGCCACTTCGCTATTGATTCCAGCg                                                       <  1:824863/150‑1 (MQ=255)
                                                                                                 agaTAGGTTATAATGACATTTGAATCCTGCGCACTTAATACAATATCTGCCAGCAAACCTTTGCGGGCAAATGCGTCATCAATACGTGAGCGCCCCGTTATCCCCTGTCGGTAAGTGATTAACGGCCACTTCGCTATTGATTCCAGCGTc                                                     <  2:292812/148‑1 (MQ=255)
                                                                                                 aCATAGGTTTTAATGACATCAGAATCCTGCGCACTTAATACAATATCTGCCAGCAAACCTTTGCGGGCAAATGCGTCATCAATACGTGAGCGCCCCGTAATCCCCTGTCGGTAAGTGATTAACGGCCACTTCGCTATTGATTCCAGCGTc                                                     <  1:1324311/150‑1 (MQ=255)
                                                                                                     aGGTTTTAATGACATCAGAATCCTGCGCACTTAATACAATATCTGCCAGCAAACCTTTGCGGGCAAATGCGTCATCAATACGTGAGCGCCCCGTAATCCCCTGTCGGTAAGTGATTAACGGCCACTTCGCTATTGATTCCAGCGTCAATg                                                 >  2:1192515/1‑150 (MQ=255)
                                                                                                      ggTTTTAATGACATCAGAATCCTGCGCACTTAATACAATATCTGCCAGCAAACCTTTGCGGGCAAATGCGTCATCAATACGTGAGCGCCCCGTAATCCCCTGTCGGTAAGTGATTAACGGCCACTTCGCTATTGATTCCAGCGTCAATgg                                                <  2:1504634/150‑1 (MQ=255)
                                                                                                      ggTTTTAATGACATCAGAATCCTGCGCACTTAATACAATATCTGCCAGCAAACCTTTGCGGGCAAATGCGTCATCAATACGTGAGCGCCCCGTAATCCCCTGTCGGTAAGTGATTAACGGCCACTTCGCTATTGATTCCAGCGTCAATgg                                                <  1:1192515/150‑1 (MQ=255)
                                                                                                      ggTTTTAATGACATCAGAATCCTGCGCACTTAATACAATATCTGCCAGCAAACCTTTGCGGGCAAATGCGTCATCAATACGTGAGCGCCCCGTAATCCCCTGTCGGTAAGTGAATAACGGCCACTACGCTATTGATTCCAGCGTCAATgg                                                >  1:1686058/1‑150 (MQ=255)
                                                                                                      ggTTTTAATGACATCAGAATCCTGAGCACTTAATACAATATCTGCCAGCAAACCTATGCGTGCAAATGCGTCATCAAAACGTGAGCGCCCAGTAATCCCCTGTCGGTAAGTGATTAACGGCCACTACGATATTGATTCCAGCGTCAATgg                                                >  1:1050235/1‑150 (MQ=255)
                                                                                                       gTTTTAATGACATCAGAAACCTGCGCACTTAATACAATATCTGCCAGCAAACCTTTGCGGGCAAATGCCACATCAATACGTTAGCGCCCCGAAATCCCCTGTCGGTAAGTGATTAACGGCCACATCGCAATTGAATCAAGCGACAATGGt                                               >  1:1662916/1‑150 (MQ=255)
                                                                                                         tttAATGACATCAGAATCCTGCGCACTTAATACAATATCTGCCAGCAAACCTTTGCGGGCAAATGCGTCATCAATACGTGAGCGCCCCGTAATCCCCTGTCGGTAAGTGATTAACGGCCACTTCGCTATTGATTCCAGCGTCAATGGTGa                                             <  2:630422/150‑1 (MQ=255)
                                                                                                              tGACATCAGAATCCTGCGCACTTAATACAATATCTGCCAGCAAACCTTTGCGGGCAAATGCGTCATCAATACGTGAGCGCCCCGTAATCCCCTGTCGGTAAGTGATTAACGGCCACTTCGCTATTGATTCCAGCGTCAATGGTGAAAttt                                        >  2:1130423/1‑150 (MQ=255)
                                                                                                               taCATCAGAATCGTGCGCACTTATTACTATATCTGCCAGCAAACCTTTGAGGTCAATTGCGTCATCAATACGTGTTAGCCCAGTAATCCGCTGTCGGTAAGTGATTTACGGCCACTTCGCTATTGATTCCTGCGTCAATGGTGAAATTTg                                       <  1:46085/149‑1 (MQ=255)
                                                                                                                aCATCAGAATCCTGCGCACTTAATACAATATCTGCCAGCAAACCTTTGCGGGCAAATGCGTCATCAATACGTGAGCGCCCCGTAATCCCCTGTCGGTAAGTGATTAACGGCCACTTCGCTATTGATTCCAGCGTCAATGGTGAAATTTGc                                      >  2:801058/1‑150 (MQ=255)
                                                                                                                         tCCTGCGCACTTAATACAATATCTGCCAGCAAACCTTTGCGGGCAAATGCGTCATCAATACGTGAGCGCCCCGTAATCCCCTGTCGGTAAGTGATTAACGGCCACTTCGCTATTGATTCCAGCGTCAATGGTGAAATTTGCGTCAAGGga                             <  1:19335/150‑1 (MQ=255)
                                                                                                                           cTGCGCACTTAATACAATATCTGCCAGCAAACCTTTGCGGGCAAATGCGTCATCAATACGTGAGCGCCCCGTAATCCCCTGTCGGTAAGTGATTAACGGCCACTTCGCTATTGATTCCAGCGTCAATGGTGAAATTTGCGTCAAGGgatg                           >  2:129257/1‑150 (MQ=255)
                                                                                                                            tGCGCACTTAATACAATATCTGCCAGCAAACCTTTGCGGGCAAATGCGTCATCAATACGTGAGCGCCCCGTAATCCCCTGTCGGTAAGTGATTAACGGCCACTTCGCTATTGATTCCAGCGTCAATGGTGAAATTTGCGTCAAGGgatga                          <  2:1066649/150‑1 (MQ=255)
                                                                                                                               gcACTTAATACAATATCTGCCAGCAAACCTTTGCGGGCAAATGCGTCATCAATACGTGAGCGCCCCGTAATCCCCTGTCGGTAAGTGATTCACGGCCACTTCGCTATTGATTGCAGCGTCAAAGGTGAAATTTTCCTCCAGGGTAGATCg                       >  2:1085245/1‑150 (MQ=255)
                                                                                                                               gcACTTAATACAATATCTGCCAGCAAACCTTTGCGGGCAAATGCGTCATCAATACGTGAGCGCCCCGTAATCCCCTGTCGGTAAGTGATTAACGGCCACTTCGCTATTGATTCCAGCGTCAATGGTGAAATTTGCGTCAAGGGATGATCg                       <  1:1677625/150‑1 (MQ=255)
                                                                                                                               gcACTTAATACAATATCTGCCAGCAAACCTTTGCGGGCAAATGCGTCATCAATACGTGAGCGCCCCGTAATCCCCTGTCGGTAAGTGATTAACGGCCACTTCGCTATTGATTCCAGCGTCAATGGTGAAATTTGCGTCAAGGGATGATCg                       <  1:597905/150‑1 (MQ=255)
                                                                                                                               acACTTAATACAATATCTGCAAGCAAACCTTTGCGGGCAAATGCGTCATCAATACGTGAGCGCCCCGTAATCCCCTGTCGGTAAGTGATTAACGGCCACTTCGCTATTGATTCCAGCGTCAATGGTGAAATTTGCGTCAAGGGATGATCg                       <  2:876728/149‑1 (MQ=255)
                                                                                                                                 aCTTAATACAATATCTGCCCGCAAACCTTAGCGGGCAAATGCGTCATCAATACGTGAGCGCCTCGTAATCCCCTGTCGGTAAGTGATTAACGGCCACTTCGCTATTGATTCCAGCGTCAATGGTGAAATTTGCGTCAAGGGATGATCgtg                     <  1:273861/150‑1 (MQ=255)
                                                                                                                                  cTTAATACAATATCTGCCAGCAAACCTTTGCGGGCAAATGCGTCATCAATACGTGAGCGCCCCGTAATCCCCTGTCGGTAAGTGATTAACGGCCACTTCGCTATTGAGTCCAGCGTCAATGGTGAAATTTGCGTCAAGGGATGATCgtgt                    >  2:1248906/1‑150 (MQ=255)
                                                                                                                                    tAATACAATATCTGCCAGCAAACCTTTGCGGGCAAATGCGTCATCAATACGTGAGCGCCCCGTAATCCCCTGTCGGTAAGTGATTAACGGCCACTTCGCTATTGATTCCAGCGTCAATGGTGAAATTTGCGTCAAGGGATGATCGTGTgg                  <  1:487502/150‑1 (MQ=255)
                                                                                                                                       tACAATATCTGCCAGCAAACCTTTGCGGGCAAATGCGTCATCATTACGTGAGCGCCCCGTAATCCCCTGTCGGTAAGTGATTAACGGCCACTTCGCTATTGATTCCAGCGTCAATGGTGAAATTTGCGGCAAGGGATGATCGTGTGGAAc               <  1:1493894/150‑1 (MQ=255)
                                                                                                                                         cAATATCTGCCAGCAAACCTTTGCGGGCAAATGCGTCATCAATACGTGAGCGCCCCGTAATCCCCTGTCGGTAAGTGATTAACGGCCACTTCGCTATTGATTCCAGCGTCAATGGTGAAATTTGCGTCAAGGGATGATCGTGTGGAACaa             <  1:801058/150‑1 (MQ=255)
                                                                                                                                          cgcatttGCCCGCAAAGGTTTGCTGGCAAATGCGTCATCAATACGTGAGCGCCCCGTAATCCCCTGTCGGTAAGTGATTAACGGCCACTTCGCTATTGATTCCAGCGTCAATGGTGAAATTTGCGTCAAGGGATGATCGTGTGGAACAAg            >  1:1714677/7‑150 (MQ=255)
                                                                                                                                          aaTATCTGCCAGCAAACCTTTGCGGGCAAATGCGTCATCAATACGTGAGCGCCCCGTAATCCCCTGTCGGTAAGTGATTAACGGCCACTTCGCTATTGATTCCAGCGTCAATGGTGAAATTTGCGTCAAGGGATGATCGTGTGGACCAAg            >  2:1171733/1‑150 (MQ=255)
                                                                                                                                           atatCTGCCAGCAAACCTTTGCGGGCAAATGCGTCATCAATACGTGAGCGCCCCGTAATCCCCTGTCGGTAAGTGATTAACGGCCACTTCGCTATTGATTCCAGCGTCAATGGTGAAATTTGCGTCAAGGGATGATCGTGTGGAACAAGc           >  2:961259/1‑150 (MQ=255)
                                                                                                                                             atttGCCCGCAAAGGTTTGCTGGCAAATGCGTCATCAATACGTGAGCGCCCCGTATCCCCTGTCGGTAAGTGATTAACGGCCACTTCGCTATTGATTCCAGCGTCAATGGTGAAATTTGCGTCAAGGGATGATCGTGTGGAACAAGCaaa        >  1:561569/4‑150 (MQ=255)
                                                                                                                                              tCTGCCAGCAAACCTTTGCGGGCAAATGCGTCATCAATACGTGAGCGCCCCGTAATCCCCTGTCGGTAAGTGATTAACGGCCACTTCGCTATTGATTCCAGCGTCAATGGTGAAATTTGCGTCAAGGGATGATCGTGTGGAACAAGCaaa        >  2:507811/1‑150 (MQ=255)
                                                                                                                                              tCTGCCAGCAAACCTTTGCGGGCAAATGCGTCATCAATACGTGAGCGCCCCGTAATCCCCGGTCGGTAAGTGATTAACGGCCACTTCGCTATTGATTCCAGCGTCAATGGTGAAATTTGCGTCAAGGGATGATCGTGTGGCACAAGCaaa        >  2:664588/1‑150 (MQ=255)
                                                                                                                                               cTGCCAGCAAACCTTTGCGGGCAAATGCGTCATCAATACGTGAGCGCCCCGTAATCCCCTGTCGGTAAGTGATTAACGGCCACATCGCTATTGATTCCAGCGTCAATGGTGAAATATTCGTCAAGGTATGATAGTGTGGAACAAGCAAAc       >  2:320559/1‑150 (MQ=255)
                                                                                                                                               cTGCCAGCAAACCTTTGCGGACAAATGCGTCATCAATACGTGAGCGCCCCGTAATCCCCTGTCGGTAAGTGATTAACGGCCACTTCGCTATTGATTCCAGCGTCAATGGTGAAATTTGCGTCAAGGGATGATCGTGTGGAACAAGCAAAc       >  2:988093/1‑150 (MQ=255)
                                                                                                                                                    aGCAAACCTTTGCGGGCAAATGCGTCATCAATACGTGAGCGCCCCGTAATCCCCTGTCGGTAAGTGATTAACGGCCACTTCGCTATTGATTCCAGCGTCAATGGTGAAATTTGCGTCAAGGGATGATCGTGTGGAACAAGCAAACTATgg  >  1:974778/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
ATCAAAAAGATGCCGCGTATCCAGGCGGATTAAATTCTCTTCCTCTTGTTCGCCACTGGATTGCTCGGCAACTAATCCGATCCCAAGCCCAAGAGCAACATAGGTTTTAATGACATCAGAATCCTGCGCACTTAATACAATATCTGCCAGCAAACCTTTGCGGGCAAATGCGTCATCAATACGTGAGCGCCCCGTAATCCCCTGTCGGTAAGTGATTAACGGCCACTTCGCTATTGATTCCAGCGTCAATGGTGAAATTTGCGTCAAGGGATGATCGTGTGGAACAAGCAAACTATGG  >  NC_000913/2060117‑2060414

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: