Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,316,904 T→G 6.3% intergenic (+77/+123) 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,9040TG 6.3% 366.8 / 13.2 127intergenic (+77/+123)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 T (73/46);  new base G (4/4);  total (77/50)
Fisher's exact test for biased strand distribution p-value = 7.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.41e-01

Only 100 of 128 total aligned reads displayed.

GGGTGTCGAGAACGATATCGCCCTGCTGAATTATCTCTCTTCAGTGACCTTAAGTCCGGCAGATAAAGACTAATCACCTGTAGGCCAGATAAGACGCGTTAGTGTCTTATCTGGCATTTGCACCGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGATTCCCATTATTTCAACAAATTACATTAAAGTAGGCCAGATAAGACGCGTCAGCGTCGCATCTGGCATTTGCACTGAATGCCGGATGCGGCGTAAACGCCTTATCCGTCCTACGAATCCCGCGATTTCCTGAC  >  NC_000913/2316755‑2317050
                                                                                                                                                     |                                                                                                                                                  
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)
          aaCGATATCGCCCTGCTGAATTATCTCTCTTCAGTGACCTTAAGTCCGGCAGATAAAGACTAATCACCTGTAGGCCAGATAAGACGCGTTAGTGTCTTATCTGGCATTTGCACCGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCt                                                                                                                                          >  2:616424/1‑150 (MQ=255)
          aaCGATATCGCCCTGCTGAATTATCTCTCTTCAGTGACCTTAAGTCCGGCAGATAAAGACTAATCACCTGTAGGCCAGATAAGACGCGTTAGTGTCTTATCTGGCATTTGCACCGATTGCCGGATAAGGCGTTTACGCCGCATCCGGcaa                                                                                                                                          <  2:146991/150‑3 (MQ=255)
           aCGATATCGCCCTGCTGAATTATCTCTCTTCAGTGACCTTAAGTCTGGCAGATAAAGAGTAATCACCTGTAGGCCAGATAAGACGCGTTAGTGTCTTATCTGGCATTTGCACCGATTGCCGGATAAGGCGTTTACGCCGCATCCGGcaat                                                                                                                                         <  1:518765/150‑4 (MQ=255)
            cGATATCGCCCTGCTGAATTATCTCTCTTCAGTGACCTTAAGTCCGGCAGATAAAGACTAATCACCTGTAGGCCAGATAAGACGCGTTAGTGTCTTATCTGGCATTTGCACCGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTAc                                                                                                                                        <  2:785917/150‑1 (MQ=255)
            cGATATCGCCCTGCTGAATTATCTCTCTTCAGTGACCTTAAGTCCGGCAGATAAAGACTAATCACCTGTAGGCCAGATAAGACGCGTTAGTGTCTTATCTGGCATTTGCACCGATTGCCGGATGCGGAGTAAACGCCTTATCCGGCCTAc                                                                                                                                        >  1:95355/1‑150 (MQ=255)
              atatCGCCCTGCTGAATTATCTCTCTTCAGTGACCTTAAGTCCGGCAGATAAAGACTAATCACCTGTAGGCCAGATAAGACGCGTTAGTGTCTTATCTGGCATTTGCACCGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGa                                                                                                                                      >  2:242532/1‑150 (MQ=255)
              atatCGCCCTGCTGAATTATCTCTCTTCAGTGACCTTAAGTCCGGCAGATAAAGACTAATCACCTGTAGGCCAGATAAGACGCGTTAGTGTCTTATCTGGCATTTGCACCGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGa                                                                                                                                      >  2:708949/1‑150 (MQ=255)
               tatCGCCCTGCTGAATTATCTCTCTTCAGTGACCTTAAGTCCGGCAGATAAAGACTAATCACCTGTAGGCCAGATAAGACGCGTTAGTGTCTTATCTGGCATTTGCACCGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGAt                                                                                                                                     <  2:499568/150‑1 (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)
                        gctgAATTATCTCTCTTCAGTGACCTTAAGTCCGGCAGATAAAGACTAATCACCTGTAGGCCAGATAAGACGCGTTAGTGTCTTATCTGGCATTTGCACCGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGATTCCCattat                                                                                                                            <  1:598613/150‑1 (MQ=255)
                          tgAATTATCTCTCTTCAGTGACCTTAAGTCCGGCAGATAAAGACTAATCACCTGTAGGCCAGATAAGACGCGTTAGTGTCTTATCTGGCATTTGCACCGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGATTCCCATTAttt                                                                                                                          <  2:379595/150‑1 (MQ=255)
                                aTCTCTCTTCAGTGACATTAAGTCCGGCAGATAAAGACTAATCACCTGTAGGCCAGATAAGACGCGTTAGTGTCTTATCTGGCATTTGCACCGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGATTCCCATTATATcaacaa                                                                                                                    >  1:870514/1‑150 (MQ=255)
                                     tctTCAGTGACCTTAAGTCCGGCAGATAAAGACTAATCACCTGTAGGCCAGATAAGACGCGTTAGTGTCTTATCTGGCATTTGCACCGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGATTCCCATTATTTCAACAAATTAc                                                                                                               >  1:1043950/1‑150 (MQ=255)
                                          aGTGACCTTTAGTCCGGCAGATAAAGACTAATCACCTGTAGGCCAGATAAGACGCGTTAGTGTCTTATCTGGCATTTGCACCGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGATTCCCATTATTTCAACAAATTACATTaa                                                                                                          <  2:296083/150‑1 (MQ=255)
                                          aGTGACCTTAAGTCCGGCAGATAAAGACTAATCACCTGTAGGCCAGATAAGACGCGTTAGTGTCTTATCTGGCATTTGCACCGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGATTCCCATTATTTCAACAAATTACATTaa                                                                                                          >  2:65554/1‑150 (MQ=255)
                                          aGTGACCTTAAGTCCGGCAGATAAAGACTAATCACCTGTAGGCCAGATAAGACGCGTTAGTGTCTTATCTGGCATTTGCACCGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGATTCCCATTATTTCAACAAATTACATTaa                                                                                                          >  1:433219/1‑150 (MQ=255)
                                           gTGACCTTAAGTCCGGCAGATAAAGACTAATCACCTGTAGGCCAGATAAGACGCGTTAGTGTCTTATCTGGCATTTGCACCGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGATTCCCATTATTTCAACAAATTACATTaaa                                                                                                         >  1:1697647/1‑150 (MQ=255)
                                             gACCTTAAGTCCGGCAGATAAAGACTAATCACCTGTAGGCCAGATAAGACGCGTTAGTGTCTTATCTGGCATTTGCACCGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGATTCCCATTATTTCAACAAATTACATTAAAGt                                                                                                       >  1:1689823/1‑150 (MQ=255)
                                                cTTAAGTCCGGCAGATAAAGACTAATCACCTGTAGGCCAGATAAGACGCGTTAGTGTCTTATCTGGCATTTGCACCGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGATTCCCATTATTTCAACAAATTACATTAAAGTAgg                                                                                                    <  1:1276355/150‑1 (MQ=255)
                                                 ttAAGTCCGGCAGATAAAGACTAATCACCTGTAGGCCAGATAAGACGCGTTAGTGTCTTATCTGGCATTTGCACCGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGATTCCCATTATTTCAACAAATTACATTAAAGTAGGc                                                                                                   >  2:1552305/1‑150 (MQ=255)
                                                   aaGTCCGGCAGATAAAGACTAATCACCTGTAGGCCAGATAAGACGCGTTAGTGTCTTATCTGGCATTTGCACCGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGATTCCCATTATTTCAACAAATTACATTAAAGTAGGCCa                                                                                                 <  2:429672/150‑1 (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: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)
                                                                              tGTAGGCCAGATAAGACGCGTTAGTGTCTTATCTGGCATTTGCACCGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGATTCCCATTATTTCAACAAATTACATTAAAGTAGGCCAGATAAGACGCGTCAGCGTCGCATCTgg                                                                      <  1:1151386/150‑1 (MQ=255)
                                                                                 aGGCCAGATAAGACGCGTTAGTGTCTTATCTGGCATTTGCACCGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGATTCCCATTATTTCAACAAATTACATTAAAGTAGGCCAGATAAGACGCGTCAGCGTCGCATCTGGCAt                                                                   >  2:1282848/1‑150 (MQ=255)
                                                                                  ggCCAGATAAGACGCGTTAGTGTCTTATCTGGCATTTGCACCGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGATTCCCATTATTTCAACAAATTACATTAAAGTAGGCCAGATAAGACGCGTCAGCGTCGCATCTGGCAtt                                                                  >  1:482995/1‑150 (MQ=255)
                                                                                    ccAGATAAGACGCGTTAGTGTCTTATCTGGCATTTGCACCGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGATTCCCATTATTTCAACAAATTACATTAAAGTAGGCCAGATAAGACGCGTCAGCGTCGCATCTGGCATTTg                                                                >  2:1437714/1‑150 (MQ=255)
                                                                                    ccAGATAAGACGCGTTAGTGTCTTATCTGGCATTTGCACCGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGATTCCCATTATTTCAACAAATTACATTAAAGTAGGCCAGATAAGACGCGTCAGCGTCGCATCTGGCATTTg                                                                >  1:1464224/1‑150 (MQ=255)
                                                                                    ccAGATAAGACGCGTTAGTGTCTTATCTGGCATTTGCACCGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGATTCCCATTATTTCAACAAATTACATTAAAGTAGGCCAGATAAGACGCGTCAGCGTCGCATCTGGCATTTg                                                                <  2:1527306/150‑1 (MQ=255)
                                                                                    ccAGATAAGACGCGTTAGTGTCTTATCTGGCATTTGCACCGATTGCCGGATGCGGCGTAAACGCCTTATCAGGCCTACGATTCCCATTATTTCAACAAATTACATTAAAGTAGGCCAGATAAGACGCGTCAGCGTCGCATCTGGCATTTg                                                                >  1:1105731/1‑150 (MQ=255)
                                                                                    ccAGATAAGACACTAACGCGTCTTATCTGGCATTTGCACCGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGATTCCCATTATTTCAACAAATTACATTAAAGTAGGCCAGATAAGACGCGTCAGCGTCGCATCTGGCATTTg                                                                >  1:1599852/1‑150 (MQ=38)
                                                                                    ccAGATAAGACACTAACGCGTCTTATCTGGCATTTGCACCGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGATTCCCATTATTTCAACAAATTACATTAAAGTAGGCCAGATAAGACGCGTCAGCGTCCCATCTGGCATTTg                                                                >  2:136062/1‑150 (MQ=37)
                                                                                     cAGATAAGACGCGTTAGTGTCTTATCTGGCATTTGCACCGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGATTCCCATTATTTCAACAAATTACATTAAAGTAGGCCAGATAAGACGCGTCAGCGACGCATCTGGCATTTga                                                               >  1:691430/1‑149 (MQ=255)
                                                                                     cAGATAAGACGCGTTAGTGTCTTATCTGGCATTTGCACCGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGATTCCCATAAATTCAACAAATTACATTAAAGTAGGCCAGATAAGACGCGTCAGAGTCCCATCTGGCAATTGc                                                               >  1:321612/1‑150 (MQ=37)
                                                                                      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)
                                                                                              cgcgTTAGTGTCTTATCTGGCATTTGCACCGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGATTCCCATTATTTCAACAAATTACATTAAAGTAGGCCAGATAAGACGCGTCAGCGTCGCATCTGGCATTTGCACTGAATGc                                                      >  1:1211271/1‑150 (MQ=255)
                                                                                              cgcgTTAGTGTCTTATCTGGCATTTGCACCGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGATTCCCATTATTTCAACAAATTACATTAAAGTAGGCCAGATAAGACGCGTCAGCGTCGCATCTGGCATTTGCACTGAATGc                                                      >  2:1109338/1‑150 (MQ=255)
                                                                                                    aGTGTCTTATCTGGCATTTGCACCGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGATTCCCATTATTTCAACAAATTACATTAAAGTAGGCCAGATAAGACGCGTCAGCGTCGCATCTGGCATTTGCACTGAATGCCGGATg                                                >  1:881089/1‑150 (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)
                                                                                                              cTGGCATTTGCACCGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGATTCCCATTATTACAACAAATTACATTAAAGAAGGCCAGATAACACGCTTCAGCGTCGCATCTGCCATTTGCACAGAATTACGGATGCGGCGTaaaa                                      >  1:875485/1‑149 (MQ=37)
                                                                                                                   aTTTGCACCGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGATTCCCATTATTTCAACAAATTACATTAAAGTAGGCCAGATAAGACGCGTCAGCGTCGCATCTGGCATTTGCACTGAATGCCGGATGCGGCGTAAACGCCtt                                 <  2:1043950/150‑1 (MQ=255)
                                                                                                                   aTTTGCACCGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGATTCCCATTATTTCAACAAATTACATTAAAGTAGGCCAGATAAGACGCGTCAGCGTCGCATCTGGCATTTGCACTGAATGCCGGATGCGGCGTAAACACCtt                                 <  2:1689823/150‑1 (MQ=255)
                                                                                                                     ttGCACCGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGATTCCCATTATTTCAACAAATTACATTAAAGTAGGCCAGATAAGACGCGTCAGCGTCGCATCTGGCATTTGCACTGAATGCCGGATGCGGCGTAAACGCCTTAt                               <  2:1632675/150‑1 (MQ=255)
                                                                                                                        cACCGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGATTCCCATTATTTCAACAAATTACATTAAAGTAGGCCAGATAAGACGCGTCAGCGTCGCATCTGGCATTTGCACTGAATGCCGGATGCGGCGTAAACGCCTTATCCg                            <  1:1379323/150‑1 (MQ=255)
                                                                                                                        cACCGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGATTCCCATTATTTCAACAAATTACATTAAAGTAGGCCAGATAAGACGCGTCAGCGTCGCATCTGGCATTTGCACTGAATGCCGGATGCGGCGTAAACGCCTTATCCg                            <  1:1137441/150‑1 (MQ=255)
                                                                                                                        cACCGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGATTCCCATTATTTCAACAAATTACATTAAAGTAGGCCAGATAAGACGCGTCAGCGTCGCATCTGGCATTTGCACTGAATGCCGGATGCGGCGTAAACGCCTTATCCg                            >  1:757095/1‑150 (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)
                                                                                                                               ggCCGGATAAGGCGTTTACGCAGCATCCGGCCTACGATTCCCATTATTTCAACAAATTACATTAAAGTAGGCCAGATAAGACGCGTCAGCGTCGCATCTGGCATTTGCACTGAATGCCGGATGCGGCGTAAACGCCTTATCCGTCCTACg                     >  1:965199/2‑150 (MQ=25)
                                                                                                                                gCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGATTCCCATTATTTCAACAAATTACATTAAAGTAGGCCAGATAAGACGCGTCAGCGTCGCATCTGGCATTTGCACTGAATGCCGGATGCGGCGTAAACGCCTTATCCGTCCTACGa                    <  2:1211271/150‑1 (MQ=255)
                                                                                                                                   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)
                                                                                                                                          gcgTAAACGCCTTATCCGGCCTACGATTCCCATTATTTCCACAAATTACATTAAAGTAGGCCAGATAAGACGCGTCAGCGTCGCATCTGGCATTTGCACTGAATGCCGGATGCGGCGTAAACGCCTTATCCGTCCTACGAATCCCGCGAt          >  2:957435/1‑150 (MQ=255)
                                                                                                                                          gcgTAAACGCCTTATCCGGCCTACGATTCCCATTATTTCAACAAATTACATTAAAGTAGGCCAGATAAGACGCGTCAGCGTCGCATCTGGCTTTTGCACTGAATGCCGGATGCGGCGTAAACGCCTTATCCGTCCTACGAATCCCGCGAt          >  1:1123282/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)
                                                                                                                                                aCGCCTTATCCGGCCTACGATTCCCATTATTTCAACAAATTACATTAAAGTAGGCCAGATAAGACGCGTCAGCGTCGCATCTGGCATTTGCACTGAATGCCGGATGCGGCGTAAACGCCTTATCCGTCCTACGAATCCCGCGATTTCCTg    >  1:959834/1‑150 (MQ=255)
                                                                                                                                                  gCCTTATCCGGCCTACGATTCCCATTATTTCAACAAATTACATTAAAGTAGGCCAGATAAGACGCGTCAGCGTCGCATCTGGCATTTGCACTGAATGCCGGATGCGGCGTAAACGCCTTATCCGTCCTACGAATCCCGCGATTTCCTGAc  <  2:299572/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
GGGTGTCGAGAACGATATCGCCCTGCTGAATTATCTCTCTTCAGTGACCTTAAGTCCGGCAGATAAAGACTAATCACCTGTAGGCCAGATAAGACGCGTTAGTGTCTTATCTGGCATTTGCACCGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGATTCCCATTATTTCAACAAATTACATTAAAGTAGGCCAGATAAGACGCGTCAGCGTCGCATCTGGCATTTGCACTGAATGCCGGATGCGGCGTAAACGCCTTATCCGTCCTACGAATCCCGCGATTTCCTGAC  >  NC_000913/2316755‑2317050

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: