Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,767,670 C→A 16.7% 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,6700CA16.7% 402.1 / 65.5 155intergenic (+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 (69/60);  new base A (17/9);  total (86/69)
Fisher's exact test for biased strand distribution p-value = 2.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.23e-01

Only 100 of 158 total aligned reads displayed.

GTGCTGTGAGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCGGTCTACCGCTAT  >  NC_000913/2767571‑2767767
                                                                                                   |                                                                                                 
gTGCTGTGAGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTcgc                                                                                                   <  1:2496825/100‑1 (MQ=255)
 tGCTGTGAGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTcgcg                                                                                                  <  2:754134/100‑1 (MQ=255)
    tgtgAGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGtt                                                                                               >  2:1315897/1‑100 (MQ=255)
    tgtgAGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGtt                                                                                               >  1:388498/1‑100 (MQ=255)
    tgtcAGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCGGAACAAGCCGGTGAACGGCAGAACCTCGCgg                                                                                                 <  2:2339655/98‑1 (MQ=39)
     gtgAGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGCACAATCCGATGAACGGCAGAACATCGCGGtt                                                                                               >  2:374423/1‑99 (MQ=255)
     gtgAGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTc                                                                                              >  2:3110676/1‑100 (MQ=255)
      tgAGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTc                                                                                              >  2:3371284/1‑99 (MQ=255)
      tgAGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCt                                                                                             >  1:2061689/1‑100 (MQ=255)
       gAGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTg                                                                                            >  2:1144706/1‑100 (MQ=255)
       gAGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCGCGAGGTTCTg                                                                                            <  1:580895/100‑1 (MQ=255)
        aGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGc                                                                                           >  1:1864508/1‑100 (MQ=255)
        aGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGc                                                                                           <  1:3764112/100‑1 (MQ=255)
        aGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCGCGAGGTTCTGc                                                                                           <  1:1774731/100‑1 (MQ=255)
          ttCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCg                                                                                         <  2:1573811/100‑1 (MQ=255)
            cAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGtt                                                                                       >  1:715360/1‑100 (MQ=255)
             aGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGttt                                                                                      >  2:253005/1‑100 (MQ=255)
             aGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGttt                                                                                      >  2:1286385/1‑100 (MQ=255)
             aGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCGCGAGGTTCTGCCGttc                                                                                      <  1:2922237/100‑2 (MQ=255)
              gAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGtttt                                                                                     >  2:1200187/1‑100 (MQ=255)
              gAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACATCGCGGTTCTGCCGtt                                                                                       >  2:2381076/1‑98 (MQ=255)
                 aGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGt                                                                                  <  1:857935/100‑1 (MQ=255)
                  gATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGtt                                                                                 <  1:543282/100‑1 (MQ=255)
                    ttCTGCAATCTATGTGCATTGAAAAAAGCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGtttt                                                                               <  2:2094520/100‑1 (MQ=255)
                      cTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTAc                                                                             >  1:3457084/1‑100 (MQ=255)
                      cTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTAc                                                                             <  1:3216964/100‑1 (MQ=255)
                        gCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACAGCAGAACCTCGCGGTTCTGCAGTTTTTTTTTTAcg                                                                            >  2:2316013/1‑98 (MQ=255)
                           aTCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTc                                                                        <  2:575718/100‑1 (MQ=255)
                            tCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCa                                                                       <  1:1842891/100‑1 (MQ=255)
                            tCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCa                                                                       <  1:2310113/100‑1 (MQ=255)
                                tgtgCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAgagg                                                                   <  1:1531067/100‑1 (MQ=255)
                                tgtgCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAgagg                                                                   <  2:462322/100‑1 (MQ=255)
                                  tgCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAgaggag                                                                 >  2:194917/1‑100 (MQ=255)
                                    cATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAgtg                                                               >  2:1188736/1‑100 (MQ=255)
                                       tGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAgtggtg                                                            >  2:2467920/1‑100 (MQ=255)
                                       tGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAgtggtg                                                            <  1:3356856/100‑1 (MQ=255)
                                        gAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGtga                                                           >  1:38073/1‑100 (MQ=255)
                                        gAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGtga                                                           >  1:2051566/1‑100 (MQ=255)
                                         aaaaaaTCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGtgat                                                          >  2:2746010/1‑100 (MQ=255)
                                          aaaaaTCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGtgatg                                                         <  2:1333557/100‑1 (MQ=255)
                                                  tCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAAttt                                                 <  1:3318581/100‑1 (MQ=255)
                                                   cAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAAtttt                                                <  1:2600977/100‑1 (MQ=255)
                                                   cAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAAtttt                                                <  2:2435147/100‑1 (MQ=255)
                                                    agatagatATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTg                                               >  2:2099458/1‑100 (MQ=255)
                                                    agatagatATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTg                                               <  1:621857/100‑1 (MQ=255)
                                                     gatagataTCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGa                                              <  2:185347/100‑1 (MQ=255)
                                                      atagataTCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGaa                                             <  1:925392/100‑1 (MQ=255)
                                                      atagataTCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGaa                                             <  2:2194802/100‑1 (MQ=255)
                                                       tagataTCATTCTCCTGAACAATCCGATGAACGGAAGAACCCCGCGGTTCTGCCGTTTTTGGTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAg                                            >  2:2080365/1‑100 (MQ=255)
                                                         gataTCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGtt                                          >  1:3666622/1‑100 (MQ=255)
                                                           tatCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCt                                        >  1:1362666/1‑100 (MQ=255)
                                                           tatCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCt                                        >  1:2255370/1‑100 (MQ=255)
                                                             tCATTCTCCTGAACTATCCGATGTACGGGAGAACCTCGCGGTTCTGACTTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGATAATTTTTCAAGTTCTGc                                      <  2:2241340/100‑1 (MQ=255)
                                                              cATTCTCCTGAGCATTCCGATGAGCGGCAGAGCCTCGCGGTTGTGCCGTTTTTGTTGTACTTGTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCa                                     <  2:3100764/100‑1 (MQ=255)
                                                              cATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCa                                     >  1:473398/1‑100 (MQ=255)
                                                                 tctcCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGcc                                  <  1:3249227/100‑1 (MQ=255)
                                                                   tcCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCAc                                >  2:1992037/1‑100 (MQ=255)
                                                                    ccTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACt                               >  1:264295/1‑100 (MQ=255)
                                                                      tGAACAATCCGATGAACGGCAGACCCTCGCGGTGCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTa                              <  2:3313550/99‑1 (MQ=255)
                                                                       gAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTAc                             >  2:3718054/1‑99 (MQ=255)
                                                                       gAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACa                            >  1:324179/1‑100 (MQ=255)
                                                                       gAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACa                            >  1:2460681/1‑100 (MQ=255)
                                                                       gAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACa                            >  2:2696789/1‑100 (MQ=255)
                                                                           aaTCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAAc                        >  2:1902776/1‑100 (MQ=255)
                                                                           aaTCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAAc                        <  2:242936/100‑1 (MQ=255)
                                                                            aTCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTTAAGTTATGCAGCCACTACAGAAc                        >  2:3640192/1‑99 (MQ=255)
                                                                            aTCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACt                       >  2:1455843/1‑100 (MQ=255)
                                                                            aTCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACt                       >  2:2592831/1‑100 (MQ=255)
                                                                             aaacaaaaaCGGCAGAACCGCGAGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACt                       >  2:1441078/8‑99 (MQ=37)
                                                                                gATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCt                   >  1:2810310/1‑100 (MQ=255)
                                                                                  tGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCttt                 <  2:736358/100‑1 (MQ=255)
                                                                                  aaaaCGGCAGAACCGCGAGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCttt                 <  1:1529451/98‑1 (MQ=255)
                                                                                   gAACGGCAGAACCTCGCGGTTCTGCCGTTTTTTTTTTTTTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTc                >  2:1858978/1‑100 (MQ=255)
                                                                                   gAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTAAGAGAAGTGGTGATAAAAAATTTTGAAGTTCTGCAGCCAATACAGAACTACCTTTc                >  2:3415102/1‑100 (MQ=255)
                                                                                   aaaCGGCAGAACCGCGAGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTc                >  1:3581863/2‑100 (MQ=255)
                                                                                    aaCGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTGCCTTTcc               >  2:2048620/1‑100 (MQ=255)
                                                                                    aaCGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTcc               <  2:948118/100‑1 (MQ=255)
                                                                                    aaCGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTcc               >  1:2346080/1‑100 (MQ=255)
                                                                                     aCGGCAGAACCGCGAGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCg              >  2:990873/1‑100 (MQ=255)
                                                                                      cGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCgg             <  2:1790770/100‑1 (MQ=255)
                                                                                      cGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCgg             <  2:890183/100‑1 (MQ=255)
                                                                                      cGGCAGAACCGCGAGGTTCTGCCGTTTTTTTTTTTCTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCgg             >  2:2129152/1‑100 (MQ=39)
                                                                                      cGGCAGAACCGCGAGGTTCTGCCGTTTTTGTTTTTTTTTTTTTTGTAGTGGTGATGAAAAATTTTGATGTTCTTCAGCCACTACAGAACTCCCTTTCCgt             >  2:2710786/1‑99 (MQ=25)
                                                                                       ggCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTTCTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCGGt            >  2:2770675/1‑100 (MQ=255)
                                                                                       ggCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCGGt            >  2:3022918/1‑100 (MQ=255)
                                                                                       ggCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCGGt            >  1:52050/1‑100 (MQ=255)
                                                                                        ncAGAACCGCGAGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCGGTc           >  1:3045901/2‑100 (MQ=255)
                                                                                        gCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCgg             <  2:683777/98‑1 (MQ=255)
                                                                                        gCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCGGt            >  2:152302/1‑99 (MQ=255)
                                                                                        gCAGAACCGCGAGGTTCTGCCGTTTTTTTTTTTTTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCGGTc           >  2:1515184/1‑100 (MQ=38)
                                                                                        gCAGAACCGCGAGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCGGTc           >  1:354685/1‑100 (MQ=255)
                                                                                        gCAGAACCGCGAGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCGGTc           >  1:1061813/1‑100 (MQ=255)
                                                                                        gCAGAACCGCGAGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCGGTc           >  1:3404872/1‑100 (MQ=255)
                                                                                         cAGAACCTCGCTGTTCTGCCGTTGTTGTTTTACTGTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCGGTCt          <  2:3301580/100‑1 (MQ=255)
                                                                                         cAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCGGTCt          >  1:1570731/1‑100 (MQ=255)
                                                                                         cAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCGGTCt          >  1:1570730/1‑100 (MQ=255)
                                                                                          aGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCGGt            <  2:1656908/97‑1 (MQ=255)
                                                                                            aaCCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCTTTTCCGGTCTaac       <  1:2592850/100‑3 (MQ=255)
                                                                                             aCCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCGGTCTACCg      <  1:3070152/100‑1 (MQ=255)
                                                                                                 ngcgGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACGACAGAACTCCCTTTCCGGTCTACCGCTAt  >  2:2291908/2‑100 (MQ=255)
                                                                                                   |                                                                                                 
GTGCTGTGAGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCGGTCTACCGCTAT  >  NC_000913/2767571‑2767767

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