Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,906,597 C→A intergenic (‑30/+256) bglG ← / ← phoU transcriptional antiterminator of the bgl operon/negative regulator of PhoR/PhoB two‑component regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,906,5970CA81.7% 302.7 / 58.4 132intergenic (‑30/+256)bglG/phoUtranscriptional antiterminator of the bgl operon/negative regulator of PhoR/PhoB two‑component regulator
Reads supporting (aligned to +/- strand):  ref base C (12/12);  new base A (60/47);  total (72/60)
Fisher's exact test for biased strand distribution p-value = 6.53e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.52e-01

Only 100 of 135 total aligned reads displayed.

CCGCTGTTGATCATCAATAACCACCACAACATTATTGTTGAGAATTTTGGTGATTTGCATGTTCATAGCAAGGACCTTTTTTATAAACAAAAAAACCCGACTTCACCAGTATTCTCTGGTTATGTCAGGTTTTGCCTGCGAATGCAGTAACAATCCAGTCATTTATTAATGGTTTTTATAACGAACATCC  >  NC_000913/3906502‑3906691
                                                                                               |                                                                                              
ccGCTGTTGATCATCAATAACCACCACAACATTATTGTTGAGAATTTTGGTGATTTGCATGTTCATAGCAAGGACCTTTTTTATAAACAAAAAAAACCGa                                                                                            <  1:2776884/100‑1 (MQ=255)
    tGTTGATCATCAATAACCACCACAACATTATTGTTGAGAATTTTGGTGATTTGCATGTTCATAGCAAGGACCTTTTTTATAAACAAAAAAAACCGACTTc                                                                                        >  1:698432/1‑100 (MQ=255)
     gTTGATCATCAATAACCACCACAACATTATTGTTGAGAATTTTGGTGATTTGCATGTTCATAGCAAGGACCTTTTTTATAAACAAAAAAACCCGACTTCa                                                                                       >  2:1677850/1‑100 (MQ=255)
     gTTGATCATCAATAACCACCACAACATTATTGTTGAGAATTTTGGTGATTTGCATGTTCATAGCAAGGACCTTTTTTATAAACAAAAAAAACCGACTTCa                                                                                       >  2:458642/1‑100 (MQ=255)
      ttGATCATCAATAACCACCACAACATTATTGTTGAGAATTTTGGTGATTTGCATGTTCATAGCAAGGACCTTTTTTATAAACAAAAAAAACCGACTTCAc                                                                                      <  2:2961116/100‑1 (MQ=255)
       tGATCATCAATAACCACCACAACATTATTGTTGAGAATTTTGGTGATTTGCATGTTCATAGCAAGGACCTTTTTTATAAACAAAAAAACCCGACTTCAcc                                                                                     >  1:2737004/1‑100 (MQ=255)
       tGATCATCAATAACCACCACAACATTATTGTTGAGAATTTTGGTGATTTGCATGTTCATAGCAAGGACCTTTTTTATAAACAAAAAAAACCGACTTCAcc                                                                                     <  2:1872948/100‑1 (MQ=255)
       tGATCATCAATAACCACCACAACATTATTGTTGAGAATTTTGGTGATTTGCATGTTCATAGCAAGGACCTTTTTTATAAACAAAAAAAACCGACTTCAcc                                                                                     <  2:3352153/100‑1 (MQ=255)
        gATCATCAATAACCACCACAACATTATTGTTGAGAATTTTGGTGATTTGCATGTTCATAGCAAGGACCTTTTTTATAAACAAAAAAAACCGACTTCACCa                                                                                    >  2:1518773/1‑100 (MQ=255)
        gATCATCAATAACCACCACAACATTATTGTTGAGAATTTTGGTGATTTGCATGTTCATAGCAAGGACCTTTTTTATAAACAAAAAAAACCGACTTCACCa                                                                                    >  1:1430170/1‑100 (MQ=255)
        gATCATCAATAACCACCACAACATTATTGTTGAGAATTTTGGTGATTTGCATGTTCATAGCAAGAACCTTTTTTATAAACAAAAAAAACCGACTTCACCa                                                                                    >  1:2644537/1‑100 (MQ=255)
         atcatcAATAACCACCACAACGTTATTGTTGAGAATTTTGGTGATTTGCATGTTCATAGCAAGGACCTTTTTTATAAACAAAAAAAACCGACTTCACCAg                                                                                   >  1:236377/1‑100 (MQ=255)
         atcatcAATAACCACCACAACATTATTGTTGAGAATTTTGGTGATTTGCATGTTCATAGCAAGGACCTTTTTTATAAACAAAAAAACCCGACTTCACCAg                                                                                   <  2:723684/100‑1 (MQ=255)
         atcatcAATAACCACCACAACATTATTGTTGAGAATTTTGGTGATTTGCATGTTCATAGCAAGGACCTTTTTTATAAACAAAAAAAACCGACTTCACCAg                                                                                   <  2:2751615/100‑1 (MQ=255)
         atcatcAATAACCACCACAACATTATTGTTGAGAATTTTGGTGATTTGCATGTTCATAGCAAGGACCTTTTTTATAAACAAAAAAAACCGACATCACCAg                                                                                   <  2:1821142/100‑1 (MQ=255)
         atcataAATAACCACAACAACATTATTGTTGAGAATTTTGATAAATTGCATGTTCATAGCAAGGACCTTAAAAAAAAAAAAAAAAAACCGa                                                                                            <  1:191476/91‑1 (MQ=255)
          tcatcaATAACCACCACAACATTATTGTTGAGAATTTTGGTGATTTGCATGTTCATAGCAAGGACCTTTTTTATAAACAAAAAAACCCGACTTCACCAGt                                                                                  >  1:2626657/1‑100 (MQ=255)
          tcatcaATAACCACCACAACATTATTGTTGAGAATTTTGGTGATTTGCATGTTCATAGCAAGGACCTTTTTTATAAACAAAAAAAACCGACTTCACCAGt                                                                                  >  1:2919174/1‑100 (MQ=255)
          tcatcaATAACCACCACAACATTATTGTTGAGAATTTTGGTGATTTGCATGTTCATAGCAAGGACCTTTTTTATAAAAAAAAAAAACCGACTTCACCAGt                                                                                  <  1:1768837/100‑1 (MQ=255)
           catcaATAACCACCACAACATTATTGTTGAGAATTTTGGTGATTTGCATGTTCATAGCAAGGACCTTTTTTATAAACAAAAAAAACCGACTTCACCAGTa                                                                                 <  2:2106436/100‑1 (MQ=255)
           catcaATAACCACCACAACATTATTGTTGAGAATTTTGGTGATTTGCATGTTCATAGCAAGGACCTTTTTTATAAACAAAAAAAACCGACTTCACCAGTa                                                                                 >  1:2225355/1‑100 (MQ=255)
            atcaATAACCACCACAACATTATTGTTGAGAATTTTGGTGATTTGCATGTTCATAGCAAGGACCTTTTTTATAAACAAAAAAAACCGACTTCACCAGTAt                                                                                <  1:858081/100‑1 (MQ=255)
             tcaATAACCACCACAACATTATTGTTGAGAATTTTGGTGATTTGCATGTTCATAGCAAGGACCTTTTTTATAAACAAAAAAAACCGACTTCACCAGTAtt                                                                               >  2:1265627/1‑100 (MQ=255)
               aaTAACCACCACAACATTATTGTTGAGAATTTTGGTGATTTGCATGTTCATAGCAAGGACCTTTTTTATAAACAAAAAAACCCGACTTCACCAGTATtct                                                                             >  1:3031054/1‑100 (MQ=255)
               aaTAACCACCACAACATTATTGTTGAGAATTTTGGTGATTTGCATGTTCATAGCAAGGACCTTTTTTATAAACAAAAAAAACCGACTTCACCAGTATtct                                                                             >  2:3409500/1‑100 (MQ=255)
               aaTAACCACCACAACATTATTGTTGAGAATTTTGGTGATTTGCATGTTCATAGCAAGGACCTTTTTTATAAACAAAAAAAACCGACTTCACCAGTATtct                                                                             >  2:1212533/1‑100 (MQ=255)
                aTAACCACCACAACATTATTGTTGAGAATTTTGGTGATTTGCATGTTCATAGCAAGGACCTTTTTTATAAACAAAAAAAACCGACTTCACCAGTATtct                                                                             >  2:2397824/1‑99 (MQ=255)
                 tAACCACCACAACATTATTGTTGAGAATTTTGGTGATTTGCATGTTCATAGCAAGGACCTTTTTTATAAACAAAAAAAACCGACTTCACCAGTATTctct                                                                           <  2:1098935/100‑1 (MQ=255)
                  aaCCACCACAACATTATTGTTGAGAATTTTGGTGATTTGCATGTTCATAGCAAGGACCTTTTTTATAAACAAAAAAAACCGACTTCACCAGTATTCTCTg                                                                          <  1:3409448/100‑1 (MQ=255)
                      accacAACATTATTGTTGAGAATTTTGGTGATTTGCATGTTCATAGCTAGGACCTTTTTTATAAACAAAAAAAACCGACTTCACCAGTATTCTCTGGTTa                                                                      >  1:658003/1‑100 (MQ=255)
                      accacAACATTATTGTTGAGAATTTTGGTGATTTGCATGTTCATAGCAAGGACCTTTTTTATAAACAAAAAAAACCGACTTCACCAGTATTCTCTGGTTa                                                                      <  1:2348998/100‑1 (MQ=255)
                      accacAACATTATTGTTGAGAATTTTGGTGATTTGCATGTTCATAGCAAGGACCTTTTTTAAAAAAAAAAAAATCCGACTTCACCAGTATTCTCTGGTTa                                                                      <  1:3137240/100‑1 (MQ=255)
                       ccacAACATTATTGTTGAGAATTTTGGTGATTTGCATGTTCATAGCAAGGACCTTTTTTATAAACAAAAAAAACCGACTTCACCAGTATTCTCTGGTTAt                                                                     >  1:2234881/1‑100 (MQ=255)
                       ccacAACATTATTGTTGAGAATTTTGGTGATTTGCATGTTCATAGCAAGGACCTTTTTTATAAACAAAAAAAACCGACTTCACCAGTATTCTCTGGTTAt                                                                     >  2:1002605/1‑100 (MQ=255)
                       ccacAACATTATTGTTGAGAATTTTGGTGATTTGCATGTTCATAGCAAGGACCTTTTTTATAAACAAAAAAAACCGACTTCACCAGTATTCTCTGGTTAt                                                                     >  2:1778662/1‑100 (MQ=255)
                       ccacAACATTATTGTTGAGAATTTTGGTGATTTGCATGTTCATAGCAAGGACCTTTTTTATAAACAAAAAAAACCGACTTCACCAGTATTCTCTGGTTAt                                                                     <  1:1435816/100‑1 (MQ=255)
                       ccaaAACATTATTGTTGAGAATTTTGGTGATTTGCATGTTCATAGCAAGGACCTTTTTTATAAACAAAAAAAACCGACTTCACCAGTATTCTCTGGTTAt                                                                     <  2:675624/100‑1 (MQ=255)
                           aacaTTATTGTTGAGAATTTTGGTGATTTGCATGTTCATAGCAAGGACCTTTTTTATAAACAAAAAAAACCGACTTCACCAGTATTCTCTGGTTATGTCa                                                                 <  1:222510/100‑1 (MQ=255)
                            acaTTATTGTTGAGAATTTTGGTGATTTGCATGTTCATAGCAAGGACCTTTTTTATAAACAAAAAAACCCGACTTCACCAGTATTCTCTGGTTATGTCAg                                                                >  2:3052450/1‑100 (MQ=255)
                                tattGTTGAGAATTTTGGTGATTTGCATGTTCATAGCAAGGACCTTTTTTATAAACAAAAAAAACCGACTTCACCAGTATTCTCTGGTTATTTCAGGttt                                                            >  2:2760670/1‑100 (MQ=255)
                                  ttgttgAGAATTTTGCTGATGTGTATGTTCATAGCAATTACCTTTTAAAAAAAAAAAAAAAACCGACTTCACCAGTATTCTCTGGTTATGTCAGGTTTTg                                                          <  1:827608/100‑1 (MQ=255)
                                    gttgAGAATTTTGGTGATTTGCATGTTCATAGCAAGGACCTTTTTTATAAACAAAAAAACCCGACTTCACCAGTATTCTCTGGTTATGTCAGGTTTTGcc                                                        <  1:549835/100‑1 (MQ=255)
                                       gagaATTTTGGTGATTTGCATGTTCATAGCAAGGACCTTTTTTATAAACAAAAAAACCCGACTTCACCAGTATTCTCTGGTTATGTCAGGTTTTGCCTGc                                                     >  1:1296288/1‑100 (MQ=255)
                                       gagaATTTTGGTGATTTGCATGTTCATAGCAAGGACCTTTTTTATAAACAAAAAAACACGACTTCACCAGTATTCTCTGGTGATGTAAGGTTTTGACTGc                                                     >  1:1426852/1‑100 (MQ=255)
                                       gagaATTTTGGTGATTTGCATGTTCATAGCAAGGACCTTTTTTATAAACAAAAAAAACCGACTTCACCAGTATTCTCTGGTTATGTCAGGTTTTGCCTGc                                                     >  1:2452664/1‑100 (MQ=255)
                                       gagaATTTTGGTGATTTGCATGTTCATAGCAAGGACCTTTTTTATAAACAAAAAAAACCGACTTCACCAGTATTCTCTAGTTATGTCAGGTTTTCCCtac                                                     >  1:1458877/1‑98 (MQ=255)
                                        agaATTTTGGTGATTTGCATGTTCATAGCAAGGACCTTTTTTATAAACAAAAAAAACCGACTTCACCAGTATTCTCTGGTTATGTCAGGTTTTGCCTGCg                                                    >  2:1381539/1‑100 (MQ=255)
                                         aaaTTTTAGTGATATGCATTTTGATAGAAAGTAACTGTATAATAAAAAAAAAAAACCGACTTCACCAGTATTCTCTGGTTATGTCAGGTTTTGCCTGCg                                                    <  1:3194233/98‑1 (MQ=255)
                                                gttaaTTTGCATGTTCATAGCAAGGACCTTTTTTATAAACAAAAAAAACCGACTTCACCAGTATTCTCTGGTTATGTCAGGTTTTGCCTGCGAATGCAGt                                            <  1:464698/96‑1 (MQ=255)
                                                ggTGATTTGCATGTTCATAGCAAGGACCTTTTTTATAAACAAAAAAACCCGACTTCACCAGTATTCTCTGGTTATGTCAGGTTTTGCCTGCGAATGCAGt                                            <  2:2268317/100‑1 (MQ=255)
                                                ggTGATTTGCATGTTCATAGCAAGGACCTTTTTTATAAACAAAAAAAACCGACTTCACCAGTATTCTCTGGTTATGTCAGGTTTTGCCTGCGAATGCAGt                                            <  2:2599761/100‑1 (MQ=255)
                                                    aTTTGCATGTTCATAGCAAGGACCTTTTTTATAAACAAAAAAAACCGACTTCACCAGTATTCTCTGGTTATGTCAGGTTTTGCCTGCGAATGCAGTAACa                                        <  2:2241875/100‑1 (MQ=255)
                                                        gCATGTTCATAGCAAGGACCTTTTTTATAAACAAAAAAAACCGACTTCACCAGTATTCTCTGGTTATGTCAGGTTTTGCCTGCGAATGCAGTAACAATcc                                    <  1:3346836/100‑1 (MQ=255)
                                                         cATGATCATAGCAAGGAACTTTTTTATAAACAAAAAAAACCGACTGCAGCAATATTTTCTAATTAAATCAGAGTATTCCTGCAAATTCATTAATAATGCa                                   >  1:1156390/1‑100 (MQ=255)
                                                            gTTCGTAGCAAGGACCTTTTTGATAAACAAAAAAAACCGACTTCACCAGGATTCTCTGGTTATGTCAGGTTTGGCCTGCGGATGGAGTGAGAGTCCAGTc                                >  1:2076239/1‑100 (MQ=255)
                                                              tCATAGCAAGGACCTTTTTTATAAACAAAAAAAACCGACTTCACCAGTATTCTCTGGTTATGTCAGGTTTTGTCTGCGAAAGCAGTAACAATTCAGTAAt                              >  1:2944568/1‑100 (MQ=255)
                                                               cATAGCAAGGACCTTTTTTATAAACAAAAAAAACCGACTTCACCAGTATTCTCTGGTTATGTCAGGTTTTGCCTGCGAATGCAGTAACAATCCAGACAtt                             <  1:547453/100‑1 (MQ=255)
                                                               cATAGCAAGGACCTTTTTAATAAACAAAAAAAACCGACTTCACCAGTATTCTCTGGTTATGTCAGGTTTTGCCTGCGAATGCAGTAACAATCCAGTCAtt                             <  1:1660439/100‑1 (MQ=255)
                                                                aTAGCAAGGACCTTTTTTATAAACAAAAAAACCCGACTTCACCAGTATTCTCTGGTTATGTCAGGTTTTGCCTGCGAATGCAGTAACAATCCAGTCAttt                            >  1:996303/1‑100 (MQ=255)
                                                                aTAGCAAGGACCTTTTTTATAAACAAAAAAACCCGACTTCACCAGTATTCTCTGGTTATGTCAGGTTTTGCCTGCGAATGCAGTAACAATCCAGTCAttt                            <  2:2737047/100‑1 (MQ=255)
                                                                 tAGCAAGGACCTTTTTTATAAACAAAAAAAACCGACTTCACCAGTATTCTCTGGTTATGTCAGGTTTTGCCTGCGAATGCAGTAACAATCCAGTCAttt                            >  2:438271/1‑99 (MQ=255)
                                                                 tAGCAAGGACCTTTTTTATAAACAAAAAAAACCGACTTCACCAGTATTCTCTGGTTATGTCAGGTTTTGCCTGCGAATGCAGTAACAATCCAGTCAttt                            >  2:894289/1‑99 (MQ=255)
                                                                  aGCAAGGACCTTTTTTATAAACAAAAAAACCCGACTTCACCAGTATTCTCTGGTTATGTCAGGTTTTGCCTGCGAATGCAGTAACAATCCAGTCATttat                          <  2:1705633/100‑1 (MQ=255)
                                                                   gCAAGGACCTTTTTTATAAACAAAAAAACCCGACTTCACCAGTATTCTCTGGTTATGTCAGGTTTTGCCTGCGAATGCAGTAACAATCCAGTCATttatt                         >  2:1721172/1‑100 (MQ=255)
                                                                   gCAAGGACCTTTTTTATAAACAAAAAAAACCGACTTCACCAGTATTCTCTGGTTATGTCAGGTTTTGCCTGCGAATGCAGTAACAATCCAGTCATttatt                         >  1:2200396/1‑100 (MQ=255)
                                                                   gCAAGGACCTTTTTTATAAACAAAAAAAACCGACTTCACCAGTATTCTCTGGTTATGTCAGGTTTTGCCTGCGAATGCAGTAACAATCCAGTCATttatt                         >  1:3289433/1‑100 (MQ=255)
                                                                   gCAAGGACCTTTTTTATAAACAAAAAAAACCGACTTCACCAGTATTCTCTGGTTATGTCAGGTTTTGCCTGCGAATGCAGTAACAATCCAGTCATttatt                         >  2:3036752/1‑100 (MQ=255)
                                                                   gCAAGGACCTTTTTTATAAACAAAAAAAACCGACTTCACCAGTATTCTCTGGTTATGTCAGGTTTTGCCTGCGAATGCAGTAACAATCCAGTCATttatt                         >  1:3087877/1‑100 (MQ=255)
                                                                    cAAGGACCTTTTTTATACACAAAAAAAACCGACTTAACCAGTATTCTGAGGTTATGTCAGGTTTTGCCAGCGAATGCAGTAACAATCCAGGCATgtatta                        >  1:1485735/1‑100 (MQ=255)
                                                                    cAAGGACCTTTTTTATAAACAAAAAAAACCGACTTCACCAGTATTCTCTGGTTATGTCAGGTTTTGCCTGCGAATGCAGTAACAATCCAGTCATttatta                        <  1:1968188/100‑1 (MQ=255)
                                                                      aGGACCTTTTTTATAAACAAAAAAACCCGACTTCACCAGTATTCTCTGGTTATGTCAGGTTTTGCCTGCGAATGCAGTAACAATCCAGTCATTTATTAAt                      <  2:2420213/100‑1 (MQ=255)
                                                                      aGGACCTTTTTTATAAACAAAAAAACCCGACTTCACCAGTATTCTCTGGTTATGTCAGGTTTTGCCTGCGAATGCAGTAACAATCCAGTCATTTATTAAt                      <  2:2413289/100‑1 (MQ=255)
                                                                       ggACCTTTTTTATAAACAAAAAAAACCGACTTCACCAGTATTCTCTGGTTATGTCAGGTTTTGCCTGCGAATGCAGTAACAATCCAGTCATTTATTAATg                     <  2:3360995/100‑1 (MQ=255)
                                                                         aCCTTTTTTATAAACAAAAAAAACCGACTTCACCAGTATTCTCTGGTTATGTCAGGTTTTGCCTGCGAATGCAGTAACAATCCAGTCATTTATTAATGGt                   <  1:2628064/100‑1 (MQ=255)
                                                                          ccTTTTTTATAAACAAAAAAAACCGACTTCACCAGTATTCTCTGGTTATGTCAGGTTTTGCCTGCGAATGCAGTAACAATCCAGTCATTTATTAATGGtt                  <  2:2824820/100‑1 (MQ=255)
                                                                           cTTTTTTATAAACAAAAAAAACCGACTTCACCAGTATTCTTTGGTTATGTCAGGTTTTGCCTGCGAATGCAGTAACAACCCAGTCATTTATTAATGGttt                 >  1:3203529/1‑100 (MQ=255)
                                                                           cTTTTTTATAAACAAAAAAAACCGACTTCACCAGTATTCTCTGGTTATGTCAGGTTTTGCCTGCGAATGCAGTAACAATCCAGTCATTTATTAATGGttt                 >  1:954867/1‑100 (MQ=255)
                                                                           cTTTTTTATAAACAAAAAAAACCGACTTCACCAGTATTCTCTGGTTATGTCAGGTTTTGCCTGCGAATGCAGTAACAATCCAGTCATTTATTAATGGttt                 >  2:2970119/1‑100 (MQ=255)
                                                                           cTTTTTTATAAACAAAAAAAACCGACTTCACCAGTATTCTCTGGTTATGTCAGGTTTTGCCTGCGAATGCAGTAACAATCCAGTCATTTATTAATGGttt                 >  2:813797/1‑100 (MQ=255)
                                                                           cTTTTTTATAAACAAAAAAAACCGACTTCACCAGTATTCTCTGGTTATGTCAGGTTTTGCCTGCGAATGCAGTAACAATCCAGTCATTTATTAATGGttt                 >  1:1137823/1‑100 (MQ=255)
                                                                           cTTTTTTATAAACAAAAAAAACCGACTTCACCAGTATTCTCTGGTTATGTCAGGTTTTGCCTGCGAATGCAGTAACAATCCAGTCATTTATTAATGGttt                 >  1:2701854/1‑100 (MQ=255)
                                                                           cTTTTTTATAAACAAAAAAAACCGACTTCACCAGTATTCTCTGGTTATGTCAGGTTTTGCCTGCGAATGCAGTAACAATCCAGTCATTTATTAATGGttt                 >  2:1568841/1‑100 (MQ=255)
                                                                            ttttttATAAACAAAAAAAACCGACTTCACCAGTATTCTCTGGTTATGTCAGGTTTTGCCTGCGAATGCAGTAACAATCCAGTCATTTATTAATGGtttt                <  2:309139/100‑1 (MQ=255)
                                                                             tttttATAAACAAAAAAAACCGACTTCACCAGTATTCTCTGGTTATGTCAGGTTTTGCCTGCGAATGCAGTAACAATCCAGTCATTTATTAATGGttttt               <  1:421160/100‑1 (MQ=255)
                                                                              ttttATAAACAAAAAAAACCGACTTCACCAGTATTCTCTGGTTATGTCAGGTTTTGCCTGCGAATGCAGTAACAATCCAGTCATTTATTAATGGTTTTta              <  1:2921470/100‑1 (MQ=255)
                                                                              ttttATAAACAAAAAAAACCGACTTCACCAGTATTCTCTGGTTATGTCAGGTTTTGCCTGCGAATGCAGTAACAATCCAGTCATTTATTAATGGTTTTta              <  2:516121/100‑1 (MQ=255)
                                                                                ttATAAACAAAAAAACCCGACTTCACCAGTATTCTCTGGTTATGTCAGGTTTTGCCTGCGAATGCAGTAACAATCCAATCATTTATTAATGGTTTTtata            >  1:1685436/1‑100 (MQ=255)
                                                                                 tataAACAAAAAAAACCGACTTCACCAGTATTCGCTGATTATGTCAGGTTTTGCCTGCGAATGCAGTCACACTCGAGTCATTTATTAATGGTTTTtata            >  1:1719755/1‑99 (MQ=255)
                                                                                  ataAACAAAAAAAACCGACTTCACCAGTATTCTCTGGTTATGTCAGGTTTTGCCTGCGAATGCAGTAACAATCCAGTCATTTATTAATGGTTTTTATAAc          >  2:1940918/1‑100 (MQ=255)
                                                                                    aaaCAAAAAAAACCGACTTCACCAGTATTCTCTGGTTATGTCAGGTTTTGCCTGCGAATGCAGTAACAATCCAGTCATTTATTAATGGTTTTTATAACGa        <  1:898015/100‑1 (MQ=255)
                                                                                     aaCAAAAAAAACCGACTTCACCAGTATTCTCTGGTTATGTCAGGTTTTGCCTGCGAATGCAGTAACAATCCAGTCATTTATTAATGGTTTTTATAACGaa       >  1:2367378/1‑100 (MQ=255)
                                                                                      agaaaaaaaaCTGACTTCACAAGTATGGGATGGTTATGTCAGGTTTTGCCTGCGAATGCAGTAACAAGCCAGTCATTTATTAATGGTTTTTATAACGAAc      <  1:696805/98‑1 (MQ=255)
                                                                                      aCAAAAAAAACCGACTTCACCAGTATTCTCTGGTTATGTCAGGTTTTGCCTGCGAATGCAGTAACAATCCAGTCATTTATTAATGGTTTTTATAACGAAc      >  2:2607185/1‑100 (MQ=255)
                                                                                      aCAAAAAAAACCGACTTCACCAGTATTCTCTGGTTATGTCAGGTTTTGCCTGCGAATGCAGTAACAATCCAGTCATTTATTAATGGTTTTTATAACGAAc      <  2:1163946/100‑1 (MQ=255)
                                                                                       cAAAAAAAACCGACTTCACCAGTATTCTCTGGTTATGTCAGGTTTTGCCTGCGAATGCAGTAACAATCCAGTCATTTATTAATGGTTTTTATAACGAACa     >  2:908358/1‑100 (MQ=255)
                                                                                         aaaaaaACCGACTTCACCAGTATTCTCTGGTTATGTCAGGTTTTGCCTGCGAATGCAGTAACAATCCAGTCATTTATTAATGGTTTTTATAACGAACATc   >  1:3349118/1‑100 (MQ=255)
                                                                                         aaaaaaACAGACTTCACCAGTATTCTCTGGTTATGTAAGGTTTTGCCTGCGAATGCAGTAACAATCCAGGCATTTATTAATGGTTTTTATAACGAAAAt    >  1:419794/1‑99 (MQ=255)
                                                                                          aaaaaaCCGACTTCACCAGTATTCTCTGGTTATGTCAGGTTTTGCCTGCGAATGCAGTAACAATCCAGTCATTTATTAATGGTTTTTATAACGAACATcc  >  2:1945229/1‑100 (MQ=255)
                                                                                          aaaaaaCCGACTTCACCAGTATTCTCTGGTTATGTCAGGTTTTGCCTGCGAATGCAGTAACAATCCAGTCATTTATTAATGGTTTTTATAACGAACATcc  >  2:1955934/1‑100 (MQ=255)
                                                                                          aaaaaCCCGACTTCACCAGTATTCTCTGGTTATGTCAGGTTTTGCCTGCGAATGCAGTAACAATCCAGTCATTTATTAATGGTTTTTATAACGAACATcc  <  2:150676/100‑1 (MQ=255)
                                                                                               |                                                                                              
CCGCTGTTGATCATCAATAACCACCACAACATTATTGTTGAGAATTTTGGTGATTTGCATGTTCATAGCAAGGACCTTTTTTATAAACAAAAAAACCCGACTTCACCAGTATTCTCTGGTTATGTCAGGTTTTGCCTGCGAATGCAGTAACAATCCAGTCATTTATTAATGGTTTTTATAACGAACATCC  >  NC_000913/3906502‑3906691

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: