Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,767,670 C→A 18.1% intergenic (+315/‑40) 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,6700CA18.1% 489.8 / 90.5 197intergenic (+315/‑40)rnlB/yfjPCP4‑57 prophage; uncharacterized protein/CP4‑57 prophage; 50S ribosome‑binding GTPase family protein
Reads supporting (aligned to +/- strand):  ref base C (74/87);  new base A (22/14);  total (96/101)
Fisher's exact test for biased strand distribution p-value = 1.39e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.44e-01

Only 100 of 204 total aligned reads displayed.

GTGCTGTGAGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCGGTCTACCGCTATGG  >  NC_000913/2767571‑2767769
                                                                                                   |                                                                                                   
  gCTGTGAGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTcgc                                                                                                     <  1:1611899/98‑1 (MQ=255)
   cTGTGAGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGt                                                                                                  <  2:1072779/100‑1 (MQ=255)
   cTGTGAGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGt                                                                                                  <  2:4361947/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)
        aGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCATCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCGCGAGGTTCTGc                                                                                             <  1:4296238/100‑1 (MQ=255)
         nttCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGcc                                                                                            >  1:3271566/2‑100 (MQ=255)
         gTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCGCGAGGTTCTGcc                                                                                            <  2:1295162/100‑1 (MQ=255)
         gTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCGCGAGGTTCTGcc                                                                                            <  1:3492621/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:2652164/1‑100 (MQ=255)
            cAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGtt                                                                                         >  2:58681/1‑100 (MQ=255)
            cAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCGCGAGGTTCTGCCGtt                                                                                         <  2:1595142/100‑1 (MQ=255)
             aGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCGCGAGGTTCTGCCGttc                                                                                        <  1:1504332/100‑2 (MQ=255)
                 aGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGt                                                                                    <  2:4260995/100‑1 (MQ=255)
                  gATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCGCGAGGTTCTGCCGttcatcgg                                                                                   <  2:2690233/100‑7 (MQ=38)
                    ttCTGCAATCTATGTGCATGGAAAAAATCGTCAGATAGATATCATTCGTCTGAACAATTCGATAAACGGCAGAACCTCGCGGTTTCGGCGTTTTTGtttt                                                                                 >  1:448421/1‑100 (MQ=255)
                      cTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTAc                                                                               <  1:1649816/100‑1 (MQ=255)
                       tGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACt                                                                              >  2:2423763/1‑100 (MQ=255)
                         cAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACttt                                                                            <  2:2733723/100‑1 (MQ=255)
                            tCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCa                                                                         <  2:840098/100‑1 (MQ=255)
                               aggtgCATTGAAAAAATCTTCAGATAGATATCAGTCTCCTGAACAATCCGATGAACCGCAGAACCTCGCGGTTCTGCCGTGTTTGATTTACTGTTCag                                                                        <  1:2543490/96‑1 (MQ=255)
                                tgtgCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAgagg                                                                     <  2:147544/100‑1 (MQ=255)
                                tgtgCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAgagg                                                                     <  1:3226291/100‑1 (MQ=255)
                                tgtgCATTGAAAAAATCTTCAGATAGATATCATTCTCCCGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAgagg                                                                     <  2:3439619/100‑1 (MQ=255)
                                 gtgCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAgagg                                                                     <  2:2451729/99‑1 (MQ=255)
                                   gCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAgt                                                                  >  2:2145584/1‑100 (MQ=255)
                                    cATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAgtg                                                                 >  2:815408/1‑100 (MQ=255)
                                    cATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAgtg                                                                 >  1:4028745/1‑100 (MQ=255)
                                     aTTGAAAAAATCTTCAGATAGATATCATTCTCCTTAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAgtgg                                                                <  1:596998/100‑1 (MQ=255)
                                     aTTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAgtgg                                                                <  2:133181/100‑1 (MQ=255)
                                     aTTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAgtgg                                                                <  2:620103/100‑1 (MQ=255)
                                        gAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGtga                                                             >  1:2479718/1‑100 (MQ=255)
                                           aaaaTCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGtgatga                                                          <  1:1455194/100‑1 (MQ=255)
                                                   cAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAAtttt                                                  >  1:3404555/1‑100 (MQ=255)
                                                   cAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAAtttt                                                  >  2:3054428/1‑100 (MQ=255)
                                                   cAGATAGATATCATTATCCTTAACAATCCGATGAACGGCAGAATCTCGCGGTTCTGCAGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAAtttt                                                  <  1:2083879/100‑1 (MQ=255)
                                                    agatagatATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTg                                                 <  1:3131159/100‑1 (MQ=255)
                                                    agatagatATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTg                                                 >  1:1923674/1‑100 (MQ=255)
                                                    agatagatATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTg                                                 >  2:3698908/1‑100 (MQ=255)
                                                      atagataTCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGaa                                               <  2:1359279/100‑1 (MQ=255)
                                                       tagataTCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAg                                              <  2:2677614/100‑1 (MQ=255)
                                                       tagataTCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAg                                              <  2:901870/100‑1 (MQ=255)
                                                       tagataTCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAg                                              >  2:3791736/1‑100 (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)
                                                              cATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTATCAGACGAGTGGTGCTGAAAAATTTTGAAGTTCTGCa                                       >  1:2968279/1‑100 (MQ=255)
                                                               aTTCTCTTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAg                                      <  1:4269260/100‑1 (MQ=255)
                                                               aTTCTCCTGAACAATCCGATGAACGGCAGAACGTCGCGGTTCTGCCGTTTTTGTTTGACTTTTAAGAGGAGTGGCAATGAAAAAATTTGAAGTTATGCAg                                      >  1:2138245/1‑100 (MQ=255)
                                                               aTTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAg                                      <  1:2300816/100‑1 (MQ=255)
                                                                 tctcCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGcc                                    <  2:2846321/100‑1 (MQ=255)
                                                                  ctcCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCa                                   >  2:3314833/1‑100 (MQ=255)
                                                                   tcCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCAc                                  <  2:2545321/100‑1 (MQ=255)
                                                                     cTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGCTTATGTTTTACTTTTCAGAGGAGAGGTGATGAAAAATTTTCAAGTTCTGCAGCCACTa                                >  1:156371/1‑100 (MQ=255)
                                                                       gAACAATCCGATGAACGGAAGAACCTGTCGGTTCTGCCGTTTTTGTTTTACTGTTCAGAGGAGTGGTGATGAAAAATATTGAAGTTATGCAGCCAGTACa                              <  1:392085/100‑1 (MQ=255)
                                                                        aaCAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAg                             <  1:3116049/100‑1 (MQ=255)
                                                                         aCAATCCGCTGAACGGCAGAACCTCGCGGTTCTGCCGTGTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTAc                               >  1:4325327/1‑97 (MQ=255)
                                                                         aCAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGa                            <  1:3971152/100‑1 (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                      >  1:3748253/1‑100 (MQ=255)
                                                                               cGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTccc                      <  2:143470/100‑1 (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)
                                                                                    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)
                                                                                         cAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACACCCTTTCCGGTCt            >  1:3198637/1‑100 (MQ=255)
                                                                                           gtaCCTCGAGGTTCTGCCGTTTTTGTTTTACTTTTAAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGTCACTACAGAACTCCCTTTCCGGTCTAc          <  1:1903044/98‑1 (MQ=255)
                                                                                           gAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCGGTCtcc          >  1:1821861/1‑98 (MQ=255)
                                                                                             aCCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGATGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCGGTCTAcc         <  1:4225370/99‑1 (MQ=255)
                                                                                              ccTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCGGTCTACCGc       <  1:2731714/100‑1 (MQ=255)
                                                                                              ccTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCGGTCTACCGc       >  2:1209857/1‑100 (MQ=255)
                                                                                                  gcgGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCGGTCTACCGCTATg   <  2:4202791/100‑1 (MQ=255)
                                                                                                  gcgGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCGGTCTACCGCTATg   <  2:4185220/100‑1 (MQ=255)
                                                                                                   cgGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCGGTCTACCGCTATgg  <  2:2479730/100‑1 (MQ=255)
                                                                                                   cgGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCGGTCTACCGCTATgg  <  1:3843629/100‑1 (MQ=255)
                                                                                                   |                                                                                                   
GTGCTGTGAGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCGGTCTACCGCTATGG  >  NC_000913/2767571‑2767769

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: