Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,767,667 T→G 17.2% intergenic (+312/‑43) rnlB → / → yfjP CP4‑57 prophage; uncharacterized protein/CP4‑57 prophage; 50S ribosome‑binding GTPase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,767,6670TG17.2% 450.4 / 97.8 201intergenic (+312/‑43)rnlB/yfjPCP4‑57 prophage; uncharacterized protein/CP4‑57 prophage; 50S ribosome‑binding GTPase family protein
Reads supporting (aligned to +/- strand):  ref base T (75/91);  new base G (22/13);  total (97/104)
Fisher's exact test for biased strand distribution p-value = 6.45e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.96e-02

Only 100 of 202 total aligned reads displayed.

ATGGTGCTGTGAGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCGGTCTACCGC  >  NC_000913/2767568‑2767764
                                                                                                   |                                                                                                 
   gTGCTGTGAGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTcgn                                                                                                <  2:3147643/100‑2 (MQ=255)
     gCTGTGAGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTcgc                                                                                                <  1:1611899/98‑1 (MQ=255)
     gCTGTGAGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCgg                                                                                              <  1:3876459/100‑1 (MQ=255)
      cTGTGAGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGt                                                                                             <  2:4361947/100‑1 (MQ=255)
      cTGTGAGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGt                                                                                             <  2:1072779/100‑1 (MQ=255)
       tgtgAGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGtt                                                                                            >  2:3258410/1‑100 (MQ=255)
       tgtgAGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGtt                                                                                            >  2:3550738/1‑100 (MQ=255)
        gtgAGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTc                                                                                           <  2:193634/100‑1 (MQ=255)
         ttaGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCt                                                                                          <  1:3637560/98‑1 (MQ=255)
         tgAGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCt                                                                                          <  1:695222/100‑1 (MQ=255)
          gAGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTg                                                                                         <  1:1329740/100‑1 (MQ=255)
           aGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGc                                                                                        >  2:1821891/1‑100 (MQ=255)
            gTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCGCGAGGTTCTGcc                                                                                       <  1:3492621/100‑1 (MQ=255)
            gTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCGCGAGGTTCTGcc                                                                                       <  2:1295162/100‑1 (MQ=255)
             ttCAGAAAGATTCTGGAATCTATGTGCATTTAAAAAATCTTCAGATAGATATCATTCTTCTGAACAATCCGATGAACGTCAGAACCGCGAGGTTCTGcc                                                                                       <  1:1190168/99‑1 (MQ=255)
             ttCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCg                                                                                      <  2:1013475/100‑1 (MQ=255)
              tCAGAAAGATTCTGCAATCTATGTGCATTGGAATAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGcc                                                                                       >  1:2906259/1‑98 (MQ=255)
              tCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGt                                                                                     >  1:2371951/1‑100 (MQ=255)
              tCAGAAAGATTCTGCAATCTATGTGCATTAAAAAAATCTTCAGATAGATATCATCCTCCTGAACAATCAGATGAACGACAGAACCTCGCGGTTCTGCCg                                                                                      >  1:1288362/1‑99 (MQ=255)
               cAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGtt                                                                                    >  2:58681/1‑100 (MQ=255)
               cAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCGCGAGGTTCTGCCGtt                                                                                    <  2:1595142/100‑1 (MQ=255)
                aGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCGCGAGGTTCTGCCGttc                                                                                   <  1:1504332/100‑2 (MQ=255)
                 gAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGtttt                                                                                  <  1:3167569/100‑1 (MQ=255)
                  aatgtttCTGTAATCTATGCGCCTGTAAAACATCCTCAGATAGTTATCGTTCTCCTTAACAATCCGATGAACGTCATACCCGCGAGCTTCTGCCGtt                                                                                    <  1:4253235/92‑1 (MQ=255)
                    aGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGt                                                                               <  2:4260995/100‑1 (MQ=255)
                     gATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCGCGAGGTTCTGCCGttcatcgn                                                                              <  2:2690202/100‑7 (MQ=38)
                       ttCTGCAATCTATGTGCATGGAAAAAATCGTCAGATAGATATCATTCGTCTGAACAATTCGATAAACGGCAGAACCTCGCGGTTTCGGCGTTTTTGtttt                                                                            >  1:448421/1‑100 (MQ=255)
                         cTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTAc                                                                          <  1:1649816/100‑1 (MQ=255)
                         cTGCAATCTATGTGCATTGAAAAAAGCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGtt                                                                              <  1:3260208/96‑1 (MQ=255)
                          tGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACt                                                                         >  2:2423763/1‑100 (MQ=255)
                            cAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACttt                                                                       <  2:2733723/100‑1 (MQ=255)
                            cAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACtt                                                                        >  1:4085581/1‑99 (MQ=255)
                            cAATCTATGTGAATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACtt                                                                        <  1:3207896/99‑1 (MQ=255)
                               tCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCa                                                                    <  2:840098/100‑1 (MQ=255)
                                  aggtgCATTGAAAAAATCTTCAGATAGATATCAGTCTCCTGAACAATCCGATGAACCGCAGAACCTCGCGGTTCTGCCGTGTTTGATTTACTGTTCag                                                                   <  1:2543490/96‑1 (MQ=255)
                                   tgtgCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAgagg                                                                <  2:147544/100‑1 (MQ=255)
                                   tgtgCATTGAAAAAATCTTCAGATAGATATCATTCTCCCGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAgagg                                                                <  2:3439619/100‑1 (MQ=255)
                                    gtgCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAgagg                                                                >  2:24809/1‑99 (MQ=255)
                                    gtgCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAgagg                                                                <  2:2451729/99‑1 (MQ=255)
                                      gCATTGAAAAAATCTTCAGATAGATATCCTTCTCCTTAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAgt                                                             <  1:3941851/100‑1 (MQ=255)
                                      gCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAgt                                                             >  2:2145584/1‑100 (MQ=255)
                                       cATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAgtg                                                            >  2:815408/1‑100 (MQ=255)
                                        aTTGAAAAAATCTTCAGATAGATATCATTCTCCTTAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAgtgg                                                           <  1:596998/100‑1 (MQ=255)
                                        aTTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAgtgg                                                           <  2:620103/100‑1 (MQ=255)
                                        aTTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAgtgg                                                           <  2:133181/100‑1 (MQ=255)
                                           gAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGtga                                                        >  1:2479718/1‑100 (MQ=255)
                                              aaaaTCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGtgatga                                                     <  1:1455194/100‑1 (MQ=255)
                                                  tCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGaaaaa                                                 >  1:3129196/1‑100 (MQ=255)
                                                      cAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAAtttt                                             >  2:3054428/1‑100 (MQ=255)
                                                      cAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAAtttt                                             >  1:3404555/1‑100 (MQ=255)
                                                      cAGATAGATATCATTATCCTTAACAATCCGATGAACGGCAGAATCTCGCGGTTCTGCAGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAAtttt                                             <  1:2083879/100‑1 (MQ=255)
                                                       agatagatATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTg                                            >  2:3698908/1‑100 (MQ=255)
                                                       agatagatATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTg                                            >  1:1923674/1‑100 (MQ=255)
                                                         atagataTCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGaa                                          <  2:1359279/100‑1 (MQ=255)
                                                          tagataTCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAg                                         >  2:3791736/1‑100 (MQ=255)
                                                          tagataTCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAg                                         <  2:901870/100‑1 (MQ=255)
                                                             atatCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGCGTGGTGATGAAAAATTTTAAAGTTc                                      >  1:834385/1‑100 (MQ=255)
                                                                tCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGc                                   <  2:688815/100‑1 (MQ=255)
                                                                tCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGc                                   <  2:2329420/100‑1 (MQ=255)
                                                                 cATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCa                                  <  2:4085611/100‑1 (MQ=255)
                                                                 cATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCa                                  <  2:271044/100‑1 (MQ=255)
                                                                 cATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTATCAGACGAGTGGTGCTGAAAAATTTTGAAGTTCTGCa                                  >  1:2968279/1‑100 (MQ=255)
                                                                 cATTCTCCTGAACAATCCAATGAGCGGCAGAACCTCGAGGTTCTGACGTTTTTGTTTTACTTTTGAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTg                                    <  1:3276671/98‑1 (MQ=255)
                                                                  aTTCTCCTGAACAATCCGATGAACGGCAGAACGTCGCGGTTCTGCCGTTTTTGTTTGACTTTTAAGAGGAGTGGCAATGAAAAAATTTGAAGTTATGCAg                                 >  1:2138245/1‑100 (MQ=255)
                                                                  aTTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAg                                 <  1:2300816/100‑1 (MQ=255)
                                                                   cgatcttGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGc                                <  1:4289973/94‑1 (MQ=255)
                                                                    tctcCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGcc                               <  2:2846321/100‑1 (MQ=255)
                                                                     ctcCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCa                              >  2:3314833/1‑100 (MQ=255)
                                                                        cTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGCTTATGTTTTACTTTTCAGAGGAGAGGTGATGAAAAATTTTCAAGTTCTGCAGCCACTa                           >  1:156371/1‑100 (MQ=255)
                                                                           aaCAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAg                        <  1:3116049/100‑1 (MQ=255)
                                                                            aCAATCCGCTGAACGGCAGAACCTCGCGGTTCTGCCGTGTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTAc                          >  1:4325327/1‑97 (MQ=255)
                                                                            aCAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGa                       <  1:1419739/100‑1 (MQ=255)
                                                                             cAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGaa                      >  2:3829637/1‑100 (MQ=255)
                                                                              taaaacaaaaaCGGCAGAACCGCGAGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAAc                     >  2:1894763/10‑100 (MQ=37)
                                                                               aaaacaaaaaCGGCAGAACCGCGAGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACt                    >  2:3627180/9‑100 (MQ=37)
                                                                                  cGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTccc                 <  2:4202110/100‑1 (MQ=255)
                                                                                  cGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTccc                 <  2:143470/100‑1 (MQ=255)
                                                                                  cGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTccc                 >  1:3748253/1‑100 (MQ=255)
                                                                                  acaaaaaCGGCAGAACCGCGAGGTTCTGCGGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGaagatcg                 >  1:1253834/6‑95 (MQ=37)
                                                                                   gATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCt                <  1:219446/100‑1 (MQ=255)
                                                                                   caaaaaCGGCAGAACCGCGAGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCt                >  2:1241055/5‑100 (MQ=38)
                                                                                   caaaaaCGGCAGAACCGCGAGGTTCTGCCGTTTTTGTTTCACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCt                >  2:504877/5‑100 (MQ=38)
                                                                                     tGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTTTTTTTTTTTTTAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCttt              >  1:1108095/1‑100 (MQ=255)
                                                                                     tGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCtt               <  1:3563972/99‑1 (MQ=255)
                                                                                       aaCGGCAGAGCCTGGCGGTTCTGCCGTTTTTGTTTTACGTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTcc            <  1:3426940/100‑1 (MQ=255)
                                                                                       aaCGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTcc            <  1:3988506/100‑1 (MQ=255)
                                                                                       aaCGGCAGAACCGCGAGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTcc            >  2:2644486/1‑100 (MQ=255)
                                                                                       aaCGGCAGAACCGCGAGGTTCTGCCGTTTTTGTTTCACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTcc            >  2:1720231/1‑100 (MQ=39)
                                                                                        aCGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCg           >  2:1484643/1‑100 (MQ=255)
                                                                                        aCGGCAGAACCGCGAGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCg           >  2:1524794/1‑100 (MQ=255)
                                                                                          ggCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCGGt         <  2:2075669/100‑1 (MQ=255)
                                                                                          ggCAGAACCGCGAGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCGGt         >  2:3069378/1‑100 (MQ=255)
                                                                                          ggCAGAACCGCGAGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCGGt         >  2:77290/1‑100 (MQ=255)
                                                                                           gCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTCCTTTTCAGAGGAGGGGTGATGAAAACTTTTGAAGTTCTGCAGCCCCTACAGAACTCCCTTCCCg           >  2:1832244/1‑97 (MQ=255)
                                                                                           gCAGAACCGCGAGGTTCTGCCGTTTTTTTTTTTTTTTTTAGAGGAGAGGTGATGAAAAATTTTGATGTTCTGCAGCCACTACAGAACTCCCTTTCCGGTc        >  1:806543/1‑100 (MQ=255)
                                                                                              gtaCCTCGAGGTTCTGCCGTTTTTGTTTTACTTTTAAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGTCACTACAGAACTCCCTTTCCGGTCTAc     <  1:1903044/98‑1 (MQ=255)
                                                                                              gAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCGGTCtcc     >  1:1821861/1‑98 (MQ=255)
                                                                                                 ccTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCGGTCTACCGc  <  1:2731714/100‑1 (MQ=255)
                                                                                                 ccTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCGGTCTACCGc  <  2:2654101/100‑1 (MQ=255)
                                                                                                 ccTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCGGTCTACCGc  >  2:1209857/1‑100 (MQ=255)
                                                                                                   |                                                                                                 
ATGGTGCTGTGAGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCGGTCTACCGC  >  NC_000913/2767568‑2767764

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: