Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,316,890 G→A 5.2% intergenic (+63/+137) 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,8900GA 5.2% 356.9 / 5.5 114intergenic (+63/+137)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 G (64/44);  new base A (2/4);  total (66/48)
Fisher's exact test for biased strand distribution p-value = 2.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.83e-01

Only 100 of 119 total aligned reads displayed.

GCGATGATGAAGCTGGGTGTCGAGAACGATATCGCCCTGCTGAATTATCTCTCTTCAGTGACCTTAAGTCCGGCAGATAAAGACTAATCACCTGTAGGCCAGATAAGACGCGTTAGTGTCTTATCTGGCATTTGCACCGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGATTCCCATTATTTCAACAAATTACATTAAAGTAGGCCAGATAAGACGCGTCAGCGTCGCATCTGGCATTTGCACTGAATGCCGGATGCGGCGTAAACGCCTTATCCGTCCTACGAATCCCGC  >  NC_000913/2316741‑2317039
                                                                                                                                                     |                                                                                                                                                     
gCGATGATGAAGCTGGGTGTCGAGAACGATATCGCCCTGCTGAATTATCTCTCTTCAGTGACCTTAAGTCCGGCAGATAAAGACTAATCACCTGTAGGCCAGATAAGACGCGTTAGTGTCTTATCTGGCATTTGCACCGATTGCCGGATg                                                                                                                                                       >  2:810254/1‑150 (MQ=255)
   atgatgAAGCTGGGTGTCGAGAACGATATCGCCCTGCTGAATTATCTCTCTTCAGTGACCTTAAGTCCGGCAGATAAAGACTAATCACCTGTAGGCCAGATAAGACGCGTTAGTGTCTTATCTGGCATTTGCACCGATTGCCGGATgcgg                                                                                                                                                    <  2:770829/150‑1 (MQ=255)
   atgatgAAGCTGGGTGTCGAGAACGATATCGCCCTGCTGAATTATCTCTCTTCAGTGACCTTAAGTCCGGCAGATAAAGACTAATCACCTGTAGGCCAGATAAGACGCGTTAGTGTCTTATCTGGCATTTGCACCGATTGCCGGATgcgg                                                                                                                                                    >  2:1156987/1‑150 (MQ=255)
            cTGGGTGTCGAGAACGATATCGCCCTGCTGAATTATCTCTCTTCAGTGACCTTAAGTCCGGCAGATAAAGACTAATCACCTGTAGGCCAGATAAGACGCGTTAGTGTCTTATCTGGCATTTGCACCGATTGCCGGATGCGGAGTAAACGc                                                                                                                                           >  2:1665297/1‑150 (MQ=255)
              gggTGTCTAGAACGATATCGCCCTGCTGAATTATCTCTCTTCAGTGACCTTAAGTCCGGCAGATAAAGACTAATCACCTGTAGGCCAGATAAGACGCGTTAGTGTCTTATCTGGCATTTGCACCGATTGCCGGATGCGGCGTTATCGCCt                                                                                                                                         <  2:1226334/150‑1 (MQ=255)
                  gtCGAGAACGATATCGCCCTGCTGAATTATCTCTCTTCAGTGACCTTAAGTCCGGCAGATAAAGACTAATCACCTGTAGGCCAGATAAGACGCGTTAGTGTCTTATCTGGCATTTGCACCGATTGCCGGATGCGGCGTAAACGCCTTATc                                                                                                                                     >  2:562417/1‑150 (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)
                      agaACGATATCGCCCTGCTGAATTATCTCTCTTCAGTGACCTTAAGTCCGGCAGATAAAGACTAATCACCTGTAGGCCAGATAAGACGCGTTAGTGTCTTATCTGGCATTTGCACCGATTGCCGGATAAGGCGTTTACGCCGCATCCGGc                                                                                                                                 <  1:323234/150‑1 (MQ=38)
                        aaCGATATCGCCCTGCTGAATTATCTCTCTTCAGTGACCTTAAGTCCGGCAGATAAAGACTAATCACCTGTAGGCCAGATAAGACGCGTTAGTGTCTTATCTGGCATTTGCACCGATTGCCGGATGCGGCGTAAACGCCTTTTCCGGCCt                                                                                                                               >  2:694322/1‑150 (MQ=255)
                        aaCGATATCGCCCTGCTGAATTATCTCTCTTCAGTGACCTTAAGTCCGGCAGATAAAGACTAATCACCTGTAGGCCAGATAAGACGCGTTAGTGTCTTATCTGGCATTTGCACCGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCt                                                                                                                               >  2:616424/1‑150 (MQ=255)
                         aCGATATCGCCCTGCTGAATTATCTCTCTTCAGTGACCTTAAGTCCGGCAGATAAAGACTAATCACCTGTAGGCCAGATAAGACGCGTTAGTGTCTTATCTGGCATTTGCACCGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTa                                                                                                                              >  1:299572/1‑150 (MQ=255)
                          cGATATCGCCCTGCTGAATTATCTCTCTTCAGTGACCTTAAGTCCGGCAGATAAAGACTAATCACCTGTAGGCCAGATAAGACGCGTTAGTGTCTTATCTGGCATTTGCACCGATTGCCGGATGCGGAGTAAACGCCTTATCCGGCCTAc                                                                                                                             >  1:95355/1‑150 (MQ=255)
                             tatCGCCCTGCTGAATTATCTCTCTTCAGTGACCTTAAGTCCGGCAGATAAAGACTAATCACCTGTAGGCCAGATAAGACGCGTTAGTGTCTTATCTGGCATTTGCACCGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGAt                                                                                                                          <  2:499568/150‑1 (MQ=255)
                                cGCCCTGCTGAATTATCTCTCTTCAGTGACCTTAAGTCCGGCAGATAAAGACTAATCACCTGTAGGCCAGATAAGACGCGTTAGTGTCTTATCTGGCATTTGCACCGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGATTcc                                                                                                                       >  2:1057063/1‑150 (MQ=255)
                                 gCCCTGCTGAATTATCTCTCTTCAGTGACCTTAAGTCCGGCAGATAAAGACTAATCACCTGTAGGCCAGATAAGACGCGTTAGTGTCTTATCTGGCATTTGCACCGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGATTccc                                                                                                                      >  1:750077/1‑150 (MQ=255)
                                  cccTGCTGAATTATCTCTCTTCAGTGACCTTAAGTCCGGCAGATAAAGACTAATCACCTGTAGGCCAGATAAGACGCGTTAGTGTCTTATCTGGCATTTGCACCGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGATTCCCa                                                                                                                     <  2:531009/150‑1 (MQ=255)
                                    ctgctgAATTATCTCTCTTCAGTGACCTTAAGTCCGGCAGATAAAGACTAATCACCTGTAGGCCAGATAAGACGCGTTAGTGTCTTATCTGGCATTTGCACCGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGATTCCCatt                                                                                                                   >  1:1688279/1‑150 (MQ=255)
                                    ctgctgAATTATCTCTCTTCAGTGACCTTAAGTCCGGCAGATAAAGACTAATCACCTGTAGGCCAGATAAGACGCGTTAGTGTCTTATCTGGCATTTGCACCGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGATTCCCatt                                                                                                                   <  2:726058/150‑1 (MQ=255)
                                     tgctgAATTATCTCTCTTCAGTGACCTTAAGTCCGGCAGATAAAGACTAATCACCTGTAGGCCAGATAAGACGCGTTAGTGTCTTATCTGGCATTTGCACCGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGATTCCCatta                                                                                                                  <  2:415202/150‑1 (MQ=255)
                                      gctgAATTATCTCTCTTCAGTGACCTTAAGTCCGGCAGATAAAGACTAATCACCTGTAGGCCAGATAAGACGCGTTAGTGTCTTATCTGGCATTTGCACCGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGATTCCCattat                                                                                                                 <  1:598613/150‑1 (MQ=255)
                                        tgAATTATCTCTCTTCAGTGACCTTAAGTCCGGCAGATAAAGACTAATCACCTGTAGGCCAGATAAGACGCGTTAGTGTCTTATCTGGCATTTGCACCGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGATTCCCATTAttt                                                                                                               <  2:379595/150‑1 (MQ=255)
                                         gAATTATCTCTCTTCAGTGACCTTAAGTCCGGCAGATAAAGACTAATCACCTGTAGGCCAGATAAGACGCGTTAGTGTCTTATCTGGCATTTGCACCGATTGCCGGATGCGGCGt                                                                                                                                                 <  2:1219271/115‑1 (MQ=255)
                                         gAATTATCTCTCTTCAGTGACCTTAAGTCCGGCAGATAAAGACTAATCACCTGTAGGCCAGATAAGACGCGTTAGTGTCTTATCTGGCATTTGCACCGATTGCCGGATGCGGCGt                                                                                                                                                 >  1:1219271/1‑115 (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)
                                                        aGTGACCTTTAGTCCGGCAGATAAAGACTAATCACCTGTAGGCCAGATAAGACGCGTTAGTGTCTTATCTGGCATTTGCACCGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGATTCCCATTATTTCAACAAATTACATTaa                                                                                               <  2:296083/150‑1 (MQ=255)
                                                        aGTGACCTTAAGTCCGGCAGATAAAGACTAATCACCTGTAGGCCAGATAAGACGCGTTAGTGTCTTATCTGGCATTTGCACCGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGATTCCCATTATTTCAACAAATTACATTaa                                                                                               >  1:433219/1‑150 (MQ=255)
                                                        aGTGACCTTAAGTCCGGCAGATAAAGACTAATCACCTGTAGGCCAGATAAGACGCGTTAGTGTCTTATCTGGCATTTGCACCGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGATTCCCATTATTTCAACAAATTACATTaa                                                                                               >  2:65554/1‑150 (MQ=255)
                                                         gTGACCTTAAGTCCGGCAGATAAAGACTAATCACCTGTAGGCCAGATAAGACGCGTTAGTGTCTTATCTGGCATTTGCACCGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGATTCCCATTATTTCAACAAATTACATTaaa                                                                                              >  1:1697647/1‑150 (MQ=255)
                                                          tGACCTTAAGTCCGGCAGATAAAGACTAATCACCTGTAGGCCAGATAAGACGCGTTAGTGTCTTATCTGGCATTTGCACCGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGATTCCCATTATTTCAACAAATTACATTAAAg                                                                                             <  1:810254/150‑1 (MQ=255)
                                                              cTTAAGTCCGGCAGATAAAGACTAATCACCTGTAGGCCAGATAAGACGCGTTAGTGTCTTATCTGGCATTTGCACCGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGATTCCCATTATTTCAACAAATTACATTAAAGTAgg                                                                                         <  1:1276355/150‑1 (MQ=255)
                                                               ttAAGTCCGGCAGATAAAGACTAATCACCTGTAGGCCAGATAAGACGCGTTAGTGTCTTATCTGGCATTTGCACCGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGATTCCCATTATTTCAACAAATTACATTAAAGTAGGc                                                                                        >  2:1552305/1‑150 (MQ=255)
                                                                   gTCCGGCAGATAAAGACTAATCACCTGTAGGCCAGATAAGACGCGTTAGTGTCTTATCTGGCATTTGCCCCGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGATTCCCATTATTTCAACAAATTACATTAATGTAGGCCAGa                                                                                    >  2:269249/1‑150 (MQ=255)
                                                                   gTCCGGCAGATAAAGACTAATCACCTGTAGGCCAGATAAGACGCGTTAGTGTCTTATCTGGCATTTGCACCGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGATTCCCATTATTTCAACAAATTACATTAAAGTAGGCCAGa                                                                                    >  2:62814/1‑150 (MQ=255)
                                                                      cGGCAGATAAAGACTAATCACCTGTAGGCCAGATAAGACGCGTTAGTGTCTTATCTGGCATTTGCACCGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGATTCCCATTATTTCAACAAATTACATTAAAGTAGGCCAGATaa                                                                                 <  2:750077/150‑1 (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:1431606/150‑1 (MQ=255)
                                                                              aaaGACTAATCACCTGTAGGCCAGATAAGACGCGTTAGTGTCTTATCTGGCATTTGCACCGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGATTCCCATTATTTCAACAAATTACATTAAAGTAGGCCAGATAAGACGCGTc                                                                         >  2:750370/1‑150 (MQ=255)
                                                                                  aCTAATCACCTGTAGGCCAGATAAGACGCGTTAGTGTCTTATCTGGCATTTGCACCGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGATTCCCATTATTTCAACAAATTACATTAAAGTAGGCCAGATAAGACGCGTCAGCg                                                                     >  1:1594907/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)
                                                                                             gTAGGCCAGATAAGACGCGTTAGTGTCTTATCTGGCATTTGCACCGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGATTCCCATTATTTCAACAAATTACATTAAAGTAGGCCAGATGCGACGCTGACGCGTCTTATCTGGc                                                          <  1:172747/150‑1 (MQ=25)
                                                                                             gTAGGCCAGATAAGACGCGTTAGTGTCTTATCTGGCATTTGCACCGATTGCCGGATGCGGCGTAAACGCCTTATCCGCCCTACGATTCCCATTATTTCAACAAATTACATTGAAGTAGGCCAGATAAGACGCGTCAGCGTCGCATCTGGc                                                          <  1:1120185/150‑1 (MQ=255)
                                                                                               aGGCCAGATAAGACGCGTTAGTGTCTTATCTGGCATTTGCACCGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGATTCCCATTATTTCAACAAATTACATTAAAGTAGGCCAGATAAGACGCGTCAGCGTCGCATCTGGCAt                                                        >  2:1282848/1‑150 (MQ=255)
                                                                                                ggCCAGATAAGACGCGTTAGTGTCTTATCTGGCATTTGCACCGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGATTCCCATTATTTCAACAAATTACATTAAAGTAGGCCAGATAAGACGCGTCAGCGTCGCATCTGGCAtt                                                       >  1:482995/1‑150 (MQ=255)
                                                                                                  ccAGATAAGACGCGTTAGTGTCTTATCTGGCATTTGCACCGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGATTCCCATTATTTCAACAAATTACATTAAAGTAGGCCAGATAAGACGCGTCAGCGTCGCATCTGGCATTTg                                                     <  2:1527306/150‑1 (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                                                    <  1:932158/150‑1 (MQ=255)
                                                                                                   cAGATAAGACGCGTTAGTGTCTTATCTGGCATTTGCACCGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGATTCCCATAAATTCAACAAATTACATTAAAGTAGGCCAGATAAGACGCGTCAGAGTCCCATCTGGCAATTGc                                                    >  1:321612/1‑150 (MQ=37)
                                                                                                    aGATAAGACGCGTTAGTGTCTTATCTGGCATTTGCACCGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGATTCCCATTATTTCAACAAATTACATTAAAGTAGGCCAGATAAGACGCGTCAGTGTCGCATCTGGCATTTGCa                                                   >  2:117073/1‑150 (MQ=255)
                                                                                                    aGATAAGACGCGTTAGTGTCTTATCTGGCATTTGCACCGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGATTCCCATTATTTCAACAAATTACATTAAAGTAGGCCAGATAAGACGCGTCAGCGTCGCATCTTGCATTTGCa                                                   >  1:622752/1‑150 (MQ=255)
                                                                                                    aGATAAGACGCGTTAGTGTCTTATCTGGCATTTGCACCGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGATTCACATTATTTCAACAAATTACATTAAAGTAGGCAAGATAAGACGCGTCAGCGACGCATCTTGCATTTgaa                                                   >  1:461884/1‑148 (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:881089/1‑150 (MQ=255)
                                                                                                                  aGTGTCTTATCTGGCATTTGCACCGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGATTCCCATTATTTCAACAAATTACATTAAAGTAGGCCAGATAAGACGCGTCAGCGTCGCATCTGGCATTTGCACTGAATGCCGGATg                                     <  1:395116/150‑1 (MQ=255)
                                                                                                                   gtgtCTTATCTGGCATTTGCACCGATTGCCGGATGCGGCGTAAACGCCTTATACGGCCTACGATTCCCATTATATCAACAAATTAAATTAAAGTAGTCCAGCTAACACGAGTCAGCGTCGCATCTGTCATTTGCACTAAATGCCGCAtac                                    >  1:148756/1‑148 (MQ=37)
                                                                                                                     gtCTTATCTGGCATTTGCACCGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGATTCCCATTATTTCAACAAATTACATTAAAGTAGGCCAGATAAGACGCGTCAGCGTCGCATCTGGCATTTGCACTGAATGCCGGATgcgg                                  <  1:511700/150‑1 (MQ=255)
                                                                                                                      tCTTATCTGGCATTTGCACCGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGATTCCCATTATTTCAACAAATTACATTAAAGTAGGACAGATAAGACGCGTCAGCGTCGCATCTGGCATTTGCACTGAATGCCGGATgcggc                                 >  1:543943/1‑150 (MQ=255)
                                                                                                                        ttATCTGGCATTTGCACCGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGATTCCCATTATTTCAACAAATTACATTAAAGTAGGCCAGATAAGACGCGTCAGCGTCGCATCTGGCATTTGCACTGAATGCCGGATGCGGCGt                               <  1:1200120/150‑1 (MQ=255)
                                                                                                                          aTCTGGCATTTGCACCGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGATTCCCATTATTTCAACAAATTACATTAAAGTAGGCCAGATAAGACGCGTCAGCGTCGCATCTGGCATTTGCACTGAATGCCGGATGCGGCGTaa                             <  2:1688279/150‑1 (MQ=255)
                                                                                                                            cTGGCATTTGCACCGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGATTCCCATTATTTCAACAAATTACATTAAAGTAGGCCAGATAAGACGCGTCAGCGTCGCATCTGGCATTTGCACTGAATGCCGGATGCGGCGTAAAc                           >  1:766490/1‑150 (MQ=255)
                                                                                                                            cTGGCATTTGCACCGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGATTCCCATTATTTCAACAAAATACATTAAAGTAGGCCAGATAAGACGCGTCAGCGTAGCATCTGGCATTTGAACTGAATGCCGGATGCGGCGTAAAc                           >  1:223578/1‑150 (MQ=255)
                                                                                                                            cTGGCATTTGCACCGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGAGTCCCAGTATTTCAACAAATTACATTAAAGTAGGCCAGATAAGAGGCGTCAGCGTCGCATCTGGCATTTGCACTGAATGCCGGATGCGGCGTAAAc                           >  2:86403/1‑150 (MQ=255)
                                                                                                                                cATTTGCACCGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGATTCCCATTATTTCAACAAATTACATTAAAGTAGGCCAGATAAGACGCGTCAGCGTCGCATCTGGCATTTGCACTGAATGCCGGATGCGGCGTAAACGCCt                       <  1:562417/150‑1 (MQ=255)
                                                                                                                                 aTTTGCACCGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGATTCCCATTATTTCAACAAATTACATTAAAGTAGGCCAGATAAGACGCGTCAGCGTCGCATCTGGCATTTGCACTGAATGCCGGATGCGGCGTAAACGCCtt                      <  2:1043950/150‑1 (MQ=255)
                                                                                                                                 aTTTGCACCGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGATTCCCATTATTTCAACAAATTACATTAAAGTAGGCCAGATAAGACGCGTCAGCGTCGCATCTGGCATTTGCACTGAATGCCGGATGCGGCGTAAACACCtt                      <  2:1689823/150‑1 (MQ=255)
                                                                                                                                    tGCACCGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGATTCCCATTATTTCAACAAATTACATTAAAGTAGGCCAGATAAGACGCGTCAGCGTCGCATCTGGCATTTGCACTGAATGCCGGATGCGGCGTAAACGCCTTATc                   <  1:1162569/150‑1 (MQ=255)
                                                                                                                                     gCACCGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGATTCCCATTATTTCAACAAAATACATAAAAGTAGGCCAGATAAGACGCGTCAGCGTCGCATCTGGCATTTGCACTGAATGCCGGATGCGGCGTAACCGCCTTATcc                  >  1:1392292/1‑150 (MQ=255)
                                                                                                                                      cACCGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGATTCCCATTATTTCAACAAATTACATTAAAGTAGGCCAGATAAGACGCGTCAGCGTCGCATCTGGCATTTGCACTGAATGCCGGATGCGGCGTAAACGCCTTATCCg                 >  1:757095/1‑150 (MQ=255)
                                                                                                                                      cACCGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGATTCCCATTATTTCAACAAATTACATTAAAGTAGGCCAGATAAGACGCGTCAGCGTCGCATCTGGCATTTGCACTGAATGCCGGATGCGGCGTAAACGCCTTATCCg                 <  1:1137441/150‑1 (MQ=255)
                                                                                                                                      cACCGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGATTCCCATTATTTCAACAAATTACATTAAAGTAGGCCAGATAAGACGCGTCAGCGTCGCATCTGGCATTTGCACTGAATGCCGGATGCGGCGTAAACGCCTTATCCg                 <  1:1379323/150‑1 (MQ=255)
                                                                                                                                         cGATTGCCGGATGAGGCTTAAACGCCTTATCCGGCCTACGATTCCCATTATTTCAACAAATTACATTAAAGTAGGCCAGATAAGACGCGTCAGCGTCGCATCTGGCATTTGCACTGAATGCCGGATGCGGCGTAAACGCCTTATCCGTcc              <  1:627885/150‑1 (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)
                                                                                                                                                    tGCGGCGTAAACGCCTTATCCGGCCTACGATTCCCATTATTTCAACAAATTACATTAAAGTAGGCCAGATAAGACGCGTCAGCGTCGCATCTGGCATTTGCACTGAATGCCGGATGCGGCGTAAACGCCTTATCCGTCCTACGAATCCcg   <  1:242532/150‑1 (MQ=255)
                                                                                                                                                     gcggcgTAAACGCCTTATCCGGCCTACGATTCCCATTATTTCAACAAATTACATTAAAGTAGGCCAGATAAGACGCGTCAGCGTCGCATCTGGCATTTGCACTGAATGCCGGATGCGGCGTAAACGCCTTATCCGTCCTACGAATCCcgc  >  1:1339684/1‑150 (MQ=255)
                                                                                                                                                     gcggcgTAAACGCCTTATCCGGCCTACGATTCCCATTATTTCAACAAATTACATTAAAGTAGGCCAGATAAGACGCGTCAGCGTCGCATCTGGCATTTGCACTGAATGCCGGATGCGGCGTAAACGCCTTATCCGTCCTACGAATCCcgc  <  2:908777/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GCGATGATGAAGCTGGGTGTCGAGAACGATATCGCCCTGCTGAATTATCTCTCTTCAGTGACCTTAAGTCCGGCAGATAAAGACTAATCACCTGTAGGCCAGATAAGACGCGTTAGTGTCTTATCTGGCATTTGCACCGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGATTCCCATTATTTCAACAAATTACATTAAAGTAGGCCAGATAAGACGCGTCAGCGTCGCATCTGGCATTTGCACTGAATGCCGGATGCGGCGTAAACGCCTTATCCGTCCTACGAATCCCGC  >  NC_000913/2316741‑2317039

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