Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,068,924 G→A 5.0% intergenic (‑116/+23) argO ← / ← mscS arginine transporter/mechanosensitive channel protein, small conductance

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,068,9240GA 5.0% 382.8 / 5.4 121intergenic (‑116/+23)argO/mscSarginine transporter/mechanosensitive channel protein, small conductance
Reads supporting (aligned to +/- strand):  ref base G (59/55);  new base A (2/4);  total (62/59)
Fisher's exact test for biased strand distribution p-value = 4.35e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.29e-01

Only 100 of 126 total aligned reads displayed.

CAAGTGCAAGACCTTGAAAGTAATAAGAAAACACGTTGGATATTCCGAATTAAGTATCTGATGCAGAGACTATATCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACGGGAAGCTGATACCGGCGGCATCAAATTCACGTTTAATACGCTCCAGCACATCCCAGTACACGTTTTGCAG  >  NC_000913/3068775‑3069072
                                                                                                                                                     |                                                                                                                                                    
cAAGTGCAAGACCTTGAAAGTAATAAGAAAACACGTTGGATATTCCGAATTAAGTATCTGATGCAGAGACTATATCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCagagag                                                                                                                                                      >  2:681210/1‑150 (MQ=255)
aaaGTGCAAGACCTTGAAAGTAATAAGAAAACACGTGGGATATTCCTAATTAAGTCTCTGATGCAGAGACTATATCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCagagag                                                                                                                                                      <  2:70362/149‑1 (MQ=255)
   gTGCAAGACCTTGAAAGTAATAAGAAAACACGTTGGATATTCCGAATTAAGTATCTGATGCAGAGACTATATCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGtat                                                                                                                                                   >  1:697118/1‑150 (MQ=255)
    tGCAAGACCTTGAAAGTAATAAGAAAACACGTTGGATATTCCGAATTAAGTATCTGATGCAGAGACTATATCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGATATAGGCCAAAGAGtatt                                                                                                                                                  >  1:996481/1‑150 (MQ=255)
     gCAAGACCTTGAAAGTAATAAGAAAACACGTTGGATATTCCGAATTAAGAATCAGATGCAGAGACTATATCGCGGGATAGTTATAAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGACGCCAGAGAGtatta                                                                                                                                                 <  2:320406/150‑1 (MQ=255)
      cAAGACCTTGAAAGTAATAAGAAAACACGTTGGATATTCCGAATTAAGTATCTGATGCAGAGACTATATCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGtattat                                                                                                                                                >  2:1031001/1‑150 (MQ=255)
            cTTGAAAGTAATAAGAAAACACGTTGGATATTCCGAATTAAGTATCTGATGCAGAGACTATATCGCGTGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGcc                                                                                                                                          <  1:1157868/150‑1 (MQ=255)
            cTTGAAAGTAATAAGAAAACACGTTGGATATTCCGAATTAAGTATCTGATGCAGAGACTATATCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGATAATACTCTCTGGcc                                                                                                                                          <  1:24538/150‑1 (MQ=255)
             ttGAAAGTAATAAGAAAACACGTTGGATATTCCGAATTAAGTATCTGATGCAGAGACTATATCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCt                                                                                                                                         <  1:1094404/150‑1 (MQ=255)
              tGAAAGTAATAAGAAAACACGTTGGATATTCCGAATTAAGTATCTGATGCAGAGACTATATCGCGGGATTGTTATTAGCGGAATTTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTc                                                                                                                                        <  2:1714937/150‑1 (MQ=255)
                aaaGTAATAAGAAAACACGTTGGATATTCCGAATTAAGTATCTGATGCAGAGACTATATCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAg                                                                                                                                      <  1:550714/150‑1 (MQ=255)
                  aGTAATAAGAAAACACGTTGGATATTCCGAATTAAGTATCTGATGCAGAGACTATATCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGATAATACTCTCTGGCCTctctc                                                                                                                                    <  1:1702830/150‑5 (MQ=255)
                    taataaGAAAACACGTTGGATATTCCGAATTAAGTATCTGATGCAGAGACTATATCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGtt                                                                                                                                  >  2:451905/1‑150 (MQ=255)
                    taataaCAAAACACGTTGGATATTCCGAATTAAGTATCTGATGCAGAGACTATATCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGtt                                                                                                                                  <  2:1746396/150‑1 (MQ=255)
                     aataaGAAAACACGTTGGATATTCCGAATTAAGTATCTGATGCAGAGACTATATCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTg                                                                                                                                 <  2:1455337/150‑1 (MQ=255)
                       taaGAAAACACGTTGGATATTCCGAATTAAGTATCTGATGCAGAGACTATATCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGAt                                                                                                                               >  1:1552019/1‑150 (MQ=255)
                        aaGAAAACACGTTGGATATTCCGAATTAAGTATCTGATGCAGAGACTATATCTCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGAtt                                                                                                                              <  1:499778/150‑1 (MQ=255)
                        aaGAAAACACGTTGGATATTCCGAATTAAGTATCTGATGCAGAGACTATATCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGAtt                                                                                                                              <  2:1276646/150‑1 (MQ=255)
                        aaGAAAACACGTTGGATATTCCGAATTAAGTATCTGATGCAGAGACTATATCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGAtt                                                                                                                              >  2:1903631/1‑150 (MQ=255)
                          gAAAACACGTTGGATATTCCGAATTAAGTATCTGATGCAGAGACTATATCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTTGAAATCAGAAAAACTAATAACAGAGAGGCCAGATAGTATTATCTGGCCTCAGCGTTCATTAc                                                                                                                            >  2:763712/1‑150 (MQ=255)
                                acGTTGGATATTCCGAATTAAGTATCTGATGCAGAGACTATATCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCt                                                                                                                      <  2:569403/150‑1 (MQ=255)
                                   ttGGATATTCCGAATTAAGTATCTGATGCAGAGACTATATCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTg                                                                                                                   <  1:80060/150‑1 (MQ=255)
                                   ttGGATATTCCGAATTAAGTATCTGATGCAGAGACTATATCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTg                                                                                                                   <  1:1066916/150‑1 (MQ=255)
                                    tGGATATTCCGAATTAAGTATCTGATGCAGAGACTATATCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGt                                                                                                                  >  1:1254390/1‑150 (MQ=255)
                                          ttCCGAATTAAGTATCTGATGCAGAGACTATATCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTcttctt                                                                                                            <  2:1339085/150‑1 (MQ=255)
                                          ttCCGAATTAAGTATCTGATGCAGAGACTATATCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTcttctt                                                                                                            >  2:1157049/1‑150 (MQ=255)
                                          ttCCGAATTAAGTATCTGATGCAGAGACTATATCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTcttctt                                                                                                            <  1:398870/150‑1 (MQ=255)
                                           tcagAATTAAGTATCTGATGCAGAGACTATATCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCttt                                                                                                           >  2:400171/4‑150 (MQ=255)
                                           tCCGAATTAAGTATCTGATGCAGAGACTATATCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCAAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCttt                                                                                                           >  1:1256399/1‑150 (MQ=255)
                                            ccGAATTAAGTATCTGATGCAGAGACTATATCGCGGGATTGTTATTAGCGGAAATTGTTAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTc                                                                                                          >  2:1389511/1‑150 (MQ=255)
                                                  tAAGTATCTGATGCAGAGACTATATCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGc                                                                                                    >  1:1340791/1‑150 (MQ=255)
                                                   aaGTATCTGATGCAGAGACTATATCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCt                                                                                                   <  2:180443/150‑1 (MQ=255)
                                                   aaGTATCTGATGCAGAGACTATATCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCt                                                                                                   >  1:444305/1‑150 (MQ=255)
                                                   aaGTATCTGATGCAGAGACTATATCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCt                                                                                                   <  1:1809925/150‑1 (MQ=255)
                                                   aaGTATCTGATGCAGAGACTATATCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCt                                                                                                   <  2:286456/150‑1 (MQ=255)
                                                   aaGTATCTGATGCAGAGACTATATCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACACGCt                                                                                                   >  1:434736/1‑150 (MQ=255)
                                                     gTATCTGATGCAGAGACTATATCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGATGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTa                                                                                                 >  1:1341136/1‑150 (MQ=255)
                                                     gTATCTGATGCAGAGACTATATCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTa                                                                                                 >  2:1094623/1‑150 (MQ=255)
                                                         cTGATGCAGAGACTATATCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGt                                                                                             >  1:287595/1‑150 (MQ=255)
                                                           gATGCAGAGACTATATCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTc                                                                                           >  2:1585980/1‑150 (MQ=255)
                                                             tGCAGAGACTATATCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTcac                                                                                         >  1:43463/1‑150 (MQ=255)
                                                             tGCAGAGACTATATCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTcac                                                                                         <  2:1378602/150‑1 (MQ=255)
                                                              gCAGAGACTATATCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTcaca                                                                                        <  1:1820553/150‑1 (MQ=255)
                                                               cAGAGACTATATCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACAt                                                                                       >  1:1041112/1‑150 (MQ=255)
                                                               cAGAGACTATATCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTGGAGTCCGACGCGTTGTCTTCTTTCACCAGCTGAAAGTTCACAt                                                                                       >  2:955680/1‑150 (MQ=255)
                                                                  agaCTATATCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCa                                                                                    >  1:1031019/1‑150 (MQ=255)
                                                                      tataTCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTg                                                                                <  2:264095/150‑1 (MQ=255)
                                                                        tatCGCTGGATTGTTATTAGAGGAAATTGATAATATTAATTGGCCATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCg                                                                              <  2:176824/150‑1 (MQ=255)
                                                                          tCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCAAGAGAGTATTATCTGGCCTCAGCGATGATTACTCAGCTTTGTCTTCTTTCACCCGATTAAAGTTCACATCCATTTGCggg                                                                            >  1:1176570/1‑150 (MQ=255)
                                                                           cgcgGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGt                                                                           >  2:386239/1‑150 (MQ=255)
                                                                           cgcgGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGt                                                                           <  2:1226641/150‑1 (MQ=255)
                                                                             cgGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTAc                                                                         <  1:1031001/150‑1 (MQ=255)
                                                                             cgGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTAc                                                                         <  2:74322/150‑1 (MQ=255)
                                                                                gATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACggg                                                                      <  1:1266826/150‑1 (MQ=255)
                                                                                    gTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTAAGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACGGGAAGc                                                                  >  1:481316/1‑150 (MQ=255)
                                                                                    gTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTATAAGAGAAAAGCCAGAGAGTAATATCAGGCCTACGGCTTGATTTCGCAGCTTTGTCTTATTTAACCCGCTTAAAGTTAAAATACATTTTCGGGTTCGGGAAACt                                                                 >  2:1743473/1‑150 (MQ=255)
                                                                                      tattaGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACGGGAAGCTg                                                                >  1:1627339/1‑150 (MQ=255)
                                                                                         taGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACGGGAAGCTGATa                                                             >  1:1828647/1‑150 (MQ=255)
                                                                                         taGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACGGGAAGCTGATa                                                             <  1:16736/150‑1 (MQ=255)
                                                                                           gCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTATTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACGGGAAGCTGAAAcc                                                           >  1:861907/1‑150 (MQ=255)
                                                                                                    gATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACGGGAAGCTGATACCGGCGGCAta                                                  >  1:1814187/1‑149 (MQ=255)
                                                                                                       aataTTAATTGGCAATCATAAAAACTAATAAGAGAGAGGCCATAGAGTATTATCTGGCCTCAGCGTTGATAAAGCAGCTTTGTCTGCTTTCAACCGCTTAAAGTTCACATCCATTTTCGGTAACTGGAAGCTGATACCCTCGGCATCaaa                                               >  1:170609/1‑150 (MQ=255)
                                                                                                       aataTTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTGCGGGGAGCTGATGCCGGCGGCAtgagg                                               >  2:1778222/1‑146 (MQ=255)
                                                                                                         tatTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACGGGAAGCTGATACCGGCGGCATCAAAtt                                             <  1:1684801/150‑1 (MQ=255)
                                                                                                         tatTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACGGGAAGCTGATACCGGCGGCATCAAAtt                                             >  2:1560880/1‑150 (MQ=255)
                                                                                                           ttAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACGGGAAGCTGATACCGGCGGCATCAAATTCa                                           >  1:1623305/1‑150 (MQ=255)
                                                                                                            tAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACGGGAAGCTGATACCGGCGGCATCAAATTCAc                                          <  2:428433/150‑1 (MQ=255)
                                                                                                            tAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACACAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACGGGAAGCTGATACCGGCGGCATCAAATTCAc                                          >  1:703583/1‑150 (MQ=255)
                                                                                                                tGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACGGGAAGCTGATACCGGCGGCATCAAATTCACGttt                                      <  1:1507179/150‑1 (MQ=255)
                                                                                                                 ggCAATCAGAAAAACTACTAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACGGGAAGCTGATACCGGCGGCATCAAATTCACGTTTa                                     <  2:773680/150‑1 (MQ=255)
                                                                                                                 ggCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGACTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACGGGAAGCTGATACCGGCGGCATCAAATTCACGTTTa                                     >  1:1116202/1‑150 (MQ=255)
                                                                                                                  gctaTCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACGGGAAGCTGATACCGGCGGCATCAAATTCACGTTTaa                                    <  2:63163/147‑1 (MQ=255)
                                                                                                                  gCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACGGGAAGCTGATACCGGCGGCATCAAATTCACGTTTaa                                    <  1:1440010/150‑1 (MQ=255)
                                                                                                                    aaTCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACGGGAAGCTGATACCGGCGGCATCAAATTCACGTTTAATa                                  >  2:1109496/1‑150 (MQ=255)
                                                                                                                    aaTCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACGGGAAGCTGATACCGGCGGCATCAAATTCACGTTTAATa                                  <  2:1486716/150‑1 (MQ=255)
                                                                                                                     aTCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACGGGAAGCTGATACCGGCGGCATCAAATTCACGTTTAATAc                                 <  2:1733980/150‑1 (MQ=255)
                                                                                                                       cAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACGGGAAGCTGATACCGGCGGCAGCAAATTCACGTTTAATACGc                               >  2:1234084/1‑150 (MQ=255)
                                                                                                                        aGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACGGGAAGCTGATACCGGCGGCATCAAATTCACGTTTAATACGCt                              >  1:1142847/1‑150 (MQ=255)
                                                                                                                           aaaaCTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACGGGAAGCTGATACCGGCGGCATCAAATTCACGTTTAATACGCTCCa                           <  1:1031016/150‑1 (MQ=255)
                                                                                                                             aaCTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACGGGAAGCTGATACCGGCGGCATCAAATTCACGTTTAATACGCTCCAGc                         <  1:935251/150‑1 (MQ=255)
                                                                                                                              aCTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACGGGAAGCTGATATCGGCGGCATCAAATTCACGTTTAATACGCTCCAGca                        <  1:400171/150‑1 (MQ=255)
                                                                                                                               ataataaGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACGGGAAGCTGATACCGGCGGCATCAAATTCACGTTTAATACGCTCCAGcac                       <  2:835712/149‑1 (MQ=255)
                                                                                                                                taataaGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACGGGAAGCTGATACCGGCGGCATCAAATTCACGTTTAATACGCTCCAGcaca                      <  1:133692/150‑1 (MQ=255)
                                                                                                                                 aataaGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACGGGAAGCTGATACCGGCGGCATCAAATTCACGTTTAATCCGCTCCAGCACAt                     >  1:1498376/1‑150 (MQ=255)
                                                                                                                                 aataaGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACGGGAAGCTGATACCGGCGGCATCAAATTCACGTTTAATACGCTCCAGCACAt                     >  2:1011127/1‑150 (MQ=255)
                                                                                                                                 aataaGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACGGGAAGCTGATACCGGCGGCATCAAATTCACGTTTATACGCTCCAGCACATc                    >  2:1848141/1‑150 (MQ=255)
                                                                                                                                   taaGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACGGGAAGCTGATACCGGCGGCATCAAATTCACGTTTAATACGCTCCAGCACATcc                   >  1:750940/1‑150 (MQ=255)
                                                                                                                                      gagagagGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACGGGAAGCTGATACCGGCGGCATCAAATTCACGTTTAATACGCTCCAGCACATccctg                >  2:1887105/1‑148 (MQ=255)
                                                                                                                                      gagagagGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACGGGAAGCTGATACCGGCGGCATCAAATTCACGTTTAATACGCTCCAGCACATCCCAg                >  2:914566/1‑150 (MQ=255)
                                                                                                                                      gagagagGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACGGGAAGCTGATACCGGCGGCATCAAATTCACGTTTAATACGCTACAGCACATCCCAg                >  2:1351638/1‑150 (MQ=255)
                                                                                                                                       gctgagGCCAGATAATACTCTCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACGGGAAGCTGATACCGGCGGCATCAAATTACCGTTTAATACGCTCCAGCACCTCCCAga               >  2:791180/4‑149 (MQ=255)
                                                                                                                                        gagagGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACGGGAAGCTGATACCGGCGGCATCAAATTAACGTTTAATACGCTCCAGCACATCCCAGTa              >  2:440226/1‑150 (MQ=255)
                                                                                                                                        ctgagGCCAGATAATACTCTCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACGGGAAGCTGATACCGGCGGCATCAAATTCACGTTTAATACGCTCCAGCACATCCCAGTa              >  2:1164239/3‑150 (MQ=255)
                                                                                                                                            ggCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACGGGAAGCTGATACCGGCGGCATCAAATTCACGTTTAATACGCTCCAGCACATCCCAGTACACg          >  2:281958/1‑150 (MQ=255)
                                                                                                                                            ggCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACGGGAAGCTGATACCGGCGGCATCAAATTCACGTTTAATACGCTCCAGCACATCCCAGTACACg          <  2:1449440/150‑1 (MQ=255)
                                                                                                                                            ggCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACGGGAAGCTGATACCGGCGGCATCAAATTCACGTTTAATACGCTCCAGCACATCCCAGTACACg          <  2:45664/150‑1 (MQ=255)
                                                                                                                                            ggCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACGGGAAGCTGATACCGGCGGCATCAAATTCACGTTTAATACGCTCCAGCACATCCCAGTACACg          <  2:39566/150‑1 (MQ=255)
                                                                                                                                                agagagTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACGGGAAGCTGATACCGGCGGCATCAAATTCACGTTTAATACGCTCCAGCACATCCCAGTACACGtttt      <  2:1468067/150‑1 (MQ=255)
                                                                                                                                                    agTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACGGGAAGCTGATACCGGCGGCATCAAATTCACGTTTAATACGCTCCAGCACGTCCCAGTACACGTTTTGCAg  >  2:130049/1‑150 (MQ=255)
                                                                                                                                                    agTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACGGGAAGCTGATACCGGCGGCATCAAATTCACGTTTAATACGCTCCAGCACATCCCAGTACACGTTTTGCAg  <  1:1644693/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CAAGTGCAAGACCTTGAAAGTAATAAGAAAACACGTTGGATATTCCGAATTAAGTATCTGATGCAGAGACTATATCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACGGGAAGCTGATACCGGCGGCATCAAATTCACGTTTAATACGCTCCAGCACATCCCAGTACACGTTTTGCAG  >  NC_000913/3068775‑3069072

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: