Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,000,244 G→C 9.4% intergenic (+99/+48) uvrD → / ← yigE DNA‑dependent ATPase I and helicase II/DUF2233 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,000,2440GC 9.4% 262.8 / 11.7 91intergenic (+99/+48)uvrD/yigEDNA‑dependent ATPase I and helicase II/DUF2233 family protein
Reads supporting (aligned to +/- strand):  ref base G (41/41);  new base C (6/3);  total (47/44)
Fisher's exact test for biased strand distribution p-value = 4.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 103 total aligned reads displayed.

CAGGGTATTAAATGGCTGGTGGCGGCATACGCCCGGCTGGAGTCGGTGTAACGTTGCCGGATGCGGTGCTGCGCACCTTATTTGGCCTAAAAAATCATTCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGACGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGCCCCGCCTTTCATATACATATGTGAAATCGTGCCATCAAGAT  >  NC_000913/4000095‑4000392
                                                                                                                                                     |                                                                                                                                                     
cAGGGTATTAACTGGCTGGTGGCGGCATACGCCCGGCTGGAGTCGGTGTAACGTTGCCGGATGCGGTGCTGCGCACCTTATTTGGCCTAAAAAATCATTCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTAcgccg                                                                                                                                                       <  2:713672/150‑1 (MQ=255)
cAGGGTATTAAATGGCTGGTGGCGGCATACGCCCGGCTGGAGTCGGTGTAACGTTGCCGGATGCGGTGCTGCGCACCTTATTTGGCCTAAAAAATCATTCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTAcgccg                                                                                                                                                       <  1:879520/150‑1 (MQ=255)
 aGGGTATTAAATGGCTGGTGGCGGCATACGCCCGGCTGGAGTCGGTGTAACGTTGCCGGATGCGGTGCTGCGCACCTTATTTGGCCTAAAAAATCATTCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTAcgccgc                                                                                                                                                      <  1:911364/150‑1 (MQ=255)
   ggTATTAAATGGCTGGTGGCGGCATACGCCCGGCTGGAGTCGGTGTAACGTTGCCGGATGCGGTGCTGCGCACCTTATTTGGCCTAAAAAATCATTCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCAt                                                                                                                                                    <  1:1325272/150‑1 (MQ=255)
   ggTATTAAATGGCTGGTGGCGGAATACGCCCGGCTGGAGTCGGTGTAACGTTGCCGGATGCGGTGCTGCGCACCTTATTTGGCCTAAAAAATCATTCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCAt                                                                                                                                                    <  1:1797317/150‑1 (MQ=255)
     tATTAAATGGCTGGTGGCGGCATACGCCCGGCTGGAGTCGGTGTAACGTTGCCGGATGCGGTGCTGCGCACCTTATTTGGCCTAAAAAATCATTCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATGCGGCGTAAACGCCCTATcc                                                                                                                                                  <  1:1620933/150‑1 (MQ=255)
     tATTAAATGGCTGGTGGCGGCATACGCCCGGCTGGAGTCGGTGTAACGTTGCCGGATGCGGTGCTGCGCACCTTATTTGGCCTAAAAAATCATTCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATcc                                                                                                                                                  <  1:390126/150‑1 (MQ=255)
     tATTAAATGGCTGGTGGCGGCATACGCCCGGCTGGAGTCGGTGTAACGTTGCCGGATGCGGTGCTGCGCACCTTATTTGGCCTAAAAAATCATTCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATcc                                                                                                                                                  <  1:1810609/150‑1 (MQ=255)
      aTTAAATGGCTGGTGGCGGCATACGCCCGGCTGGAGTCGGTGTAACGTTGCCGGATGCGGTGCTGCGCACCTTATTTGGCCTAAAAAATCATTCAGATTCAATAAATTGCAACGTCATGTAGGCCCGATAGGGCGTTTACGCCGCATCCg                                                                                                                                                 >  1:54729/1‑150 (MQ=255)
       ttAAATGGCTGGGGCCGGCATACGCCCGGCTGGAGTCGGGGTAACGTTGCCGGATGCGGTGCTGCGCACCTTATTTGGCCTAAAAAATCATTCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATGCGGCGTAAACGCCCTATCCgg                                                                                                                                                <  1:185504/150‑1 (MQ=255)
         aaaTGGCTGGTGGCGGCATACGCCCGGCTGGAGTCGGTGTAACGTTGCCGGATGCGGTGCTGCGCCCCTTATTTGGCCTAAAAAATCATTCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCa                                                                                                                                              >  2:668367/1‑150 (MQ=255)
         aaaTGGCTGGTGGCGGCATACGCCCGGCTGGAGTCGGTGTAACGTTGCCGGATGCGGTGCTGCGCACCTTATTTGGCCTAAAAAATCATTCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCa                                                                                                                                              <  1:72250/150‑1 (MQ=255)
         aaaTGGCTGGTGGCGGCATACGCCCGGCTGGAGGCGGTGTAACGTTGCCGGATGCGGTGCTGCGCACCTTATTTGGCCTAAAAAATCATTCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCa                                                                                                                                              <  2:1798426/150‑1 (MQ=255)
            tGGCTGGTGGCGGCATACGCCCGGCTGGAGTCGGTGTAACGTTGCCGGATGCGGTGCTGCGCACCTTATTTGGCCTAAAAAATCATTCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATGCGGCGTAAACGCCCTATCCGGcctac                                                                                                                                           <  2:150374/150‑5 (MQ=255)
              gCTGGTGGCGGCATACGCCCGGCTCGAGTCGGTGTAACGTTGCCGGATGCGGTGCTGCGCACCTTATTTGGCCTAAAAAATCATTCAGATTCAATAAATTCCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGAATCTgc                                                                                                                                         >  2:777698/1‑150 (MQ=255)
               cTGGTGGCGGCATACGCCCGGCTGGAGTCGGTGTAACGTTGCCGGATGCGGTGCTGCGCACCTTATTTGGCCTAAAAAATCATTCAGATGCAATAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTgcg                                                                                                                                        >  2:942810/1‑150 (MQ=255)
                  gtggCGGCATACGCCCGGCTGGAGTCGGTGTAACGTTGCCGGATGCGGTGCTGCGCACCTTATTTGGCCTAAAAAATCATTCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCa                                                                                                                                     <  1:495181/150‑1 (MQ=255)
                       ggcATACGCCCGGCTGGAGTCGGTGTAACGTTGCCGGATGCGGTGCTGCGCACCTTATTTGGCCTAAAAAATCATTCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTc                                                                                                                                <  1:1547613/150‑1 (MQ=255)
                       ggcATACGCCCGGCTGGAGTCGGTGTAACGTTGCCGGATGCGGTGCTGCGCACCTTATTTGGCCTAAAAAATCATTCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTc                                                                                                                                <  2:1047794/150‑1 (MQ=255)
                         cATACGCCCGTCTGGAGTCGTTGTAACGTTGCCGGATGCGGTGCTGCGCACCTTATTTGGCCTAAAAAATCATTCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAg                                                                                                                              <  1:1433559/150‑1 (MQ=255)
                          aTACGCCCGGCTGGAGTCGGTGTAACGTTGCCGGATGCGGTGCTGCGCACCTTATTTGGCCTAAAAAATCAGTCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGt                                                                                                                             >  2:218905/1‑150 (MQ=255)
                             cGCCCGGCTGGAGTCGGTGTAACGTTGCCGGATGCGGTGCTGCGCACCTTATTTGGCCTAAAAAATCATTCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATc                                                                                                                          >  2:1088381/1‑150 (MQ=255)
                               cccGGCTGGAGTCGGTGTAACGTTGCCGGATGCGGTGCTGCGCACCTTATTTGGCCTAAAAAATCATTCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTg                                                                                                                        <  1:1660240/150‑1 (MQ=255)
                               cccGGCTGGAGTCGGTGTAACGTTGCCGGATGCGGTGCTGCGCACCTTATTTGGCCTAAAAAATCATTCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTg                                                                                                                        <  2:1495710/150‑1 (MQ=255)
                                    cTGGAGTCGGTGTAACGTTGCCGGATGCGGTGCTGCGCACCTTATTTGGCCTAAAAAATCATTCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGAcaca                                                                                                                   >  1:36320/1‑150 (MQ=255)
                                            ggTGTAACGTTGCCGGATGCGGTGCTGCGCACCTTATTTGGCCTAAAAAATCATTCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATATGACACAAAACTATc                                                                                                           >  1:1343357/1‑150 (MQ=255)
                                             gtgtAACGTTGCCGGATGCGGTGCTGCGCACCTTATTTGGCCTAAAAAATCATTCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCg                                                                                                          <  2:1235000/150‑1 (MQ=255)
                                               gtAACGTTGCCGGATGCGGTGCTGCGCACCTTATTTGGCCTAAAAAATCATTCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGtt                                                                                                        <  2:870772/150‑1 (MQ=255)
                                               gtAACGTTGCCGGATGCGGTGCTGCGCACCTTATTTGGCCTAAAAAATCATTCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGtt                                                                                                        <  2:862924/150‑1 (MQ=255)
                                                 aaCGTTGCCGGATGCGGTGCTGCGCACCTTATTTGGCCTAAAAAATCATTCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGtttt                                                                                                      <  2:1092282/150‑1 (MQ=255)
                                                  aCGTTGCCGGATGCGGTGCTGCGCACCTTATTTGGCCTAAAAAATCATTCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTa                                                                                                     <  2:1234796/150‑1 (MQ=255)
                                                  aCGTTGCCGGATGCGGTGCTGCGCACCTTATTTGGCCTAAAAAATCATTCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTa                                                                                                     <  2:81586/150‑1 (MQ=255)
                                                   cGTTGCCGGATGCGGTGCTGCGCACCTTATTTGGCCTAAAAAATCATTCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTaa                                                                                                    >  2:231734/1‑150 (MQ=255)
                                                    gTTGCCGGATGCGGTGCTGCGCACCTTATTTGGCCTAAAAAATCATTCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAAc                                                                                                   >  2:1757953/1‑150 (MQ=255)
                                                    gTTGCCGGATGCGGTGCTGCGCACCTTATTTGGCCTAAAAAATCATTCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAAc                                                                                                   <  2:333879/150‑1 (MQ=255)
                                                    gTTGCCAGATGCGGTGCTGCGCACCTTATTTGGCCTAAAAAATCATTCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAAc                                                                                                   <  2:1209139/150‑1 (MQ=255)
                                                      tGCCGGATGCGGTGCTGCGCACCTTATTTGGCCTAAAAAATCATTCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCt                                                                                                 >  1:46898/1‑150 (MQ=255)
                                                          ggATGCGGTGCTGCGCACCTTATTTGGCCTAAAAAATCATTCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATAGGTCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTcc                                                                                             <  2:1399786/150‑1 (MQ=255)
                                                               cGGTGCTGCGCACCTTATTTGGCCTAAAAAATCATTCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCg                                                                                        <  2:1086619/150‑1 (MQ=255)
                                                               cGGTGCTGCGCACCTTATTTGGCCTAAAAAATCATTCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCg                                                                                        <  2:926136/150‑1 (MQ=255)
                                                                ggTGCTGCGCACCTTATTTGGCCTAAAAAATCATTCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGa                                                                                       >  2:1579951/1‑150 (MQ=255)
                                                                         cACCTTATTTGGCCTAAAAAATCATTCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGACGGAAATCa                                                                             >  1:639273/1‑150 (MQ=255)
                                                                          aCCTTATTTGGCCTAAAAAATCATTCAGATTCAATAAATTGTAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGTCGAGCGGaaaggtt                                                                            <  1:1380672/150‑5 (MQ=39)
                                                                           ccTTATTTGGCCTAAAAAATCATTCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGACGGAAATCATg                                                                           >  2:1485875/1‑150 (MQ=255)
                                                                           ccTTATTTGGCCTAAAAAATCATTCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGACGGAAATCATg                                                                           >  2:1695851/1‑150 (MQ=255)
                                                                           ccTTATTTGGCCTAAAAAATCATTCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATAGGGAGTTTACGCCGCATCCGCCATCTGCGGCATCTTCAGTATCTGGCACAAAACTATCGATTTAACCTTTCCGATCGACGGAAATCATg                                                                           >  2:1775609/1‑150 (MQ=255)
                                                                                    gCCTAAAAAATCATTCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGACGGAAATCATGGTGACaaaa                                                                  >  1:1830221/1‑150 (MQ=255)
                                                                                     ccTAAAAAATCATTCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGACGGAAATCATGGTGACAAAAg                                                                 >  2:1068013/1‑150 (MQ=255)
                                                                                     ccTAAAAAATCATTCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGACGGAAATCATGGTGACAAAAg                                                                 >  1:1367937/1‑150 (MQ=255)
                                                                                     ccTAAAAAATCATTCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGACGGAAATCATGGTGACAAAAg                                                                 >  1:1536249/1‑150 (MQ=255)
                                                                                       tAAAAAATCATTCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGACGGAAATAATGGTGACAAAAGGa                                                               >  1:1813994/1‑150 (MQ=255)
                                                                                       tAAAAAATCATTCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAAATATCGTTTTAACCTTTCCGCTCGACGGAAATCATGGTGACAAAAggt                                                               >  1:180992/1‑149 (MQ=255)
                                                                                        aaaaaaTCATTCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGACGGAAATCATGGTGACAAAAGGAt                                                              <  2:143442/150‑1 (MQ=255)
                                                                                                 ttCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGACGGAAATCATGGTGACAAAAGGATAGCGTTGcc                                                     >  1:1637666/1‑150 (MQ=255)
                                                                                                  tCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGACGGAAATCATGGTTACAAAAGGATAGCGTTGccc                                                    >  2:73074/1‑149 (MQ=255)
                                                                                                      aTTCAATAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGACGGAAATCATGGTGACAAAAGGATAGCGTTGCCACggg                                                <  1:616160/150‑1 (MQ=255)
                                                                                                          aaTAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAACACTATCGTTTTAACCTTTCCGCTCGACGGAAATCATGGTGACAAAAGGATAGCGTTGCCCCGGGATTg                                            >  2:847535/1‑150 (MQ=255)
                                                                                                          aaTAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGACGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTg                                            >  2:1098304/1‑150 (MQ=255)
                                                                                                          aaTAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGACGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTg                                            >  1:518435/1‑150 (MQ=255)
                                                                                                           aTAAATTGCAACGTCATGTAGGACGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGACGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGc                                           <  1:208434/150‑1 (MQ=255)
                                                                                                             aaaTTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGACGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGccc                                         >  2:420289/1‑150 (MQ=255)
                                                                                                             aaaTTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGACGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGacc                                         >  1:551854/1‑150 (MQ=255)
                                                                                                              aaTTTCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGACGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGcccc                                        <  2:778136/150‑1 (MQ=255)
                                                                                                              aaTTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGACGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGcccc                                        <  2:778133/150‑1 (MQ=255)
                                                                                                                   cAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGACGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGCCCCGCCtt                                   >  2:867472/1‑150 (MQ=255)
                                                                                                                     aCGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGACGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGCCCCGCCTTTc                                 <  2:1536249/150‑1 (MQ=255)
                                                                                                                     aCGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGACGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGCCCCGCCTTTc                                 <  1:1251158/150‑1 (MQ=255)
                                                                                                                      cGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGACGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGCCCCGCCTTTCa                                <  2:1418048/150‑1 (MQ=255)
                                                                                                                      cGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCAGCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGACGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGCCCCGCCTTTCa                                <  1:305744/150‑1 (MQ=255)
                                                                                                                       gTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGACGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGCCCCGCCTTTCat                               >  2:848769/1‑150 (MQ=255)
                                                                                                                        tCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGACGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGCCCCGCCTTTCata                              <  2:645109/150‑1 (MQ=255)
                                                                                                                         cATGTAGGCCGGATAGGGCGTTTACGCCGCATCAGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGACGGAAGTCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGCCCCGCCTTTCatat                             >  1:548624/1‑150 (MQ=255)
                                                                                                                          gcagatgCCGGATGCGGCGTAAACGCCCTATCCGGCATCTGCGCAATCTTCAGTATCTGACACAAAACTATCGTTTTATCCTTTCCGCTCGACGGAATTCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGCCCCGCCTTTCAtata                            >  1:1098981/7‑150 (MQ=38)
                                                                                                                           cagatgCCGGATGCGGCGTAAACGCCCTATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGACGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGCCCCGCCTTTCATATAc                           >  1:116344/6‑150 (MQ=255)
                                                                                                                             tAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGACGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGCCCCGCCTTTCATATACat                         >  1:455874/1‑150 (MQ=255)
                                                                                                                             gatgCCGGATGCGGCGTAAACGCCCTATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGACGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGCCCCGCCTTTCATATACat                         >  1:1628669/4‑150 (MQ=255)
                                                                                                                              cggCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCATCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGACGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGCCCCGCCTTTCATATACata                        <  2:1577908/149‑1 (MQ=255)
                                                                                                                              aGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGACGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGCCCCGCCTTTCATATACata                        >  2:669876/1‑150 (MQ=255)
                                                                                                                               tgCCGGATGCGGCGTAAACGCCCTATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGACGGAAATCATGGTGCCAAAAGGATAGCGTTGCCACGGGATTGCCCCGCCGTTCATATACatat                       >  2:1737029/2‑150 (MQ=255)
                                                                                                                               tgCCGGATGCGGCGTAAACGCCCTATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGACGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGCCCCGCCTATCATATACattt                       >  2:1324546/2‑148 (MQ=255)
                                                                                                                               tgCCGGATGCGGCGTAAACGCCCTATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACGATCGTGTTACCCTTTCCGCTCGACGGAAATCATGGGGACAACAGGATAGCGTTGCCACGGGATTGCCCCGCCTTGCATATACatat                       >  2:1318213/2‑150 (MQ=255)
                                                                                                                                gCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGACGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGCCCCGCCTTTCATATACATAtg                      <  1:1741157/150‑1 (MQ=255)
                                                                                                                                gCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGACGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGCCCCGCCTTTCATATACATAtg                      >  1:1488200/1‑150 (MQ=255)
                                                                                                                                 ccGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGCCGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGCCCCGCCTTTCATATACATAtgt                     >  1:1588615/1‑150 (MQ=255)
                                                                                                                                 ccGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGACGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGCCCCGCCTTTCATATACATAtgt                     >  1:1749430/1‑150 (MQ=255)
                                                                                                                                      tAGGGCGGTTACGCCGCATGCGGCATCTGCGCCATCTTCAGTATGTGACACAAAACTATCGTTTTAACCTTTGCGCTCGACGGAAATCATTGTGACAAAAGGATAGGGTTGCCACGGGATTGCCCCGCCTTTCATATACATATGTGaatt                <  2:1478030/150‑3 (MQ=255)
                                                                                                                                         ggCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGACGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGCCCCGCCTTTCATATACATATGTGAAATCGt             >  1:1606757/1‑150 (MQ=255)
                                                                                                                                         ggCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGACGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGCCCCGCCTTTCATATACATATGTGAAATCGt             >  1:387532/1‑150 (MQ=255)
                                                                                                                                           cGTTTACGCCGCATCCGGCAGCTGCGCCATCTTCAGTATCTGACACAATACTATCGTTTTAACCTTTCCGATCGACGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGCCCCGCCTTTCATATACATATGTGAAATCGTGc           <  1:319502/150‑1 (MQ=255)
                                                                                                                                             tttACGCCGCATCCGGCATCTGCGACATCTTCAGTATATGACACAAAACTATCGTTTTAACATTTCCGCTCGACGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGTATTGCCCCGCATTTCAAATACATATGTGAAATCGTGCCa         <  1:1225649/150‑1 (MQ=255)
                                                                                                                                                aCGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGACGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGCCCCGCCTTTCATATACATATGTGAAATCGTGCCATCa      <  1:1088381/150‑1 (MQ=255)
                                                                                                                                                 cgccgcATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGACGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGCCCCGCCTTTCATATACATATGTGAAATCGTGCCATCaa     >  1:880966/1‑150 (MQ=255)
                                                                                                                                                 cgccgcATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGACGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGCCCCGCCTTTCATATACATATGTGAAATCGTGCCATCaa     <  1:867472/150‑1 (MQ=255)
                                                                                                                                                 cgccgcATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGACGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGCCCCGCCTTTCATATACATATGTGAAATCGTGCCATCaa     >  2:1460456/1‑150 (MQ=255)
                                                                                                                                                 cgccgcATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGACGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGCCCCGCCTTTCATATACATATGTGAAATCGTGCCATCaa     <  2:36320/150‑1 (MQ=255)
                                                                                                                                                 cgccgcATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGACGGAAATCATGGTGACAAAAGGATAGCGTAGCCACGGGATTGCCCCGCCTTTCATATACATATGTGAAATCGTGGCATCaa     >  1:182111/1‑150 (MQ=255)
                                                                                                                                                  gccgcATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGACGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGCCCCGCCTTTCATATACATATGTGAAATCGTGCCATCAAg    >  1:473957/1‑150 (MQ=255)
                                                                                                                                                  gccgcAGACGGCATCTGCGCAATCTTCAGTATCTGACTCAAAACTATCGTGTAAACCTTTCCGCTCGACGGAATTCATGGTGTCAAAAGGATAGCGTGGCCACGGGAGTGCCCCGCCTTTCATATACATATGTGAAATCGTGCCATCAAg    <  2:279247/150‑1 (MQ=255)
                                                                                                                                                   ccgcATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTAGCGTTTTAACCTTTCCGCTCGCCGGAAATCATGGTGACAAAAGGATAGCCTTGCCACGGGCTTGCCCCGCCTGTCTGATCCATATGTGATATCGTGCCATAAAGa   >  2:1126428/1‑150 (MQ=255)
                                                                                                                                                    cgcATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGACGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGCCCCGCCTTTCATATACATATGTGAAATCGTGCCATCAAGAt  >  2:855843/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CAGGGTATTAAATGGCTGGTGGCGGCATACGCCCGGCTGGAGTCGGTGTAACGTTGCCGGATGCGGTGCTGCGCACCTTATTTGGCCTAAAAAATCATTCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGACGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGCCCCGCCTTTCATATACATATGTGAAATCGTGCCATCAAGAT  >  NC_000913/4000095‑4000392

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