Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,068,924 G→A 6.8% 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 6.8% 419.7 / 11.4 134intergenic (‑116/+23)argO/mscSarginine transporter/mechanosensitive channel protein, small conductance
Reads supporting (aligned to +/- strand):  ref base G (56/67);  new base A (2/7);  total (58/76)
Fisher's exact test for biased strand distribution p-value = 2.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.28e-01

Only 100 of 136 total aligned reads displayed.

CAAGTGCAAGACCTTGAAAGTAATAAGAAAACACGTTGGATATTCCGAATTAAGTATCTGATGCAGAGACTATATCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACGGGAAGCTGATACCGGCGGCATCAAATTCACGTTTAATACGCTCCAGCACATCCCAGTACACGTTTTGCA  >  NC_000913/3068775‑3069071
                                                                                                                                                     |                                                                                                                                                   
  aGTGCAAGACCTTGAAAGTAATAAGAAAACACGTTGGATATTCCGAATTAAGTATCTGATGCAGAGACTATATCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGta                                                                                                                                                   >  2:705619/1‑150 (MQ=255)
     gCAAGACCTTGAAAGTAATAAGAAAACACGTTGGATATTCCGAATTAAGTATCTGATGCAGAGACTATATCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGtatta                                                                                                                                                >  1:347634/1‑150 (MQ=255)
     gCAAGACCTTGAAAGTAATAAGAAAACACGTTGGATATTCCGAATTAAGTATCTGATGCAGAGACTATATCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGtatta                                                                                                                                                <  2:439512/150‑1 (MQ=255)
     gCAAGACCTTGAAAGTAATAAGAAAACACGTTGGATATTCCGAATTAAGTATCTGATGCAGAGACTATATCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGtatta                                                                                                                                                >  1:1302846/1‑150 (MQ=255)
     gCAAGACCTTGAAAGTAATAAGAAAACACGTTGGATATTCCGAATTAAGTATCTGATGCAGAGACTATATCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGtatta                                                                                                                                                >  2:800558/1‑150 (MQ=255)
     gCAAGACCTTGAAAGTAATAAGAAAACACGTTGGATATTCCGAATTAAGTATCTGATGCAGAGACTATATCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGtatta                                                                                                                                                >  2:705763/1‑150 (MQ=255)
       aaGACCTTGAAAGTAATAAGAAAACACGTTGGATATTCCGAATTAAGTATCTGATGCAGAGACTATATCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATc                                                                                                                                              <  1:1633985/150‑1 (MQ=255)
        aGACCTTGAAAGTAATAAGAAAACACGTTGGATATTCCGAATTAAGTATCTGATGCAGAGACTATATCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCt                                                                                                                                             <  2:219871/150‑1 (MQ=255)
         gACCTTGAAAGTAATAAGAAAACACGTTGGATATTCCGAATTAAGTATCTGATGCAGAGACTATATCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTg                                                                                                                                            <  1:414425/150‑1 (MQ=255)
          aCCTTGAAAGTAATAAGAAAACACGTTGGATATTCCGAATTAAGTATCTGATGCAGAGACTATATCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTgg                                                                                                                                           <  1:1590159/150‑1 (MQ=255)
           ccTTGAAAGTAATAAGAAAACACGTTGGATATTCCGAATTAAGTATCTGATGCAGAGACTATATCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGc                                                                                                                                          >  2:1302164/1‑150 (MQ=255)
            cTTGAAAGTAATAAGAAAACACGTTGGATATTCCGAATTAAGTATCTGATGCAGAGACTATCTCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGcc                                                                                                                                         <  2:97183/150‑1 (MQ=255)
            cTTGAAAGTAATAAGAAAACACGTTGGATATTCCGAATTAAGTATCTGATGCAGAGACTATATCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATTAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGcc                                                                                                                                         >  1:1618461/1‑150 (MQ=255)
            cTTGAAAGTAATAAGAAAACACGTTGGATATTCCGAATTAAGTATCTGATGCAGAGACTATATCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGATAATACTCTCTGGcc                                                                                                                                         <  1:340994/150‑1 (MQ=255)
            cTTGAAAGTAATAAGAAAACACGTTGGATATTCCGAATTAAGTATCTGATGCAGAGACTATATCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGATAATACTCTCTGGcc                                                                                                                                         <  2:829292/150‑1 (MQ=255)
              tGAAAGTAATAAGAAAACACGTTGGATATTCCGAATTAAGTATCTGATGCAGAGACTATATCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGATAATACTCTCTGGCCTc                                                                                                                                       <  1:238054/150‑1 (MQ=255)
              tGAAAGTAATAAGAAAACACGTTGGATATTCCGAATTAAGTATCTGATGCAGAGACTATATCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGATAATACTCTCTGGCCTc                                                                                                                                       <  2:27485/150‑1 (MQ=255)
               gAAAGTAATAAGAAAACACGTTGGATATTCCGAATTAAGTATCTGATGCAGAGACTATATCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGATAATACTCTCTGGCCTct                                                                                                                                      <  1:787973/150‑2 (MQ=255)
                aaaGTAATAAGAAAACACGTTGGATATTCCGAATTAAGTATCTGATGCAGAGACTATATCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAg                                                                                                                                     >  2:269904/1‑150 (MQ=255)
                    taataaGAAAACACGTTGGATATTCCGAATTAAGTATCTGATGCAGAGACTATATCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGtt                                                                                                                                 >  1:763735/1‑150 (MQ=255)
                     aataaGAAAACACGTTGGATATTCCGAATTAAGTATCTGATGCAGAGACTATATCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTg                                                                                                                                <  2:1410883/150‑1 (MQ=255)
                      ataaGAAAACAGGTTGGATATTCCGAATTAAGTATCTGATGCAGAGACTATATCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCTATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGa                                                                                                                               <  1:1229723/150‑1 (MQ=255)
                        aaGAAAACACGTTGGATATTCAGAATTAAGTATCTGATGCAGAGACTATATCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGAtt                                                                                                                             >  1:1121558/1‑150 (MQ=255)
                         aGAAAACACGTTGGATATTCCGAATTAAGTATCTGATGCAGAGACTATATCGCGGGATTGTTAATAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTa                                                                                                                            <  2:1619445/150‑1 (MQ=255)
                          gAAAACACGTTGGATATTCGGAATTAAGTATCTGATGCAGAGACTATATCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTAc                                                                                                                           <  2:206919/150‑1 (MQ=255)
                          gAAAACACGTTGGATATTCCGAATTAAGTATCTGATGCAGAGACTATATCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTAc                                                                                                                           >  2:733178/1‑150 (MQ=255)
                          gAAAACACGTTGGATATTCCGAATTAAGTATCTGATGCAGAGACTATATCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGAGTAGATTaa                                                                                                                           >  1:1069706/1‑149 (MQ=255)
                           aaaaCACGTTGGATATTCCGAATTAAGTATCTGATGCAGAGACTATATCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACg                                                                                                                          <  2:317808/150‑1 (MQ=255)
                             aaCACGTTGGATATTCCGAATTAAGTATCTGATGCAGAGACTATATCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCa                                                                                                                        <  1:384442/150‑1 (MQ=255)
                             aaCACGTTGGATATTCCGAATTAAGTATCTGATGCAGAGACTATATCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCa                                                                                                                        <  2:317465/150‑1 (MQ=255)
                               cacGTTGGATATTCCGAATTAAGTATCTGATGCAGAGACTATATCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGc                                                                                                                      >  1:1110685/1‑150 (MQ=255)
                                acGTTGGATATTCCGAATTAAGTATCTGATGCAGAGACTATATCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCt                                                                                                                     <  1:191464/150‑1 (MQ=255)
                                   tttgATATTCCGAATTAAGTATCTGATTCAGAGACTATATCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTg                                                                                                                  <  2:820382/147‑1 (MQ=255)
                                   ttGGATATTCCGAATTAAGTATCTGATGCAGAGACTATATCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTg                                                                                                                  >  1:330481/1‑150 (MQ=255)
                                   ttGGATATTCCGAATTAAGTATCTGATGCAGAGACTATATCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTg                                                                                                                  <  2:1647558/150‑1 (MQ=255)
                                     ggATATTCCGAATTAAGTATCTGATGCAGAGACTATATCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGtc                                                                                                                >  1:724983/1‑150 (MQ=255)
                                      gATATTCCGAATTAAGTATCTGATGCAGAGACTATATCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGtct                                                                                                               <  1:800558/150‑1 (MQ=255)
                                      gATATTCCGAATTAAGTATCTGATGCAGAGACTATATCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGtct                                                                                                               <  2:505877/150‑1 (MQ=255)
                                        tatTCCGAATTAAGTATCTGATGCAGAGACTATATCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGtcttc                                                                                                             <  1:1549268/150‑1 (MQ=255)
                                          ttCCGAATTAAGTATCTGATGCAGAGACTATATCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTcttctt                                                                                                           >  2:725947/1‑150 (MQ=255)
                                           tCCGAATTAAGTATCTGCTGCAGAGACTATATCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCTGAGAGTATTATCTGGCCTCAGCGTGGATTCCGCAGCTTTGTcttctta                                                                                                          >  1:872658/1‑149 (MQ=255)
                                            ccGAATTAAGTATCTTATGCAGAGACTATATCTGGGGATTGTTATTAGCGGAAATTGATAATATTAATTGTCATTAAGAAAAACTATTTTGATATAGGCCAGAGATTATTATCTTGCCTGAGCGTTGATTACGCTGCTTTGTCTTCtttg                                                                                                         <  2:1610190/150‑2 (MQ=255)
                                            acGAATTAAGTATCTGATGCAGAGACTATATCGCGTGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTTATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTc                                                                                                         <  1:1130061/149‑1 (MQ=255)
                                             cGAATTAAGTATCTGATGCAGAGACTATATCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCa                                                                                                        <  2:724983/150‑1 (MQ=255)
                                                aTTAAGTATCTGATGCAGAGACTATATCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCAccc                                                                                                     <  2:1256060/150‑1 (MQ=255)
                                                 ttAAGTATCTGATGCAGAGCCTATATCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTTGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCg                                                                                                    <  2:349985/150‑1 (MQ=255)
                                                 ttAAGTATCTGATGCAGAGACTATATCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCg                                                                                                    <  1:733178/150‑1 (MQ=255)
                                                     gTATCTGATGCAGAGACTATATCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTa                                                                                                <  2:194482/150‑1 (MQ=255)
                                                        tCTGATGCAGAGACTATATCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGATAATACTCTCTGGCCTctctct                                                                                                                                  <  1:1411021/113‑6 (MQ=38)
                                                        tCTGATGCAGAGACTATATCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGATAATACTCTCTGGCCTctctct                                                                                                                                  >  2:1411021/1‑108 (MQ=38)
                                                        tCTGATGCAGAGACTATATCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAg                                                                                             <  2:1310895/150‑1 (MQ=255)
                                                         cTGATGCAGAGACTATATCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGt                                                                                            >  1:754925/1‑150 (MQ=255)
                                                              gCAGAGACTATATCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTcaca                                                                                       <  1:499910/150‑1 (MQ=255)
                                                                agagACTATATCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATc                                                                                     <  2:564502/150‑1 (MQ=255)
                                                                     cTATATCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCAttt                                                                                >  2:1200304/1‑150 (MQ=255)
                                                                      tataTCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTg                                                                               <  2:499527/150‑1 (MQ=255)
                                                                           cgcgGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGt                                                                          >  2:1344525/1‑150 (MQ=255)
                                                                             cgGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTAc                                                                        >  1:1505856/1‑150 (MQ=255)
                                                                              gggATTGTTTTTAGCGGAAATTGATAATATTAATGGGCAATCAGAAAAACTAATAAGCGCGCGGCCAGCGACGATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACg                                                                       <  2:1084564/150‑1 (MQ=255)
                                                                                gATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACggg                                                                     <  1:1496323/150‑1 (MQ=255)
                                                                                gATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACggg                                                                     >  2:499877/1‑150 (MQ=255)
                                                                                    gTTAATAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACGGGAAGc                                                                 <  2:60856/150‑1 (MQ=255)
                                                                                       attaGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACGGGAAGCTGa                                                              <  1:275948/150‑1 (MQ=255)
                                                                                         taGCGGAAATTGATAATATTAATTTGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACGGGAAGCTGATa                                                            >  2:1394677/1‑150 (MQ=255)
                                                                                         taGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACGGGAAGCTGATa                                                            <  2:1110685/150‑1 (MQ=255)
                                                                                          aGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACGGGAAGCTGATAc                                                           >  2:372819/1‑150 (MQ=255)
                                                                                          aGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACGGGAAGCTGATAc                                                           >  1:1081168/1‑150 (MQ=255)
                                                                                              gAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACGGGAAGCTGATACcggc                                                       >  2:846681/1‑150 (MQ=255)
                                                                                                 aTTGATAATATTAATTTGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGATTTGTCATATTTAACCCGCTTAAAGTTCAAATCCATTTGCGGGTACGGGAAGCTGATACCgtcggc                                                    >  2:785572/1‑150 (MQ=255)
                                                                                                    gATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACGGGAAGCTGATACCGGCGGCATc                                                 >  1:1283325/1‑150 (MQ=255)
                                                                                                     ataataTTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACGGGAAGCTGATACCGGCGGCATCa                                                <  1:899369/150‑1 (MQ=255)
                                                                                                     ataataTTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACGGGAAGCTGATACCGGCGGCATCa                                                <  2:560719/150‑1 (MQ=255)
                                                                                                               ttGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTCTTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACGGGAAGCTGATACCGGCGGCATCAAATTCACGtt                                      <  1:269904/150‑1 (MQ=255)
                                                                                                                tGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACGGGAAGCTGATACCGGCGGCATCAAATTCACGttt                                     >  1:20246/1‑150 (MQ=255)
                                                                                                                tGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACGGGAAGCTGATACCGGCGGCATCAAATTCACGttt                                     >  2:1545214/1‑150 (MQ=255)
                                                                                                                   cAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACGGGAAGCTGATACCGGCGGCATCAAATTCACGTTTAAt                                  >  1:922277/1‑150 (MQ=255)
                                                                                                                   cAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACGGGAAGCTGATACCGGCGGCATCAAATTCACGTTTAAt                                  >  2:1586838/1‑150 (MQ=255)
                                                                                                                    aaTCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACGGGAAGCTGATACCGGCGGCATCAAATTCACGTTTAATa                                 >  2:882882/1‑150 (MQ=255)
                                                                                                                    aaTCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACGGGAAGCTGATACCGGCGGCATCAAATTCACGTTTAATa                                 <  1:1419172/150‑1 (MQ=255)
                                                                                                                     aTCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACGTGAATTTGATGCCGCCGGTATCAGCTTCCCGtacccgca                                <  2:763735/150‑8 (MQ=37)
                                                                                                                     aTCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACGGGAAGCTGATACCGGCGGCATCAAATTCACGTTTAATAc                                >  2:1579232/1‑150 (MQ=255)
                                                                                                                     aTCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACGGGAAGCTGATACCGGCGGCATCAAATTCACGTTAAATAc                                >  1:391562/1‑150 (MQ=255)
                                                                                                                       cAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACGGGAAGCTGATACCGGCGGCATCAAATTCACGTTTAATACGc                              <  1:1629109/150‑1 (MQ=255)
                                                                                                                       cAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACGGGAAGCTGATACCGGCGGCATCAAATTCACGTTTAATACGc                              >  2:179537/1‑150 (MQ=255)
                                                                                                                       cAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACGGGAAGCTGATACCGGCGGCATCAAATTCACGTTTAATACGc                              <  1:950454/150‑1 (MQ=255)
                                                                                                                        aGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACGGGAAGCTGATACCGGCGGCATCAAATTCACGTTTAATACGCt                             <  1:705763/150‑1 (MQ=255)
                                                                                                                        aGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACGGGAAGCTGATACCGGCGGCATCAAATTCACGTTTAATACGCt                             <  2:1595213/150‑1 (MQ=255)
                                                                                                                        aGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCATTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACGGGAAGCTGATACCGGCGGCATCAAATTCACGTTTAATACGCt                             <  2:330481/150‑1 (MQ=255)
                                                                                                                         gAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACGGGAAGCTGATACCGGCGGCATCAAATTCACGTTTAATACGCTc                            >  1:1579762/1‑150 (MQ=255)
                                                                                                                            aaaCTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCTGGTACTGGAAGCTTATACCGTCGGCATCAAATTCACGTTTAATACGCTCCAg                         <  2:1121558/150‑1 (MQ=255)
                                                                                                                            aaaCTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACGGGAAGCTGATACCGGCGGCATCAAATTCACGTTTAATACGCTCCAg                         <  2:1283325/150‑1 (MQ=255)
                                                                                                                                    aaGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACGGGAAGCTGATACCGGCGGCATCAAATTCACGTTTAATACGCTCCAGCACATccc                 <  1:1586838/150‑1 (MQ=255)
                                                                                                                                      gagagagGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACGGGAAGCTGATACCGGCGGCATCAAATTCACGTTTAATACGCTCCAGCACATCCCAg               <  1:811824/150‑1 (MQ=255)
                                                                                                                                        gagagGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACGGGAAGCTGATACCGGCGGCATCAAATTCACGTTTAATACGCTCCAGCACATCCCAGTa             <  1:1394677/150‑1 (MQ=255)
                                                                                                                                        ctgagGCCAGATAATACTCTCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACGGGAAGCTGATACCGGCGGCATCAAATTCACGTTTAATACGCTCCAGCACATCCCAGTa             >  1:1076886/3‑150 (MQ=255)
                                                                                                                                          gagGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACGGGAAGCTGATACCGGCGGCATCAAATTCACGTTTAATACGCTCCAGCACATCCCAGTaca           >  1:580390/1‑150 (MQ=255)
                                                                                                                                             gCCAGAGAGTATTATCTGGCCTCAGCCTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACGGGAAGCTGATACCTGCGGCATAAAATTCACGTTTAATTCGCTCCAGCACATCCCAGTACACGt        >  1:1571876/1‑150 (MQ=255)
                                                                                                                                               cAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACGGGAAGCTGATACCGGCGGCATCAAATTCACGTTTAATACGCTCCAGCACATCCCAGTACACGttt      >  1:1585557/1‑150 (MQ=255)
                                                                                                                                               cAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACGGGAAGCTGATACCGGCGGCATCAAATTCACGTTTAATACGCTCCAGCACATCCCAGTACACGttt      <  1:882882/150‑1 (MQ=255)
                                                                                                                                                   gagTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACGGGAAGCTGATACCGGCGGCATCAAATTCACGTTTAATACGCTCCAGCACATCCCAGTACACGTTTTGCa  <  1:521338/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
CAAGTGCAAGACCTTGAAAGTAATAAGAAAACACGTTGGATATTCCGAATTAAGTATCTGATGCAGAGACTATATCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACGGGAAGCTGATACCGGCGGCATCAAATTCACGTTTAATACGCTCCAGCACATCCCAGTACACGTTTTGCA  >  NC_000913/3068775‑3069071

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