Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,767,667 T→G 16.8% 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,6670TG16.8% 394.3 / 75.1 164intergenic (+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 (69/67);  new base G (17/11);  total (86/78)
Fisher's exact test for biased strand distribution p-value = 4.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.32e-01

Only 100 of 166 total aligned reads displayed.

ATGGTGCTGTGAGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCGGTCTACCGCTA  >  NC_000913/2767568‑2767766
                                                                                                   |                                                                                                   
 tGGTGCTGTGAGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTc                                                                                                    >  2:1620263/1‑100 (MQ=255)
 tGGTGCTGTGAGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTc                                                                                                    >  1:1360929/1‑100 (MQ=255)
 tGGTGCTGTGAGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTc                                                                                                    >  1:252806/1‑100 (MQ=255)
   gTGCTGTGAGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGACCCTcgc                                                                                                  <  2:1044836/100‑1 (MQ=255)
   gTGCTGTGAGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTcgc                                                                                                  <  2:1967494/100‑1 (MQ=255)
   gTGCTGTGAGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTcgc                                                                                                  >  2:3353255/1‑100 (MQ=255)
     gCTGTGAGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCgg                                                                                                >  2:1714410/1‑100 (MQ=255)
     gCTGTGAGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCgg                                                                                                <  1:1846384/100‑1 (MQ=255)
       tgtgAGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGtt                                                                                              >  1:1172016/1‑100 (MQ=255)
       tgtgAGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGt                                                                                               <  2:252446/99‑1 (MQ=255)
       tctgAGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCGCGAGGtt                                                                                              <  1:2329685/98‑1 (MQ=255)
        gtgAGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTc                                                                                             >  2:1295156/1‑100 (MQ=255)
        gtgAGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTc                                                                                             >  1:2266806/1‑100 (MQ=255)
         tgAGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCt                                                                                            >  2:1925675/1‑100 (MQ=255)
           aGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGc                                                                                          <  2:1253761/100‑1 (MQ=255)
           aGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGc                                                                                          >  1:1903201/1‑100 (MQ=255)
            gTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGcc                                                                                         >  1:2477028/1‑100 (MQ=255)
            gTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCGCGAGGTTCTGcc                                                                                         <  1:3129688/100‑1 (MQ=255)
            gTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCGCGAGGTTCTGcc                                                                                         <  1:1034585/100‑1 (MQ=255)
              tCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGt                                                                                       <  2:2094128/100‑1 (MQ=255)
              tCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCGCGAGGTTCTGCCGt                                                                                       <  2:2364865/100‑1 (MQ=255)
               ctggaaGATTCTGCAATCTATGTGCATTGAAAAAATCTGCAGATAGATATCATTCTCCTGAACAATCCTATGAACGGCAGAACCTCGCGGTTCTGCCGtt                                                                                      >  1:495875/5‑100 (MQ=255)
               cAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGtt                                                                                      >  2:2281821/1‑100 (MQ=255)
               cAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGt                                                                                       >  1:2400291/1‑99 (MQ=255)
                 gAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGtttt                                                                                    >  1:2703509/1‑100 (MQ=255)
                   aaGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGATCAATCCGATGAACGGCAGAACCGCGAGGTTCTGCCGt                                                                                       <  1:575131/95‑1 (MQ=255)
                    aGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTg                                                                                  >  1:3143160/1‑99 (MQ=255)
                    aGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGt                                                                                 >  1:3329587/1‑100 (MQ=255)
                      aTTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGttt                                                                               <  1:2574502/100‑1 (MQ=255)
                      aTTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGttt                                                                               <  2:2321332/100‑1 (MQ=255)
                          tGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTTTTTTACt                                                                           >  1:1096106/1‑100 (MQ=255)
                           gCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACtt                                                                          >  1:2417046/1‑100 (MQ=255)
                           gCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACtt                                                                          <  2:1706931/100‑1 (MQ=255)
                            cAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACttt                                                                         >  1:1743269/1‑100 (MQ=255)
                             aaTCTATGGGCATGGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTGTACtttt                                                                        >  1:1511770/1‑100 (MQ=255)
                              aTCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTc                                                                       <  1:228537/100‑1 (MQ=255)
                              aTCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTGACTTTTc                                                                       <  2:1773452/100‑1 (MQ=255)
                                cTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCag                                                                     <  1:1546269/100‑1 (MQ=255)
                                  aTGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCagag                                                                   <  2:2399761/100‑1 (MQ=255)
                                   tgtgCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAgagg                                                                  <  2:2840723/100‑1 (MQ=255)
                                     tgCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAgaggag                                                                <  2:2987529/100‑1 (MQ=255)
                                       cATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAAGCCGATGAACGGCAGAACCGCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAgtg                                                              <  1:552858/100‑1 (MQ=255)
                                        aTTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAgtgg                                                             >  1:1809712/1‑100 (MQ=255)
                                         ttGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAgtggt                                                            >  1:2879386/1‑100 (MQ=255)
                                          tGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAgtggtg                                                           >  1:1889930/1‑100 (MQ=255)
                                           gAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGtga                                                          >  1:315670/1‑100 (MQ=255)
                                            aaaaaaTCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAgtggtg                                                           <  1:2025313/98‑1 (MQ=255)
                                            aaaaaaTCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAgtggt                                                            <  1:14406/97‑1 (MQ=255)
                                            aaaaaaTCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGtgat                                                         >  2:643544/1‑100 (MQ=255)
                                             tccgaTCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGtgatg                                                        <  1:678338/96‑1 (MQ=255)
                                             aaaaaTCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGtgatg                                                        >  2:2150020/1‑100 (MQ=255)
                                              aaaaTCTTCAGATAGATATCATTCTCCTTCACATTCCGATGAACTGCAGAACCTCGCGGTTCTGGCGTTTTTGTTTTACTTTTCAGAGTAgtgg                                                             <  1:2610740/94‑1 (MQ=255)
                                                aaTCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGaaa                                                     <  2:2555303/100‑1 (MQ=255)
                                                     tCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGtgatgagatcgga                                                >  2:678341/1‑93 (MQ=255)
                                                      cAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAAtttt                                               >  2:1106230/1‑100 (MQ=255)
                                                       agatagatATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTg                                              >  2:353604/1‑100 (MQ=255)
                                                       agatagatATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTg                                              <  2:3127384/100‑1 (MQ=255)
                                                       agatagatATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGATTTACTTTTCAGAGGAGTGGTGATGAAAAAtttt                                               >  1:604116/1‑99 (MQ=255)
                                                          tagataTTGTTCTCCTGAACATTCCGATGTTCTGCTTAACAGAGCGGTTCTGACGTTTGTGTTTTACTTTTCAGACGAGTGGTGATGAAAAATTTTGAAg                                           <  1:769506/100‑1 (MQ=255)
                                                          tagataTCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAg                                           <  2:1586088/100‑1 (MQ=255)
                                                            gataTCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTCTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGtt                                         >  1:806695/1‑100 (MQ=255)
                                                             atatCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTc                                        >  2:864989/1‑100 (MQ=255)
                                                               attaTTCTCCTGAACTATCTGATGTAGGGCAAAGGCTAGCGGTTCTGCCGTTTTTTTTTTCATTTTCATAGTTGTGGTGATGTAAAATTTTGAAGTGCTg                                      <  1:2201531/97‑1 (MQ=255)
                                                                tCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGc                                     >  2:1376331/1‑100 (MQ=255)
                                                                tCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGc                                     >  1:473593/1‑100 (MQ=255)
                                                                tCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGc                                     <  2:1860105/100‑1 (MQ=255)
                                                                 cATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCa                                    <  1:3239453/100‑1 (MQ=255)
                                                                        cTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAGTTTTGAAGTTCTGCAGCCACTa                             >  2:2397388/1‑100 (MQ=255)
                                                                        cTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTa                             <  1:1246215/100‑1 (MQ=255)
                                                                         tGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGACGTTCTGCAGCCACTAc                            >  1:366259/1‑100 (MQ=255)
                                                                            aCAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGa                         <  2:2944380/100‑1 (MQ=255)
                                                                              aaTCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAAc                       <  1:891681/100‑1 (MQ=255)
                                                                              aaTCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTGTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCAATACAGa                         <  1:2326467/98‑1 (MQ=255)
                                                                                  cGATGAACTGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTccc                   <  1:215488/100‑1 (MQ=255)
                                                                                  cGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTccc                   >  2:546899/1‑100 (MQ=255)
                                                                                  cGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTccc                   >  1:189467/1‑100 (MQ=255)
                                                                                  acaaaaaCGGCAGAACCGCGAGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTccc                   >  1:517219/6‑100 (MQ=38)
                                                                                   gATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCt                  <  2:1080371/100‑1 (MQ=255)
                                                                                    aTGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCt                  >  2:2555228/1‑99 (MQ=255)
                                                                                     tGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGCAGTGGTGATGAAAAATTTTGAAATTCTGCAGCCACTACAGAACTCCCtt                 <  1:95626/99‑1 (MQ=255)
                                                                                       aaCGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTcc              <  2:1448049/100‑1 (MQ=255)
                                                                                        aCGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTTTGCAGCCACTACAGAACTCCCTTTCCg             >  1:3374120/1‑100 (MQ=255)
                                                                                         cGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCgg            <  2:1815718/100‑1 (MQ=255)
                                                                                         cGGCAGAACCGCGAGGTTCTGCCGTTTTTTTTTTTTTTTTCAGAGGAGGGGAGGGGTAAACTATTGTAGTTCTCCTGACCCTACAGAAc                       >  1:2086292/1‑89 (MQ=9)
                                                                                          ggCAGAACCGCGAGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCGGGt           >  2:153215/1‑100 (MQ=255)
                                                                                          ggCAGAACCGCGAGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCGGt           >  2:423021/1‑100 (MQ=255)
                                                                                           gCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCGGTc          >  1:259828/1‑100 (MQ=255)
                                                                                           gCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCGGTc          >  2:454381/1‑100 (MQ=255)
                                                                                           gCAGAACCGCGAGGTTCTGCCGTTTTTTTTTTTTTTTTTATAGGAGCGCTGATCTTATATTTTGTCGTTCTTCAGCCGCTACAGAGCTGCCt                  >  1:902301/1‑92 (MQ=9)
                                                                                           gCAGAACCGCGAGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCGGTc          >  2:1458390/1‑100 (MQ=255)
                                                                                           gCAGAACCGCGAGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCGGTc          >  2:2476820/1‑100 (MQ=255)
                                                                                           gCAGAACCGCGAGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCGGTc          >  2:840099/1‑100 (MQ=255)
                                                                                            cAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCTCTACAGAACTCCCTTTCCGGTCt         <  1:422373/100‑1 (MQ=255)
                                                                                            cAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCGGTCt         >  2:776715/1‑100 (MQ=255)
                                                                                                 ccTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCGGTCTAcc      <  2:894992/98‑1 (MQ=255)
                                                                                                 ccTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCGGTCTACCGc    <  2:2947186/100‑1 (MQ=255)
                                                                                                 ccTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTAATGAAAAATTTCGAAGTTCCGCAGCCACTACAGAACTCCCTTTCCCGTCTACCGc    <  1:1555936/100‑1 (MQ=255)
                                                                                                  cTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCGGTCTACCGCt   <  1:2887283/100‑1 (MQ=255)
                                                                                                  cTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCGGTCTACCGCt   >  1:3054358/1‑100 (MQ=255)
                                                                                                   tCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCGGTCTACCGCTa  >  1:1288544/1‑100 (MQ=255)
                                                                                                   |                                                                                                   
ATGGTGCTGTGAGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCGGTCTACCGCTA  >  NC_000913/2767568‑2767766

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: