Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,388,284 A→T 6.9% A281A (GCT→GCA) ‡ aaeA ← p‑hydroxybenzoic acid efflux system component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,388,2840AT 6.9% 294.1 / 8.4 85A281A (GCT→GCA) ‡aaeAp‑hydroxybenzoic acid efflux system component
Reads supporting (aligned to +/- strand):  ref base A (37/42);  new base T (2/4);  total (39/46)
Fisher's exact test for biased strand distribution p-value = 6.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.51e-01

CCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATACGAACCGGAACACGTTGCGCAAGACGCACCCATTCAAGGTTAGAGTCTATAGTCGCCATCCCTTTGTCGTCACGCGTGCTGCTGGCGTTGG  >  NC_000913/3388137‑3388426
                                                                                                                                                    |                                                                                                                                              
ccAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAg                                                                                                                                               >  2:136196/1‑149 (MQ=255)
ccAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAg                                                                                                                                               >  3:101130/1‑149 (MQ=255)
 cAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGc                                                                                                                                              >  1:385722/1‑149 (MQ=255)
  aGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCa                                                                                                                                             <  1:673482/149‑1 (MQ=255)
     tttACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAgtg                                                                                                                                          <  2:668865/149‑1 (MQ=255)
     tttACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAgtg                                                                                                                                          <  4:2391/149‑1 (MQ=255)
        aCCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGCTACAgtggtg                                                                                                                                       <  3:314176/149‑1 (MQ=255)
        aCCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGCTACAgtggtg                                                                                                                                       <  2:536167/149‑1 (MQ=255)
         ccGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGCTACAgtggtgg                                                                                                                                      <  3:460338/149‑2 (MQ=255)
         ccGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGCTACAgtggtgg                                                                                                                                      <  3:505239/149‑2 (MQ=255)
               aGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAg                                                                                                                                <  4:163422/149‑1 (MQ=255)
                  gAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGcc                                                                                                                             >  3:198306/1‑149 (MQ=255)
                  gAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGcc                                                                                                                             >  1:120423/1‑149 (MQ=255)
                   aaTATGTTGGTTAGCAATCGAGAAAATACCCAGCGTGATTAACCAAAGTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAAACCTGGCTTTCGTCGCGGTCTTGTTTGCCAGTGACCACGGCTGTAGGAGTGGTGCCCGCAGGCCa                                                                                                                            <  4:629944/149‑1 (MQ=255)
                     tatGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGa                                                                                                                          <  4:198306/149‑1 (MQ=255)
                       tGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATg                                                                                                                        >  4:263087/1‑149 (MQ=255)
                               aGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTg                                                                                                                <  2:345384/149‑1 (MQ=255)
                                    tGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGt                                                                                                           >  3:131314/1‑149 (MQ=255)
                                     ggAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGtt                                                                                                          <  2:435154/149‑1 (MQ=255)
                                     ggAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGtt                                                                                                          <  1:823802/149‑1 (MQ=255)
                                      gagaAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTg                                                                                                         >  2:62852/1‑149 (MQ=255)
                                        gaAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTc                                                                                                       <  3:231019/149‑1 (MQ=255)
                                          aaaTACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCga                                                                                                     <  3:630474/149‑1 (MQ=255)
                                             tACCCACCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGAc                                                                                                  <  2:390602/149‑1 (MQ=255)
                                              aCCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACg                                                                                                 >  4:328826/1‑149 (MQ=255)
                                                ccATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGaa                                                                                               <  2:747077/149‑1 (MQ=255)
                                                ccATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGaa                                                                                               >  1:650375/1‑149 (MQ=255)
                                                 cATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTCCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAAt                                                                                              >  1:160553/1‑149 (MQ=255)
                                                   tCGTGCTTAACCAAACTCACGCAGGCGATGAGCCATGTTACGGAAGCACGAATCCTTCCTTTCGCCGCGATCTTGTTTGTCAGTGATAACCACTGTAGCAGGGGCGTTCGCAAGCCAGATGTTCTCTTGCTCGTTGTCGGGACGAATAc                                                                                            >  4:837128/1‑149 (MQ=255)
                                                   tCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGTTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATAc                                                                                            <  4:694844/149‑1 (MQ=255)
                                                            aCCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATACGAACCGGaa                                                                                   >  2:141970/1‑149 (MQ=255)
                                                             ccAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATACGAACCGGAac                                                                                  <  2:264093/149‑1 (MQ=255)
                                                             ccAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATACGAACCGGAac                                                                                  >  4:200050/1‑149 (MQ=255)
                                                                    cACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATACGAACCGGAACACGTTgc                                                                           <  4:77408/149‑1 (MQ=255)
                                                                        cAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATACGAACCGGAACACGTTGCGCaa                                                                       <  4:640389/149‑1 (MQ=255)
                                                                         aGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATACGAAc                                                                                        <  1:285746/131‑1 (MQ=255)
                                                                         aGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATACGAACCGGAACACGTTGCGCAAg                                                                      >  1:428908/1‑149 (MQ=255)
                                                                         aGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTAGTTGTCGAGACGAATACGAAc                                                                                        >  2:285746/1‑131 (MQ=255)
                                                                          ggCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATACGAACCGGAACACGTTGCGCAAGa                                                                     >  1:242529/1‑149 (MQ=255)
                                                                             gATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATACGAACCGGAACACGTTGCGCAAGACGc                                                                  >  1:565070/1‑149 (MQ=255)
                                                                             gATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATACGAACCGGAACACGTTGCGCAAGACGc                                                                  >  4:445919/1‑149 (MQ=255)
                                                                                  gCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGa                                                                                                                          <  4:222844/88‑1 (MQ=255)
                                                                                  gCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGa                                                                                                                          >  3:222844/1‑88 (MQ=255)
                                                                                      ttttACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATACGAACCGGAACACGTTGCGCAAGACGCACCCATTCa                                                         >  1:512692/1‑149 (MQ=255)
                                                                                      ttttACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATACGAACCGGAACACGTTGCGCAAGACGCACCCATTCa                                                         <  3:285837/149‑1 (MQ=255)
                                                                                       tttACGGAAGCACGAATCCTGGCTTTTGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATACGAACCGGAACACGTTGCGCAAGACGCACCCATTCaa                                                        <  1:641894/149‑1 (MQ=255)
                                                                                        ttACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATACGAACCGGAACACGTTGCGCAAGTCGCACCCATTCAAg                                                       <  1:671013/149‑1 (MQ=255)
                                                                                        ttACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATACGAACCGGAACACGTTGCGCAAGACGCACCCATTCAAg                                                       <  3:121885/149‑1 (MQ=255)
                                                                                        ttACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATACGAACCGGAACACGTTGCGCAAGACGCACCCATTCAAg                                                       <  2:120423/149‑1 (MQ=255)
                                                                                            ggAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATACGAACCGGAACACGTTGCGCAAGACGCACCCATTCAAGGTTa                                                   >  3:86057/1‑149 (MQ=255)
                                                                                            ggAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATACGAACCGGAACACGTTGCGCAAGACGCACCCATTCAAGGTTa                                                   >  3:477172/1‑149 (MQ=255)
                                                                                              acgaaCGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATACGAACCGGAACACGTTGCGCAAGACGCACCCATTCAAGGTTaga                                                 <  1:436350/147‑1 (MQ=255)
                                                                                              aagaaCGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATACGAACCGGAACACGTTGCGCAAGACGCACCCATTCAAGGTTaga                                                 <  2:497042/149‑1 (MQ=255)
                                                                                              aagaaCGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATACGAACCGGAACACGTTGCGCAAGACGCACCCATTCAAGGTTaga                                                 <  4:532248/149‑1 (MQ=255)
                                                                                                gaaCGAATCCTTGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATACGAACCGGAACACGTTGCGCAAGACGCACCCATTCAAGGTTAGAGt                                               >  3:704957/1‑149 (MQ=255)
                                                                                                gaaCGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATACGAACCGGAACACGTTGCGCAAGACGCACCCATTCAAGGTTAGAGt                                               >  1:314412/1‑149 (MQ=255)
                                                                                                     aaTCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATACGAACCGGAACACGTTGCGCAAGACGCACCCATTCAAGGTTAGAGTCtata                                          <  3:660407/149‑1 (MQ=255)
                                                                                                     aaTCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATACGAACCGGAACACGTTGCGCAAGACGCACCCATTCAAGGTTAGAGTCtata                                          <  2:650375/149‑1 (MQ=255)
                                                                                                       tCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGaa                                                                                               >  2:691267/1‑94 (MQ=255)
                                                                                                       tCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGaa                                                                                               <  1:691267/94‑1 (MQ=255)
                                                                                                       tCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATACGAACCGGAACACGTTGCGCAAGACGCACCCATTCAAGGTTAGAGTCTATAGt                                        >  4:661812/1‑149 (MQ=255)
                                                                                                        ccTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATACGAACCGGAACACGTTGCGCAAGACGCACCCATTCAAGGTTAGAGTCTATAGTc                                       <  3:661812/149‑1 (MQ=255)
                                                                                                        ccTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATACGAACCGGAACACGTTGCGCAAGACGCACCCATTCAAGGTTAGAGTCTATAGTc                                       <  3:263087/149‑1 (MQ=255)
                                                                                                              tttCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTtc                                                                                                                     <  3:158055/65‑1 (MQ=255)
                                                                                                              tttCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTtc                                                                                                                     >  4:158055/1‑65 (MQ=255)
                                                                                                              tttCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATACGAACCGGAACACGTTGCGCAAGACGCACCCATTCAAGGTTAGAGTCTATAGTCGCCATc                                 >  4:510213/1‑149 (MQ=255)
                                                                                                              tttCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATACGAACCGGAACACGTTGCGCAAGACGCACCCATTCAAGGTTAGAGTCTATAGTCGCCATc                                 >  1:843602/1‑149 (MQ=255)
                                                                                                                  gtcgCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATACGAACCGGAACACGTTGCGCAAGACGCACCCATTCAAGGTTAGAGTCTATAGTCGCCATCCCtt                             <  3:445919/149‑1 (MQ=255)
                                                                                                                  gtcgCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATACGAACCGGAACACGTTGCGCAAGACGCACCCATTCAAGGTTAGAGTCTATAGTCGCCATCCCtt                             <  4:337993/149‑1 (MQ=255)
                                                                                                                   tcgCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATACGAACCGGAACACGTTGCGCAAGACGCACCCATTCAAGGTTAGAGTCTATAGTCGCCATCCCttt                            <  1:578310/149‑1 (MQ=255)
                                                                                                                   tcgCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATACGAACCGGAACACGTTGCGCAAGACGCACCCATTCAAGGTTAGAGTCTATAGTCGCCATCCCttt                            <  4:392484/149‑1 (MQ=255)
                                                                                                                   tcgCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATACGAACCGGAACACGTTGCGCAAGACGCACCCATTCAAGGTTAGAGTCTATAGTCGCCATCCCttt                            <  4:86057/149‑1 (MQ=255)
                                                                                                                      cgATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATACGAACCGGAACACGTTGCGCAAGACGCACCCATTCAAGGTTAGAGTCTATAGTCGCCATCCCTTTgtc                         <  4:84543/149‑1 (MQ=255)
                                                                                                                      cgATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATACGAACCGGAACACGTTGCGCAAGACGCACCCATTCAAGGTTAGAGTCTATAGTCGCCATCCCTTTgtc                         >  1:271110/1‑149 (MQ=255)
                                                                                                                       gATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATACGAACCGGAACACGTTGCGCAAGACGCACCCATTCAAGGTTAGAGTCTATAGTCGCCATCCCTTTgtcg                        >  4:531804/1‑149 (MQ=255)
                                                                                                                        aTCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATACGAACCGGAACACGTTGCGCAAGACGCACCCATTCAAGGTTAGAGTCTATAGTCGCCATCCCTTTgtcgt                       >  4:691566/1‑149 (MQ=255)
                                                                                                                         tCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATACGAACCGGAACACGTTGCGCAAGACGCACCCATTCAAGGTTAGAGTCTATAGTCGCCATCCCTTTgtcgtc                      <  3:604317/149‑1 (MQ=255)
                                                                                                                                 gCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATACGAACCGGAACACGTTGCGCAAGACGCACCCATTCAAGGTTAGAGTCTATAGTCGCCATCCCTTTGTCGTCACGCGtgc              >  1:823709/1‑149 (MQ=255)
                                                                                                                                        ggcaccacTGCTACAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATACGAACCGGAACACGTTGCGCAAGACGCACCCATTCAAGGTTAGAGTCTATAGTCGCCATCCCTTTGTCGTCACGCGTGCTGCTGGc       >  1:862796/3‑149 (MQ=255)
                                                                                                                                        ggcaccacTGCTACAGTGGTGCCCGCAGGCCAGAGGTTCTCTTGCTGGTTGGCGAGACGAATACGAACCGGAACACGTTGCGCAAGACGCACCCATTCAAGGTTAGAGTCTATAGTCGCCATCCCTTTGTCGTCACGCGTGCTGCTGGc       >  2:703992/3‑149 (MQ=255)
                                                                                                                                          caccacTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATACGAACCGGAACACGTTGCGCAAGACGCACCCATTCAAGGTTAGAGTCTATAGTCGCCATCCCTTTGTCGTCACGCGTGCTGCTGGCGt     >  4:398185/1‑149 (MQ=255)
                                                                                                                                           accacTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATACGAACCGGAACACGTTGCGCAAGACGCACCCATTCAAGGTTAGAGTCTATAGTCGCCATCCCTTTGTCGTCACGCGTGCTGCTGGCGtt    >  3:59428/1‑149 (MQ=255)
                                                                                                                                           accacTGCTACAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATACGAACCGGAACACGTTGCGCAAGACGCACCCATTCAAGGTTAGAGTCTATAGTCGCCATCCCTTTGTCGTCACGCGTGCTGCTGGCGtt    <  2:862796/149‑1 (MQ=255)
                                                                                                                                            ccacTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATACGAACCGGAACACGTTGCGCAAGACGCACCCATTCAAGGTTAGAGTCTATAGTCGCCATCCCTTTGTCGTCACGCGTGCTGCTGGCGTTg   <  4:450118/149‑1 (MQ=255)
                                                                                                                                             cacTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATACGAACCGGAACACGTTGCGCAAGACGCACCCATTCAAGGTTAGAGTCTATAGTCGCCATCCCTTTGTCGTCACGCGTGCTGCTGGCGTTgg  >  3:210439/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                              
CCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATACGAACCGGAACACGTTGCGCAAGACGCACCCATTCAAGGTTAGAGTCTATAGTCGCCATCCCTTTGTCGTCACGCGTGCTGCTGGCGTTGG  >  NC_000913/3388137‑3388426

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