Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,597,539 T→C intergenic (‑250/+21) b3460 ← / ← rpoH LivJ/sigma(32) factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,597,5390TC99.6% 795.8 / NA 246intergenic (‑250/+21)b3460/rpoHLivJ/sigma(32) factor
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (123/122);  minor base A (1/0);  total (124/122)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 252 total aligned reads displayed.

ACAATCTGCGTGGGGATTGGCGTTTTGCCGGAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGTTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCATCGCGTTCTTTTCCAGCTGGCGTACACGCTCAGCGGA  >  NC_000913/3597449‑3597628
                                                                                          |                                                                                         
aCAATCTGCGTGGGGATTGGCGTTTTGCCGGAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCtc                                                                                         >  1:6067946/1‑93 (MQ=255)
aCAATCTGCGTGGGGATTGGCGTTTTGCCGGAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCtc                                                                                         <  1:2021433/93‑1 (MQ=255)
 cAATCTGCGTGGGGATTGGCGTTTTGCCGGAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGctct                                                                                        <  2:193653/93‑1 (MQ=255)
  aaTCTGCGTGGGGATTGGCGTTTTGCCGGAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCtctc                                                                                       >  2:2919927/1‑93 (MQ=255)
  aaTCTGCGTGGGGATTGGCGTTTTGCCGGAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCtctc                                                                                       >  1:4720465/1‑93 (MQ=255)
  aaTCTGCGTGGGGATTGGCGTTTTGCCGGAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCtctc                                                                                       >  2:3762205/1‑93 (MQ=255)
  aaTCTGCGTGGGGATTGGCGTTTTGCCGGAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCtctc                                                                                       >  2:6491836/1‑93 (MQ=255)
      tGCGTGGGGATTGGCGTTTTGCCGGAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCt                                                                                   <  2:1185097/93‑1 (MQ=255)
      tGCGTGGGGATTGGCGTTTTGCCGGAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCt                                                                                   >  2:12503/1‑93 (MQ=255)
       gCGTGGGGATTGGCGTTTTGCCGGAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCtt                                                                                  <  1:3267762/93‑1 (MQ=255)
       gCGTGGGGATTGGCGTTTTGCCGGAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCtt                                                                                  >  1:629244/1‑93 (MQ=255)
       gCGTGGGGATTGGCGTTTTGCCGGAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCtt                                                                                  >  2:6036260/1‑93 (MQ=255)
       gCGTGGGGATTGGCGTTTTGCCGGAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCtt                                                                                  >  2:6622280/1‑93 (MQ=255)
       gCGTGGGGATTGGCGTTTTGCCGGAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCtt                                                                                  <  2:2788808/93‑1 (MQ=255)
         gTGGGGATTGGCGTTTTGCCGGAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTaa                                                                                >  2:3617773/1‑93 (MQ=255)
          tGGGGATTGGCGTTTTGCCGGAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTAAt                                                                               >  2:5891439/1‑93 (MQ=255)
          tGGGGATTGGCGTTTTGCCGGAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTAAt                                                                               >  2:5324013/1‑93 (MQ=255)
           ggggATTGGCGTTTTGCCGGAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTAATa                                                                              <  1:1630199/93‑1 (MQ=255)
            gggATTGGCGTTTTGCCGGAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAg                                                                             <  2:3891520/93‑1 (MQ=255)
                ttGGCGTTTTGCCGGAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGa                                                                         >  2:5181621/1‑93 (MQ=255)
                ttGGCGTTTTGCCGGAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGa                                                                         <  1:1507931/93‑1 (MQ=255)
                  ggCGTTTTGCCGGAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGaaa                                                                       <  2:3630127/93‑1 (MQ=255)
                  ggCGTTTTGCCGGAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGaaa                                                                       >  1:6056129/1‑93 (MQ=255)
                   gCGTTTTGCCGGAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAAt                                                                      <  1:4424670/93‑1 (MQ=255)
                   gCGTTTTGCCGGAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAAt                                                                      <  1:658507/93‑1 (MQ=255)
                     gTTTTGCCGGAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTa                                                                    <  1:3464181/93‑1 (MQ=255)
                     gTTTTGCCGGAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTa                                                                    >  2:4071710/1‑93 (MQ=255)
                       tttGCCGGAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACg                                                                  >  1:4647499/1‑93 (MQ=255)
                        ttGCCGGAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGc                                                                 >  1:4535162/1‑93 (MQ=255)
                        ttGCCGGAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGc                                                                 <  1:5883263/93‑1 (MQ=255)
                         tGCCGGAGGGGAAATTGATTATTACAGAGGCCCAAACAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGGAATTACGCt                                                                >  2:3268433/1‑93 (MQ=255)
                         tGCCGGAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCt                                                                <  2:1083462/93‑1 (MQ=255)
                          gCCGGAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCtt                                                               >  2:6685082/1‑93 (MQ=255)
                            cGGAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGCAATTACGCTTCa                                                             >  2:468355/1‑93 (MQ=255)
                             ggAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCaa                                                            >  1:5354792/1‑93 (MQ=255)
                              gAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAAt                                                           <  1:535706/93‑1 (MQ=255)
                              gAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAAt                                                           <  2:3689031/93‑1 (MQ=255)
                              gAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAAt                                                           <  2:6241126/93‑1 (MQ=255)
                                  ggAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGgca                                                       >  2:5701219/1‑93 (MQ=255)
                                   gAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGgcag                                                      >  2:802766/1‑93 (MQ=255)
                                   gAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGgcac                                                      >  1:413412/1‑92 (MQ=255)
                                     aaTTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGgcagca                                                    >  2:1682312/1‑93 (MQ=255)
                                      aTTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCAc                                                   >  1:6518211/1‑93 (MQ=255)
                                      aTTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCAc                                                   >  1:6730094/1‑93 (MQ=255)
                                         gATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCa                                                >  2:1367215/1‑93 (MQ=255)
                                           ttattaCAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAAt                                              >  1:788234/1‑93 (MQ=255)
                                            tattaCAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAAtt                                             <  2:2022495/93‑1 (MQ=255)
                                             attaCAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAAttt                                            <  2:346141/93‑1 (MQ=255)
                                             attaCAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAAttt                                            <  1:5767969/93‑1 (MQ=255)
                                              ttaCAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAAtttt                                           <  1:775807/93‑1 (MQ=255)
                                              ttaCAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAAtttt                                           >  2:2002184/1‑93 (MQ=255)
                                              ttaCAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAAtttt                                           <  1:4773083/93‑1 (MQ=255)
                                                aCAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAAtttttt                                         <  1:6677519/93‑1 (MQ=255)
                                                  agagGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCa                                       >  2:4942392/1‑93 (MQ=255)
                                                   gagGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCAt                                      <  1:1363725/93‑1 (MQ=255)
                                                    agGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCATc                                     <  2:5298296/93‑1 (MQ=255)
                                                     ggCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCATcg                                    >  1:3843308/1‑93 (MQ=255)
                                                     ggCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCATcg                                    <  1:1967954/93‑1 (MQ=255)
                                                     ggCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCATcg                                    >  2:4292542/1‑93 (MQ=255)
                                                       cccAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCATcgcg                                  >  1:5454913/1‑93 (MQ=255)
                                                        ccAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCATCGCGt                                 >  1:1017652/1‑93 (MQ=255)
                                                         cAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCATCGCGtt                                >  1:2768137/1‑93 (MQ=255)
                                                          aaaaaaCAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCATCGCGTTc                               <  1:6161986/93‑1 (MQ=255)
                                                           aaaaacaaaaacCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCATCGCGTTCt                              >  2:618318/1‑93 (MQ=255)
                                                           aaaaacaaaaacCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCATCGCGTTCt                              <  2:4414943/93‑1 (MQ=255)
                                                             aaacgaaaacCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCATCGCGTTCttt                            >  1:1236262/1‑93 (MQ=255)
                                                             aaacaaaaacCCCGGGCTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTCCGCTTCAATGGCAGCACGCAATTTTTTCATCGCGTTCttt                            <  2:3754283/93‑1 (MQ=255)
                                                               acaaaaacCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCATCGCGTTCTTTTc                          <  1:3428277/93‑1 (MQ=255)
                                                                caaaaacCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCATCGCGTTCTTTTcc                         >  1:2872920/1‑93 (MQ=255)
                                                                 aaaaacCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCATCGCGTTCTTTTCCa                        <  2:3313825/93‑1 (MQ=255)
                                                                 aaaaacCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCATCGCGTTCTTTTCCa                        <  1:4952662/93‑1 (MQ=255)
                                                                  aaaacCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCATCGCGTTCTTTTCCAg                       <  2:681128/93‑1 (MQ=255)
                                                                  aaaacCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCATCGCGTTCTTTTCCAg                       >  2:5477596/1‑93 (MQ=255)
                                                                   aaacCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCATCGCGTTCTTTTCCAGc                      <  2:5822377/93‑1 (MQ=255)
                                                                   aaacCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCATCGCGTTCTCTTCCAGc                      >  1:2460731/1‑93 (MQ=255)
                                                                    aacCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCATCGCGTTCTTTTCCAGCt                     >  2:3201152/1‑93 (MQ=255)
                                                                       cccGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCATCGCGTTCTTTTCCAGCTGGc                  <  2:5041321/93‑1 (MQ=255)
                                                                        ccGGACTCTCATCCAGGGATCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGAAATTTTTTCATCGCGTTCTTTTCCAGCTGGCg                 >  1:174007/1‑93 (MQ=255)
                                                                          ggACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCATCGCGTTCTTTTCCAGCTGGCGTa               >  2:2720293/1‑93 (MQ=255)
                                                                          ggACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCATCGCGTTCTTTTCCAGCTGGCGTa               <  1:1922122/93‑1 (MQ=255)
                                                                          ggACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCATCGCGTTCTTTTCCAGCTGGCGTa               <  1:4292267/93‑1 (MQ=255)
                                                                          ggACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCATCGCGTTCTTTTCCAGCTGGCGTa               <  2:1003079/93‑1 (MQ=255)
                                                                           gACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCATCGCGTTCTTTTCCAGCTGGCGTac              <  1:2688788/93‑1 (MQ=255)
                                                                           gACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCATCGCGTTCTTTTCCAGCTGGCGTac              <  2:1645860/93‑1 (MQ=255)
                                                                            aCTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCATCGCGTTCTTTTCCAGCTGGCGTaca             >  1:89693/1‑93 (MQ=255)
                                                                             ctctCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCATCGCGTTCTTTTCCAGCTGGCGTacac            <  1:5312347/93‑1 (MQ=255)
                                                                              tctcATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCATCGCGTTCTTTTCCAGCTGGCGTACACg           >  2:2723633/1‑93 (MQ=255)
                                                                               ctcATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCATCGCGTTCTTTTCCAGCTGGCGTACACGc          <  1:3990547/93‑1 (MQ=255)
                                                                                tcATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCATCGCGTTCTTTTCCAGCTGGCGTACACGCt         <  2:462785/93‑1 (MQ=255)
                                                                                   tCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCATCGCGTTCTTTTCCAGCTGGCGTACACGCTCAg      <  2:2206061/93‑1 (MQ=255)
                                                                                   tCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCATCGCGTTCTTTTCCAGCTGGCGTACACGCTCAg      >  2:237764/1‑93 (MQ=255)
                                                                                    ccAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCATCGCGTTCTTTTCCAGCTGGCGTACACGCTCAGc     >  1:2116602/1‑93 (MQ=255)
                                                                                    ccAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCATCGCGTTCTTTTCCAGCTAGCCTACACGCTCAGc     >  2:4830054/1‑93 (MQ=255)
                                                                                     cAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCATCGCGTTCTTTTCCAGCTGGCGTACACGCTCAGCg    >  1:2306632/1‑93 (MQ=255)
                                                                                     cAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCATCGCGTTCTTTTCCAGCTGGCGTACACGCTCAGCg    >  1:4888280/1‑93 (MQ=255)
                                                                                      aGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCATCGCGTTCTTTTCCAGCTGGCGTACACGCTCAGCgg   >  2:5562387/1‑93 (MQ=255)
                                                                                      aGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCATCGCGTTCTTTTCCAGCTGGCGTACACGCTCAGCgg   <  1:5485724/93‑1 (MQ=255)
                                                                                      aGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCATCGCGTTCTTTTCCAGCTGGCGTACACGCTCAGCgg   <  1:1906614/93‑1 (MQ=255)
                                                                                       gggCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCATCGCGTTCTTTTCCAGCTGGCGTACACGCTCAGCGGa  <  1:3033590/93‑1 (MQ=255)
                                                                                       gggCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCATCGCGTTCTTTTCCAGCTGGCGTACACGCTCAGCGGa  >  1:5230286/1‑93 (MQ=255)
                                                                                          |                                                                                         
ACAATCTGCGTGGGGATTGGCGTTTTGCCGGAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGTTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCATCGCGTTCTTTTCCAGCTGGCGTACACGCTCAGCGGA  >  NC_000913/3597449‑3597628

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 23 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: