Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,708,595 G→A 5.4% intergenic (‑890/+21) dppA ← / ← proK dipeptide/heme ABC transporter periplasmic binding protein; dipeptide chemotaxis receptor/tRNA‑Pro

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,708,5950GA 5.4% 348.1 / 6.0 112intergenic (‑890/+21)dppA/proKdipeptide/heme ABC transporter periplasmic binding protein; dipeptide chemotaxis receptor/tRNA‑Pro
Reads supporting (aligned to +/- strand):  ref base G (56/49);  new base A (3/3);  total (59/53)
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.02e-01

Only 100 of 114 total aligned reads displayed.

GTGATAGAGGATTACTCGCCACTTCGTGGCTCGCCCTGCGGGCCGTTGCTGGCGCAACGTTCTCTCGCTTCGCTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCAGGCTGCGCCAATCACCGAATGCGGGGCGCATCTTACTGCGCAGATACGCCCTCGTCAATCCCTTAATAG  >  NC_000913/3708454‑3708744
                                                                                                                                             |                                                                                                                                                     
gTGATTGAGGATTATTGGACAATTCGTGGCTCGAGCTGCGTTACGTTGGTGGCGCATAGTTATCTGTCTTCTGTCGAGTCTAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAg                                                                                                                                               <  2:1733385/150‑1 (MQ=14)
 tGATAGAGGATTACTCGCCACTTCGTGGCTCGCCCTGCGGGCCGTTGCTGGCGCAACGTTCTCTCGCTTCGCTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAgg                                                                                                                                              >  2:1709526/1‑150 (MQ=255)
   aTAGATGATTACTCGCCACTTCGTGGCTCTCCCTGCGGGCCGTGGCTGGCGCAACGTTCTCTCGCTTCGCTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCAACGAGCAGGCt                                                                                                                                            <  1:751967/150‑1 (MQ=25)
     agagGATTACTCGCCACTTCGTGGCTCGCCCTGCGGGCCGTTGCTGGCGCAACGTTCTCTCGCTTCGCTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCttt                                                                                                                                          <  2:1361643/150‑1 (MQ=255)
     acagGTTTACTCGCCACTACGTGGCTCGCCCTGCGGTCCGTTGCTGGCGCAACGTTCTCTCGCTTCGCTCTAGTCCAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCAACGAGCAGGCttt                                                                                                                                          <  1:968620/148‑1 (MQ=25)
      gagGATTACTCGCCACTTCGTGGCTCGCCCTGCGGGCCGTTGCTGGCGCAACGTTCTCTCGCTTCGCTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAATCCTGCGCGGTGAGCAGGCtttt                                                                                                                                         >  2:1257865/1‑150 (MQ=255)
      gagGATTACTCGCCACTTCGTGGCTCGCCCTGCGGGCCGTTGCTGGCGCAACGTTCTCTCGCTTCGCTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCtttt                                                                                                                                         >  2:427301/1‑150 (MQ=255)
         taTTACTCGCCACTTCGTGGCTCGCTCGGCGGGCGGGTGCTGGCGCATCGTTCTCTCGCTTCGCTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCAACGAGCAGGCttttttg                                                                                                                                      <  2:1090023/149‑4 (MQ=21)
         gATTACTCGCCACTTCGTGGCTCGCCCTGCGGGCCGTTGCTGGCGCAACGTTCTCTCGCTTCGCTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCGCGTTGAGCAGGCTTTTCGa                                                                                                                                      >  2:609417/1‑150 (MQ=255)
             aCTCGCCACTTCGTGGCTCGCCCGGCGGGCCGGTGCTGGCGCAACGTTCTCTCGCTTCGCTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAAttt                                                                                                                                  <  1:1432129/150‑1 (MQ=255)
               tCGCCACTTCGTGGCTCGCGCTGCGGGCCGTTGCTGGAGCAACGTTCTCTCGCTTCGCTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTgg                                                                                                                                <  2:757526/150‑1 (MQ=255)
                  ccACTTCGTGGCTCGCCCTGCGGGCCGTTGCTGGAGCAACGTTCTCTCGCTTCGCTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCg                                                                                                                             <  1:1026626/150‑1 (MQ=255)
                     cTTCGTGGCTCGCCCTGCGGGCCGTTGCTGGCGCAACGTTCTCTCGCTTCGCTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTg                                                                                                                          <  2:843497/150‑1 (MQ=255)
                     cTTCGTGGCTCGCCCTGCGGGCCGTTGCTGGCGCAACGTTCTCTCGCTTCGCTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTg                                                                                                                          >  2:336008/1‑150 (MQ=255)
                           ggCTCGCCCTGCGGGCCGTTGCTGGCGCAACGTTCTCTCGCTTCGCTCGAGTCGACCCTCCCTGCCCCCGGAGGTTGGCGTCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGGGATagag                                                                                                                    >  1:1067853/1‑150 (MQ=21)
                            gCTCGCCCTGCGGGCCGTTGCTGGCGCAACGTTCTCTCGCTTCGCTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAgg                                                                                                                   >  1:906638/1‑150 (MQ=255)
                            gCTCGCCCTGCGGGCCGTTGCTGGCGCAACGTTCTCTCGCTTCGCTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAgg                                                                                                                   >  1:1001327/1‑150 (MQ=255)
                                gCCCTGCGGGCCGTTGCTGGCGCAACGTTCTCTCGCTTCGCTCGAGTCGAACCTCCCTGCCCCCGGCGGTTCTCATCCGCATGAGCTGCAGCTGCAAAAAAGCCTGCTCGGTGAGCAGGCTGTTCGAATTTGGTCGGTGATAGAGGAttg                                                                                                               >  1:385363/1‑149 (MQ=18)
                                gCCCTGCGGGCCGTTGCTGGCGCAACGTTCTCTCGCTTCGCTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTc                                                                                                               >  1:397705/1‑150 (MQ=255)
                                 cccTGCGGGCCGTTGCTGGCGCAACGTTCTCTCGCTTCGCTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCg                                                                                                              >  1:629843/1‑150 (MQ=255)
                                 cccTGCGGGCCGTTGCTGGCGCAACGTTCTCTCGCTTCGCTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCg                                                                                                              >  2:456294/1‑150 (MQ=255)
                                 cacTGCGGGCCGTTGCTGGCGGATCGTTATCTCGCTTAGCTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCg                                                                                                              <  2:1562443/148‑1 (MQ=255)
                                    ttcGGGCCGTTGATGGCGCAACGTTCTATCTCTTCGCGCGTGTCGATCATACCTGCCCCCTGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCGAAc                                                                                                           <  1:1598926/148‑1 (MQ=17)
                                       gggCCGTTGCTGGCGCAACGTTCTCTCGCTTCGCTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCGAACCTc                                                                                                        >  1:440693/1‑150 (MQ=255)
                                        ggCCGTTGCTGGCGCAACGTTCTCTCGCTTCGCTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCGAACCTcc                                                                                                       >  1:225333/1‑150 (MQ=255)
                                        ggCCGTTGCTGGCGCAACGTTCTCTCGCGTCGCTCGCGTCGAACCTCCCTGCCCCCGGAGGTTCTCCTCCTGATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGGGGCGGCTTTTCGAATTTGGTGGGTGATTGAGGATTCGAACCTcc                                                                                                       >  1:377022/1‑150 (MQ=18)
                                         gCCGTTGCTGGCGCAACGTTCTCTCGCTTCGCTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAATCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTAGAACCTccc                                                                                                      >  2:1498084/1‑149 (MQ=255)
                                         gCCGTTGCTGGCGCAACGTTCTCTCGCTTCGCTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGGTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCg                                                                                                      >  1:1591543/1‑150 (MQ=255)
                                           ggTTGCTGGCGCAACGTTCTCTCGCTTCGCTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGAc                                                                                                    <  1:732697/149‑1 (MQ=255)
                                              tGCTGGCGAAACGTTCTCTCGCTTCGCTCGAGTCGAACATCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCTGATGCAAAAAAGCCGGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGAcccc                                                                                                 <  2:82219/150‑1 (MQ=255)
                                                    cgcAACGTTCTCTCGCTTCGCTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTc                                                                                           <  2:861542/150‑1 (MQ=255)
                                                     gcAACGTTCTCTCGCTTCTCTCTAGTCGAACCTCCCTGCCCCCTGAGGTTCTCATCCTCATTAGCTTCAGATCCAAAAAAGCCTGCTCTTTTAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTcc                                                                                          <  2:1202884/150‑1 (MQ=18)
                                                     gcAACGTTCTCTCGCTTCGCTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGGTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTcc                                                                                          <  1:609417/150‑1 (MQ=255)
                                                      cAACGTTCTCTCGCTTCGCTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGGTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTccc                                                                                         >  1:366224/1‑150 (MQ=255)
                                                         cGTTCTCTCGCTTCGCTCGAGTCTAACATACCTGCCCCCGGAGGTTCTCATCCTCATTAGCTTCAGATGCAAAAAAGCCTTCTCGTTGATCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCcgaa                                                                                      <  2:1237257/150‑1 (MQ=18)
                                                           ttCTCTCGCTTCGCTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCCTCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATGTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTGCGTCCggatcg                                                                                    >  1:276833/1‑150 (MQ=17)
                                                              tctcGCTTCGCTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAg                                                                                 >  1:65864/1‑150 (MQ=255)
                                                                 cgcttcgctCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTgc                                                                              <  1:1199588/150‑1 (MQ=255)
                                                                 cgcttcgctCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTgc                                                                              <  1:1377086/150‑1 (MQ=255)
                                                                   cttcgctCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGACTTTGGTCGGTGATAGAGGATTCGAACCTCCGCCCCCTTCGTCCCGAACGAAGTgcgc                                                                            >  1:271140/1‑150 (MQ=255)
                                                                        ctCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCa                                                                       >  1:519666/1‑150 (MQ=255)
                                                                        ctCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCGTTTCGAATTTGGTCGGTGATTGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCGAccc                                                                       >  1:363422/1‑149 (MQ=255)
                                                                         tCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCAg                                                                      <  2:1368719/150‑1 (MQ=255)
                                                                           gAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCAGGc                                                                    >  1:1122237/1‑150 (MQ=255)
                                                                           gAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCAGGc                                                                    <  2:1140362/150‑1 (MQ=255)
                                                                            aGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCAGGCt                                                                   >  2:1648252/1‑150 (MQ=255)
                                                                             gTCTAAGCTCCCTGCCCCCGTAGGTTCTCATCCTCTTGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGCGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCGACCCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCAGGCTg                                                                  <  1:1132168/150‑1 (MQ=18)
                                                                             gTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGCCCCCTTCGTCCCGAACGAAGTGCGCTACCAGGCTg                                                                  >  1:1023543/1‑150 (MQ=255)
                                                                               cGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGCATTCGAACCTCCGACCCCTACGTCCCGAACGAAGTGCGCTACCAGGCTgcg                                                                >  2:1244108/1‑150 (MQ=255)
                                                                                gAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCCGGCTgcgc                                                               >  1:844936/1‑150 (MQ=255)
                                                                                 aaCCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCAGGCTGCGcc                                                              <  1:945593/150‑1 (MQ=255)
                                                                                 aaCCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCAGGCTGCGcc                                                              <  1:1648252/150‑1 (MQ=255)
                                                                                  aCCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCAGGCTGCGCCa                                                             >  2:652766/1‑150 (MQ=255)
                                                                                   ccTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCAGGCTGCGCCaa                                                            >  1:1581721/1‑150 (MQ=255)
                                                                                     gcccTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCAGGCTGCGCCAATc                                                          <  1:427301/149‑1 (MQ=255)
                                                                                      cccTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATTGAGGGTTCGAACCTCCGACCCCTTCGTCCCGAACGCAGTGCGCTACCAGGCTGCGCCAATcc                                                         >  2:1103411/1‑149 (MQ=255)
                                                                                        cTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGACCGAAGTGCGCTACCAGGCTGCGCCAATCAcc                                                       >  2:59916/1‑150 (MQ=255)
                                                                                         tGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCAGGCTGCGCCAATCACCg                                                      <  2:1122381/150‑1 (MQ=255)
                                                                                         tGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCAGGCTGCGCCAATCACCg                                                      <  2:1128237/150‑1 (MQ=255)
                                                                                                ggAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCAGGCTGCGCCAATCACCGAATGCgg                                               >  2:1126198/1‑150 (MQ=255)
                                                                                                 gAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCAGGCTGCGCCAATCACCGAATGCggg                                              >  2:1167756/1‑150 (MQ=255)
                                                                                                 gAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTAGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCAGGCTGCGCCAATCACCGAATGCggg                                              <  1:642097/150‑1 (MQ=255)
                                                                                                        ccATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCAGGCTGCGCCAATCACCGAATGCGGGGCGCATc                                       <  1:1177352/149‑1 (MQ=255)
                                                                                                         cATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCAGGCTGCGCCAATCACCGAATGCGGGGCGCATCt                                      <  1:1212989/150‑1 (MQ=255)
                                                                                                         cATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCAGGCTGCGCCAATCACCGAATGCGGGGCGCATCt                                      >  2:695235/1‑150 (MQ=255)
                                                                                                          aTCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCAGGCTGCGCCAATCACCGAATGCGGGGCGCATCtt                                     <  1:83654/150‑1 (MQ=255)
                                                                                                          aTCCTCATGAGCGGCAGATGCAAAAAAGCCGGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCAGGCTGCGCCAATCACCGAATGCGGGGCGCATCtt                                     <  2:1320564/150‑1 (MQ=255)
                                                                                                           tCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCAGGCTGCGCCAATCACCGAATGCGGGGCGCATCTTa                                    >  2:1716625/1‑150 (MQ=255)
                                                                                                             cTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCAGGCTGCGCCAATCACCGAATGCGGGGCGCATCTTACt                                  >  2:1711282/1‑150 (MQ=255)
                                                                                                             cTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCAGGCTGCGCCAATCACCGAATGCGGGGCGCATCTTACt                                  >  2:1028449/1‑150 (MQ=255)
                                                                                                              tCATTAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCAGGCTGCGCCAATCACCGAATGCGGGGCGCATCTTACTg                                 <  2:234436/150‑1 (MQ=255)
                                                                                                                 tGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGGGATAGAGGATTCGATCCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCCGGCTGCGCCAATCACCGAATGCTGGGCGCATCTTACTgcgc                              >  1:947912/1‑150 (MQ=255)
                                                                                                                   aGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCAGGCTGCGCCAATCACCGAATGCGGGGCGCATCTTACTGCGCAg                            >  2:346048/1‑150 (MQ=255)
                                                                                                                     cTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCAGGCTGCGCCAATCACCGAATGCGGGGCGCATCTTACTGCGCAGAt                          >  1:1019502/1‑150 (MQ=255)
                                                                                                                      tGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCAGGCTGCGCCAATCACCGAATGCGGGGCGCATCTTACTGCGCAGATa                         <  1:142814/150‑1 (MQ=255)
                                                                                                                          gATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTTGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCAGGCTGCGCCAATCACCGAATGCGGGGCGCATCTTACTGCGCAGATACGcc                     <  1:1069422/150‑1 (MQ=255)
                                                                                                                            tGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCAGGCTGCGCCAATCACCGAATGCGGGGCGCATCTTACTGCGCAGATACGCCCt                   <  2:1067853/150‑1 (MQ=255)
                                                                                                                             gCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCAGGCTGCGCCAATCACCGAATGCGGGGCGCATCTTACTGCGCAGATACGCCCTc                  <  2:50063/150‑1 (MQ=255)
                                                                                                                              cAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCAGGCTGCGCCAATCACCGAATGCGGGGCGCATCTTACTGCGCAGATACGCCCTCg                 <  1:230988/150‑1 (MQ=255)
                                                                                                                              cAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCAGGCTGCGCCAATCACCGAATGCGGGGCGCATCTTACTGCGCAGATACGCCCTCg                 <  2:440693/150‑1 (MQ=255)
                                                                                                                               cgaaaaGCCTGCTCAACGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCAGGCTGCGCCAATCACCGAATGCGGGGCGCATCTTACTGCGCAGATACGCCCTCgg                >  2:344945/3‑149 (MQ=255)
                                                                                                                               cgaaaaGCCTGCTCAAAGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCAGGCTGCGCCAATCACCGAATGCGGGGCGCATCTTACTGCGCAGATACGCCCTCGt                >  2:790130/3‑150 (MQ=255)
                                                                                                                               aaaaaaGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTAGGTGATAGAGGATTCGAAACTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCAGGCTGCGCCAATCACCGAATGCGGGGCGCATCTTACTGCGCAGATACGCCCTCGt                >  2:376371/1‑150 (MQ=255)
                                                                                                                                aaaaaGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCAGGCTGCGCCAATCACCGAATGCGGGGCGCATCTTACTGCGCAGATACGCCCTCGTc               >  2:5799/1‑150 (MQ=255)
                                                                                                                                aaaaaGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCAGGCTGCGCCAATCACCGAATGCGGGGCGCATCTTACTGCGCAGATACGCCCTCGTc               <  1:1709526/150‑1 (MQ=255)
                                                                                                                                 aaaaGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCAGGCTGCGCCAATCACCGAATGCGGGGCGCATCTTACTGCGCAGATACGCCCTCGTCa              <  2:1591543/150‑1 (MQ=255)
                                                                                                                                  aaaGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCAGGCTGCGCCAATCACCGAATGCGGGGCGCATCTTACTGCGCAGATACGCCCTCGTCaa             >  1:1507775/1‑150 (MQ=255)
                                                                                                                                     gCCTGCTCAACGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCAGGCTGCGCCAATCACCGAATGCGGGGCGCATCTTACTGCGCAGATACGCCCTCGTCAATcc          >  1:851146/1‑150 (MQ=255)
                                                                                                                                      ccTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTcc                                                                                                              >  1:137157/1‑48 (MQ=17)
                                                                                                                                      ccTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCg                                                                                                              <  2:137157/49‑1 (MQ=35)
                                                                                                                                      ccTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCAGGCTGCGCCAATCACCGAATGCGGGGCGCATCTTACTGCGCAGATACGCCCTCGTCAATccc         >  2:1542954/1‑150 (MQ=255)
                                                                                                                                      ccTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCAGGCTGCGCCAATCACCGAATGCGGGGCGCATCTTACTGCGCAGATACGCCCTCGTCAATccc         >  1:497800/1‑150 (MQ=255)
                                                                                                                                         gCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCAGGCTGCGCCAATCACCGAATGCGGGGCGCATCTTACTGCGCAGATACGCCCTCGTCAATCCCTTa      <  2:1045237/150‑1 (MQ=255)
                                                                                                                                           tCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCAGGCTGCGCCAATCACCGAATGCGGGGCGCATCTTACTGCGCAGATACGCCCTCGTCAATCCCTTAAt    <  1:1599663/150‑1 (MQ=255)
                                                                                                                                           tCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCAGGCTGCGCCAATCACCGAATGCGGGGCGCATCTTACTGCGCAGATACGCCCTCGTCAATCCCTTAAt    <  1:5799/150‑1 (MQ=255)
                                                                                                                                           tCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCAGGCTGCGCCAATCACCGAATGCGGGGCGCATCTTACTGCGCAGATACGCCCTCGTCAATCCCTTAAt    <  2:1480724/150‑1 (MQ=255)
                                                                                                                                           tCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCAGGCTGCGCCAATCACCGAATGCGGGGCGCATCTTACTGCGCAGATACGCCCTCGTCAATCCCTTAAt    <  2:371979/150‑1 (MQ=255)
                                                                                                                                            cGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCAGGCTGCGCCAATCACCGAATGCGGGGCGCATCTTACTGCGCAGATACGCCCTCGTCAATCCCTTAATa   >  1:670389/1‑150 (MQ=255)
                                                                                                                                             gTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCAGGCTGCGCCAATCACCGAATGCGGGGCGCATCTTACTGCGCAGATACGCCCTCGTCAATCCCTTAATAg  >  2:1353958/1‑150 (MQ=255)
                                                                                                                                             gTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCAGGCTGCGCCAATCACCGAATGCGGGGCGCATCTTACTGCGCAGATACGCCCTCGTCAATCCCTTAATAg  >  2:23403/1‑150 (MQ=255)
                                                                                                                                             |                                                                                                                                                     
GTGATAGAGGATTACTCGCCACTTCGTGGCTCGCCCTGCGGGCCGTTGCTGGCGCAACGTTCTCTCGCTTCGCTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCAGGCTGCGCCAATCACCGAATGCGGGGCGCATCTTACTGCGCAGATACGCCCTCGTCAATCCCTTAATAG  >  NC_000913/3708454‑3708744

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: