Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 298,399 C→T 5.5% V110I (GTT→ATT)  paoD ← moco insertion factor for PaoABC aldehyde oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913298,3990CT 5.5% 239.1 / 3.4 72V110I (GTT→ATT) paoDmoco insertion factor for PaoABC aldehyde oxidoreductase
Reads supporting (aligned to +/- strand):  ref base C (39/29);  new base T (2/2);  total (41/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.54e-01

GTGGGCAGGCACACCAGCGATTGTGCTTGCGGATCGTAGCGCAGCCCCACCGGTTTTCTCTGTTCCAGGCGGTTCAGCACGGCGAGCAGAGGCTGTGCCGAGCGTAGTTTATGGAGCGTCAGCGTGATCCCACCGCCGCAGGGCAGAACGATGTCAAACCACGGCGAACCTTCGCCATAGCGAATTTCGCGATCGCGGCCTGAGCCCATCATCTCCAGCGCTTCAAAAGCGGCAGCGGCCTCCACGCAGCCGCCAGAGACAAAACCGCAGTAACGACCATCTTCGCGCACCAC  >  NC_000913/298251‑298543
                                                                                                                                                    |                                                                                                                                                
gTGGGCAGGCACACCAGCGATTGTGCTTGCGGATCGTAGCGCAGCCCCACCGGTTTTCTCTGTTCCAGGCGGTTCAGCACGGCGAGCAGAGGCTGTGCCGAGCGTAGTTTATGGAGCGTCAGCGTGATCCCACCGCCGCAGGGCAGAAc                                                                                                                                                  <  1:172341/149‑1 (MQ=255)
  gggCAGGCACACCAGCGATTGTGCTTGCGGATCGTAGCGCAGCCCCACCGGTTTTCTCTGTTCCAGGCGGTTCAGCACGGCGAGCAGAGGCTGTGCCGAGCGTAGTTTATGGAGCGTCAGCGTGATCCCACCGCCGCAGGGCAGAACGa                                                                                                                                                >  4:735549/1‑149 (MQ=255)
  gggCAGGCACACCAGCGATTGTGCTTGCGGATCGTAGCGCAGCCCCACCGGTTTTCTCTGTTCCAGGCGGTTCAGCACGGCGAGCAGAGGCTGTGCCGAGCGTAGTTTATGGAGCGTCAGCGTGATCCCACCGCCGCAGGGCAGAACGa                                                                                                                                                <  1:228386/149‑1 (MQ=255)
       ggcaCACCAGCGATTGTGCTTGCGGATCGTAGCGCAGCCCCACCGGTTTTCTCTGTTCCAGGCGGTTCAGCACGGCGAGCAGAGGCTGTGCCGAGCGTAGTTTATGGAGCGTCAGCGTGATCCCACCGCCGCAGGGCAGAACGATGTCa                                                                                                                                           >  3:145957/1‑149 (MQ=255)
        gcaCACCAGCGCTTGTGCTTGCGGATCGTAGCGCAGCCCCACCGGTTTTCTCTGTTCCAGGCGGTTCAGCACGGCGAGCAGAGGCTGTGCCGAGCGTAGTTTATGGAGCGTCAGCGTGATCCCACCGCCGCAGGGCAGAACGATGTCaa                                                                                                                                          >  3:655500/1‑149 (MQ=255)
           cacCAGCGATTGTGCTTGCGGATCGTAGCGCAGCCCCACCGGTTTTCTCTGTTCCAGGCGGTTCAGCACGGCGAGCAGAGGCTGTGCCGAGCGTAGTTTATGGAGCGTCAGCGTGATCCCACCGCCGCAGGGCAGAACGATGTCAAAcc                                                                                                                                       >  2:621358/1‑149 (MQ=255)
            acCAGCGATTGTGCTTGCGGATCGTAGCGCAGCCCCACCGGTTTTCTCTGTTCCAGGCGGTTCAGCACGGCGAGCAGAGGCTGTGCCGAGCGTAGTTTATGGAGCGTCAGCGTGATCCCACCGCCGCAGGGCAGAACGATGTCAAACCa                                                                                                                                      >  4:703430/1‑149 (MQ=255)
            acCAGCGATTGTGCTTGCGGATCGTAGCGCAGCCCCACCGGTTTTCTCTGTTCCAGGCGGTTCAGCACGGCGAGCAGAGGCTGTGCCGAGCGTAGTTTATGGAGCGTCAGCGTGATCCCACCGCCGCAGGGCAGAACGATGTCAAACCa                                                                                                                                      >  2:406999/1‑149 (MQ=255)
             ccAGCGATTGTGCTTGCGGATCGTAGCGCAGCCCCACCGGTTTTCTCTGTTCCAGGCGGTTCAGCACGGCGAGCAGAGGCTGTGCCGAGCGTAGTTTATGGAGCGTCAGCGTGATCCCACCGCCGCAGGGCAGAACGATGTCAAACCAc                                                                                                                                     >  2:666200/1‑149 (MQ=255)
                  gATTGTGCTTGCGGATCGTAGCGCAGCCCCACCGGTTTTCTCTGTTCCAGGCGGTTCAGCACGGCGAGCAGAGGCTGTGCCGAGCGTAGTTTATGGAGCGTCAGCGTGATCCCACCGCCGCAGGGCAGAATGATGTCAAACCACGGCGa                                                                                                                                >  1:182226/1‑149 (MQ=255)
                  gATTGTGCTTGCGGATCGTAGCGCAGCCCCACCGGTTTTCTCTGTTCCAGGCGGTTCAGCACGGCGAGCAGAGGCTGTGCCGAGCGTAGTTTATGGAGCGTCAGCGTGATCCCACCGCCGCAGGGCAGAATGATGTCAAACCACGGCGa                                                                                                                                >  1:753884/1‑149 (MQ=255)
                   aTTGTGCTTGCGGATCGTAGCGCAGCCCCACCGGTTTTCTCTGTTCCAGGCGGTTCAGCGCGGCGAGCAGAGGCTGTGCCGAGCGTAGTTTATGGAGCGTCAGCGTGATCCCACCGCCGCAGGGCAGAACGATGTCAAACCACGGCGaa                                                                                                                               <  2:247236/149‑1 (MQ=255)
                    ttGTGCTTGCGGATCGTAGCGCAGCCCCACCGGTTTTCTCTGTTCCAGGCGGTTCAGCACGGCGAGCAGAGGCTGTGCCGAGCGTAGTTTATGGAGCGTCAGCGTGATCCCACCGCCGCAGGGCAGAACGATGTCAAACCACGGCGAAc                                                                                                                              <  2:602038/149‑1 (MQ=255)
                          ttGCGGATCGTAGCGCAGCCCCACCGGTTTTCTCTGTTCCAGGCGGTTCAGCACGGCGAGCAGAGGCTGTGCCGAGCGTAGTTTATGGAGCGTCAGCGTGATCCCACCGCCGCAGGGCAGAACGATGTCAAACCACGGCGAACCTTCGc                                                                                                                        >  3:593561/1‑149 (MQ=255)
                          ttGCGGATCGTAGCGCAGCCCCACCGGTTTTCTCTGTTCCAGGCGGTTCAGCACGGCGAGCAGAGGCTGTGCCGAGCGTAGTTTATGGAGCGTCAGCGTGATCCCACCGCCGCAGGGCAGAACGATGTCAAACCACGGCGAACCTTCGc                                                                                                                        <  1:533889/149‑1 (MQ=255)
                              ggATCGTAGCGCAGCCCCACCGGTTTTCTCTGTTCCAGGCGGTTCAGCACGGCGAGCAGAGGCTGTGCCGAGCGTAGTTTATGGAGCGTCAGCGTGATCCCACCGCCGCAGGGCAGAACGATGTCAAACCACGGCGAACCTTCGCCATa                                                                                                                    <  2:270271/149‑1 (MQ=255)
                                  cGTAGCGCAGCCCCACCGGTTTTCTCTGTTCCAGGCGGTTCAGCACGGCGAGCAGAGGCTGTGCCGAGCGTAGTTTATGGAGCGTCAGCGTGATCCCACCGCCGCAGGGCAGAACGATGTCAAACCACGGCGAAGGTTCGccgtggttt                                                                                                                <  1:503523/149‑8 (MQ=255)
                                   gTAGCGCAGCCCCACCGGTTTTCTCTGTTCCAGGCGGTTCAGCACGGCGAGCAGAGGCTGTGCCGAGCGTAGTTTATGGAGCGTCAGCGTGATCCCACCGCCGCAGGGCAGAACGATGTCAAACCACGGCGAACCTTCGCCATAGCGaa                                                                                                               >  2:536220/1‑149 (MQ=255)
                                     aGCGCAGCCCCACCGGTTTTCTCTGTTCCAGGCGGTTCAGCACGGCGAGCAGAGGCTGTGCCGAGCGTAGTTTATGGAGCGTCAGCGTGATCCCACCGCCGCAGGGCAGAACGATGTCAAACCACGGCGAACCTTCGCTATGGCGaagg                                                                                                             <  2:569883/149‑3 (MQ=255)
                                      gcgcAGCCCCACCGGTTTTCTCTGTTCCAGGCGGTTCAGCACGGCGAGCAGAGGCTGTGCCGAGCGTAGTTTATGGAGCGTCAGCGTGATCCCACCGCCGCAGGGCAGAACGATGTCAAACCACGGCGAACCTTCGCCATAGCGAAttt                                                                                                            <  2:17376/149‑1 (MQ=255)
                                         cAGCCCCACCGGTTTTCTCTGTTCCAGGCGGTTCAGCACGGCGAGCAGAGGCTGTGCCGAGCGTAGTTTATGGAGCGTCAGCGTGATCCCACCGCCGCAGGGCAGAACGATGTCAAACCACGGCGAACCTTCGCCATAGCGAATTTcgc                                                                                                         <  3:735549/149‑1 (MQ=255)
                                         cAGCCCCACCGGTTTTCTCTGTTCCAGGCGGTTCAGCACGGCGAGCAGAGGCTGTGCCGAGCGTAGTTTATGGAGCGTCAGCGTGATCCCACCGCCGCAGGGCAGAACGATGTCAAACCACGGCGAACCTTCGCCATAGCGAATTTcgc                                                                                                         <  4:704692/149‑1 (MQ=255)
                                           gCCCCACCGGTTTTCTCTGTTCCAGGCGGTTCAGCACGGCGAGCAGAGGCTGTGCCGAGCGTAGTTTATGGAGCGTCAGCGTGATCCCACCGCCGCAGGGCAGAACGATGTCAAACCACGGCGAACCTTCGCCATAGCGAATTTCGCGa                                                                                                       <  3:282320/149‑1 (MQ=255)
                                              ccACCGGTTTTCTCTGTTCCAGGCGGTTCAGCACGGCGAGCAGAGGCTGTGCCGAGCGTAGTTTATGGAGCGTCAGCGTGATCCCACCGCCGCAGGGCAGAACGATGTCAAACCACGGCGAACCTTCGCCATAGCGAATTTCGCGATcg                                                                                                    >  4:257234/1‑149 (MQ=255)
                                                aCCGGTTTTCTCTGTTCCAGGCGGTTCAGCACGGCGAGCAGAGGCTGTGCCGAGCGTAGTTTATGGAGCGTCAGCGTGATCCCACCGCCGCAGGGCAGAACGATGTCAAACCACGGCGAACCTTCGCCATAGCGAATTTCGCGATcgcg                                                                                                  >  4:324561/1‑149 (MQ=255)
                                                 ccGGTTTTCTCTGTTCCAGGCGGTTCAGCACGGCGAGCAGAGGCTGTGCCGAGCGTAGTTTATGGAGCGTCAGCGTGATCCCACCGCCGCAGGGCAGAACGATGTCAAACCACGGCGAACCTTCGCCATAGCGAATTTCGCGATCGCgg                                                                                                 <  3:267714/149‑1 (MQ=255)
                                                  cGGTTTTCTCTGTTCCAGGCGGTTCAGCACGGCGAGCAGAGGCTGTGCCGAGCGTAGTTTATGGAGCGTCAGCGTGATCCCACCGCCGCAGGGCAGAACGATGTCAAACCACGGCGAACCTTCGCCATAGCGAATTTCGCGATCGCGGc                                                                                                >  1:615105/1‑149 (MQ=255)
                                                  cGGTTTTCTCTGTTCCAGGCGGTTCAGCACGGCGAGCAGAGGCTGTGCCGAGCGTAGTTTATGGAGCGTCAGCGTGATCCCACCGCCGCAGGGCAGAACGATGTCAAACCACGGCGAACCTTCGCCATAGCGAATTTCGCGATCGCGGc                                                                                                >  2:373943/1‑149 (MQ=255)
                                                  cGGTTTTCTCTGTTCCAGGCGGTTCAGCACGGCGAGCAGAGGCTGTGCCGAGCGTAGTTTAGGGAGCGTCAGCGTGATCCCACCGCCGCAGGGCAGAACGATGTCAAACCACGGCGAACCTTCGTCATAGCGAATTTCGCGATcgcgac                                                                                                <  2:8736/149‑3 (MQ=255)
                                                     ttttCTCTGTTCCAGGCGGTTCAGCACGGCGAGCAGAGGCTGTGCCGAGCGTAGTTTATGGAGCGTCAGCGTGATCCCACCGCCGCAGGGCAGAACGATGTCAAACCACGGCGAACCTTCGCCATAGCGAATTTCGCGATCGCGGCCTg                                                                                             <  4:145957/149‑1 (MQ=255)
                                                          tctGTTCCAGGCGGTTCAGCACGGCGAGCAGAGGCTGTGCCGAGCGTAGTTTATGGAGCGTCAGCGTGATCCCACCGCCGCAGGGCAGAACGATGTCAAACCACGGCGAACCTTCGCCATAGCGAATTTCGCGATCGCGGCCTGAGccc                                                                                        >  2:217061/1‑149 (MQ=255)
                                                          tctGTTCCAGGCGGTTCAGCACGGCGAGCAGAGGCTGTGCCGAGCGTAGTTTATGGAGCGTCAGCGTGATCCCACCGCCGCAGGGCAGAACGATGTCAAACCACGGCGAACCTTCGCCATAGCGAATTTCGCGATCGCGGCCTGAGccc                                                                                        >  2:285096/1‑149 (MQ=255)
                                                          tctGTTCCAGGCGGTTCAGCACGGCGAGCAGAGGCTGTGCCGAGCGTAGTTTATGGAGCGTCAGCGTGATCCCACCGCCGCAGGGCAGAACGATGTCAAACCACGGCGAACCTTCGCCATAGCGAATTTCGCGATCGCGGCCTGAGccc                                                                                        >  3:691472/1‑149 (MQ=255)
                                                          tctGTTCCAGGCGGTTCAGCACGGCGAGCAGAGGCTGTGCCGAGCGTAGTTTATGGAGCGTCAGCGTGATCCCACCGCCGCAGGGCAGAACGATGTCAAACCACGGCGAACCTTCGCCATAGCGAATTTCGCGATCGCGGCCTGAGccc                                                                                        >  1:8928/1‑149 (MQ=255)
                                                                  aGGCGGTTCAGCACGGCGAGCAGAGGCTGTGCCGAGCGTAGTTTATGGAGCGTCAGCGTGATCCCACCGCCGCAGGGCAGAACGATGTCAAACCACGGCGAACCTTCGCCATAGCGAATTTCGCGATCGCGGCCTGAGCCCATCAtctc                                                                                >  1:235662/1‑149 (MQ=255)
                                                                  aGGCGGTTCAGCACGGCGAGCAGAGGCTGTGCCGAGCGTAGTTTATGGAGCGTCAGCGTGATCCCACCGCCGCAGGGCAGAACGATGTCAAACCACGGCGAACCTTCGCCATAGCGAATTTCGCGATCGCGGCCTGAGCCCATCAtctc                                                                                >  2:483009/1‑149 (MQ=255)
                                                                    gCGGTTCAGCACGGCGAGCAGAGGCTGTGCCGAGCGTAGTTTATGGAGCGTCAGCGTGATCCCACCGCCGCAGGGCAGAACGATGTCAAACCACGGCGAACCTTCGCCATAGCGAATTTCGCGATCGCGGCCTGAGCCCATCATCTCCa                                                                              >  1:618327/1‑149 (MQ=255)
                                                                          cAGCACGGCGAGCAGAGGCTGTGCCGAGCGTAGTTTATGGAGCGTCAGCGTGATCCCACCGCCGCAGGGCAGAACGATGTCAAACCACGGCGAACCTTCGCCATAGCGAATTTCGCGATCGCGGCCTGAGCCCATCATCTCCAGCGCtt                                                                        >  2:104607/1‑149 (MQ=255)
                                                                           aGCACGGCGAGCAGAGGCTGTGCCGAGCGTAGTTTATGGAGCGTCAGCGTGATCCCACCGCCGCAGGGCAGAACGATGTCAAACCACGGCGAACCTTCGCCATAGCGAATTTCGCGATCGCGGCCTGAGCCCATCATCTCCAGCGCTTc                                                                       >  1:78715/1‑149 (MQ=255)
                                                                            gCACGGCGAGCAGAGGCTGTGCCGAGCGTAGTTTATGGAGCGTCAGCGTGATCCCACCGCCGCAGGGCAGAACGATGTCAAACCACGGCGAACCTTCGCCATAGCGAATTTCGCGATCGCGGCCTGAGCCCATCATCTCCAGCGCTTCa                                                                      >  4:487304/1‑149 (MQ=255)
                                                                             cACGGCGAGCAGAGGCTGTGCCGAGCGTAGTTTATGGAGCGTCAGCGTGATCCCACCGCCGCAGGGCAGAACGATGTCAAACCACGGCGAACCTTCGCCATAGCGAATTTCGCGATCGCGGCCTGAGCCCATCATCTCCAGCGCTTCaa                                                                     >  4:470290/1‑149 (MQ=255)
                                                                                ggCGGGCAGAGGCTGTGCCGAGCGTAGTTTATGGAGCGTCAGCGTGATCCCACCGCCGCAGGGCAGAACGATGTCAAACCACGGCGAACCTTCGCCATAGCGAATTTCGCGATCGCGGCCTGAGCCCATCATCTCCAGCGCTTCAAAag                                                                  >  2:702831/1‑149 (MQ=255)
                                                                                ggCGAGCAGAGGCTGTGCCGAGCGTAGTTTATGGAGCGTCAGCGTGATCCCACCGCCGCAGGGCAGAACGATGTCAAACCACGGCGAACCTTCGCCATAGCGAATTTCGCGATCGCGGCCTGAGCCCATCATCTCCAGCGCTTCAAAag                                                                  <  1:647642/149‑1 (MQ=255)
                                                                                ggCGAGCAGAGGCTGTGCCGAGCGTAGTTTATGGAGCGTCAGCGTGATCCCACCGCCGCAGGGCAGAACGATGTCAAACCACGGCGAACCTTCGCCATAGCGAATTTCGCGATCGCGGCCTGAGCCCATCATCTCCAGCGCTTCAAAag                                                                  <  4:536731/149‑1 (MQ=255)
                                                                                 gcgtgtagagGCTGTGCCGAGCGTAGTCTATGGAGCGTCAGCGTGATCCCACCGCCGCAGGGCAGAACGATGTCAAACCACGGCGAACCTTCGCCATAGCGAATTTCGCGATCGCGGCCTGAGCCCATCATCTCCAGCGCTTCAAAagc                                                                 <  4:655500/143‑1 (MQ=255)
                                                                                 gCGAGCAGAGGCTGTGCCGAGCGTAGTTTATGGAGCGTCAGCGTGATCCCACCGCCGCAGGGCAGAACGATGTCAAACCACGGCGAACCTTCGCCATAGCGAATTTCGCGATCGCGGCCTGAGCCCATCATCTCCAGCGCTTCAAAagc                                                                 <  4:117881/149‑1 (MQ=255)
                                                                                         agGCTGTGCCGAGCGTAGTTTATGGAGCGTCAGCGTGATCCCACCGCCGCAGGGCAGAACGATGTCAAACCACGGCGAACCTTCGCCATAGCGAATTTCGCGATCGCGGCCTGAGCCCATCATCTCCAGCGCTTCAAAagcggcagcgg                                                         >  4:123156/1‑149 (MQ=255)
                                                                                          ggCTGTGCCGAGCGTAGTTTATGGAGCGTCAGCGTGATCCCACCGCCGCAGGGCAGAACGATGTCAAACCACGGCGAACCTTCGCCATAGCGAATTTCGCGATCGCGGCCTGAGCCCATCATCTCCAGCGCTTCAAAagcggcagcggc                                                        <  1:406999/149‑1 (MQ=255)
                                                                                           gCTGTGCCGAGCGTAGTTTATGGAGCGTCAGCGTGATCCCACCGCCGCAGGGCAGAACGATGTCAAACCACGGCGAACCTTCGCCATAGCGAATTTCGCGATCGCGGCCTGAGCCCATCATCTCCAGCGCTTCAAAAGCGGCAGCGGcc                                                       >  4:457207/1‑149 (MQ=255)
                                                                                                 ccGAGCGTAGTTTATGGAGCGTCAGCGTGATCCCACCGCCGCAGGGCAGAACGATGTCAAACCACGGCGAACCTTCGCCATAGCGAATTTCGCGATCGCGGCCTGAGCCCATCATCTCCAGCGCTTCAAAAGCGGCAGCGGCCTCCACg                                                 >  3:262904/1‑149 (MQ=255)
                                                                                                   gAGCGTAGTTTATGGAGCGTCAGCGTGATCCCACCGCCGCAGGGCAGAACGATGTCAAACCACGGCGAACCTTCGCCATAGCGAATTTCGCGATCGCGGCCTGAGCCCATCATCTCCAGCGCTTCAAAagcggc                                                              <  3:50745/134‑1 (MQ=255)
                                                                                                   gAGCGTAGTTTATGGAGCGTCAGCGTGATCCCACCGCCGCAGGGCAGAACGATGTCAAACCACGGCGAACCTTCGCCATAGCGAATTTCGCGATCGCGGCCTGAGCCCATCATCTCCAGCGCTTCAAAagcggc                                                              >  4:50745/1‑134 (MQ=255)
                                                                                                               tGGAGCGTCAGCGTGATCCCACCGCCGCAGGGCAGAATGATGTCAAACCACGGCGAACCTTCGCCATAGCGAATTTCGCGATCGCGGCCTGAGCCCATCATCTCCAGCGCTTCAAAAGCGGCAGCGGCCTCCACGCAGCCGCCAGAGAc                                   <  2:182226/149‑1 (MQ=255)
                                                                                                               tGGAGCGTCAGCGTGATCCCACCGCCGCAGGGCAGAATGATGTCAAACCACGGCGAACCTTCGCCATAGCGAATTTCGCGATCGCGGCCTGAGCCCATCATCTCCAGCGCTTCAAAAGCGGCAGCGGCCTCCACGCAGCCGCCAGAGAc                                   <  2:753884/149‑1 (MQ=255)
                                                                                                               tGGAGCGTCAGCGTGATCCCACCGCCGCAGGGCAGAACGATGTCAAACCACGGCGAACCTTCGCCATAGCGAATTTCGCGATCGCGGCCTGAGCCCATCATCTCCAGCGCTTCAAAAGCGGCAGCGGCCTCCACGCAGCCGCCAGAGAc                                   <  4:308526/149‑1 (MQ=255)
                                                                                                                  aGCGTCAGCGTGATCCCACCGCCGCAGGGCAGAACGATGTCAAACCACGGCGAACCTTCGCCATAGCGAATTTCGCGATCGCGGCCTGAGCCCATCATCTCCAGCGCTTCAAAAGCGGCAGCGGCCTCCACGCAGCCGCCAGAGACaaa                                >  1:24691/1‑149 (MQ=255)
                                                                                                                   gCGTCAGCGTGATCCCACCGCCGCAGGGCAGAACGATGTCAAACCACGGCGAACCTTCGCCATAGCGAATTTCGCGATCGCGGCCTGAGCCCATCATCTCCAGCGCTTCAAAAGCGGCAGCGGCCTCCACGCAGCCGCCAGAGACaaaa                               >  4:65826/1‑149 (MQ=255)
                                                                                                                   gCGTCAGCGTGATCCCACCGCCGCAGGGCAGAACGATGTCAAACCACGGCGAACCTTCGCCATAGCGAATTTCGCGATCGCGGCCTGAGCCCATCATCTCCAGCGCTTCAAAAGCGGCAGCGGCCTCCACGCAGCCGCCAGAGACaaaa                               <  1:621358/149‑1 (MQ=255)
                                                                                                                   gCGTCAGCGTGATCCCACCGCCGCAGGGCAGAACGATGTCAAACCACGGCGAACCTTCGCCATAGCGAATTTCGCGATCGCGGCCTGAGCCCATCATCTCCAGCGCTTCAAAACCGGCAGCGGCCTCCACGCAGCCGCCAGAGACaaat                               <  4:593561/149‑2 (MQ=255)
                                                                                                                      tCAGCGTGATCCCACCGCCGCAGGGCAGAACGATGTCAAACCACGGCGAACCTTCGCCATAGCGAATTTCGCGATCGCGGCCTGAGCCCATCATCTCCAGCGCTTCAAAAGCGGCAGCGGCCTCCACGCAGCCGCCAGAGACAAAACCg                            >  3:500108/1‑149 (MQ=255)
                                                                                                                        aGCGTGATCCCACCGCCGCAGGGCAGAACGATGTCAAACCACGGCGAACCTTCGCCATAGCGAATTTCGCGATCGCGGCCTGAGCCCATCATCTCCAGCGCTTCAAAAGCGGCAGCGGCCTCCACGCAGCCGCCAGAGACAAAACCGCa                          <  3:703803/149‑1 (MQ=255)
                                                                                                                        aGCGTGATCCCACCGCCGCAGGGCAGAACGATGTCAAACCACGGCGAACCTTCGCCATAGCGAATTTCGCGATCGCGGCCTGAGCCCATCATCTCCAGCGCTTCAAAAGCGGCAGCGGCCTCCACGCAGCCGCCAGAGACAAAACCGCa                          >  2:495923/1‑149 (MQ=255)
                                                                                                                             gATCCCACCGCCGCAGGGCAGAACGATGTCAAACCACGGCGAACCTTCGCCATAGCGAATTTCGCGATCGCGGCCTGAGCCCATCATCTCCAGCGCTTCAAAAGCGGCAGCGGCCTCCACGCAGCCGCCAGAGACAAAACCGCAGTAAc                     >  4:414587/1‑149 (MQ=255)
                                                                                                                             gATCCCACCGCCGCAGGGCAGAACGATGTCAAACCACGGCGAACCTTCGCCATAGCGAATTTCGCGATCGCGGCCTGAGCCCATCATCTCCAGCGCTTCAAAAGCGGCAGCGGCCTCCACGCAGCCGCCAGAGACAAAACCGCAGTAAc                     <  1:172750/149‑1 (MQ=255)
                                                                                                                               tCCCACCGCCGCAGGGCAGAACGATGTCAAACCACGGCGAACCTTCGCCATAGCGAATTTCGCGATCGCGGCCTGAGCCCATCATCTCCAGCGCTTCAAAAGCGGCAGCGGCCTCCACGCAGCCGCCAGAGACAAAACCGCAGTAACGa                   <  2:368522/149‑1 (MQ=255)
                                                                                                                               tCCCACCGCCGCAGGGCAGAACGATGTCAAACCACGGCGAACCTTCGCCATAGCGAATTTCGCGATCGCGGCCTGAGCCCATCATCTCCAGCGCTTCAAAAGCGGCAGCGGCCTCCACGCAGCCGCCAGAGACAAAACCGCAGTAACGa                   <  1:666200/149‑1 (MQ=255)
                                                                                                                                cccACCGCCGCAGGGCAGAACGATGTCAAACCACGGCGAACCTTCGCCATAGCGAATTTCGCGATCGCGGCCTGAGCCCATCATCTCCAGCGCTTCAAAAGCGGCAGCGGCCTCCACGCAGCCGCCAGAGACAAAACCGCAGTAACGAc                  <  4:262904/149‑1 (MQ=255)
                                                                                                                                 ccACCGCCGCAGGGCAGAACGATGTCAAACCACGGCGAACCTTCGCCATAGCGAATTTCGCGATCGCGGCCTGAGCCCATCATCTCCAGCGCTTCAAAAGCGGCAGCGGCCTCCACGCAGCCGCCAGAGACAAAACCGCAGTAACGAcc                 >  1:414218/1‑149 (MQ=255)
                                                                                                                                 ccACCGCCGCAGGGCAGAACGATGTCAAACCACGGCGAACCTTCGCCATAGCGAATTTCGCGATCGCGGCCTGAGCCCATCATCTCCAGCGCTTCAAAAGCGGCAGCGGCCTCCACGCAGCCGCCAGAGACAAAACCGCAGTAACGAcc                 <  2:284034/149‑1 (MQ=255)
                                                                                                                                 ccACCGCCGCAGGGCAGAACGATGTCAAACCACGGCGAACCTTCGCCATAGCGAATTTCGCGATCGCGGCCTGAGCCCATCATCTCCAGCGCTTCAAAAGCGGCAGCGGCCTCCACGCAGCCGCCAGAGACAAAACCGCAGTAACGAcc                 >  2:217269/1‑149 (MQ=255)
                                                                                                                                  cACCGCCGCAGGGCAGAACGATGTCAAACCACGGCGAACCTTCGCCATAGCGAATTTCGCGATCGCGGCCTGAGCCCATCATCTCCAGCGCTTCAAAAGCGGCAGCGGCCTCCACGCAGCCGCCAGAGACAAAACCGCAGTAACGACCa                >  1:360776/1‑149 (MQ=255)
                                                                                                                                   aCCGCCGCAGGGCAGAACGATGTCAAACCACGGCGAACCTTCGCCATAGCGAATTTCGCGATCGGGGCCTGAGCCCATCATCTCCAGCGCTTCAAAAGCGGCAGCGGCCTCCACGCAGCCGCCAGAGACAAAACCGCAGTAACGACCAt               >  1:528258/1‑149 (MQ=255)
                                                                                                                                                aGAACGATGTCAAACCACGGCGAACCTTCGCCATAGCGAATTTCGCGATCGCGGCCTGAGCCCATCATCTCCAGCGCTTCAAAAGCGGCAGCGGCCTCCACGCAGCCGCCAGAGACAAAACCGCAGTAACGACCATCTTCGCGcaccac  >  4:412475/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                
GTGGGCAGGCACACCAGCGATTGTGCTTGCGGATCGTAGCGCAGCCCCACCGGTTTTCTCTGTTCCAGGCGGTTCAGCACGGCGAGCAGAGGCTGTGCCGAGCGTAGTTTATGGAGCGTCAGCGTGATCCCACCGCCGCAGGGCAGAACGATGTCAAACCACGGCGAACCTTCGCCATAGCGAATTTCGCGATCGCGGCCTGAGCCCATCATCTCCAGCGCTTCAAAAGCGGCAGCGGCCTCCACGCAGCCGCCAGAGACAAAACCGCAGTAACGACCATCTTCGCGCACCAC  >  NC_000913/298251‑298543

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑