Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 361,895 A→C 6.0% intergenic (‑35/+31) lacA ← / ← lacY thiogalactoside acetyltransferase/lactose permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913361,8950AC 6.0% 272.1 / 4.3 101intergenic (‑35/+31)lacA/lacYthiogalactoside acetyltransferase/lactose permease
Reads supporting (aligned to +/- strand):  ref base A (45/49);  new base C (2/5);  total (47/54)
Fisher's exact test for biased strand distribution p-value = 4.45e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01

Only 100 of 104 total aligned reads displayed.

TGATTAAACTCATACATTAACGTTTTCCCACGAAGTCTTTTTTCCGGTAAGCCTTCGCACATATCGGTAAATAGCTTGCCTGCTCTTATTCTTTCGGTCATTGGCATGTTCAATGCGATCACTCCGTTATGATATGTTGGTCGGATAAGGCGCTCGCGCCGCATCCGACATTGATTGCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTAAGGTGAAGCCCAGCGCCACCAGACCCAGCACCAGATAAGCGCCCT  >  NC_000913/361748‑362044
                                                                                                                                                   |                                                                                                                                                     
tGATTAAACTCATACATTAACGTTTTCCCACGAAGTCTTTTTTCCGGTAAGCCTTCGCACATATCGGTAAATAGCTTGCCTGCTCTTATTCTTTCGGTCATTGGCATGTTCAATGCGATCACTCCGTTATGATATGTTGGTCGGATAAgg                                                                                                                                                     <  1:1646026/150‑1 (MQ=255)
     aaaCTCATACATTAACGTTTTCCCACGAAGTCTTTTTTCCGATAAGCCTTCGCACATATCGGTAAATAGCTTGCCTGCTCTTATTCTTTCGGTCATTGGCATGTTCAATGCGATCACTCCGTTATGATATGTTGGTCGGATAAGGCGCTc                                                                                                                                                >  1:1717595/1‑150 (MQ=255)
          aaTACATTAACGTTTTCCCACGAAGTCTTTTTTCCGGTAAGCCTTCGCACATATCGGTAAATAGCTTGCCTGCTCTTATTCTTTCGGTCATTGGCATGTTCAATGCGATCACTCCGTTATGATATGTTGGTCGGATAAGGCGCTCGcgcc                                                                                                                                           <  1:1675405/149‑1 (MQ=255)
            tACATTAACGTTTTCCCACGAAGTCTTTTTTCCGGTAAGCCTTCGCACATATCGGTAAATAGCTTGCCTGCTCTTATTCTTTCGGTCATTGGCATGTTCAATGCGATCACTCCGTTATGATATGTTGGTCGGATAAGGCGCTCGcgccgc                                                                                                                                         >  2:1455347/1‑150 (MQ=255)
             aCATTAACGTTTTCCCACGAAGTCTTTTTTCCGGTAAGCCTTCGCACATATCGGTAAATAGCTTGCCTGCTCTTATTCTTTCGGTCATTGGCATGTTCAATGCGATCACTCCGTTATGATATGTTGGTCGGATAAGGCGCTCGCGCCGCa                                                                                                                                        <  1:1271605/150‑1 (MQ=255)
             aCATTAACGTTTTCCCACGAAGTCTTTTTTCCGGTAAGCCTTCGCACATATCGGTAAATAGCTTGCCTGCTCTTATTCTTTCGGTCATTGGCATGTTCAATGCGATCACTCCGTTATGATATGTTGGTCGGATAAGGCGCTCGCGCCGCa                                                                                                                                        >  2:1546107/1‑150 (MQ=255)
               aTTAACGTTTTCCCACGAAGTCTTTTTTCCGGTTAGCCTTCGCACATATCGGTAAATAGCTTGCCTGCTCTTATTCTTTCGGTCATTGGCATGTTCAATGCGATCACTCCGTTATGATATGTTGGTCGGATGCGGCGCGAGcgccttatc                                                                                                                                      <  1:1096572/150‑6 (MQ=255)
               aTTAACGTTTTCCCACGAAGTCTTTTTTCCGGTAAGCCTTCGCACATATCGGTAAATAGCTTGCCTGCTCTTATTCTTTCGGTCATTGGCATGTTCAATGCGATCACTCCGTTATGATATGTTGGTCGGATGCGGCGCGAGcgccttatc                                                                                                                                      <  1:1447219/150‑6 (MQ=255)
               aTTAACGTTTTCCCACGAAGTCTTTTTTCCGGTAAGCCTTCGCACATATCGGTAAATAGCTTGCCTGCTCTTATTCTTTCGGTCATTGGCATGTTCAATGCGATCACTCCGTTATGATATGTTGGTCGGATAAGGCGCGAGcgccttatc                                                                                                                                      <  1:404039/150‑6 (MQ=255)
                   aCGTTTTCCCACGAAGTCTTTTTTCCGGTAAGCCTTCGCACATATCGGTAAATAGCTTGCCTGCTCTTATTCTTTCGGTCATTGGCATGTTCAATGCGATCACTCCGTTATGATATGTTGGTCGGATGCGGCGCGAGCGCCTTATCCGAc                                                                                                                                  <  1:109466/150‑1 (MQ=255)
                   aCGTTTTCCCACGAAGTCTTTTTTCCGGTAAGCCTTCGCACATATCGGTAAATAGCTTGCCTGCTCTTATTCTTTCGGTCATTGGCATGTTCAATGCGATCACTCCGTTATGATATGTTGGTCGGATAAGGCGCTCGCGCCGCATCCGAc                                                                                                                                  >  2:705053/1‑150 (MQ=255)
                    cGTTTTCCCACGAAGTCTTTTTTCCGGTAAGCCTTCGCACATATCGGTAAATAGCTTGCCTGCTCTTATTCTTTCGGTCATTGGCATGTTCAATGCGATCACTCCGTTATGATATGTTGGTCGGATAAGGCGCTCGCGCCGCATCCGACa                                                                                                                                 >  2:1362214/1‑150 (MQ=255)
                    cGTTTTCCCACGAAGTCTTTTTTCCGGTAAGCCTTCGCACATATCGGTAAATAGCTTGACTGCTCTTATTCTTTCGGTCATTGGCATGTTCAATGCGATCACTCCGTTATGATATGTTGGTCGGATAAGGCGCTCGCGCCGCATCCGACa                                                                                                                                 >  2:102235/1‑150 (MQ=255)
                          cccACGAAGTCTTTTTTCCGGTAAGCCTTCGCACATATCGGTAAATAGCTTGCCTGCTCTTATTCTTTCGGTCATTGGCATGTTCAATGCGATCACTCCGTTATGATATGTTGGTCGGATAAGGCGCTCGCGCCGCATCCGACattgatt                                                                                                                           <  2:59008/150‑1 (MQ=255)
                           ccACGAAGTCTTTTTTCCGGTAAGCCTTCGCACATATCGGTAAATAGCTTGCCTGCTCTTATTCTTTCGGTCATTGGCATGTTCAATGCGATCACTCCGTTATGATATGTTGGTCGGATAAGGCGCTCGCGCCGCATCCGACattgattg                                                                                                                          >  2:921511/1‑150 (MQ=255)
                           ccACGAACTCTATTTTCCGGTAAGCCTTCGAACATAACGGTAAATAGCTTGCCTGCTATTATTCTTTCGGTCATTGGCATGTTCAATGCGATCACTCCGTTATGATATGTTGGTCGGATAAGGCGCTCGCGCCGCATCCGACattgattg                                                                                                                          >  2:200176/1‑150 (MQ=255)
                            cACGAAGTCTTTTTTCCGGTAAGCCTTCGCACATTTCGGTAAATAGCTTGCCTGCTCTTATTCTTTCGGTCATTGGCATGTTCAATGCGATCACTCCGTTTTGATATGTTGGTCGGATAAGGCGCTCGCGCCGCATCCGCCATTGATTGc                                                                                                                         >  2:1096744/1‑150 (MQ=255)
                            cACGAAGTCTTTTTTCCGGTAAGCCTTCGCACATATCGGTAAATAGCTTGCCTGCTCTTATTCTTTCGGTCATTGGCATGTTCAATTCGATCACTCCGTTATGATATGTTGGTCGGATAAGGCGCTCGCGCCGCATCCGACATAGATTGc                                                                                                                         >  2:888717/1‑150 (MQ=255)
                            cACGAAGTCTTTTTTCCGGTAAGCCTTCGCACATATCGGTAAATAGCTTGCCTGCTCTTATTCTTTCGGTCATTGGCATGTTCAATGCGATCACTCCGTTATGATATGTTGGTCGGATAAGGCGCTCGCGCCGCATCCGACATTGATTGc                                                                                                                         >  2:675234/1‑150 (MQ=255)
                            cACGAAGTCTTTTTTCCGGTAAGCCTTCGCACATATCGGTAAATAGCTTGCCTGCTCTTATTCTTTCGGTCATTGGCATGTTCAATGCGATCACTCCGTTATGATATGTTGGTCGGATAAGGCGCTCGCGCCGCATCCGACATTGATTGc                                                                                                                         >  1:666801/1‑150 (MQ=255)
                             aCGAAGTCTTTTTTCCGGTAAGCCTTCGCACATATCGGTAAATAGCTTGCCTGCTCTTATTCTTTCGGTCATTGGCATGTTCAATGCGATCACTCCGTTATGATATGTTGGTCGGATAAGGCGCTCGCGCCGCATCCGACATTGATTGCt                                                                                                                        <  2:1291533/150‑1 (MQ=255)
                             aCGAAGTCTTTTTTCCGGTAAGCCTTCGCACATATCGGTAAATAGCTTGCCTGCTCTTATTCTTTCGGTCATTGGCATGTTCAATGCGATCACTCCGTTATGATATGTTGGTCGGATAAGGCGCTCGCGCCGCATCCGACATTGATTGCt                                                                                                                        <  2:833088/150‑1 (MQ=255)
                             aCGAAGTCGTTATTCCGGTAAGCCTTCGCACATATCCGGAAATAGCTTGCCTGATCTTATTCTTTCGGTCATTGGCATGTTCAATGCGATCACTCCGTTATGATATGTTGGTCGGATAAGTCGCTCGCGCCGCATCCGACATTGATTGCt                                                                                                                        <  1:310584/150‑1 (MQ=255)
                                 aGTCTTTTTTCCGGTAAGCCTTCGCACATATCGGTAAATAGCTTGCCTGCTCTTATTCTTTCGGTCATTGGCATGTTCAATGCGATCACTCCGTTATGATATGTTGGTCGGATAAGGCGCTCGCGCCGCATCCGACATTGATTGCTTAAg                                                                                                                    <  2:1151105/150‑1 (MQ=255)
                                    cTTTTTTCCGGTAAGCCTTCGCACATATCGGTAAATAGCTTGCCTGCTCTTATTCTTTCGGTCATTGGCATGTTCAATGCGATCACTCCGTTATGATATGTTGGTCGGATAAGGCGCTCGCGCCGCATCCGACATTGATTGCTTAAGCGa                                                                                                                 >  2:1698694/1‑150 (MQ=255)
                                    cTTTTTTCCGGTAAGCCTTCGCACATATCGGTAAATAGCTTGCCTGCTCTTATTCTTTCGGTCATTGGCATGTTCAATGCGATCACTCCGTTATGATATGTTGGTCGGATAAGGCGCTCGCGCCGCATCCGACATTGATTGCTTAAGCGa                                                                                                                 >  2:1207781/1‑150 (MQ=255)
                                    cTTTTTTCCGGTAAGCCTTCGCACATATCGGTAAATAGCTTGCCTGCTCTTATTCTTTCGGTCATTGGCATGTTCAATACGATCACTCCGTTATGATATGTTGGTCGGATAAGGCGCTCGCGCCGCATCAGACATTGATTGATTAAGCGa                                                                                                                 >  2:102102/1‑150 (MQ=255)
                                    cTTTTTTCCGGTAAGCCTTCGCACATATCGGTAAATAGCTTGCCTGCTCTTATTCTTTCGGTCATTGGCATGATCAATGCGAACACACCGATATGATATGTTGGTCGTATAAGGCGCTCCAGCCGCAACCGACAATGATTCCTTAAGCGa                                                                                                                 >  2:334229/1‑150 (MQ=255)
                                      tttttCCGGTAAGCCTTCGCACATATCGGTAAATAGCTTGCCTGCTCTTATTCTTTCGGTCATTGGCATGTTCAATGCGATCACTCCGTTATGATATGTTGGTCGGATAAGGCGCTCGCGCCGCATCCGACATTGATTGCTTAAGCGACt                                                                                                               <  2:1341735/150‑1 (MQ=255)
                                         ttCCGGTAAGCCTTCGCACATATCGGTAAATAGCTTGCCTGCTCTTATTCTTTCGGTCATTGGCATGTTCAATGCGATCACTCCGTTATGATATGTTGGTCGGATAAGGCGCTCGCGCCGCATCCGACATTGATTGCTTAAGCGACttca                                                                                                            <  2:1229268/150‑1 (MQ=255)
                                              gTAAGCCTTCGCACATATCGGTAAATAGCTTGCCTGCTCTTATTCTTTCGGTCATTGGCATGTTCAATGCGATCACTCCGTTATGATATGTTGGTCGGATAAGGCGCTCGCGCCGCATCCGACATTGATTGCTTAAGCGACTTCATTCAc                                                                                                       <  1:624154/150‑1 (MQ=255)
                                                aaGCCTTCGCACATCTCGGTAAATAGCTTGCCTGCTCTTATTCTTTCGGTCATTGGCATGTGCAATGCGAGCACTCCGTTATGATATGTTGGTCGGATAAGGCGCTCGCGCCGCATCCGACATTGATTGCTTAAGCGACTTCATTCACCt                                                                                                     <  2:1094311/150‑1 (MQ=255)
                                                  gCCTTCGCACATATCGGTAAATAGCTTGCCTGCTCTTATTCTTTCGGTCATTGGCATGTTCAATGCGATCACTCCGTTATGATATGTTGGTCGGATAAGGCGCTCGCGCCGCATCCGACATTGATTGCTTAAGCGACTTCATTCACCTga                                                                                                   >  1:769776/1‑150 (MQ=255)
                                                    cTTCGCACATATCGGTAAATAGCTTGCCTGCTCTTATTCTTTCGGTCATTGGCATGTTCAATGCGATCACTCCGTTATGATATGTTGGTCGGATAAGGCGCTCGCGCCGCATCCGACATTGATTGCTTAAGCGACTTCATTCACCTgacg                                                                                                 >  1:419011/1‑150 (MQ=255)
                                                    cTTCGCACATATCGGTAAATAGCTTGCCTGCTCTTATTCTTTCGGTCATTGGCATGTTCAATGCGATCACTCCGTTATGATATGTTGGTCGGATAAGGCGCTCGCGCCGCATCAGACATTGATTGCTTAAGCGACTTCATTCACCTgacg                                                                                                 >  2:629383/1‑150 (MQ=255)
                                                        gCACATATCGGTAAATAGCTTGCCTGCTCTTATTCTTTCGGTCATTGGCATGTTCAATGCGATCACTCCGTTATGATATGTTGGTCGGATAAGGCGCTCGCGCCGCATCCGACATTGATTGCTTAAGCGACTTCATTCACCTGACGACgc                                                                                             <  1:818822/150‑1 (MQ=255)
                                                         cacaTATCGGTAAATAGCTTGCCTGCTCTTATTCTTTCGGTCATTGGCATGTTCAATGCGATCACTCCGTTATGATATGTTGGTCGGATAAGGCGCTCGCGCCGCATCCGACATTGATTGCTTAAGCGACTTCATTCACCTGACGACgca                                                                                            >  2:115494/1‑150 (MQ=255)
                                                           caTATCGGTAAATAGCTTGCCTGCTCTTATTCTTTCGGTCATTGGCATGTTCAATGCGATCACTCCGTTATGATATGTTGGTCGGATAAGGCGCTCGCGCCGCATCCGACATTGATTGCTTAAGCGACTTCATTCACCTGACGACgcagc                                                                                          <  1:1116725/150‑1 (MQ=255)
                                                             tatCGGTAAATAGCTTGCCTGCTCTTATTCTTTCGGTCATTGGCATGTTCAATGCGATCACTCCGTTATGATATGTTGGTCGGATAAGGCGCTCGCGCCGCATCCGACATTGATTGCTTAAGCGACTTCATTCACCTGACGACGcagcag                                                                                        >  1:448717/1‑150 (MQ=255)
                                                               tCGGTAAATAGCTTGCCTGCTCTTATTCTTTCCGTCATTGGCATGTTCAATGCGATCACTCCGTTATGATATGTTGGTCGGATAAGGCGCTCGCGCAGCATCCGACATTGATTGCTTAAGCGACATCATTCACCTGACGACGCAGCAggg                                                                                      >  1:1723926/1‑150 (MQ=255)
                                                                    aaaTAGCTTGCCTGCTCTTATTCTTTCGGTCATTGGCATGTTCAATGCGATCACTCCGTTATGATATGTTGGTCGGATAAGGCGCTCGCGCCGCATCCGACATTGATTGCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGc                                                                                 >  2:633886/1‑150 (MQ=255)
                                                                      aTAGCTTGCCTGCTCTTATTCTTTCGGTCATTGGCATGTTCAATGCGATCACACCGTTATGATATGTTGGTCGGATAAGGCGCTCGCGCCGCATCCGACATTGATTGCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCgg                                                                               <  1:273128/150‑1 (MQ=255)
                                                                              ccTGCTATTAATCTTTCGGTCATTGGCAGGTTCACTGCGAGCACTCCGTTATGATTTGTGGATCGGATAAGGCGCTCGCGCCGCATCCGACATTGATTGCTTAAGCGACGTCATTCACCTGACGACGCAGCAGGGAAAGCgggccggggc                                                                       <  2:1086818/150‑1 (MQ=255)
                                                                               cTGCTCTTATTCTTTCGGTCATTGGCATGTTCACTGCGATCACTCCGTTATGATATGTTGGTCGGATAAGGCGCTCGCGCCGCATCCGACATTGATTGCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCgggccggggcc                                                                      <  1:1681392/150‑1 (MQ=255)
                                                                               cTGCTCTTATTCTTTCGGTCATTGGCATGTTCAATGCGATCACTCCGTTATGATATGTTGGTCGGATAAGGCGCTCGCGCCGCATCCGACATTGATTGCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCgggccggggcc                                                                      <  1:72907/150‑1 (MQ=255)
                                                                               cTGCTCTTATTCTTTCGGTCATTGGCATGTTCAATGCGATCACTCCGTTATGATATGTTGGTCGGATAAGGCGCTCGCGCCGCATCCGACATTGATTGCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCgggccggggcc                                                                      >  2:1185912/1‑150 (MQ=255)
                                                                                 gCTCTTATTCTTTCGGTCATTGGCATGTTCAATGCGATCACTCCGTTATGATATGTTGGTCGGATAAGGCGCTCGCGCCGCATCCGACATTGATTGCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGc                                                                    <  2:769776/150‑1 (MQ=255)
                                                                                  ctctTATTCTTTCGGTCATTGGCATGTTCAATGCGATCACTCCGTTATGATATGTTGGTCGGATAAGGCGCTCGCGCCGCATCCGACATTGATTGCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCt                                                                   >  2:145528/1‑150 (MQ=255)
                                                                                  ctctTATTCTTTCGGTCATTGGCATGTTCAATGCGAGCACTCCGTTATGATATGTTGGTCGGATACGGCGCTCGCGCCGCATCCGACATTGATTGCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCt                                                                   <  1:1734089/150‑1 (MQ=255)
                                                                                   tctTATTCTTTCGGTCATTGGCATGTTCAATGCGATCACTCCGTTATGATATGTTGGTCGGATAAGGCGCTCGAGCCGCATCCGACATTGATTGCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTa                                                                  >  1:54554/1‑150 (MQ=255)
                                                                                    ctTATTCTTTCGGTCATTGGCATGTTCAATGCGATCACTCCGTTATGATATGTTGGTCGGATAAGGCGCTCGCGCCGCATCCGACATTGATTGCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTaa                                                                 >  2:540100/1‑150 (MQ=255)
                                                                                     ttATTCTTTCGGTCATTGGCATGTTCAATGCGATCACTCCGTTATGATATGTTGGTCGGATAAGGCGCTCGCGCCGCATCCGACATTGATTGCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGACGGGGCCGCTAAg                                                                <  1:1584612/150‑1 (MQ=255)
                                                                                       cttCTTTCGGTCATTGTCAGGTTCAATGCGCTCACTCCGTTATGATCTGTTGGTCGGATACGGCGATCGCGCCGCATCCGACATTGATTGCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAACGCGGGCCGGGGCCGCTAAGCg                                                              <  1:1060788/149‑1 (MQ=255)
                                                                                       aTTCTTTCGGTCATTGGCATGTTCAATGCGATCACTCCGTTATGATATGTTGGTCGGATAAGGCGCTCGCGCCGCATCCGACATTGATTGCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCg                                                              <  2:1210735/150‑1 (MQ=255)
                                                                                       aTTCTTTCGGTCATTGGCATGTTCAATGCGATCACTCCGTTATGATATGTAGGTCGGATAAGGCGCTCGCGCCGCATCCGACATTGATTGCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGCCCCGGCCCGCtttccc                                                              <  2:848063/150‑6 (MQ=255)
                                                                                        ttCTTTCGGTCATTGGCATGTTCAATGCGATCACTCCGTTATGATATGTTGGTCGGATAAGGCGCTCGCGCCGCATCCGACATTGATTGCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGt                                                             <  2:787085/150‑1 (MQ=255)
                                                                                         tCTTTCGGTCATTGGCATGTTCAATGCGATCACTCCGTTATGATATGTTGGTCGGATAAGGCGCTCGCGCCGCATCCGACATTGATTGCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTg                                                            <  2:154982/150‑1 (MQ=255)
                                                                                          cTTTCGGTCATTGGCATGTTCAATGCGATCACTCCGTTATGATATGTTGGTCGGATAAGGCGCTCGCGCCGCATCCGACATTGATTGCTTAAGCGCCTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGa                                                           >  1:73058/1‑150 (MQ=255)
                                                                                          cTTTCGGTCATTGGCATGTTCAATGCGATCACTCCGTTATCATATGTTGGTCGGATAAGGCGCTCGCGCCGCATCCGACATTGATTGCTTAAGCGACTTCATTCACCTGACGACGCAGAAGGGAAAGCGCTCCCGGGCCGCTAAGCgaga                                                           >  1:220366/1‑150 (MQ=255)
                                                                                           tttCGGTCATTGGCATGTTCAATGCGATCACTCCGTTATGATATGTTGGTCGGATAAGGCGCTCGCGCCACATCCGACATTGATTGCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGaa                                                          <  1:707352/150‑1 (MQ=255)
                                                                                           tttCGGTCATTGGCATGGGCAATGCGATCACTCCGTTATGATATGGTGGTCGGATAAGGCGCTCGCGCCGCATCCGACATTGATTGCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGaa                                                          <  1:240314/150‑1 (MQ=255)
                                                                                            ttCGGTCATTGGCATGTTCAATGCGATCACTCCGTTATGATATGTTGGTCGGATAAGGCGCTCGCGCCGCATCCGACATTGATTGCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAac                                                         <  1:758208/150‑1 (MQ=255)
                                                                                               ggTCATTGGCATGTTCAATGCGATCACTCCGTTATGATATGTTGGTCGGATAAGGCGCTCGCGCCGCATCCGACATTGATTGGTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGATAAGCGTGAacacc                                                      >  1:1609347/1‑149 (MQ=255)
                                                                                               ggTCATTGGCATGTTCAATGCGATCACTCCGTTATGATATGTTGGTCGGATAAGGCGCTCGCGCCGCATCCGACATTGATTGCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACg                                                      <  1:533297/150‑1 (MQ=255)
                                                                                               ggTCATTGGCATGTTCAATGCGATCACTCCGTTATGATATGTTGGTCGGATAAGGCGCTCGCGCCGCATCCGACATTGATTGCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACg                                                      <  1:526540/150‑1 (MQ=255)
                                                                                                   aTTGGCATGTTCAATGCGATCACTCCGTTATGATATGTTGGTCGGATAAGGCGCTCGCGCCGCATCCGACATTGATTGCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGaaa                                                  >  1:328801/1‑150 (MQ=255)
                                                                                                      ggCATGTTCAATGCGATCACTCCGTTATGATATGTTGGTCGGATAAGGCGCTCGCGCCGCATCCGACATTGATTGCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTa                                               >  2:425201/1‑150 (MQ=255)
                                                                                                      ggCATGTTCAATGCGATCACTCCGTTATGATATGTTGGTCGGATAAGGCGCTCGCGCCGCATCCGACATTGATTGCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTa                                               >  1:1696530/1‑150 (MQ=255)
                                                                                                          tgatCAATGCGATCACTCCGTTATGATATGTTGGTCGGATAAGGCGCTCGCGCCGCATCCGACATTGATTGCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTAAGGt                                           <  1:629383/147‑1 (MQ=255)
                                                                                                          tGTTCAATGCGATCACTCCGTTATGATATGTTGGTCGGATAAGGCGCTCGCGCCGCATCCGACATTGATTGCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTAAGGt                                           <  2:1717595/150‑1 (MQ=255)
                                                                                                          tGTTCAATGCGATCACTCCGTTATGATATGTTGGTCGGATAAGGCGCTCGCGCCGCATCCGACATTGATTGCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTAAGGt                                           >  1:1686947/1‑150 (MQ=255)
                                                                                                           gTTCAATGCGATCACTCCGTTATGATATGTTGGTCGGATAAGGCGCTCGCGCCGCATCCGACATTGATTGCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTAAGGTg                                          >  1:873534/1‑150 (MQ=255)
                                                                                                           gTTCAATGCGATCACTCCGTTATGATATGTTGGTCGGATAAGGCGCTCGCGCCGCATCCGACATTGATTGCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTAAGGTg                                          >  1:1347687/1‑150 (MQ=255)
                                                                                                            ttCAATGCGATCACTCCGTTATGATATGTTGGTCGGATAAGGCGCTCGCGCCGCATCCGACATTGATTGCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTAAGGTGa                                         <  1:158565/150‑1 (MQ=255)
                                                                                                              cAATGCGATCACTCCGTTATGATATGTTGGTCGGATAAGGCGCTCGCGCCGCATCCGACATTGATTGCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTAAGGTGAAg                                       >  2:250677/1‑150 (MQ=255)
                                                                                                               aaTGCGATCACTCCGTTATGATATGTTGGTCGGATAAGGCGCGCGCGCCGCATCCGACATTGATTGCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTAAGGTGAAGc                                      <  2:668982/150‑1 (MQ=255)
                                                                                                                aTGCGATCACTCCGTTATGATATGTTGGTCGGATAAGGCGCTCGCGCCGCATCCGACATTGATTGCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATCAAGGTGAAGcc                                     >  2:1068577/1‑150 (MQ=255)
                                                                                                                aTGCGATCACTCCGTTATGATATGTTGGTCGGATAAGGCGCTCGCGCCGCATCCGACATTGATTGCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATCAAGGTGAAGcc                                     >  2:683829/1‑150 (MQ=255)
                                                                                                                 tGCGATCACTCCGTTATGATATGTTGGTCGGATAAGGCGCTCGCGCCGCATCCGACATTGATTGCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTAAGGTGAAGccc                                    >  2:1129206/1‑150 (MQ=255)
                                                                                                                   cGATCACTCCGTTATGATATGTTGGTCGGATAAGGCGCTCGCGCCGCATCCGACATTGATTGCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTAAGGTGAAGCCCAg                                  >  2:1639503/1‑150 (MQ=255)
                                                                                                                         cTCCGTTATGATATGTTGGTCGGATAAGGCGCTCGCGCCGCATCCGACATTGATTGCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTAAGGTGAAGCCCAGCGccac                            <  1:153715/150‑1 (MQ=255)
                                                                                                                          tCCGTTATGATATGTTGGTCGGATAAGGCGCTCGCGCCGCATCCGACATTGATTGCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTAAGGTGAAGCCCAGCGccacc                           <  1:1455347/150‑1 (MQ=255)
                                                                                                                          tCCGTTATGATATGTTGGTCGGATAAGGCGCTCGCGCCGCATCCGACATTGATTGCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTAAGGTGAAGCCCAGCGccacc                           <  2:1285770/150‑1 (MQ=255)
                                                                                                                           ccGTTATGATATGTTGGTCGGATAAGGCGCTCGCGCCGCATCCGACATTGATTGCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTAAGGTGAAGCCCAGCGccacca                          <  1:675234/150‑1 (MQ=255)
                                                                                                                            cGTTATGATATGTTGGTCGGATAAGGCGCTCGCGCCGCATCCGACATTGATTGCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTAAGGTGAAGCCCAGCGCCACCAg                         <  1:34161/150‑1 (MQ=255)
                                                                                                                            cGTTATGATATGTTGGTCGGATAAGGCGCTCGCGCCGCAGCCGACATTGATTGCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTAAGGTGAAGCCCAGCGCCACCAg                         <  1:317279/150‑1 (MQ=255)
                                                                                                                                aTGATATGTTGGTCGGATAAGGCGCTCGCGCCGCATCCGACATTGATTGCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTAAGGTGAAGCCCAGCGCCACCAGAccc                     <  1:1533041/150‑1 (MQ=255)
                                                                                                                                 tGATATGTTGGTCGGATAAGGCGCTCGCGCCGCATCCGACATTGATTGCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTAAGGTGAAGCCCCGCGCCACCAGCCCCa                    >  1:1014681/1‑150 (MQ=255)
                                                                                                                                 tGATATGTTGGGCGGATAAGGCTCGCGCGCCGCATCCGACATTGATTGCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTAAGGTGAAGCCCAGCGCCACCAGACCCa                    <  2:830150/150‑1 (MQ=255)
                                                                                                                                   atatGTTGGTCGGATAAGGCGCTCGCGCCGCATCCGACATTGATTGCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTAAGGTGAAGCCCAGCGCCACCAGACCCAGc                  <  2:698475/150‑1 (MQ=255)
                                                                                                                                    tatGTTGGTCGGATAAGGCGCTCGCGCCGCATCCGACATTGATTGCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTAAGGTGAAGCCCAGCGCCACCAGACCCAGCa                 <  1:79127/150‑1 (MQ=255)
                                                                                                                                      tcaatgTCGGATGCGGCGCGAGCGCCTTATCCGACATTGATTGCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTAAGGTGAAGCCCAGCGCCACCAGACCCAGCAcc               >  1:1642462/6‑150 (MQ=255)
                                                                                                                                      tGTTGGGCGGATCAGGCGCTCGCGCCGCATCCGACATTGATTGCTTAAGCGACTTCATACCCCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTAAGGTGAAGCCCAGCGCCACCAGACCCAGCAcc               <  2:47756/150‑1 (MQ=255)
                                                                                                                                        aatgTCGGATGCGGCGCGAGCGCCTTATCCGACGTAGATTGCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATAAAGGTGAAGCCCAGCGCCACCAGACCCAGCACCAg             >  1:1160218/4‑150 (MQ=255)
                                                                                                                                          ggTCGGATAAGGCGCTCGCGCCGCATCCGACATTGATTGCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGc                                                                    <  2:1130580/93‑1 (MQ=255)
                                                                                                                                          ggTCGGATAAGGCGCTCGCGCCGCATCCGACATTGATTGCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGc                                                                    >  1:1130580/1‑93 (MQ=255)
                                                                                                                                           gTCGGATAAGGCGCTCGCGCCGCATACGACATTGATTGCTTAAGCGACTTCATTCTCCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTAAGGTGAAGCCCAGCTCCACCAGACCCAGCACCAGATa          >  2:610482/1‑150 (MQ=255)
                                                                                                                                                   aGGCGCTCGCGCCGCATCCGACATTGATTGCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGTAATTAAGGTGAAGCCCAGCGCCCCCAGACCCAGCACCAGATAAGCGCCCt  >  2:1076782/1‑150 (MQ=255)
                                                                                                                                                   aGGCGCTCGCGCCGCATCCGACATTGATTGCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTAAGGTGAAGCCCAGCGCCACCAGACCCAGCACCAGATAAGCGCCCt  >  2:1177835/1‑150 (MQ=255)
                                                                                                                                                   aGGCGCTCGCGCCGCATCCGACATTGATTGCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTAAGGTGAAGCCCAGCGCCACCAGACCCAGCACCAGATAAGCGCCCt  >  2:134993/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
TGATTAAACTCATACATTAACGTTTTCCCACGAAGTCTTTTTTCCGGTAAGCCTTCGCACATATCGGTAAATAGCTTGCCTGCTCTTATTCTTTCGGTCATTGGCATGTTCAATGCGATCACTCCGTTATGATATGTTGGTCGGATAAGGCGCTCGCGCCGCATCCGACATTGATTGCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTAAGGTGAAGCCCAGCGCCACCAGACCCAGCACCAGATAAGCGCCCT  >  NC_000913/361748‑362044

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: