Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 70,371 T→C 100% intergenic (‑323/‑16) araB ← / → araC ribulokinase/DNA‑binding transcriptional dual regulator AraC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091370,3710TC96.8% 6796.9 / 121.6 1883intergenic (‑323/‑16)araB/araCribulokinase/DNA‑binding transcriptional dual regulator AraC
Reads supporting (aligned to +/- strand):  ref base T (35/26);  new base C (947/871);  total (985/898)
Fisher's exact test for biased strand distribution p-value = 4.37e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.66e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 1929 total aligned reads displayed.

AGCGGGGTTACCGGTTGGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGTGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAACGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGC  >  NC_000913/70273‑70469
                                                                                                  |                                                                                                  
  cGGGGTTACCGGTTTGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAAtggcgg                                                                                                  <  2:1192619/99‑1 (MQ=255)
    gggTTACCGGTTTGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGa                                                                                                >  1:1655884/1‑99 (MQ=255)
    gggTTACCGGTTTGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGa                                                                                                <  2:779877/99‑1 (MQ=255)
     ggTTACCGGTTTGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAg                                                                                               <  1:630551/99‑1 (MQ=255)
     ggTTACCGGTTTGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAg                                                                                               <  2:374716/99‑1 (MQ=255)
     ggTTACCGGTTTGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAg                                                                                               >  2:568485/1‑99 (MQ=255)
      gTTACCGGTTTGGTTAGAGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGt                                                                                              >  1:1815154/1‑99 (MQ=255)
       ttACCGGTTTGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTa                                                                                             <  1:774672/99‑1 (MQ=255)
       ttACCGGTTTGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTa                                                                                             <  1:501030/99‑1 (MQ=255)
       ttACCGGTTGGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGTGGGAGTa                                                                                             <  1:429379/99‑1 (MQ=255)
        tACCGGTTTGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTAt                                                                                            <  1:123531/99‑1 (MQ=255)
           cGGTTTGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGaa                                                                                         >  2:1913007/1‑99 (MQ=255)
             gTTTGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGaaaa                                                                                       >  2:239870/1‑99 (MQ=255)
              tttggTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAg                                                                                      >  2:1762656/4‑99 (MQ=255)
              tttggTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAg                                                                                      >  2:1885318/4‑99 (MQ=255)
               ttggTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGt                                                                                     <  1:343238/97‑1 (MQ=255)
                tggTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTa                                                                                    <  1:673034/98‑1 (MQ=255)
                tggTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTa                                                                                    <  1:531462/98‑1 (MQ=255)
                   ttAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATgg                                                                                 <  1:1628533/99‑1 (MQ=255)
                   ttAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATgg                                                                                 >  1:1295864/1‑99 (MQ=255)
                   ttAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATgg                                                                                 >  2:329176/1‑99 (MQ=255)
                     aGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCt                                                                               <  1:1696202/99‑1 (MQ=255)
                      gCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTg                                                                              >  1:752195/1‑99 (MQ=255)
                        gagaAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGaa                                                                            <  1:1595819/99‑1 (MQ=255)
                         agaagaGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAg                                                                           <  1:8906/99‑1 (MQ=255)
                           aagaGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAgcg                                                                         <  1:248058/99‑1 (MQ=255)
                            agagCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAgcgc                                                                        >  2:1698751/1‑99 (MQ=255)
                             gagCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCa                                                                       >  2:1027808/1‑99 (MQ=255)
                              agCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCaa                                                                      >  2:1666281/1‑99 (MQ=255)
                               gCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCaaa                                                                     >  2:1940125/1‑99 (MQ=255)
                               gCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCaaa                                                                     <  2:669882/99‑1 (MQ=255)
                                ccAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCaaaa                                                                    <  1:1406061/99‑1 (MQ=255)
                                ccAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCaaaa                                                                    <  2:632890/99‑1 (MQ=255)
                                 cAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAAt                                                                   >  2:1155173/1‑99 (MQ=255)
                                 cAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAAt                                                                   <  1:126215/99‑1 (MQ=255)
                                      aaaGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATcc                                                              >  1:1513006/1‑99 (MQ=255)
                                        aGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATcccc                                                            <  1:1445714/99‑1 (MQ=255)
                                         gACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCct                                                           <  1:1730586/99‑1 (MQ=255)
                                         gACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGAGCAAAATGATCCCct                                                           >  1:1035566/1‑99 (MQ=255)
                                          aCGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCctg                                                          <  2:1471829/99‑1 (MQ=255)
                                           cGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCctgc                                                         <  1:1195379/99‑1 (MQ=255)
                                           cGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCctgc                                                         <  1:707496/99‑1 (MQ=255)
                                              aGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCtgctgc                                                      <  2:27917/99‑1 (MQ=255)
                                               gTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGcc                                                     >  1:1379046/1‑99 (MQ=255)
                                                tGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCg                                                    >  1:818773/1‑99 (MQ=255)
                                                 gACGGCAATGTCTGATGCAATATGGGCAATTGGTTTCTTCTCTGAATGGTGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCgg                                                   >  2:1540264/1‑99 (MQ=255)
                                                 gACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCgg                                                   >  2:1266149/1‑99 (MQ=255)
                                                  aCGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCggg                                                  >  1:276239/1‑99 (MQ=255)
                                                   cGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGa                                                 >  2:1412304/1‑99 (MQ=255)
                                                     gCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTTAATGGCGGGAGTATGAAAAGTATGGATGAAGCGCAAAATGATCCCCTGCTGCCGGGATa                                               >  1:927005/1‑99 (MQ=255)
                                                     gCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATa                                               >  1:461578/1‑99 (MQ=255)
                                                     gCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATa                                               >  2:1327351/1‑99 (MQ=255)
                                                     gCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATa                                               >  1:1056084/1‑99 (MQ=255)
                                                      cAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATAc                                              <  1:1559528/99‑1 (MQ=255)
                                                      cAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATAc                                              >  2:1440309/1‑99 (MQ=255)
                                                       aaTGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACt                                             >  2:970156/1‑99 (MQ=255)
                                                        aTGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTc                                            >  2:1300252/1‑99 (MQ=255)
                                                        aTGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTc                                            >  1:403155/1‑99 (MQ=255)
                                                           tCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGtt                                         <  1:1342814/99‑1 (MQ=255)
                                                           tCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGtt                                         >  2:482533/1‑99 (MQ=255)
                                                            cTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGttt                                        >  2:1052705/1‑99 (MQ=255)
                                                              gATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTaa                                      <  2:1734468/99‑1 (MQ=255)
                                                              gATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTaa                                      >  1:1093604/1‑99 (MQ=255)
                                                              gATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTaa                                      <  1:307179/99‑1 (MQ=255)
                                                              gATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTaa                                      >  2:842000/1‑99 (MQ=255)
                                                               aTGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTaat                                     >  2:1500799/1‑98 (MQ=255)
                                                               aTGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTaat                                     >  2:1794324/1‑98 (MQ=255)
                                                                   aaTATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGccc                                 <  2:815611/99‑1 (MQ=255)
                                                                    atatGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCa                                <  2:455772/99‑1 (MQ=255)
                                                                    atatGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCa                                <  1:1135354/99‑1 (MQ=255)
                                                                    atatGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCa                                <  1:594437/99‑1 (MQ=255)
                                                                     tatGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCAt                               <  2:1598702/99‑1 (MQ=255)
                                                                      atGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATc                              <  1:853364/99‑1 (MQ=255)
                                                                       tGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCt                             >  2:925406/1‑99 (MQ=255)
                                                                        ggACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCtg                            >  2:529466/1‑99 (MQ=255)
                                                                         gACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCtgg                           <  2:702303/99‑1 (MQ=255)
                                                                         gACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCtgg                           >  2:1374694/1‑99 (MQ=255)
                                                                           cAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCtggtg                         <  1:1860954/99‑1 (MQ=255)
                                                                           cAATTGGTTTCTTCTCTGAATGGCGGGAGTATCAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCtggtg                         <  1:1768055/99‑1 (MQ=255)
                                                                            aaTTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCtggtgg                        >  2:1354936/1‑99 (MQ=255)
                                                                            aaTTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCtggtgg                        <  1:1479044/99‑1 (MQ=255)
                                                                            aaTTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCtggtgg                        >  2:1227086/1‑99 (MQ=255)
                                                                               tGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCgg                     <  1:375761/99‑1 (MQ=255)
                                                                                 gTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGt                   <  1:1886254/99‑1 (MQ=255)
                                                                                     cttctCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTaa               <  2:736733/99‑1 (MQ=255)
                                                                                     cttctCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTaa               <  1:1167958/99‑1 (MQ=255)
                                                                                     cttctCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTaa               <  2:604558/99‑1 (MQ=255)
                                                                                     cttctCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTaa               <  1:1004429/99‑1 (MQ=255)
                                                                                     cttctCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTaa               >  2:883801/1‑99 (MQ=255)
                                                                                     cttctCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTaa               <  1:985246/99‑1 (MQ=255)
                                                                                      ttctCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAAc              >  2:419682/1‑99 (MQ=255)
                                                                                       tctcTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACg             >  1:221253/1‑99 (MQ=255)
                                                                                       tctcTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACg             <  2:1102685/99‑1 (MQ=255)
                                                                                        ctctGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGc            <  2:213616/99‑1 (MQ=255)
                                                                                        ctctGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGc            <  1:1906751/99‑1 (MQ=255)
                                                                                         tctGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGcc           >  1:561371/1‑99 (MQ=255)
                                                                                            gAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGAt        >  2:1633964/1‑99 (MQ=255)
                                                                                             aaTGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGAtt       >  2:1836120/1‑99 (MQ=255)
                                                                                             aaTGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGAtt       <  1:945397/99‑1 (MQ=255)
                                                                                               tggcggGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGATTGa     >  1:1937629/1‑99 (MQ=255)
                                                                                                  |                                                                                                  
AGCGGGGTTACCGGTTGGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGTGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAACGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGC  >  NC_000913/70273‑70469

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