Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,316,897 A→T 7.0% intergenic (+70/+130) rcsB → / ← rcsC response regulator in two‑component regulatory system with RcsC and YojN/hybrid sensory kinase in two‑component regulatory system with RcsB and YojN

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,316,8970AT 7.0% 363.0 / 11.1 124intergenic (+70/+130)rcsB/rcsCresponse regulator in two‑component regulatory system with RcsC and YojN/hybrid sensory kinase in two‑component regulatory system with RcsB and YojN
Reads supporting (aligned to +/- strand):  ref base A (71/44);  new base T (4/5);  total (75/49)
Fisher's exact test for biased strand distribution p-value = 3.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.76e-01

Only 100 of 126 total aligned reads displayed.

CTGGGTGTCGAGAACGATATCGCCCTGCTGAATTATCTCTCTTCAGTGACCTTAAGTCCGGCAGATAAAGACTAATCACCTGTAGGCCAGATAAGACGCGTTAGTGTCTTATCTGGCATTTGCACCGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGATTCCCATTATTTCAACAAATTACATTAAAGTAGGCCAGATAAGACGCGTCAGCGTCGCATCTGGCATTTGCACTGAATGCCGGATGCGGCGTAAACGCCTTATCCGTCCTACGAATCCCGCGATTTCC  >  NC_000913/2316753‑2317046
                                                                                                                                                |                                                                                                                                                     
cTGGGTGTCGAGAACGATATCGCCCTGCTGAATTATCTCTCTTCAGTGACCTTAAGTCCGGCAGATAAAGACTAATCACCTGTAGGCCAGATAAGACGCGTTAGTGTCTTATCTGGCATTTGCACCGATTGCCGGATGCGGAGTAAACGc                                                                                                                                                  >  2:1665297/1‑150 (MQ=255)
  gggTGTCTAGAACGATATCGCCCTGCTGAATTATCTCTCTTCAGTGACCTTAAGTCCGGCAGATAAAGACTAATCACCTGTAGGCCAGATAAGACGCGTTAGTGTCTTATCTGGCATTTGCACCGATTGCCGGATGCGGCGTTATCGCCt                                                                                                                                                <  2:1226334/150‑1 (MQ=255)
      gtCGAGAACGATATCGCCCTGCTGAATTATCTCTCTTCAGTGACCTTAAGTCCGGCAGATAAAGACTATCACCTGTAGGCCAGATAAGACGCGTTAGTGTCTTATCTGGCATTTGCACCGATTGCCGGATGCGGCGTAAACGCCTTATcc                                                                                                                                           <  1:201352/150‑1 (MQ=255)
        cGAGAACGATATCGCCCTGCTGAATTATCTCTCTTCAGTGACCTTAAGTCCGGCAGATAAAGACTAATCACCTGTAGGCCAGATAAGACGCGTTAGTGTCTTATCTGGCATTTGCACCGATTGCCGGATGCGGCGTAAACGCCTTATCCg                                                                                                                                          >  1:550990/1‑150 (MQ=255)
        cGAGAACGATATCGCCCTGCTGAATTATCTCTCTTCAGTGACCTTAAGTCCGGCAGATAAAGACTAATCACCTGTAGGCCAGATAAGACGCGTTAGTGTCTTATCTGGCATTTGCACCGAATGCCGGATGCGGCGTAAACGCCTTATCCg                                                                                                                                          >  2:447876/1‑150 (MQ=255)
         gagaACGATATCGCCCTGCTGAATTATCTCTCTTCAGTGACCTTAAGTCCGGCAGATAAAGACTAATCACCTGTAGGCCAGATAAGACGCGTTAGTGTCTTATCTGGCATTTGCACCGATTGCCGGATGCGGCGTAAACGCCTTATCCgg                                                                                                                                         >  2:919129/1‑150 (MQ=255)
          agaACGATATCGCCCTGCTGAATTATCTCTCTTCAGTGACCTTTAGTCCGGCAGATAAAGACTAATCACCTGTAGGCCAGATAAGACGCGTTAGTGTCTTATCTGGCATTTGCACCGATTGCCGGATAAGGCGTTTACGCCGCATCCGGc                                                                                                                                        <  1:1639681/150‑1 (MQ=38)
          agaACGATATCGCCCTGCTGAATTATCTCTCTTCAGTGACCTTAAGTCCGGCAGATAAAGACTAATCACCTGTAGGCCAGATAAGACGCGTTAGTGTCTTATCTGGCATTTGCACCGATTGCCGGATGCGGCGTAAACGCCTTATCCGGc                                                                                                                                        >  2:643589/1‑150 (MQ=255)
            aaCGATATCGCCCTGCTGAATTATCTCTCTTCAGTGACCTTAAGTCCGGCAGATAAAGACTAATCACCTGTAGGCCAGATAAGACGCGTTAGTGTCTTATCTGGCATTTGCACCGATTGCCGGATGCGGCGTAAACGCCTTTTCCGGCCt                                                                                                                                      >  2:694322/1‑150 (MQ=255)
            aaCGATATCGCCCTGCTGAATTATCTCTCTTCAGTGACCTTAAGTCCGGCAGATAAAGACTAATCACCTGTAGGCCAGATAAGACGCGTTAGTGTCTTATCTGGCATTTGCACCGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCt                                                                                                                                      >  2:616424/1‑150 (MQ=255)
            aaCGATATCGCCCTGCTGAATTATCTCTCTTCAGTGACCTTAAGTCCGGCAGATAAAGACTAATCACCTGTAGGCCAGATAAGACGCGTTAGTGTCTTATCTGGCATTTGCACCGATTGCCGGATAAGGCGTTTACGCCGCATCCGGcaa                                                                                                                                      <  2:146991/150‑3 (MQ=255)
             aCGATATCGCCCTGCTGAATTATCTCTCTTCAGTGACCTTAAGTCTGGCAGATAAAGAGTAATCACCTGTAGGCCAGATAAGACGCGTTAGTGTCTTATCTGGCATTTGCACCGATTGCCGGATAAGGCGTTTACGCCGCATCCGGcaat                                                                                                                                     <  1:518765/150‑4 (MQ=255)
             aCGATATCGCCCTGCTGAATTATCTCTCTTCAGTGACCTTAAGTCCGGCAGATAAAGACTAATCACCTGTAGGCCAGATAAGACGCGTTAGTGTCTTATCTGGCATTTGCACCGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTa                                                                                                                                     >  1:299572/1‑150 (MQ=255)
              cGATATCGCCCTGCTGAATTATCTCTCTTCAGTGACCTTAAGTCCGGCAGATAAAGACTAATCACCTGTAGGCCAGATAAGACGCGTTAGTGTCTTATCTGGCATTTGCACCGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTAc                                                                                                                                    <  2:785917/150‑1 (MQ=255)
              cGATATCGCCCTGCTGAATTATCTCTCTTCAGTGACCTTAAGTCCGGCAGATAAAGACTAATCACCTGTAGGCCAGATAAGACGCGTTAGTGTCTTATCTGGCATTTGCACCGATTGCCGGATGCGGAGTAAACGCCTTATCCGGCCTAc                                                                                                                                    >  1:95355/1‑150 (MQ=255)
                atatCGCCCTGCTGAATTATCTCTCTTCAGTGACCTTAAGTCCGGCAGATAAAGACTAATCACCTGTAGGCCAGATAAGACGCGTTAGTGTCTTATCTGGCATTTGCACCGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGa                                                                                                                                  >  2:708949/1‑150 (MQ=255)
                   acGCCCTTCTGAATTATCTCTCTTCTGTGAACTTAAGTCCTGGAGATATAGTCTATTCACGTGTAGGCCAGATTTTACGCGTTAGTGTCTTAGCTGGCATTTGCACCGATAGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGATTc                                                                                                                               <  1:238954/149‑1 (MQ=25)
                    cGCCCTGCTGAATTATCTCTCTTCAGTGACCTTAAGTCCGGCAGATAAAGACTAATCACCTGTAGGCCAGATAAGACGCGTTAGTGTCTTATCTGGCATTTGCACCGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGATTcc                                                                                                                              >  2:1057063/1‑150 (MQ=255)
                     gCCCTGCTGAATTATCTCTCTTCAGTGACCTTAAGTCCGGCAGATAAAGACTAATCACCTGTAGGCCAGATAAGACGCGTTAGTGTCTTATCTGGCATTTGCACCGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGATTccc                                                                                                                             >  1:750077/1‑150 (MQ=255)
                      cccTGCTGAATTATCTCTCTTCAGTGACCTTAAGTCCGGCAGATAAAGACTAATCACCTGTAGGCCAGATAAGACGCGTTAGTGTCTTATCTGGCATTTGCACCGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGATTCCCa                                                                                                                            <  2:531009/150‑1 (MQ=255)
                        ctgctgAATTATCTCTCTTCAGTGACCTTAAGTCCGGCAGATAAAGACTAATCACCTGTAGGCCAGATAAGACGCGTTAGTGTCTTATCTGGCATTTGCACCGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGATTCCCatt                                                                                                                          <  2:726058/150‑1 (MQ=255)
                        ctgctgAATTATCTCTCTTCAGTGACCTTAAGTCCGGCAGATAAAGACTAATCACCTGTAGGCCAGATAAGACGCGTTAGTGTCTTATCTGGCATTTGCACCGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGATTCCCatt                                                                                                                          >  1:1688279/1‑150 (MQ=255)
                         tgctgAATTATCTCTCTTCAGTGACCTTAAGTCCGGCAGATAAAGACTAATCACCTGTAGGCCAGATAAGACGCGTTAGTGTCTTATCTGGCATTTGCACCGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGATTCCCatta                                                                                                                         <  2:415202/150‑1 (MQ=255)
                          gctgAATTATCTCTCTTCAGTGACCTTAAGTCCGGCAGATAAAGACTAATCACCTGTAGGCCAGATAAGACGCGTTAGTGTCTTATCTGGCATTTGCACCGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGATTCCCattat                                                                                                                        <  1:598613/150‑1 (MQ=255)
                            tgAATTATCTCTCTTCAGTGACCTTAAGTCCGGCAGATAAAGACTAATCACCTGTAGGCCAGATAAGACGCGTTAGTGTCTTATCTGGCATTTGCACCGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGATTCCCATTAttt                                                                                                                      <  2:379595/150‑1 (MQ=255)
                              aaTTATCTCTCTTCAGTGACCTTAAGTCCGGCAGATAAAGACTAATCACCTGTAGGCCAGATAAGACGCGTTAGTGTCTTATCTGGCATTTGCACCGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGATTCCCATTATTTca                                                                                                                    >  2:172747/1‑150 (MQ=255)
                                 tATCTCTCTTCAGTGACCTTAAGTCCGGCAGATAAAGACTAATCACCTGTAGGCCAGATAAGACGCGTTAGTGTCTTATCTGGCATTTGCACCGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGATTCCCATTATTTcaaca                                                                                                                 >  2:1470735/1‑150 (MQ=255)
                                  aTCTCTCTTCAGTGACATTAAGTCCGGCAGATAAAGACTAATCACCTGTAGGCCAGATAAGACGCGTTAGTGTCTTATCTGGCATTTGCACCGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGATTCCCATTATATcaacaa                                                                                                                >  1:870514/1‑150 (MQ=255)
                                       tctTCAGTGACCTTAAGTCCGGCAGATAAAGACTAATCACCTGTAGGCCAGATAAGACGCGTTAGTGTCTTATCTGGCATTTGCACCGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGATTCCCATTATTTCAACAAATTAc                                                                                                           >  1:1043950/1‑150 (MQ=255)
                                            aGTGACCTTAAGTCCGGCAGATAAAGACTAATCACCTGTAGGCCAGATAAGACGCGTTAGTGTCTTATCTGGCATTTGCACCGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGATTCCCATTATTTCAACAAATTACATTaa                                                                                                      >  1:433219/1‑150 (MQ=255)
                                               gACCTTAAGTCCGGCAGATAAAGACTAATCACCTGTAGGCCAGATAAGACGCGTTAGTGTCTTATCTGGCATTTGCACCGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGATTCCCATTATTTCAACAAATTACATTAAAGt                                                                                                   >  1:1689823/1‑150 (MQ=255)
                                                   ttAAGTCCGGCAGATAAAGACTAATCACCTGTAGGCCAGATAAGACGCGTTAGTGTCTTATCTGGCATTTGCACCGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGATTCCCATTATTTCAACAAATTACATTAAAGTAGGc                                                                                               >  2:1552305/1‑150 (MQ=255)
                                                     aaGTCCGGCAGATAAAGACTAATCACCTGTAGGCCAGATAAGACGCGTTAGTGTCTTATCTGGCATTTGCACCGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGATTCCCATTATTTCAACAAATTACATTAAAGTAGGCCa                                                                                             <  2:429672/150‑1 (MQ=255)
                                                       gTCCGGCAGATAAAGACTAATCACCTGTAGGCCAGATAAGACGCGTTAGTGTCTTATCTGGCATTTGCCCCGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGATTCCCATTATTTCAACAAATTACATTAATGTAGGCCAGa                                                                                           >  2:269249/1‑150 (MQ=255)
                                                       gTCCGGCAGATAAAGACTAATCACCTGTAGGCCAGATAAGACGCGTTAGTGTCTTATCTGGCATTTGCACCGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGATTCCCATTATTTCAACAAATTACATTAAAGTAGGCCAGa                                                                                           >  1:384957/1‑150 (MQ=255)
                                                       gTCCGGCAGATAAAGACTAATCACCTGTAGGCCAGATAAGACGCGTTAGTGTCTTATCTGGCATTTGCACCGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGATTCCCATTATTTCAACAAATTACATTAAAGTAGGCCAGa                                                                                           >  2:62814/1‑150 (MQ=255)
                                                           ggCAGATAAAGACTAATCACCTGTAGGCCAGATAAGACGCGTTAGTGTCTTATCTGGCATTTGCACCGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGATTCCCATTATTTCAACAAATTACATTAAAGTAGGCCAGATAAg                                                                                       <  1:170624/150‑1 (MQ=255)
                                                              aGATAAAGACTAATCACCTGTAGGCCAGATAAGACGCGTTAGTGTCTTATCTGGCATTTGCACCGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGATTCCCATTATTTCAACAAATTACATTAAAGTAGGCCAGATAAGAcg                                                                                    >  2:237959/1‑150 (MQ=255)
                                                                aTAACGACTAATCACCTGTAGGCCAGATAAGACGCGTTAGTGTCTTATCTGGCATTTGCACCGATTGCCGGATGCGGCGTAAACGCCTTATCCAGCCTACGATTCCCATTATTTCAACAAATTACATTAAAGTAGGCCAGATAAGAcgcg                                                                                  <  2:858994/150‑1 (MQ=255)
                                                                aTAAAGACTAATCACCTGTAGGCCAGATAAGACGCGTTAGTGTCTTATCTGGCATTTGCACCGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGATTCCCATTATTTCAACAAATTACATTAAAGTAGGCCAGATAAGAcgcg                                                                                  <  1:1068133/150‑1 (MQ=255)
                                                                  aaaGACTAATCACCTGTAGGCCAGATAAGACGCGTTAGTGTCTTATCTGGCATTTGCACCGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGATTCCCATTATTTCAACAAATTACATTAAAGTAGGCCAGATAAGACGCGTc                                                                                >  2:750370/1‑150 (MQ=255)
                                                                       cTAATCACCTGTAGGCCAGATAAGACGCGTTAGTGTCTTATCTGGCATTTGCACCGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGATTCCCATTATTTCAACAAATTACATTAAAGTAGGCCAGATAAGACGCGTCAGCGt                                                                           >  2:586925/1‑150 (MQ=255)
                                                                          aTCACCTGTAGGCCAGATAAGACGCGTTAGTGTCTTATCTGGCATTTGCACCGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGATTCCCATTATTTCAACAAATTACATTAAAGTAGGCCAGATAAGACGCGTCAGCGTCGc                                                                        >  1:15484/1‑150 (MQ=255)
                                                                            cACCTGTAGGCCAGATAAGACGCGTTAGTGTCTTATCTGGCATTTGCACCGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGATTCCCATTATTTCAACAAATTACATTAAAGTAGGCCAGATAAGACGCGTCAGCGTCGCAt                                                                      >  1:1159632/1‑150 (MQ=255)
                                                                              ccTGTAGGCCAGATAAGACGCGTTAGTGTCTTATCTGGCATTTGCACCGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGATTCCCATTATTTCAACAAATTACATTAAAGTAGGCCAGATAAGACGCGTCAGCGTCGCATCt                                                                    >  2:511700/1‑150 (MQ=255)
                                                                                tGTAGGCCAGATAAGACGCGTTAGTGTCTTATCTGGCATTTGCACCGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGATTCCCATTATTTCAACAAATTACATTAAAGTAGGCCAGATAAGACGCGTCAGCGTCGCATCTgg                                                                  <  1:1151386/150‑1 (MQ=255)
                                                                                 gTAGGCCAGATAAGACGCGTTAGTGTCTTATCTGGCATTTGCACCGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGATTCCCATTATTTCAACAAATTACATTAAAGTAGGCCAGATGCGACGCTGACGCGTCTTATCTGGc                                                                 <  1:172747/150‑1 (MQ=25)
                                                                                   aGGCCAGATAAGACGCGTTAGTGTCTTATCTGGCATTTGCACCGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGATTCCCATTATTTCAACAAATTACATTAAAGTAGGCCAGATAAGACGCGTCAGCGTCGCATCTGGCAt                                                               >  2:1282848/1‑150 (MQ=255)
                                                                                    ggCCAGATAAGACGCGTTAGTGTCTTATCTGGCATTTGCACCGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGATTCCCATTATTTCAACAAATTACATTAAAGTAGGCCAGATAAGACGCGTCAGCGTCGCATCTGGCAtt                                                              >  1:482995/1‑150 (MQ=255)
                                                                                      ccAGATAAGACGCGTTAGTGTCTTATCTGGCATTTGCACCGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGATTCCCATTATTTCAACAAATTACATTAAAGTAGGCCAGATAAGACGCGTCAGCGTCGCATCTGGCATTTg                                                            >  2:1437714/1‑150 (MQ=255)
                                                                                      ccAGATAAGACGCGTTAGTGTCTTATCTGGCATTTGCACCGATTGCCGGATGCGGCGTAAACGCCTTATCAGGCCTACGATTCCCATTATTTCAACAAATTACATTAAAGTAGGCCAGATAAGACGCGTCAGCGTCGCATCTGGCATTTg                                                            >  1:1105731/1‑150 (MQ=255)
                                                                                      ccAGATAAGACACTAACGCGTCTTATCTGGCATTTGCACCGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGATTCCCATTATTTCAACAAATTACATTAAAGTAGGCCAGATAAGACGCGTCAGCGTCGCATCTGGCATTTg                                                            >  1:1599852/1‑150 (MQ=38)
                                                                                      ccAGATAAGACACTAACGCGTCTTATCTGGCATTTGCACCGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGATTCCCATTATTTCAACAAATTACATTAAAGTAGGCCAGATAAGACGCGTCAGCGTCCCATCTGGCATTTg                                                            >  2:136062/1‑150 (MQ=37)
                                                                                       cAGATAAGACGCGTTAGTGTCTTATCTGGCATTTGCACCGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGATTCCCATTATTTCAACAAATTACATTAAAGTAGGCCAGATAAGACGCGTCAGCGTCGCATCTGGCATTTGc                                                           >  2:1029941/1‑150 (MQ=255)
                                                                                       cAGATAAGACGCGTTAGTGTCTTATCTGGCATTTGCACCGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGATTCCCATTATTTCAACAAATTACATTAAAGTAGGCCAGATAAGACGCGTCAGCGACGCATCTGGCATTTga                                                           >  1:691430/1‑149 (MQ=255)
                                                                                       cAGATAAGACGCGTTAGTGTCTTATCTGGCATTTGCACCGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGATTCCCATAAATTCAACAAATTACATTAAAGTAGGCCAGATAAGACGCGTCAGAGTCCCATCTGGCAATTGc                                                           >  1:321612/1‑150 (MQ=37)
                                                                                        aGATAAGACGCGTTAGTGTCTTATCTGGCATTTGCACCGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGATTCCCATTATTTCAACAAATTACATTAAAGTAGGCCAGATAAGACGCGTCAGTGTCGCATCTGGCATTTGCa                                                          >  2:117073/1‑150 (MQ=255)
                                                                                        aGATAAGACGCGTTAGTGTCTTATCTGGCATTTGCACCGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGATTCCCATTATTTCAACAAATTACATTAAAGTAGGCCAGATAAGACGCGTCAGCGTCGCATCTTGCATTTGCa                                                          >  1:622752/1‑150 (MQ=255)
                                                                                         gATAAGACGCGTTAGTGTTTTATCTGGCATTTGCACCTATTGCCGGATGCGGCGTAAACGCCTTATCCGGCATACGATTCCCATTATTTCAACAAATTACATTAAAGTAGGCCAGATAAGACGCGTCAGCGTCGCATCTGGCATTTGCAc                                                         <  2:273481/150‑1 (MQ=255)
                                                                                         gATAAGACGCGTTAGTGTCTTATCTGGCATTTGCACCGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGATTCCCATTATTTCAACAAATTACATTAAAGTAGGCCAGATAAGACGCGTCAGCGTCGCATCTGGCATTTGCAc                                                         >  1:382221/1‑150 (MQ=255)
                                                                                           tAAGACGCGTTAGTGTCTTATCTGGCATTTGCACCGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGATTCCCATTATTTCAACAAATTACATTAAAGTAGGCCAGATAAGACGCGTCAGCGTCGCATCCGGCAttcagtgca                                                       <  1:1065362/150‑8 (MQ=255)
                                                                                            aaGACGCGTTAGTGTCTTATCTGGCATTTGCACCGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGATTCCCATTATTTCAACAAATTACATTAAAGTAGGCCAGATAAGACGCGTCAGCGTCGCATCTGGCATTTGCACTGa                                                      >  2:1585294/1‑150 (MQ=255)
                                                                                                cgcgTTAGTGTCTTATCTGGCATTTGCACCGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGATTCCCATTATTTCAACAAATTACATTAAAGTAGGCCAGATAAGACGCGTCAGCGTCGCATCTGGCATTTGCACTGAATGc                                                  >  1:1211271/1‑150 (MQ=255)
                                                                                                cgcgTTAGTGTCTTATCTGGCATTTGCACCGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGATTCCCATTATTTCAACAAATTACATTAAAGTAGGCCAGATAAGACGCGTCAGCGTCGCATCTGGCATTTGCACTGAATGc                                                  >  2:1109338/1‑150 (MQ=255)
                                                                                                      aGTGTCTTATCTGGCATTTGCACCGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGATTCCCATTATTTCAACAAATTACATTAAAGTAGGCCAGATAAGACGCGTCAGCGTCGCATCTGGCATTTGCACTGAATGCCGGATg                                            <  1:395116/150‑1 (MQ=255)
                                                                                                      aGTGTCTTATCTGGCATTTGCACCGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGATTCCCATTATTTCAACAAATTACATTAAAGTAGGCCAGATAAGACGCGTCAGCGTCGCATCTGGCATTTGCACTGAATGCCGGATg                                            >  1:881089/1‑150 (MQ=255)
                                                                                                       gtgtCTTATCTGGCATTTGCACCGATTGCCGGATGCGGCGTAAACGCCTTATACGGCCTACGATTCCCATTATATCAACAAATTAAATTAAAGTAGTCCAGCTAACACGAGTCAGCGTCGCATCTGTCATTTGCACTAAATGCCGCAtac                                           >  1:148756/1‑148 (MQ=37)
                                                                                                         gtCTTATCTGGCATTTGCACCGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGATTCCCATTATTTCAACAAATTACATTAAAGTAGGCCAGATAAGACGCGTCAGCGTCGCATCTGGCATTTGCACTGAATGCCGGATgcgg                                         <  1:511700/150‑1 (MQ=255)
                                                                                                            ttATCTGGCATTTGCACCGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGATTCCCATTATTTCAACAAATTACATTAAAGTAGGCCAGATAAGACGCGTCAGCGTCGCATCTGGCATTTGCACTGAATGCCGGATGCGGCGt                                      <  1:1200120/150‑1 (MQ=255)
                                                                                                                cTGGCATTTGCACCGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGATTCCCATTATTTCAACAAATTACATTAAAGTAGGCCAGATAAGACGCGTCAGCGTCGCATCTGGCATTTGCACTGAATGCCGGATGCGGCGTAAAc                                  >  1:766490/1‑150 (MQ=255)
                                                                                                                cTGGCATTTGCACCGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGATTCCCATTATTACAACAAATTACATTAAAGAAGGCCAGATAACACGCTTCAGCGTCGCATCTGCCATTTGCACAGAATTACGGATGCGGCGTaaaa                                  >  1:875485/1‑149 (MQ=37)
                                                                                                                    cATTTGCACCGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGATTCCCATTATTTCAACAAATTACATTAAAGTAGGCCAGATAAGACGCGTCAGCGTCGCATCTGGCATTTGCACTGAATGCCGGATGCGGCGTAAACGCCt                              <  1:562417/150‑1 (MQ=255)
                                                                                                                     aTTTGCACCGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGATTCCCATTATTTCAACAAATTACATTAAAGTAGGCCAGATAAGACGCGTCAGCGTCGCATCTGGCATTTGCACTGAATGCCGGATGCGGCGTAAACACCtt                             <  2:1689823/150‑1 (MQ=255)
                                                                                                                       ttGCACCGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGATTCCCATTATTTCAACAAATTACATTAAAGTAGGCCAGATAAGACGCGTCAGCGTCGCATCTGGCATTTGCACTGAATGCCGGATGCGGCGTAAACGCCTTAt                           <  2:1632675/150‑1 (MQ=255)
                                                                                                                         gCACCGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGATTCCCATTATTTCAACAAAATACATAAAAGTAGGCCAGATAAGACGCGTCAGCGTCGCATCTGGCATTTGCACTGAATGCCGGATGCGGCGTAACCGCCTTATcc                         >  1:1392292/1‑150 (MQ=255)
                                                                                                                          cACCGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGATTCCCATTATTTCAACAAATTACATTAAAGTAGGCCAGATAAGACGCGTCAGCGTCGCATCTGGCATTTGCACTGAATGCCGGATGCGGCGTAAACGCCTTATCCg                        <  1:1379323/150‑1 (MQ=255)
                                                                                                                          cACCGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGATTCCCATTATTTCAACAAATTACATTAAAGTAGGCCAGATAAGACGCGTCAGCGTCGCATCTGGCATTTGCACTGAATGCCGGATGCGGCGTAAACGCCTTATCCg                        >  1:757095/1‑150 (MQ=255)
                                                                                                                          cACCGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGATTCCCATTATTTCAACAAATTACATTAAAGTAGGCCAGATAAGACGCGTCAGCGTCGCATCTGGCATTTGCACTGAATGCCGGATGCGGCGTAAACGCCTTATCCg                        <  1:1137441/150‑1 (MQ=255)
                                                                                                                          cACCGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGATTCCCATTATTTCAACAAATTACATTAAAGTAGGCCAGATAAGACGCGTCAGCGTCGCATCTGGCATTTGCACTGAATGCCGGATGCGGCGTAAACGCCTTATCCg                        >  1:1226717/1‑150 (MQ=255)
                                                                                                                              gATTGCCTGATGCGGCGTAAACGCCTTATCCGGCCTACGATTCCCATTATTTCAACAAATTACATTAAAGTAGGCCAGATAAGACGCGTCAGCGTCGCATCTGGCATTTGCACTGAATGCCGGATGCGGCGTAAACGCCTTATCCGTCCt                    <  2:101303/150‑1 (MQ=255)
                                                                                                                              gATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGAGTCCCATTATTTCAACAAATTACATTAAAGTAGGCCAGATAAGACGCGTCAGCGTCGCATCTGGCATTTGCACTGAATGCCGGATGCGGCGTAAACGCCTTATCCGTCCt                    >  2:1519021/1‑150 (MQ=255)
                                                                                                                                ttGCCAGATGCGGCGTAAACGCCTTATCCGGCCTACGATTCCCATTATTTCAACAAATTACATAAAAGTAGGCCAGATAAGACGCGTCAGCGTCGCATCTGGCATTTGCACTGAATGCCGGATGCGGCGTAAACGCCTTATCCGTCCTAc                  <  1:1351118/150‑1 (MQ=255)
                                                                                                                                aggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACGATTCCCATTATTTCAACAAATTACATTAAAGTAGGCCAGATAAGACGCGTCAGCGTCGCATCTGGCATTTGCACTGAATGCCGGATGCGGCGTAAACGCCTTATCCGTCCTAc                  >  2:775666/3‑150 (MQ=25)
                                                                                                                                 tGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGATTCCCATTATTTCAACAAATTACATTAAAGTAGGCCAGATAAGACGCGTCAGCGTCGCATCTGGCATTTGCACTGAATGCCGGATGCGGCGTAAACGCCTTATCCGTCCTACg                 <  2:513143/150‑1 (MQ=255)
                                                                                                                                 ggCCGGATAAGGCGTTTACGCAGCATCCGGCCTACGATTCCCATTATTTCAACAAATTACATTAAAGTAGGCCAGATAAGACGCGTCAGCGTCGCATCTGGCATTTGCACTGAATGCCGGATGCGGCGTAAACGCCTTATCCGTCCTACg                 >  1:965199/2‑150 (MQ=25)
                                                                                                                                     ggataaggcgTTTACGCCGCATCCGGCCTACGATTCCCATTATTTCAACAAATTACATTAAAGTAGGCCAGATAAGACGCGTCAGCGTCGCATCTGGCATTTGCACTGAATGCCGGATGCGGCGTAAACGCCTTATCCGTCCTACGAATc             >  2:703938/7‑150 (MQ=25)
                                                                                                                                     ggataaggcgTTTACGCCGCATCCGGCCTACGATTCCCAATATTTCAACAAATTACATTAAAGTAGGCCAGATAAGACGCGTCAGCGTCGCATCTGGCATTTGCACTGAATGCCGGATGCGGCGTAACCGCCTTATCCGTCCTACGAATc             >  2:267558/7‑150 (MQ=38)
                                                                                                                                        tGCGGCGTAAACGCCTTATCCGGCCTACGATTCCCATTATTTCAACAAATTACATTAAAGTAGGCCAGATAAGACGCGTCAGCGTCGCATCTGGCATTTGCACTGAATGCCGGATGCGGCGTAAACGCCTTATCCGTCCTACGAATCCcg          <  1:242532/150‑1 (MQ=255)
                                                                                                                                         gcggcgTAAACGCCTTATCCGGCCTACGATTCCCATTATTTCAACAAATTACATTAAAGTAGGCCAGATAAGACGCGTCAGCGTCGCATCTGGCATTTGCACTGAATGCCGGATGCGGCGTAAACGCCTTATCCGTCCTACGAATCCcgc         <  2:908777/150‑1 (MQ=255)
                                                                                                                                         gcggcgTAAACGCCTTATCCGGCCTACGATTCCCATTATTTCAACAAATTACATTAAAGTAGGCCAGATAAGACGCGTCAGCGTCGCATCTGGCATTTGCACTGAATGCCGGATGCGGCGTAAACGCCTTATCCGTCCTACGAATCCcgc         >  1:1339684/1‑150 (MQ=255)
                                                                                                                                          cggcgTAAACGCCTTATCCGGCCTACGATTCCCATTATTTCAACAAATTACATTAAAGTAGGCCAGATAAGACGCGTCAGCGTCGCATCTGGCATTTGCACTGAATGCCGGATGCGGCGTAAACGCCTTATCCGTCCTACGAATCCcgcg        >  1:1571211/1‑150 (MQ=255)
                                                                                                                                          cggcgTAAACGCCTTATCCGGCCTACGATTCCCATTATTTCAACAAATTACATTAAAGTAGGCCAGATAAGACGCGTCAGCGTCGCATCTGGCATTTGCACTGAATGCCGGATGCGGCGTAAACGCCTTATCCGTCCTACGAATCCcgcg        <  2:1699839/150‑1 (MQ=255)
                                                                                                                                          cggcgTAAACGCCTTATCCGGCCTACGATTCCCATTATTTCAACAAATTACATTAAAGTAGGCCAGATAAGACGCGTCAGCGTCGCATCTGGCATTTGCACTGAATGCCGGATGCGGCGTAAACGCCTTATCCGTCCTACGAAACCcgcg        >  1:1566527/1‑150 (MQ=255)
                                                                                                                                            gcgTAAACGCCTTATCCGGCCTACGATTCCCATTATTTCCACAAATTACATTAAAGTAGGCCAGATAAGACGCGTCAGCGTCGCATCTGGCATTTGCACTGAATGCCGGATGCGGCGTAAACGCCTTATCCGTCCTACGAATCCCGCGAt      >  2:957435/1‑150 (MQ=255)
                                                                                                                                            gcgTAAACGCCTTATCCGGCCTACGATTCCCATTATTTCAACAAATTACATTAAAGTAGGCCAGATAAGACGCGTCAGCGTCGCATCTGGCTTTTGCACTGAATGCCGGATGCGGCGTAAACGCCTTATCCGTCCTACGAATCCCGCGAt      >  1:1123282/1‑150 (MQ=255)
                                                                                                                                            gcgTAAACGCCTTATCCGGCCTACGATTCCCATTATTTCAACAAATTACATTAAAGTAGGCCAGATAAGACGCGTCAGCGTCGCATCTGGCATTTGCACTGAATGCCGGATGCGGCGTAAACGCCTTATCCGTCCTACGAATCCCGCGAt      >  1:486150/1‑150 (MQ=255)
                                                                                                                                            gcgTAAACGCCTTATCCGGCCTACGATTCCCATTATTTCAACAAATTACATTAAAGTAGGCCAGATAAGACGCGTCAGCGTCGCATCTGGCATTTGCACTGAATGCCGGATGCGGCGTAAACGCCTTATCCGTCCTACGAATCCCGCGAt      >  2:1298933/1‑150 (MQ=255)
                                                                                                                                             cgTAAACGCCTTATCCGGCCTACGATTCCCATTATTTCAACAAATTACATTAAAGTAGGCCAGATAAGACGCGTCAGCGTCGCATCTGGCATTTGCACTGAATGCCGGATGCGGCGTAAACGCCTTATCCGTCCTACGAATCCCGCGAtt     >  1:1226493/1‑150 (MQ=255)
                                                                                                                                                aaaCGCCTTATCCGGCCTACGATTCCCATTATTTCAACAAATTACATTAAAGTAGGCCAGATAAGACGCGTCAGCGTCGCATCTGGCATTTGCACTGAATGCCGGATGCGGCGTAAACGCCTTATCCGTCCTACGAATCCCGCGATTTcc  >  2:555928/1‑150 (MQ=255)
                                                                                                                                                aaaCGCCTTATCCGGCCTACGATTCCCATTATTTCAACAAATTACATTAAAGTAGGCCAGATAAGACGCGTCAGCGTCGCATCTGGCATTTGCACTGAATGCCGGATGCGGCGTAAACGCCTTATCCGTCCTACGAATCCCGCGATTTcc  >  2:111971/1‑150 (MQ=255)
                                                                                                                                                |                                                                                                                                                     
CTGGGTGTCGAGAACGATATCGCCCTGCTGAATTATCTCTCTTCAGTGACCTTAAGTCCGGCAGATAAAGACTAATCACCTGTAGGCCAGATAAGACGCGTTAGTGTCTTATCTGGCATTTGCACCGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGATTCCCATTATTTCAACAAATTACATTAAAGTAGGCCAGATAAGACGCGTCAGCGTCGCATCTGGCATTTGCACTGAATGCCGGATGCGGCGTAAACGCCTTATCCGTCCTACGAATCCCGCGATTTCC  >  NC_000913/2316753‑2317046

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: