Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,767,670 C→A 15.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,6700CA15.1% 459.2 / 61.9 172intergenic (+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 (71/74);  new base A (15/11);  total (86/86)
Fisher's exact test for biased strand distribution p-value = 5.24e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.12e-01

Only 100 of 179 total aligned reads displayed.

GTGCTGTGAGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCGGTCTACCGCTATGG  >  NC_000913/2767571‑2767769
                                                                                                   |                                                                                                   
  gCTGTGAGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCgg                                                                                                   <  2:3889401/100‑1 (MQ=255)
  gCTGTGAGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCgg                                                                                                   >  2:1012558/1‑100 (MQ=255)
  gCTGTGAGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCgg                                                                                                   <  2:1666705/100‑1 (MQ=255)
   cTGTGAGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGt                                                                                                  >  2:2950064/1‑100 (MQ=255)
    tgtgAGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGtt                                                                                                 <  2:1137291/100‑1 (MQ=255)
    tgtgAGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGtt                                                                                                 <  2:1233209/100‑1 (MQ=255)
     gtgAGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTc                                                                                                >  2:316848/1‑100 (MQ=255)
     gtgAGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTc                                                                                                >  2:3069161/1‑100 (MQ=255)
       gAGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTg                                                                                              <  2:610989/100‑1 (MQ=255)
       gAGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTg                                                                                              <  2:2610235/100‑1 (MQ=255)
       gAGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTg                                                                                              <  2:1448112/100‑1 (MQ=255)
       gAGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTg                                                                                              <  2:1336943/100‑1 (MQ=255)
        aGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGc                                                                                             <  2:1429115/100‑1 (MQ=255)
         gTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGcc                                                                                            >  2:2676685/1‑100 (MQ=255)
         gTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGcc                                                                                            >  1:466335/1‑100 (MQ=255)
         gTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGcc                                                                                            >  2:3357256/1‑100 (MQ=255)
         gTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGcc                                                                                            <  2:850341/100‑1 (MQ=255)
          ttCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCGCGAGGTTCTGCCg                                                                                           <  2:3622693/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)
            cAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCGCGAGGTTCTGCCGtt                                                                                         <  2:3153261/100‑1 (MQ=255)
             aGAAAGATTCTGCAATCTATGTGCATTGAAAATATCTTCAGCTAGATATCATTGTCCTGAACAATCCGATGAACGGCAGAACCGCGAGGTTCTGCCGt                                                                                          <  1:1715965/98‑1 (MQ=38)
              gAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGTCGtttt                                                                                       <  1:741396/100‑1 (MQ=255)
              gAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTAGCCGTTCTGCCGtttt                                                                                       >  1:3626789/1‑100 (MQ=255)
                aaGATTCTGCAATCTATGTGCATTGAAAACATCTTCAGACAGATATCATTCTCCTGCACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTg                                                                                     <  1:3127182/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)
                     tCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCGCGAGGTTCTGCCGttcatcggatt                                                                                <  2:2544370/100‑10 (MQ=37)
                      cTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTAc                                                                               <  2:2084776/100‑1 (MQ=255)
                      cTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTAc                                                                               <  2:262698/100‑1 (MQ=255)
                           aTCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTc                                                                          >  1:942119/1‑100 (MQ=255)
                            tCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCa                                                                         <  2:1544174/100‑1 (MQ=255)
                               aTGTGAAATGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCag                                                                        <  1:2465987/98‑1 (MQ=255)
                                tgtgCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAgagg                                                                     >  1:539248/1‑100 (MQ=255)
                                tgtgCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAgagg                                                                     <  2:3690433/100‑1 (MQ=255)
                                  tgCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAgaggag                                                                   >  2:34070/1‑100 (MQ=255)
                                    cATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAgtg                                                                 >  2:3259247/1‑100 (MQ=255)
                                     aTTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAgtgg                                                                <  2:799437/100‑1 (MQ=255)
                                         aaaaaaTCTTTAGATAGATAACCTTCTCCTGACCAATCGTATGTACTGCAGAACCTTGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGTAGTGGtgat                                                            <  1:810170/100‑1 (MQ=255)
                                         aaaaaaTCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGtgat                                                            >  1:1438999/1‑100 (MQ=255)
                                               tCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCATAGGAGTGGTGATGaaaaa                                                      >  1:2092571/1‑100 (MQ=255)
                                                    agatagatATCATTCTCCTGAACAATCCGATGATAGGCAGCAGCTCGGGGTTCTGCCGTTTTTGTTTTACTTTTCATAGGAGTGGTGATGAAAAAtttt                                                  <  1:2486709/99‑1 (MQ=255)
                                                    agatagatATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTg                                                 <  2:226880/100‑1 (MQ=255)
                                                    agatagatATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTg                                                 >  2:1718629/1‑100 (MQ=255)
                                                     gatagataTCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGa                                                >  2:3571771/1‑100 (MQ=255)
                                                     gatagataTCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGa                                                >  2:2342205/1‑100 (MQ=255)
                                                     gatagataTAATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGa                                                <  1:4246462/100‑1 (MQ=255)
                                                       tagataTCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAg                                              <  1:2518057/100‑1 (MQ=255)
                                                             tCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGc                                        >  2:2492494/1‑100 (MQ=255)
                                                              cATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGGGGTGATGAAAAATTTTGAAGTTc                                           >  1:3662543/1‑96 (MQ=255)
                                                               aTTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAg                                      <  2:1096575/100‑1 (MQ=255)
                                                               aTTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAg                                      >  2:2451800/1‑100 (MQ=255)
                                                                  ctcCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTATTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCa                                       >  1:1796730/1‑96 (MQ=255)
                                                                    ccTTAACAAGCCGATGTACTGCAGAACCTCGCGGTTCTTCCGTTTTTGTTTTACTTTTTAGAGCAGTGGTGATGAAAAATTTTGAATTTCTGCAGCCACt                                 <  2:3715122/100‑1 (MQ=255)
                                                                    ccTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACt                                 >  2:4120225/1‑100 (MQ=255)
                                                                    ccTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACt                                 <  2:2124971/100‑1 (MQ=255)
                                                                     cTGAACAATCCGATGACCGGCCGAACCTCACCGTTCAACAGTTGTTGTTTTACTTTTCAGAGGAGTCGTGATGAAAAATTTTGAAGTTCTGCAGc                                     <  1:1069879/95‑1 (MQ=255)
                                                                      tGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTTTTTTTCAAAGGAGAGGTGATGACACATTTTTAAGTTTTGCAGTCACTAc                               >  1:1593217/1‑100 (MQ=255)
                                                                      tGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGCGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTa                                >  2:1911610/1‑99 (MQ=255)
                                                                      tGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTAc                               <  2:4141844/100‑1 (MQ=255)
                                                                       gAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGGGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACa                              >  1:2490853/1‑100 (MQ=255)
                                                                       gAACAATCCGAGGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACa                              <  1:3540031/100‑1 (MQ=255)
                                                                          cAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGaa                           >  1:1327382/1‑100 (MQ=255)
                                                                          cAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGaa                           >  1:890113/1‑100 (MQ=255)
                                                                            aTCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACt                         <  2:401404/100‑1 (MQ=255)
                                                                             tCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTc                        >  2:4253318/1‑100 (MQ=255)
                                                                              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)
                                                                                 aaaaaCGGCAGAACCGCGAGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTACCt                     >  1:265642/4‑99 (MQ=38)
                                                                                 aaaaaCGGCAGAACCGCGAGGTTCTGCCGTTTTTGTTGTACTTTTCAGAGGGGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCtt                    >  1:2268949/4‑100 (MQ=38)
                                                                                 aTGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCtt                    <  2:2393528/100‑1 (MQ=255)
                                                                                     aCGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCg                <  2:550262/100‑1 (MQ=255)
                                                                                      ccgCAGAACCGCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCgg               <  1:1864180/98‑1 (MQ=255)
                                                                                      cGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCgg               <  1:1519256/100‑1 (MQ=255)
                                                                                      cGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCgg               <  1:1740028/100‑1 (MQ=255)
                                                                                      cGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCgg               >  2:1913562/1‑100 (MQ=255)
                                                                                      cGGCAGAACCGCGAGGTTCTGCCGTTTTTGTTTTTTTTTTTTTTGGAGTGGTGATGACAAATTTTGAAGTTTTGCAGCCACTCCAGAACTCCCt                     >  1:250477/1‑94 (MQ=21)
                                                                                       ggCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCGGt              >  2:3297250/1‑100 (MQ=255)
                                                                                       ggCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCGGt              >  1:3069932/1‑100 (MQ=255)
                                                                                       ggCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCGGt              >  2:1391850/1‑100 (MQ=255)
                                                                                       ggCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCGGt              <  2:1153056/100‑1 (MQ=255)
                                                                                       ggCAGAACCGCGAGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCGGt              >  2:3029164/1‑100 (MQ=255)
                                                                                       ggCAGAACCGCGAGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCGGt              >  2:1280771/1‑100 (MQ=255)
                                                                                        gCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCGGTc             <  1:3259235/100‑1 (MQ=255)
                                                                                        gCAGAACCGCGAGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGGTGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCGGTc             >  1:4188120/1‑100 (MQ=255)
                                                                                         cAGAACCTCGCGGTTCTGCCGTTTTTTTTTTTTTTTTTAGAGGAGTGGTGAAGAAAAATTTTTAATTTTTGTTGCATCTACAGAAGTACTTTTCCGGTCt            >  1:4121933/1‑100 (MQ=255)
                                                                                             aCCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCGGTCTACCg        <  2:1523526/100‑1 (MQ=255)
                                                                                             aCCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCGGTCTACCg        >  2:3489697/1‑100 (MQ=255)
                                                                                             aCCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAAATCCCTTTCCGGTCTACCg        >  1:2395711/1‑100 (MQ=255)
                                                                                              ccTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCGGTCTACCGc       <  1:173223/100‑1 (MQ=255)
                                                                                               cTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCGGTCTACCGc       >  1:979859/1‑99 (MQ=255)
                                                                                                  gcgGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCGGTCTACCGCTATg   <  1:599485/100‑1 (MQ=255)
                                                                                                  gcgGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCGGTCTACCGCTATg   <  2:3069944/100‑1 (MQ=255)
                                                                                                   cgGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCGGTCTACCGCTATgg  <  1:1357428/100‑1 (MQ=255)
                                                                                                   cgGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCGGTCTACCGCTATgg  <  1:63762/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: