Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,474,971 A→G H37R (CAC→CGC)  tabA → biofilm modulator regulated by toxins; DUF386 family protein, cupin superfamily protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,474,9710AG98.2% 762.3 / 0.9 219H37R (CAC→CGC) tabAbiofilm modulator regulated by toxins; DUF386 family protein, cupin superfamily protein
Reads supporting (aligned to +/- strand):  ref base A (3/1);  new base G (96/118);  total (99/120)
Fisher's exact test for biased strand distribution p-value = 3.32e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.63e-01

Only 100 of 220 total aligned reads displayed.

AATATTCATAATCTTCAGCCGTGGCTACCCCAGGAGTTACGCCAGGCGATTGAGCATATCAAAGCACACGTTACGGCAGAAACGCCAAAGGGCAAGCACGATATCGAAGGCAATCGACTGTTTTATCTTATCTCGGAAGATATGACCGAGCCGTACGAAGCTCGCCGTGCGGAGTACCATGCCCGCTATCTTGAC  >  NC_000913/4474874‑4475068
                                                                                                 |                                                                                                 
  tatTCATAATCTTCAGCCGTGGCTACCCCAGGAGTTACGCCAGGCGATTGAGCATATCAAAGCACACGTTACGGCAGAAACGCCAAAGGGCAAGCGCGat                                                                                               <  1:3032930/100‑1 (MQ=255)
  tatTCATAATCTTCAGCCGTGGCTACCCCAGGAGTTACGCCAGGCGATTGAGCATATCAAAGCACACGTTACGGCAGAAACGCCAAAGGGCAAGCGCGat                                                                                               >  2:2859872/1‑100 (MQ=255)
    ttCATAATCTTCAGCCGTGGCTACCCCAGGAGTTACGCCAGGCGATTGAGCATATCAAAGCACACGTTACGGCAGAAACGCCAAAGGGCAAGCGCGatat                                                                                             <  2:2361043/100‑1 (MQ=255)
     tCATAATCTTCAGCCGTGGCAACCCCAAGAGATAAGCCAGGCGATTGAGCATATCAAAGCACACGTTACGGCAGAAACGCCAAAGGGAAAGCGCGATATc                                                                                            <  2:3309808/100‑1 (MQ=255)
      cATAATCTTCAGCCGTGGCTACCCCAGGAGTTACGCCAGGCGATTGAGCATATCAAAGCACACGTTACGGCAGAAACGCCAAAGGGCAAGCGCGATATCg                                                                                           <  2:716687/100‑1 (MQ=255)
      cATAATCTTCAGCCGTGGCTACCCCAGGAGTTACGCCAGGCGATTGAGCATATCAAAGCACACGTTACGGCAGAAACGCCAAAGGGCAAGCGCGATATCg                                                                                           >  2:706385/1‑100 (MQ=255)
       aTAATCTTCAGCCGTGGCTACCCCAGGAGTTACGCCAGGCGATTGAGCATATCAAAGCACACGTTACGGCAGAAACGCCAAAGGGCAAGCGCGATATCGa                                                                                          >  2:1917446/1‑100 (MQ=255)
       aTAATCTTCAGCCGTGGCTACCCCAGGAGTTACGCCAGGCGATTGAGCATATCAAAGCACACGTTACGGCAGAAACGCCAAAGGGCAAGCGCGATATCGa                                                                                          >  2:2617829/1‑100 (MQ=255)
       aTAATCTTCAGCCGTGGCTACCCCAGGAGTTACGCCAGGCGATTGAGCATATCAAAGCACACGTTACGGCAGAAACGCCAAAGGGCAAGCGCGATATCGa                                                                                          >  2:3087295/1‑100 (MQ=255)
        tAATCTTCAGCCGTGGCTACCCCAGGAGTTACGCCAGGCGATTGAGCATATCAAAGCACACGTTACGGCAGAAACGCCAAAGGGCAAGCGCGATATCGaa                                                                                         >  1:1527934/1‑100 (MQ=255)
         aaTCTTCAGCCGTGGCTACCCCAGGAGTTACGCCAGGCGATTGAGCATATCAAAGCACACGTTACGGCAGAAACGCCAAAGGGCAAGCGCGATATCGAAg                                                                                        <  1:2994615/100‑1 (MQ=255)
         aaTCTTCAGCCGTGGCTACCCCAGGAGTTACGCCAGGCGATTGAGCATATCAAAGCACACGTTACGGCAGAAACGCCAAAGGGCAAGCGCGATATCGAAg                                                                                        <  1:4815210/100‑1 (MQ=255)
          aTCTTCAGCCGTGGCTACCCCAGGAGTTACGCCAGGCGATTGAGCATATCAAAGCACACGTTACGGCAGAAACGCCAAAGGGCAAGCGCGATATCGAAgg                                                                                       <  1:3787851/100‑1 (MQ=255)
           tCTTCAGCCGTGGCTCCCCCAGGAGTTACGCCAGGCGATTGAGCATATCAAAGCACACGTTACGGCAGAAACGCCAAAGGGCAAGCGCGATATCGAAGGc                                                                                      <  1:2152266/100‑1 (MQ=255)
           tCTTCAGCCGTGGCTACCCCAGGAGTTACGCCAGGCGATTGAGCATATCAAAGCACACGTTACGGCAGAAACGCCAAAGGGCAAGCGCGATATCGAAGGc                                                                                      <  2:583309/100‑1 (MQ=255)
           tCTTCAGCCGTGGCTACCCCAGGAGTTACGCCAGGCGATTGAGCATATCAAAGCACACGTTACGGCAGAAACGCCAAAGGGCAAGCGCGATATCGAAGGc                                                                                      <  2:2276331/100‑1 (MQ=255)
             ttCAGCCGTGGCTACCCCAGGAGTTACGCCAGGCGATTGAGCATATCAAAGCACACGTTACGGCAGAAACGCCAAAGGGCAAGCGCGATATCGAAGGCaa                                                                                    >  1:2537866/1‑100 (MQ=255)
              tCAGCCGTGGCTACCCCAGGAGTTACGCCAGGCGATTGAGCATATCAAAGCACACGTTACGGCAGAAACGCCAAAGGGCAAGCGCGATATCGAAGGCAAt                                                                                   <  1:3116218/100‑1 (MQ=255)
                agtcttGGCTTCCCCAGGATTTACGTCAGGGGTTTTAGCAGATTATAGCAGACGTTACGGCAGAGACGCCAAAGGGCAAGCGCGATATCGAAGGCAATCg                                                                                 <  2:3465402/95‑1 (MQ=255)
                aGCCGTGGCTACCCCAGGAGTTACGCCAGGCGATTGAGCATATCAAAGCACACGTTACGGCAGAAACGCCAAAGGGCAAGCGCGATATCGAAGGCAATCg                                                                                 >  1:2705330/1‑100 (MQ=255)
                aGCCGTGGCTACCCCAGGAGTTACGCCAGGCGATTGAGCATATCAAAGCACACGTTACGGCAGAAACGCCAAAGGGCAAGCGCGATATCGAAGGCAATCg                                                                                 >  1:1712265/1‑100 (MQ=255)
                 gCCGTGGCTACCCCACGAGTTACGCCAGGCAATTGAGCATATCAATGCACACGTTACTGCAGAAACGCCAAATGGCAAGCGCGATATCGAAGGCAATCGa                                                                                >  1:4841143/1‑100 (MQ=255)
                    gTGGCTACCCCAGGAGTTACGCCAGGCGATTGAGCATATCAAAGCACACGTTACGGCAGAAACGCCAAAGGGCAAGCGCGATATCGAAGGCAATCGACTg                                                                             >  2:1162684/1‑100 (MQ=255)
                     tGGCTACCCCAGGAGTTACGCCAGGCGATTGAGCATATCAAAGCACACGTTACGGCAGAAACGCCAAAGGGCAAGCGCGATATCGAAGGCAATCGACTGt                                                                            <  1:4031736/100‑1 (MQ=255)
                       gCTACCCCAGGAGTTACGCCAGGCGATTGAGCATATCAAAGCACACGTTACGGCAGAAACGCCAAAGGGCAAGCGCGATATCGAAGGCAATCGACTGttt                                                                          <  1:1120880/100‑1 (MQ=255)
                       gCTACCCCAGGAGTTACGCCAGGCGATTGAGCATATCAAAGCACACGTTACGGCAGAAACGCCAAAGGGCAAGCGCGATATCGAAGGCAATCGACTGttt                                                                          <  1:3102651/100‑1 (MQ=255)
                          aCCCCAGGAGTTACGCCAGGCGATTGAGCATATCAAAGCACACGTTACGGCAGAAACGCCAAAGGGCAAGCGCGATATCGAAGGCAATCGACTGTTttat                                                                       <  2:2878021/100‑1 (MQ=255)
                          aCCCCAGGAGTTACGCCAGGCGATTGAGCATATCAAAGCACACGTTACGGCAGAAACGCCAAAGGGCAAGCGCGATATCGAAGGCAATCGACTGTTttat                                                                       <  2:2087999/100‑1 (MQ=255)
                          aCCCCAGGAGTTACGCCAGGCGATTGAGCATATCAAAGCACACGTTACGGCAGAAACGCCAAAGGGCAAGCGCGATATCGAAGGCAATCGACTGTTttat                                                                       >  1:1369368/1‑100 (MQ=255)
                           gcccAGGAGTTACGCCAGGCGATTGAGCATATCAAAGCACACGTTACGGTAGAAACGCCAAAGGGCAAGCGCGATATCGAAGGCAATCGACTGTTttatc                                                                      <  2:526355/99‑1 (MQ=255)
                           ccccAGGAGTTACGCCAGGCGATTGAGCATATCAAAGCACACGTTACGGCAGAAACGCCAAAGGGTAAGCGCGATATCGAAGGCAATCGACTGTTttatc                                                                      <  2:5082236/100‑1 (MQ=255)
                           ccccAGGAGTTACGCCAGGCGATTGAGCATATCAAAGCACACGTTACGGCAGAAACGCCAAAGGGCAAGCGCGATATCGAAGGCAATCGACTGTTttatc                                                                      <  1:198502/100‑1 (MQ=255)
                            cccAGGAGTTACGCCAGGCGATTGAGCATATCAAAGCACACGTTACGGCAGAAACGCCAAAGGGCAAGCGCGATATCGAAGGCAATCGACTGTTttatct                                                                     <  1:3902097/100‑1 (MQ=255)
                              cAGGAGTTACGCCAGGCGATTGAGCATATCAAAGCACACGTTACGGCAGAAACGCCAAAGGGCAAGCGCGATATCGAAGGCAATCGACTGTTttatctta                                                                   <  2:3776358/100‑1 (MQ=255)
                               aGGAGTTACGCCAGGCGATTGAGCATATCAAAGCACACGTTACGGCAGAAACGCCAAAGGGCAAGCGCGATATCGAAGGCAATCGACTGTTttatcttat                                                                  <  2:820723/100‑1 (MQ=255)
                                ggAGTTACGCCAGGCGATTGAGCATATCAAAGCACACGTTACGGCAGAAACGCCAAAGGGCAAGCGCGATATCGAAGGCAATCGACTGTTttatcttatc                                                                 >  2:2948246/1‑100 (MQ=255)
                                 gAGTTACGCCAGGCGATGGAGCATATCAAAGCACGCGTTACGGCAGAAACGCCAAAGGGCAAGCGCGATATCGAAGGCAATCGACTGTTTtatcttatct                                                                <  1:456609/100‑1 (MQ=255)
                                   tttACGCCAGGCGATTGAGCATATCAAAGCACACGTTACGGCAGAAACGCCAAAGGGCAAGCGCGATATCGAAGGCAATCGACTGTTTTATCTTATCTCg                                                              <  1:4249225/99‑1 (MQ=255)
                                   gTTACGCCAGGCGATTGAGCATATCAAAGCACACGTTACGGCAGAAACGCCAAAGGGCAAGCGCGATATCGAAGGCAATCGACTGTTTTATCTTATCTCg                                                              >  1:2678465/1‑100 (MQ=255)
                                   gTTACGCCAGGCGATTGAGCATATCAAAGCACACGTTACGGCAGAAACGCCAAAGGGCAAGCGCGATATCGAAGGCAATCGACTGTTTTATCTTATCTCg                                                              <  2:1465899/100‑1 (MQ=255)
                                    ttACGCCAGGCGATTGAGCATATCAAAGCACACGTTACGGCAGAAACGCCAAAGGGCAAGCGCGATATCGAAGGCAATCGACTGTTTTATCTTATCTCgg                                                             <  1:1069418/100‑1 (MQ=255)
                                     tACGCCAGGCGATTGAGCATATCAAAGCACACGTTACGGCAGAAACGCCAAAGGGCAAGCGCGATATCGAAGGCAATCGACTGTTTTATCTTATCTCGGa                                                            <  2:1462757/100‑1 (MQ=255)
                                       cGCCAGGCGATTGAGCATATCAAAGCACACGTTACGGCAGAAACGCCAAAGGGCAAGCGCGATATCGAAGGCAATCGACTGTTTTATCTTATCTCGGAAg                                                          >  1:5160760/1‑100 (MQ=255)
                                       cGCCAGGCGATTGAGCATATCAAAGCACACGTTACGGCAGAAACGCCAAAGGGCAAGCGCGATATCGAAGGCAATCGACTGTTTTATCTTATCTCGGAAg                                                          >  1:689436/1‑100 (MQ=255)
                                        gCCAGGCGATTGAGCATATCAAAGCACACGTTACGGCAGAAACGCCAAAGGGCAAGCGCGATATCGAAGGCAATCGACTGTTTTATCTTATCTCGGAAGa                                                         <  1:511934/100‑1 (MQ=255)
                                         ccAGGCGATTGAGCATATCAAAGCACACGTTACGGCAGAAACGCCAAAGGGCAAGCGCGATATCGAAGGCAATCGACTGTTTTATCTTATCTCGGAAGat                                                        >  1:2650314/1‑100 (MQ=255)
                                          cAGGCGATTGAGCATATCAAAGCACACGTTACGGCAGAAACGCCAAAGGGCAAGCGCGATATCGAAGGCAATCGACTGTTTTATCTTATCTCGGAAGata                                                       <  2:4428087/100‑1 (MQ=255)
                                          cAGGCGATTGAGCATATCAAAGCACACGTTACGGCAGAAACGCCAAAGGGCAAGCGCGATATCGAAGGCAATCGACTGTTTTATCTTATCTCGGAAGata                                                       <  1:975931/100‑1 (MQ=255)
                                          cAGGCGATTGAGCATATCAAAGCACACGTTACGGCAGAAACGCCAAAGGGCAAGCGCGATATCGAAGGCAATCGACTGTTTTATCTTATCTCGGAAGata                                                       <  1:1958908/100‑1 (MQ=255)
                                           aGGCGATTGAGCATATCAAAGCACACGTTACGGCAGAAACGCCAAAGGGCAAGCGCGATATCGAAGGCAATCGACTGTTTTATCTTATCTCGGAAGatat                                                      >  1:1774125/1‑100 (MQ=255)
                                             gCGATTGAGCATATCAAAGCACACGTTACGGCATAAACGCCAAAGGGCAAGCGCGATATCGAAGGCAATCGACTGTTTTATCTTATCTCGGAAGATATGa                                                    <  2:3004403/100‑1 (MQ=255)
                                              cGATTGAGCATATCAAAGCACACGTTACGGCAGAAACGCCAAAGGGCAAGCGCGATATCGAAGGCAATCGACTGTTTTATCTTATCTCGGAAGATATGAc                                                   <  2:3640927/100‑1 (MQ=255)
                                              cGATTGAGCATAACAAAGCACACGTTACGGCAGAAACGCCAAAGGGCAAGCGCGATATCGAAGGCAATCGACTGTTTTATCTTATCTCGGAAGATATGAc                                                   <  2:4987861/100‑1 (MQ=255)
                                               gATTGAGCATATCAAAGCACACGTTACGGCAGAAACGCCAAAGGGCAAGCGCGATATCGAAGGCAATCGACTGTTTTATCTTATCTCGGAAGATATGAct                                                  <  1:1884860/100‑2 (MQ=255)
                                               gATTGAGCATATCAAAGCACACGTTACGGCAGAAACGCCAAAGGGCAAGCGCGATATCGAAGGCAATCGACTGTTTTATCTTATCTCGGAAGATATGAcc                                                  >  2:4275903/1‑100 (MQ=255)
                                               gATTGAGCATATCAAAGCACACGTTACGGCAGAAACGCCAAAGGGCAAGCGCGATATCGAAGGCAATCGACTGTTTTATCTTATCTCGGAAGATATGAcc                                                  <  2:1853912/100‑1 (MQ=255)
                                                  tGAGCATATCAAAGCACACGTTACGGCAGAAACGCCAAAGGGCAAGCGCGATATCGAAGGCAATCGACTGTTTTATCTTATCTCGGAAGATATGACCGAg                                               <  1:674284/100‑1 (MQ=255)
                                                   gAGCATATCAAAGCACACGTTACGGCAGAAACGCCAAAGGGCAAGCGCGATATCGAAGGCAATCGACTGTTTTATCTTATCTCGGAAGATATGACCGAGc                                              >  2:3991865/1‑100 (MQ=255)
                                                   gAGCATATCAAAGCACACGTTACGGCAGAAACGCCAAAGGGCAAGCGCGATATCGAAGGCAATCGACTGTTTTATCTTATCTCGGAAGATATGACCGAGc                                              >  1:461651/1‑100 (MQ=255)
                                                   gAGCATATCAAAGCACACGTTACGGCAGAAACGCCAAAGGGCAAGCGCGATATCGAAGGCAATCGACTGTTTTATCTTATCTCGGAAGATATGACCGAGc                                              >  1:2802523/1‑100 (MQ=255)
                                                       atatCAAAGCACACGTTACGGCAGAAACGCCAAAGGGCAAGCGCGATATCGAAGGCAATCGACTGTTTTATCTTATCTCGGAAGATATGACCGAGCCGTa                                          <  2:265871/100‑1 (MQ=255)
                                                       atatCAAAGCACACGTTACGGCAGAAACGCCAAAGGGCAAGCGCGATATCGAAGGCAATCGACTGTTTTATCTTATCTCGGAAGATATGACCGAGCCGTa                                          >  2:5274336/1‑100 (MQ=255)
                                                        tatCAAAGCACACGTTACGGCAGAAACGCCAAAGGGCAAGCGCGATATCGAAGGCAATCGACTGTTTTATCTTATCTCGGAAGATATGACCGAGCCGTAc                                         <  1:2244404/100‑1 (MQ=255)
                                                        tatCAAAGCACACGTTACGGCAGAAACGCCAAAGGGCAAGCGCGATATCGAAGGCAATCGACTGTTTTATCTTATCTCGGAAGATATGACCGAGCCGTAc                                         >  1:414424/1‑100 (MQ=255)
                                                          tCAAAGCACACGTTACGGCAGAAACGCCAAAGGGCAAGCGCGATATCGAAGGCAATCGACTGTTTTATCTTATCTCGGAAGATATGACCGAGCCGTACGa                                       >  1:5465509/1‑100 (MQ=255)
                                                          tCAAAGCACACGTTACGGCAGAAACGCCAAAGGGCAAGCGCGATATCGAAGGCAATCGACTGTTTTATCTTATCTCGGAAGATATGACCGAGCCGTACGa                                       <  1:596223/100‑1 (MQ=255)
                                                          tCAAAGCACACGTGACGGCAGAAACGCCAAACGGCAAGCGCGATATCGAAGGCAATCGACTGTTTTATCTTATCTCGGAAGATATGACCGAGCCGTACGa                                       <  2:4307227/100‑1 (MQ=255)
                                                           cAAAGCACACGTTACGGCAGAAACGCCAAAGGGCAAGCGCGATATCGAAGGCAATCGACTGTTTTATCTTATCTCGGAAGATATGACCGAGCCGTACGaa                                      >  1:4047373/1‑100 (MQ=255)
                                                           cAAAGCACACGTTACGGCAGAAACGCCAAAGGGCAAGCGCGATATCGAAGGCAATCGACTGTTTTATCTTATCTCGGAAGATATGACCGAGCCGTACGaa                                      >  1:2321755/1‑100 (MQ=255)
                                                            aaaGCACACGTTACGGCAGAAACGCCAAAGGGCAAGCGCGATATCGAAGGCAATCGACTGTTTTATCTTATCTCGGAAGATATGACCGAGCCGTACGAAg                                     >  2:5021486/1‑100 (MQ=255)
                                                             aaGCACACGTTACGGCAGAAACGCCAAAGGGCAAGCGCGATATCGAAGGCAATCGACTGTTTTATCTTATCTCGGAAGATATGACCGAGCCGTACGAAGc                                    <  2:3802866/100‑1 (MQ=255)
                                                                cacaCGTTACGGCAGAAACGCCAAAGGGCAAGCGCGATATCGAAGGCAATCGACTGTTTTATCTTATCTCGGAAGATATGACCGAGCCGTACGAAGCTCg                                 >  2:4782808/1‑100 (MQ=255)
                                                                    cGTTACGGCAGAAACGCCAAAGGGCAAGCGCGATATCGAAGGCAATCGACTGTTTTATCTTATCTCGGAAGATATGACCGAGCCGTACGAAGCTCGCCGt                             >  1:771002/1‑100 (MQ=255)
                                                                    cGTTACGGCAGAAACGCCAAAGGGCAAGCGCGATATCGAAGGCAATCGACTGTTTTATCTTATCTCGGAAGATATGACCGAGCCGTACGAAGCTCGCCGt                             >  1:505231/1‑100 (MQ=255)
                                                                    cGTTACGGCAGAAACGCCAAAGGGCAAGCGCGATATCGAAGGCAATCGACTGTTTTATCTTATCTCGGAAGATATGACCGAGCCGTACGAAGCTCGCCGt                             <  1:3177434/100‑1 (MQ=255)
                                                                       tACGGCGGAAACGCCAAAGGGCAAGCGCGATATCGAAGGCAATCGACTGTTTTATCTTATCTCGGAAGATATGACCGAGCCGTACGAAGCTCGCCGTGCg                          <  1:5200113/100‑1 (MQ=255)
                                                                       tACGGCAGAAACGCCAAAGGGCAAGCGCGATATCGAAGGCAATCGACTGTTTTATCTTATCTCGGAAGATATGACCGAGCCGTACGAAGCTCGCCGTGCg                          <  2:1543510/100‑1 (MQ=255)
                                                                         cGGCAGAAACGCCAAAGGGCAAGCGCTATATCGAAGGCAATCGACTGTTTTATCTTGTCTCGGAAGATATGACCGAGCCGTACGAAGCTCGCCGTGCGGa                        <  2:690884/100‑1 (MQ=255)
                                                                         cGGCAGAAACGCCAAAGGGCAAGCGCGATATCGAAGGCAATCGACTGTTTTATCTTATCTCGGAAGATATGACCGAGCCGTACGAAGCTCGCCGTGCggc                        >  1:5568011/1‑99 (MQ=255)
                                                                         cGGCAGAAACGCCAAAGGGCAAGCGCGATATCGAAGGCAATCGACTGTTTTATCTTATCTCGGAAGATATGACCGAGCCGTACGAAGCTCGCCGTGCGGa                        <  1:2059952/100‑1 (MQ=255)
                                                                          ggCAGAAACGCCAAAGGGCAAGCGCGATATCGAAGGCAATCGACTGTTTTATCTTATCTCGGAAGATATGACCGAGCCGTACGAAGCTCGCCGTGCGGAg                       >  2:1262281/1‑100 (MQ=255)
                                                                           gCAGAAACGCCAAAGGGCAAGCGCGATATCGAAGGCAATCGACTGTTTTATCTTATCTCGGAAGATATGACCGAGCCGTACGAAGCTCGCCGTGCGGAGt                      >  2:4188970/1‑100 (MQ=255)
                                                                           gCAGAAACGCCAAAGGGCAAGCGCGATATCGAAGGCAATCGACTGTTTTATCTTATCTCGGAAGATATGACCGAGCCGTACGAAGCTCGCCGTGCGGAGt                      <  1:1178487/100‑1 (MQ=255)
                                                                           gCAGAAACGCCAAAGGGCAAGCGCGATATCGAAGGCAATCGACTGTTTTATCTTATCTCGGAAGATATGACCGAGCCGTACGAAGCTCGCCGTGCGGAGt                      >  1:366288/1‑100 (MQ=255)
                                                                            cAGAAACGCCACACGGCAACCGCGATATCGAAGGCAATCGACTGTTTTACCTTATCTCGGAAGATAGGGCCGTGCCGTACGAAGCTCGCCGTGCGGAGTa                     <  2:3382737/100‑1 (MQ=255)
                                                                            cAGAAACGCCAAAGGGCAAGCGCGATATCGAAGGCAATCGACTGTTTTATCTTATCTCGGAAGATATGACCGAGCCGTACGAAGCTCGCCGTGCGGAGTa                     <  1:3269319/100‑1 (MQ=255)
                                                                             aGAAACGCCAAAGGGCAAGCGCGATATCGAAGGCAATCGACTGTTTTATCTTATCTCGGAAGATATGACCGAGCCGTACGAAGCTCGCCGTGCGGAGTAc                    <  2:614090/100‑1 (MQ=255)
                                                                              gAAACGCCAAAGGGCAAGCGCGATATCGAAGGCAATCGACTGTTTTATCTTATCTCGGAAGATATGACCGAGCCGTACGAAGCTCGCCGTGCGGAGTAcc                   >  1:3448032/1‑100 (MQ=255)
                                                                              gAAACGCCAAAGGGCAAGCGCGATATCGAAGGCAATCGACTGTTTTATCTTATCTCGGAAGATATGACCGAGCCGTACGAAGCTCGCCGTGCGGAGTAcc                   >  1:3537672/1‑100 (MQ=255)
                                                                              gAAACGCCAAAGGGCAAGCGCGATATCGAAGGCAATCGACTGTTTTATCTTATCTCGGAAGATATGACCGAGCCGTACGAAGCTCGCCGTGCGGAGTAcc                   >  2:1092085/1‑100 (MQ=255)
                                                                                 aCGCCAAAGGGCAAGCGCGATATCGAAGGCAATCGACTGTTTTATCTTATCTCGGAAGATATGACCGAGCCGTACGAAGCTCGCCGTGCGGAGTACCATg                <  1:2387969/100‑1 (MQ=255)
                                                                                 aCGCCAAAGGGCAAGCGCGATATCGAAGGCAATCGACTGTTTTATCTTATCTCGGAAGATATGACCGAGCCGTACGAAGCTCGCCGTGCGGAGTACCATg                <  2:3848565/100‑1 (MQ=255)
                                                                                   gCCAAAGGGCAAGCGCGATATCGAAGGCAATCGACTGTTTTATCTTATCTCGGAAGATATGACCGAGCCGTACGAAGCTCGCCGTGCGGAGTACCATGcc              >  2:5176128/1‑100 (MQ=255)
                                                                                    ccAAAGGGCAAGCGCGATATCGAAGGCAATCGACTGTTTTATCTTATCTCGGAAGATATGACCGAGCCGTACGAAGCTCGCCGTGCGGAGTACCATGccc             >  1:530724/1‑100 (MQ=255)
                                                                                    ccAAAGGGCAAGCGCGATATCGAAGGCAATCGACTGTTTTATCTTATCTCGGAAGATATGACCCAGCCGTGCGAAGCTCGCCGTGCGGCGTACCATGccc             >  2:1078471/1‑100 (MQ=255)
                                                                                       aaGGGCAAGCGCGTGATCGAAGGCAATTGACTGTTTTATGTGAGCTCGGAAGAGATGACCGAGCCGTACGAAGCTCGCCGTGCGTAGTACCATGCCCGCt          <  2:3915781/100‑1 (MQ=255)
                                                                                        aGGGCAAGCGCGATATCGAAGGCAATCGACTGTTTTATCTTATCTCGGAAGATATGACCGAGCCGTACGAAGCTCGCCGTGCGGAGTACCATGCCCGCTa         <  1:3940369/100‑1 (MQ=255)
                                                                                         gggCAAGCGCGATATCGAAGGCAATCGACTGTTTTATCTTATCTCGGAAGATATGACCGAGCCGTACGAAGCTCGCCGTGCGGAGTACCATGCCCGCTAt        >  2:5500284/1‑100 (MQ=255)
                                                                                          ggCAAGCGCGATATCGAAGGCAATCGACTGTTTTATCTTATCTCGGAAGATATGACCGAGCCGTACGAAGCTCGCCGTGCGGAGTACCATGCCCGCTATc       <  2:5073487/100‑1 (MQ=255)
                                                                                              aGCGCGATATCGAAGGCAATCGACTGTTTTATCTTATCTCGGAAGATATGACCGAGCCGTACGAAGCTCGCCGTGCGGAGTACCATGCCCGCTATCTTGa   >  2:4503009/1‑100 (MQ=255)
                                                                                                 |                                                                                                 
AATATTCATAATCTTCAGCCGTGGCTACCCCAGGAGTTACGCCAGGCGATTGAGCATATCAAAGCACACGTTACGGCAGAAACGCCAAAGGGCAAGCACGATATCGAAGGCAATCGACTGTTTTATCTTATCTCGGAAGATATGACCGAGCCGTACGAAGCTCGCCGTGCGGAGTACCATGCCCGCTATCTTGAC  >  NC_000913/4474874‑4475068

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