Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,767,667 T→G 15.1% 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,6670TG15.1% 441.6 / 71.1 178intergenic (+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/75);  new base G (15/12);  total (90/88)
Fisher's exact test for biased strand distribution p-value = 6.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.58e-01

Only 100 of 179 total aligned reads displayed.

ATGGTGCTGTGAGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCGGTCTACCGC  >  NC_000913/2767568‑2767764
                                                                                                   |                                                                                                 
  ggTGCTGTGAGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTcg                                                                                                 >  1:573552/1‑100 (MQ=255)
  ggTGCTGTGAGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTcg                                                                                                 >  1:78684/1‑100 (MQ=255)
   gTGCTGTGAGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTcgc                                                                                                <  2:3339033/100‑1 (MQ=255)
     gCTGTGAGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCgg                                                                                              >  2:1713859/1‑100 (MQ=255)
      cTGTGAGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGt                                                                                             >  2:2950064/1‑100 (MQ=255)
        gtgAGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTc                                                                                           >  2:3069161/1‑100 (MQ=255)
        gtgAGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGAGAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTc                                                                                           >  2:2674078/1‑100 (MQ=255)
         tgAGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCt                                                                                          <  2:3978679/100‑1 (MQ=255)
          gAGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTg                                                                                         <  2:610989/100‑1 (MQ=255)
           aGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCGCGAGGTTCTGc                                                                                        <  2:3397005/100‑1 (MQ=255)
            gTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGcc                                                                                       <  2:850341/100‑1 (MQ=255)
            gTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGcc                                                                                       >  1:466335/1‑100 (MQ=255)
            gTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGcc                                                                                       >  1:850339/1‑100 (MQ=255)
             ttCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCg                                                                                      <  2:1273026/100‑1 (MQ=255)
             ttCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCg                                                                                      <  2:1525142/100‑1 (MQ=255)
             ttCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCGCGAGGTTCTGCCg                                                                                      <  2:3622693/100‑1 (MQ=255)
             ttCAGAAAGATTCTGCAATAAAAGAAAAATGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGcc                                                                                       <  1:3677931/99‑1 (MQ=255)
              tCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGt                                                                                     <  2:2404024/100‑1 (MQ=255)
              tCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGt                                                                                     >  1:1098403/1‑100 (MQ=255)
              tCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCGCGAGGTTCTGCCg                                                                                      <  2:683299/99‑1 (MQ=255)
              tCAGAAAGATTCTGCAATCAAAATGAATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGt                                                                                     <  1:2650921/100‑1 (MQ=255)
                aGAAAGATTCTGCAATCTATGTGCATTGAAAATATCTTCAGCTAGATATCATTGTCCTGAACAATCCGATGAACGGCAGAACCGCGAGGTTCTGCCGt                                                                                     <  1:1715965/98‑1 (MQ=38)
                aGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGttt                                                                                   >  2:355185/1‑100 (MQ=255)
                aGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAAACCGATGAACGGCAGAACCTCGCGGTTCTGCCGttt                                                                                   >  2:2099057/1‑100 (MQ=255)
                 gAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGtttt                                                                                  >  1:1482927/1‑100 (MQ=255)
                   aaGATTCTGCAATCTATGTGCATTGAAAACATCTTCAGACAGATATCATTCTCCTGCACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTg                                                                                <  1:3127182/100‑1 (MQ=255)
                   aaGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTg                                                                                <  2:1227487/100‑1 (MQ=255)
                   aaGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACAGCAGAAAATCGCGGTTCTGCCGTTTTTg                                                                                >  1:2960549/1‑100 (MQ=255)
                   aaGATTCTGCAATCTATATGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTg                                                                                <  1:1166536/100‑1 (MQ=255)
                    aGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGt                                                                               >  2:1185764/1‑100 (MQ=255)
                       ttCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGtttt                                                                            <  2:1402718/100‑1 (MQ=255)
                         cTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTAc                                                                          <  2:262698/100‑1 (MQ=255)
                          tGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACt                                                                         >  2:1371537/1‑100 (MQ=255)
                              aTCTATGTGCATTGAAAAAATCTTCAGATAGAGATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTc                                                                     <  2:2520973/100‑1 (MQ=255)
                                  aTGTGAAATGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCag                                                                   <  1:2465987/98‑1 (MQ=255)
                                   tgtgCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAgagg                                                                <  2:161714/100‑1 (MQ=255)
                                     tgtagtGAAAAAATCTTCAGATATATATCATGCTCCTGAACAATCCGATGAACGGCAGATCCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAgaggag                                                              <  1:3638966/95‑1 (MQ=38)
                                       cATTGAAACAATCTTCAGATACATATCATTCTCCTGAACAATCCGATGAACGGCAGAACATCGCGGTTATGCCGTTATTGTTTTACTTTTCAgaggag                                                              >  1:1748964/1‑98 (MQ=255)
                                       cATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAgtg                                                            <  2:3158463/100‑1 (MQ=255)
                                        aTTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAgtgg                                                           <  2:799437/100‑1 (MQ=255)
                                           gAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGtga                                                        >  2:425882/1‑100 (MQ=255)
                                           gAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGtga                                                        >  1:249413/1‑100 (MQ=255)
                                             aaaaaTCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGtgatg                                                      <  1:905515/100‑1 (MQ=255)
                                             aaaaaTCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGtgatg                                                      >  1:763343/1‑100 (MQ=255)
                                                aaTCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGaaa                                                   <  2:1058900/100‑1 (MQ=255)
                                                  tCTTCATCTAGATATCATTCTCCTCACCATTTCGATGAACGGCATCACTTCGCGTTTCTGCCGTTTTTGTTTTACTTTTCAGAGTATTGGTGATGaa                                                    <  1:989722/97‑1 (MQ=255)
                                                  tCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCATAGGAGTGGTGATGaaaaa                                                 >  1:2092571/1‑100 (MQ=255)
                                                   cTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAAt                                                >  1:61413/1‑100 (MQ=255)
                                                     tCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAAttt                                              <  2:1445207/100‑1 (MQ=255)
                                                      cAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAAtttt                                             <  2:3435652/100‑1 (MQ=255)
                                                       agatagatATCATTCTCCTGAACAATCCGATGATAGGCAGCAGCTCGGGGTTCTGCCGTTTTTGTTTTACTTTTCATAGGAGTGGTGATGAAAAAtttt                                             <  1:2486709/99‑1 (MQ=255)
                                                        gatagataTAATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGa                                           <  1:4246462/100‑1 (MQ=255)
                                                         atagataTCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGATGATGAAAAATTTTGaa                                          <  1:2858231/100‑1 (MQ=255)
                                                          tagataTCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAg                                         <  1:2518057/100‑1 (MQ=255)
                                                          tagataTCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAg                                         <  2:1141779/100‑1 (MQ=255)
                                                          tagataTCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAg                                         <  2:3713514/100‑1 (MQ=255)
                                                            gataTCATTCTCCTGAACAATCCGATGTACGGCAGAACTTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGATGAGCTGTGCTGAAAACTTTTGa                                           >  1:231445/1‑96 (MQ=255)
                                                                 cATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTCACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCa                                  <  2:129440/100‑1 (MQ=255)
                                                                   ttCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAg                                 <  2:4088776/99‑1 (MQ=255)
                                                                     ctcCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTATTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCa                                  >  1:1796730/1‑96 (MQ=255)
                                                                       ccTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACt                            >  2:3104589/1‑100 (MQ=255)
                                                                        cTGAACAATCCGATGACCGGCCGAACCTCACCGTTCAACAGTTGTTGTTTTACTTTTCAGAGGAGTCGTGATGAAAAATTTTGAAGTTCTGCAGc                                <  1:1069879/95‑1 (MQ=255)
                                                                        cTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTGACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACt                            <  2:4122694/99‑1 (MQ=255)
                                                                         tGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTTTTTTTCAAAGGAGAGGTGATGACACATTTTTAAGTTTTGCAGTCACTAc                          >  1:1593217/1‑100 (MQ=255)
                                                                         tGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTAc                          <  2:4141861/100‑1 (MQ=255)
                                                                         tGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTAc                          <  2:1411241/100‑1 (MQ=255)
                                                                          gAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACa                         >  2:2194013/1‑100 (MQ=255)
                                                                          gAACAATCCGAGGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACa                         <  1:3540031/100‑1 (MQ=255)
                                                                             cAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGaa                      >  1:1327382/1‑100 (MQ=255)
                                                                              aaTCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAAc                     >  2:1478200/1‑100 (MQ=255)
                                                                                aaacaaaaaCGGCAGAACCGCGAGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTc                   >  2:1985530/8‑100 (MQ=37)
                                                                                 ccGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGATTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTcc                  <  1:2420863/100‑1 (MQ=255)
                                                                                 ccGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTcc                  <  2:1564574/100‑1 (MQ=255)
                                                                                 ccGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTcc                  >  2:4088217/1‑100 (MQ=255)
                                                                                 ccGATGAACGGCAGAACCTCGCGGTTCTGCCGTGTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTcc                  <  1:2033158/100‑1 (MQ=255)
                                                                                  cGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTccc                 <  2:2452817/100‑1 (MQ=255)
                                                                                    aTGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCtt               <  2:2393528/100‑1 (MQ=255)
                                                                                     tGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCttt              >  2:2594638/1‑100 (MQ=255)
                                                                                        aCGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCg           <  2:550262/100‑1 (MQ=255)
                                                                                         ccgCAGAACCGCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCgg          <  1:1864180/98‑1 (MQ=255)
                                                                                         cGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCgg          >  2:3256977/1‑100 (MQ=255)
                                                                                         cGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCgg          <  1:1519256/100‑1 (MQ=255)
                                                                                         cGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCgg          >  2:1913562/1‑100 (MQ=255)
                                                                                         cGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCgg          >  2:2869560/1‑100 (MQ=255)
                                                                                         cGGCAGAACCGCGAGGTTCTGCCGTTTTTGTTTTTTTTTTTTTTGGAGTGGTGATGACAAATTTTGAAGTTTTGCAGCCACTCCAGAACTCCCt                >  1:250477/1‑94 (MQ=21)
                                                                                          ggCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCGGt         <  2:2313125/100‑1 (MQ=255)
                                                                                          ggCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCGGt         >  1:3069932/1‑100 (MQ=255)
                                                                                          ggCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCGGt         >  2:3297250/1‑100 (MQ=255)
                                                                                          ggCAGAACCGCGAGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCgg          >  2:1113500/1‑99 (MQ=255)
                                                                                          ggCAGAACCGCGAGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCGGt         >  2:3029164/1‑100 (MQ=255)
                                                                                          ggCAGAACCGCGAGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCGGt         >  2:1901736/1‑100 (MQ=255)
                                                                                           gCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCGGTc        <  2:3887066/100‑1 (MQ=255)
                                                                                           gCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCGGTc        <  1:3259235/100‑1 (MQ=255)
                                                                                           gCAGAACCGCGAGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGGTGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCGGTc        >  1:4188120/1‑100 (MQ=255)
                                                                                           gCAGAACCGCGAGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCGGTc        >  2:36821/1‑100 (MQ=255)
                                                                                            cAGAACCTCGCGGTTCTGCCGTTTTTTTTTTTTTTTTTAGAGGAGTGGTGAAGAAAAATTTTTAATTTTTGTTGCATCTACAGAAGTACTTTTCCGGTCt       >  1:4121933/1‑100 (MQ=255)
                                                                                                aCCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCGGTCTACCg   <  1:3571758/100‑1 (MQ=255)
                                                                                                aCCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAAATCCCTTTCCGGTCTACCg   >  1:2395711/1‑100 (MQ=255)
                                                                                                 ccTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCGGTCTACCGc  <  1:173223/100‑1 (MQ=255)
                                                                                                  cTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCGGTCTACCGc  >  1:979859/1‑99 (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: