Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,388,284 A→T 5.0% 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 5.0% 348.6 / 6.9 99A281A (GCT→GCA) ‡aaeAp‑hydroxybenzoic acid efflux system component
Reads supporting (aligned to +/- strand):  ref base A (50/44);  new base T (2/3);  total (52/47)
Fisher's exact test for biased strand distribution p-value = 6.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.11e-01

GCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATACGAACCGGAACACGTTGCGCAAGACGCACCCATTCAAGGTTAGAGTCTATAGTCGCCATCCCTTTGTCGTCACGCGTGCTGCTGGCGTTGGTGA  >  NC_000913/3388136‑3388429
                                                                                                                                                    |                                                                                                                                                 
gCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTa                                                                                                                                                   <  3:662416/149‑1 (MQ=255)
 ccAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAg                                                                                                                                                  >  4:567411/1‑149 (MQ=255)
   aGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCa                                                                                                                                                <  4:250261/149‑1 (MQ=255)
    gTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAg                                                                                                                                               <  2:376175/149‑1 (MQ=255)
         aCCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGCTACAgtggtg                                                                                                                                          <  1:621016/149‑1 (MQ=255)
         aCCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGCTACAgtggtg                                                                                                                                          <  1:466350/149‑1 (MQ=255)
         aCCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGCTACAgtggtg                                                                                                                                          <  3:128455/149‑1 (MQ=255)
             cAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAACCTCACGCAGGCGATGGGCCATTTACGGAAGAACGAATCCTGGCCTTCGTCGCGATCTTGTTTGCCAGTGACCACTACTGTAGCAGTGGTGCCCGc                                                                                                                                     >  3:605087/1‑149 (MQ=255)
                 gCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAgg                                                                                                                                  >  1:623627/1‑149 (MQ=255)
                  cGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGc                                                                                                                                 >  3:686546/1‑149 (MQ=255)
                              ttAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTctct                                                                                                                     <  3:547582/149‑1 (MQ=255)
                              ttAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTctct                                                                                                                     >  4:153324/1‑149 (MQ=255)
                              ttAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTctct                                                                                                                     >  2:685179/1‑149 (MQ=255)
                                  cAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCt                                                                                                                 <  3:153324/149‑1 (MQ=255)
                                   aaTGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTg                                                                                                                <  1:488751/149‑1 (MQ=255)
                                         gaAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTc                                                                                                          <  3:567411/149‑1 (MQ=255)
                                           aaaTACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCga                                                                                                        <  4:319462/149‑1 (MQ=255)
                                             aTACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCgaga                                                                                                      <  3:250055/149‑1 (MQ=255)
                                             aTACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCgaga                                                                                                      <  4:674326/149‑1 (MQ=255)
                                             aTACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCgaga                                                                                                      <  1:581449/149‑1 (MQ=255)
                                             aTACCCATCGTGATTAAACAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCgaga                                                                                                      <  2:680491/149‑1 (MQ=255)
                                              tACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGAc                                                                                                     >  4:356585/1‑149 (MQ=255)
                                              tACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGAc                                                                                                     <  1:561021/149‑1 (MQ=255)
                                                 ccATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGaa                                                                                                  >  1:8706/1‑149 (MQ=255)
                                                  cATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAAt                                                                                                 <  4:612878/149‑1 (MQ=255)
                                                   aTCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGaa                                                                                                  >  1:731297/1‑147 (MQ=255)
                                                   aTCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGaa                                                                                                  <  2:731297/147‑1 (MQ=255)
                                                   aTCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATa                                                                                                >  4:655388/1‑149 (MQ=255)
                                                   aTCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAAGCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATa                                                                                                <  4:481340/149‑1 (MQ=255)
                                                    tCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATGTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGAt                                                                                                                            >  3:301360/1‑120 (MQ=255)
                                                    tCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGAt                                                                                                                            <  4:301360/120‑1 (MQ=255)
                                                    tCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCgaga                                                                                                      <  2:434060/142‑1 (MQ=255)
                                                    tCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCgaga                                                                                                      >  1:434060/1‑142 (MQ=255)
                                                    tCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATAc                                                                                               <  1:102448/149‑1 (MQ=255)
                                                    tCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATAc                                                                                               >  4:360942/1‑149 (MQ=255)
                                                       tGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATACGaa                                                                                            <  4:308197/149‑1 (MQ=255)
                                                              ccAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATACGAACCGGAac                                                                                     >  2:728684/1‑149 (MQ=255)
                                                                aaaCTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATa                                                                                                >  4:133666/1‑136 (MQ=255)
                                                                aaaCTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATa                                                                                                <  3:133666/136‑1 (MQ=255)
                                                                aaaCTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATACGAACCGGAacac                                                                                   >  4:191546/1‑149 (MQ=255)
                                                                aaaCTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATACGAACCGGAacac                                                                                   >  4:249886/1‑149 (MQ=255)
                                                                   cTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATACGAACCGGAACACGtt                                                                                <  3:408081/149‑1 (MQ=255)
                                                                   cTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTATAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATACGAACCGGAACACGtt                                                                                <  2:107904/149‑1 (MQ=255)
                                                                      aCGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATACGAACCGGAACACGTTgcg                                                                             >  4:221959/1‑149 (MQ=255)
                                                                        gCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATACGAACCGGAACACGTTGCGCa                                                                           >  3:572468/1‑149 (MQ=255)
                                                                           ggCGGTGGGCCATTGTACGGACGAACGAATCCTGCCTTTCGTGGCGATCTCGGTTGCCAGTGACCAACACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGTTTGTCGAGACGAATACGAACCGGAACACGTTGCGCAAgg                                                                        >  2:201836/1‑148 (MQ=255)
                                                                            gCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATACGAACCGGAACACGTTGCGCAAGAc                                                                       >  1:504411/1‑149 (MQ=255)
                                                                                   gCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAg                                                                                                                                   <  1:237497/82‑1 (MQ=255)
                                                                                   gCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAg                                                                                                                                   >  2:237497/1‑82 (MQ=255)
                                                                                    ccATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATACGAACCGGAACACGTTGCGCAAGACGCACCCAt                                                               >  2:46207/1‑149 (MQ=255)
                                                                                    ccATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATACGAACCGGAACACGTTGCGCAAGACGCACCCAt                                                               >  3:577846/1‑149 (MQ=255)
                                                                                    ccATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATACGAACCGGAACACGTTGCGCAAGACGCACCCAt                                                               <  3:693455/149‑1 (MQ=255)
                                                                                      aTTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATACGAACCGGAACACGTTGCGCAAGACGCACCCATTc                                                             >  4:28835/1‑149 (MQ=255)
                                                                                         ttACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATACGAACCGGATCACGTTGCGCAAGACGCACCCATTCAAg                                                          <  3:51550/149‑1 (MQ=255)
                                                                                         ttACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATACGAACCGGAACACGTTGCGCAAGACGCACCCATTCAAg                                                          <  1:284467/149‑1 (MQ=255)
                                                                                          tACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATACGAACCGGAACACGTTGCGCAAGACGCACCCATTCAAgg                                                         >  3:370267/1‑149 (MQ=255)
                                                                                             ggAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATACGAACCGGAACACGTTGCGCAAGACGCACCCATTCAAGGTTa                                                      <  3:250182/149‑1 (MQ=255)
                                                                                             ggAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATACGAACCGGAACACGTTGCGCAAGACGCACCCATTCAAGGTTa                                                      >  4:401320/1‑149 (MQ=255)
                                                                                             ggAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATACGAACCGGAACACGTTGCGCAAGACGCACCCATTCAAGGTTa                                                      >  1:309603/1‑149 (MQ=255)
                                                                                             ggAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATACGAACCGGAACACGTTGCGCAAGACGCACCCATTCAAGGTTa                                                      <  1:728684/149‑1 (MQ=255)
                                                                                             ggAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATACGAACCGGAACACGTTGCGCAAGACGCACCCATTCAAGGTTa                                                      <  3:249886/149‑1 (MQ=255)
                                                                                               aagaaCGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGGCCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATACGAACCGGAACACGTTGCGCAAGACGCACCCATTCAAGGTTaga                                                    <  2:27267/149‑1 (MQ=255)
                                                                                               aagaaCGAATCCTGGCTTTCGTCGCGATCTTGTGTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATACGAACCGGAACACGTTGCGCAAGACGCACCCATTCAAGGTTaga                                                    <  4:605087/149‑1 (MQ=255)
                                                                                                agaaCGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATACGAACCGGAACACGTTGCGCAAGACGCACCCATTCAAGGTTagag                                                   >  3:603537/1‑149 (MQ=255)
                                                                                                agaaCGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTATTGCTGGTTGTCGAGACGAATACGAACCGGAACACGTTTCGCAAGAGGCACCCATTCAAGGTTagag                                                   >  2:17042/1‑149 (MQ=255)
                                                                                                 gaaCGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATACGAACCGGAACACGTTGCGCAAGACGCACCCATTCAAGGTTAGAGt                                                  >  2:230672/1‑149 (MQ=255)
                                                                                                   aCGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATACGAACCGGAACACGTTGCGCAAGACGCACCCATTCAAGGTTAGAGTCt                                                >  1:80690/1‑149 (MQ=255)
                                                                                                          cTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGCGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATACGAACCGGAACACGTTGCGCAAGACGCACCCATTCAAGGTTAGAGTCTATAGTCg                                         >  3:205187/1‑149 (MQ=255)
                                                                                                           tGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATACGAACCGGAACACGTTgcgc                                                                            >  3:311028/1‑113 (MQ=255)
                                                                                                           tGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATACGAACCGGAACACGTTgcgc                                                                            <  4:311028/113‑1 (MQ=255)
                                                                                                           tGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATACGAACCGGAACACGTTGCGCAAGACGCACCCATTCAAGGTTAGAGTCTATAGTCGc                                        >  3:28843/1‑149 (MQ=255)
                                                                                                           tGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATACGAACCGGAACACGTTGCGCAAGACGCACCCATTCAAGGTTAGAGTCTATAGTCGc                                        >  4:435902/1‑149 (MQ=255)
                                                                                                           tGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATACGAACCGGAACACGTTGCGCAAGACGCACCCATTCAAGGTTAGAGTCTATAGTCGc                                        >  2:548278/1‑149 (MQ=255)
                                                                                                           tGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATACGAACCGGAACACGTTGCGCAAGACGCACCCATTCAAGGTTAGAGTCTATAGTCGc                                        >  3:360048/1‑149 (MQ=255)
                                                                                                             gCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATACGAACCGGAACACGTTGCGCAAGACGCACCCATTCAAGCTTAGAGTCTATAGTCGCCa                                      >  4:517637/1‑149 (MQ=255)
                                                                                                               tttCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGa                                                                                                   >  1:569682/1‑86 (MQ=255)
                                                                                                               tttCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGa                                                                                                   <  2:569682/86‑1 (MQ=255)
                                                                                                                ttCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATACGAACCGGAACACGTTGCGCAAGACGCACCCATTCAAGGTTAGAGTCTATAGTCGCCATcc                                   <  4:603537/149‑1 (MQ=255)
                                                                                                                ttCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATACGAACCGGAACACGTTGCGCAAGACGCACCCATTCAAGGTTAGAGTCTATAGTCGCCATcc                                   <  2:587357/149‑1 (MQ=255)
                                                                                                                    tcgCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATACGAACCGGAACACGTTGCGCAAGACGCACCCATTCAAGGTTAGAGTCTATAGTCGCCATCCCttt                               <  4:625739/149‑1 (MQ=255)
                                                                                                                       cgATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATACGAACCGGAACACGTTGCGCAAGACGCACCCATTCAAGGTTAGAGTCTATAGTCGCCATCCCTTTgtc                            <  3:401320/149‑1 (MQ=255)
                                                                                                                        gATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATACGAACCGGAACACGTTGCGCAAGACGCACCCATTCAAGGTTAGAGTCTATAGTCGCCATCCCTTTgtcg                           >  4:220097/1‑149 (MQ=255)
                                                                                                                        gATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATACGAACCGGAACACGTTGCGCAAGACGCACCCATTCAAGGTTAGAGTCTATAGTCGCCATCCCTTTgtcg                           >  2:666314/1‑149 (MQ=255)
                                                                                                                          tCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATACGAACCGGAACACGTTGCGCAAGACGCACCCATTCAAGGTTAGAGTCTATAGTCGCCATCCCTTTgtcgtc                         >  4:486216/1‑149 (MQ=255)
                                                                                                                           cTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATACGAACCGGAACACGTTGCGCAAGACGCACCCATTCAAGGTTAGAGTCTATAGTCGCCATCCCTTTGTCGTCa                        >  3:30206/1‑149 (MQ=255)
                                                                                                                            ttGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATACGAACCGGAACACGTTGCGCAAGACGCACCCATTCAAGGTTAGAGTCTATAGTCGCCATCCCTTTGTCGTCAc                       >  1:136109/1‑149 (MQ=255)
                                                                                                                             atctggCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATACGAACCGGAACACGTTGCGCAAGACGCACCCATTCAAGGTTAGAGGCTATAGTCGCCATCCCTTTGTCGTCAcg                      >  2:705522/6‑149 (MQ=255)
                                                                                                                              gTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATACGAACCGGAACACGTTGCGCAAGACGCACCCATTCAAGGTTAGAGTCTATAGTCGCCATCCCTTTGTCGTCAcgc                     <  2:623627/149‑1 (MQ=255)
                                                                                                                              gTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATACGAACCGGAACACGTTGCGCAAGACGCACCCATTCAAGGTTAGAGTCTATAGTCGCCATCCCTTTGTCGTCAcgc                     <  3:221959/149‑1 (MQ=255)
                                                                                                                                    cAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATACGAACCGGAACACGTTGCGCAAGACGCACCCATTCAAGGTTAGAGTCTATAGTCGCCATCCCTTTGTCGTCACGCGtgctg               >  1:98917/1‑149 (MQ=255)
                                                                                                                                       tGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATACGAACCGGAACACGTTGCGCAAGACGCACCCATTCAAGGTTAGAGTCTATAGTCGCCATCCCTTTGTCGTCACGCGTGctgctg            <  1:76518/149‑1 (MQ=255)
                                                                                                                                       tGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATACGAACCGGAACACGTTGCGCAAGACGCACCCATTCAAGGTTAGAGTCTATAGTCGCCATCCCTTTGTCGTCACGCGTGctgctg            <  4:572468/149‑1 (MQ=255)
                                                                                                                                          gcaccacTGCTACAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATACGAACCGGAACACGTTGCGCAAGACGCACCCATTCAAGGTTAGAGTCTATAGTCGCCATCCCTTTGTCGTCACGCGTGCTGCTGGCg         >  2:676332/2‑149 (MQ=255)
                                                                                                                                           caccacTGCTACAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATACGAACCGGAACACGTTGCGCAAGACGCACCCATTCAAGGTTAGAGTCTATAGTCGCCATCCCTTTGTCGTCACGCGTGCTGCTGGCGt        >  2:348580/1‑149 (MQ=255)
                                                                                                                                            accacTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATACGAACCGGAACACGTTGCGCAAGACGCACCCATTCAAGGTTAGAGTCTATAGTCGCCATCCCTTTGTCGTCACGCGTGCTGCTGGCGtt       >  1:50948/1‑149 (MQ=255)
                                                                                                                                             ccacTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATACGAACCGGAACACGTTGCGCAAGACGCACCCATTCAAGGTTAGAGTCTATAGTCGCCATCCCTTTGTCGTCACGCGTGCTGCTGGCGTTg      <  1:46207/149‑1 (MQ=255)
                                                                                                                                              cacTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATACGAACCGGAACACGTTGCGCAAGACGCACCCATTCAAGGTTAGAGTCTATAGTCGCCATCCCTTTGTCGTCACGCGTGCTGCTGGCGTTgg     <  4:370267/149‑1 (MQ=255)
                                                                                                                                              cacTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATACGAACCGGAACACGTTGCGCAAGACGCACCCATTCAAGGTTAGAGTCTATAGTCGCCATCCCTTTGTCGTCACGCGTGCTGCTGGCGTTgg     >  3:547359/1‑149 (MQ=255)
                                                                                                                                              cacTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATACGAACCGGAACACGTTGCGCAAGACGCACCCATTCAAGGTTAGAGTCTATAGTCGCCATCCCTTTGTCGTCACGCGTGCTGCTGGCGTTgg     >  3:380663/1‑149 (MQ=255)
                                                                                                                                                 tGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATACGAACCGGAACACGTTGCGCAAGACGCACCCATTCAAGGTTAGAGTCTATAGTCGCCATCCCTTTGTCGTCACGCGTGCTGCTGGCGTTGGTGa  <  3:356585/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
GCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATACGAACCGGAACACGTTGCGCAAGACGCACCCATTCAAGGTTAGAGTCTATAGTCGCCATCCCTTTGTCGTCACGCGTGCTGCTGGCGTTGGTGA  >  NC_000913/3388136‑3388429

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: