Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,638,584 A→T 8.3% F23L (TTT→TTA)  bamB ← BamABCDE complex OM biogenesis lipoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,638,5840AT 8.3% 283.4 / 13.8 84F23L (TTT→TTA) bamBBamABCDE complex OM biogenesis lipoprotein
Reads supporting (aligned to +/- strand):  ref base A (37/40);  new base T (4/3);  total (41/43)
Fisher's exact test for biased strand distribution p-value = 7.10e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.99e-01

GACAACGTTGTCCGCCAGTGCCGGATGAAGATTGGAATAGAAGTTGCCAATGCCGCTACCAACGGAAGTGCTCCACGCCGTGGTCGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTAATTTACGCAATTGCATCGGGTCCCTCTCAGATGGACAAATTATTAATTTTCATCTGCATCATTTCGCTCAGTGCCGGAGTAACATCGCTTTTCACG  >  NC_000913/2638437‑2638732
                                                                                                                                                   |                                                                                                                                                    
gACAACGTTGTCCGCCAGTGCCGGATGAAGATTGGAATAGAAGTTGCCAATGCCGCTACCAACGGAAGTGCTCCACGCCGTGGTCGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTaa                                                                                                                                                     >  3:237660/1‑149 (MQ=255)
 aCAACGTTGTCCGCCAGTGCCGGATGAAGATTGGAATAGAAGTTGCCAATGCCGCTACCAACGGAAGTGCTCCACTCCGTGGTCGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTaaa                                                                                                                                                    <  4:529576/149‑1 (MQ=255)
 aCAACGTTGTCCGCCAGTGCCGGATGAAGATTGGAATAGAAGTTGCCAATGCCGCTACCAACGGAAGTGCTCCACGCCGTGGTCGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTaaa                                                                                                                                                    <  4:260915/149‑1 (MQ=255)
 aCAACGTTGTCCGCCAGTGCCGGATGAAGATTGGAATAGAAGTTGCCAATGCCGCTACCAACGGAAGTGCTCCACGCCGTGGTCGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTaaa                                                                                                                                                    >  1:284213/1‑149 (MQ=255)
 aCAACGTTGTCCGCCAGTGCCGGATGAAGATTAGAATAGAAGTTGCCAATGCCGCTACCAACGGAAGTGCTCCACGCCGTGGTCGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTaaa                                                                                                                                                    <  1:646874/149‑1 (MQ=255)
   aaCGTTGTCCGCCAGTGCCGGATGAAGATTGGAATAGAAGTTGCCAATGCCGCTACCAACGGAAGTGCTCCACGCCGTGGTCGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTTAACa                                                                                                                                                  >  2:238714/1‑149 (MQ=255)
   aaCGTTGTCCGCCAGTGCCGGATGAAGATTGGAATAGAAGTTGCCAATGCCGCTACCAACGGAAGTGCTCCACGCCGTGGTCGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACa                                                                                                                                                  >  3:33509/1‑149 (MQ=255)
    aCGTTGTCCGCCAGTGCCGGATGAAGATTGGAATAGAAGTTGCCAATGCCGCTACCAACGGAAGTGCTCCACGCCGTGGTCGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAg                                                                                                                                                 >  1:382822/1‑149 (MQ=255)
        tGTCCGCCAGTGCCGGATGAAGATTGGAATAGAAGTTGCCAATGCCGCTACCAACGGAAGTGCTCCACGCCGTGGTCGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTTAACAGCGaa                                                                                                                                             <  1:238714/149‑1 (MQ=255)
         gTCCGCCAGTGCCGGATGAAGATTGGAATAGAAGTTGCCAATGCCGCTACCAACGGAAGTGCTCCACGCCGTGGTCGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGGCCATCTTTACCACATCTTCTTCGCTGTTTAACAGCGaag                                                                                                                                            <  1:179337/149‑2 (MQ=255)
          tCCGCCAGTGCCGGATGAAGATTGGAATAGAAGTTGCCAATGCCGCTACCAACGGAAGTGCTCCACGCCGTGGTCGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTaaa                                                                                                                                                    <  3:138066/140‑1 (MQ=255)
          tCCGCCAGTGCCGGATGAAGATTGGAATAGAAGTTGCCAATGCCGCTACCAACGGAAGTGCTCCACGCCGTGGTCGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTaaa                                                                                                                                                    >  4:138066/1‑140 (MQ=255)
           ccGCCAGTGCCGGATGAAGATTGGAATAGAAGTTGCCAATGCCGCTACCAACGGAAGTGCTCCACGCCGTGGTCGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTTAACAGCGaagaa                                                                                                                                          <  1:289088/149‑4 (MQ=255)
             gCCAGTGCCGGATGAAGATTGGAATAGAAGTTGCCAATGCCGCTACCAACGGAAGTGCTCCACGCCGTGGTCGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGCACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGcc                                                                                                                                        >  3:67902/1‑149 (MQ=255)
              ccAGTGCCGGATGAAGATTGGAATAGAAGTTGCCAATGCCGCTACCAACGGAAGTGCTCCACGCCGTGGTCGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCg                                                                                                                                       <  2:428204/149‑1 (MQ=255)
                aGTGCCGGATGAAGTTTGGACTAGAAGTCGCCAATGCCGCTACCAACGGAAGTGCTCCACGCCGTGGTCGGCGTAACCTGGTTTTCAACGGTTGCCAATGGGGACATCTTTACCACCTCTTCTTCGCCGTTAAACAGCGCAAAGGCGCt                                                                                                                                     <  2:763589/149‑1 (MQ=255)
                    ccGGATGAAGATTGGAATAGAAGTTGCCAATGCCGCTACCAACGGAAGTGCTCCACGCCGTGGTCGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTaaa                                                                                                                                 >  1:347556/1‑149 (MQ=255)
                    ccGGATGAAGATTGGAATAGAAGTTGCCAATGCCGCTACCAACGGAAGTGCTCCACGCCGTGGTCGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTaaa                                                                                                                                 >  1:532714/1‑149 (MQ=255)
                     cGGATGAAGATTGGAATAGAAGTTGCCAATGCCGCTACCAACGGAAGTGCTCCACGCCGTGGTCGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTaaaa                                                                                                                                <  2:187507/149‑1 (MQ=255)
                       gATGAAGATTGGAATAGAAGTTGCCAATGCCGCTACCAACGGAAGTGCTCCACGCCGTGGTCGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAgg                                                                                                                              <  1:751147/149‑1 (MQ=255)
                                tGGAATAGAAGTTGCCAATGCCGCTACCAACGGAAGTGCTCCACGCCGTGGTCGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGaa                                                                                                                     <  4:67902/149‑1 (MQ=255)
                                  gAATAGAAGTTGCCAATGCCGCTACCAACGGAAGTGCTCCACGCCGTGGTCGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAg                                                                                                                   <  4:194978/149‑1 (MQ=255)
                                    aTAGAAGTTGCCAATGCCGCTACCAACGGAAGTGCTCCACGCCGTGGTCGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCa                                                                                                                 <  3:251080/149‑1 (MQ=255)
                                      aGAAGTTGCCAATGCCGCTACCAACGGAAGTGCTCCACGCCGTGGTCGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGt                                                                                                               <  3:74473/149‑1 (MQ=255)
                                       gAAGTTGCCAATGCCGCTACCAACGGAAGTGCTCCACGCCGTGGTCGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTc                                                                                                              <  2:151119/149‑1 (MQ=255)
                                          gTTGCCAATGCCGCTACCAACGGAAGTGCTCCACGCCGTGGTCGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTg                                                                                                           <  2:382822/149‑1 (MQ=255)
                                          gTTGCCAATGCCGCTACCAACGGAAGTGCTCCACGCCGTGGTCGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTg                                                                                                           >  1:256025/1‑149 (MQ=255)
                                           ttGCCAATGCCGCTACCAACGGAAGTGCTCCACGCCGTTGTCGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTCAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTgg                                                                                                          <  2:87546/149‑1 (MQ=255)
                                           ttGCCAATGCCGCTACCAACGGAAGTGCTCCACGCCGTGGTCGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTgg                                                                                                          <  4:33509/149‑1 (MQ=255)
                                            tGCCAATGCCGCTACCAACGGAAGTGCTCCACGCCGTGGTCGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGgc                                                                                                         <  3:729949/149‑1 (MQ=255)
                                              ccAATGCCGCTACCAACGGAAGTGCTCCACGCCGTGGTCGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGgcag                                                                                                       <  2:347556/149‑1 (MQ=255)
                                                aaTGCCGCTACCAACGGAAGTGCTCCACGCCGTGGTCGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGgcagca                                                                                                     <  3:569719/149‑1 (MQ=255)
                                                            aaCGGGAGTGCTCCACGCCGTGGTCGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTAATTTACGCa                                                                                         >  2:601414/1‑149 (MQ=255)
                                                               ggAAGTGCTCCACGCCGTGGTCGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTAATTTACGCAAtt                                                                                      <  1:519593/149‑1 (MQ=255)
                                                                    tGCTCCACGCCGTGGTCGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTAATTTACGCAATTGCATc                                                                                 >  3:196690/1‑149 (MQ=255)
                                                                       tCCACGCCGTGGTCGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTAATTTACGCAATTGCATCggg                                                                              >  3:754637/1‑149 (MQ=255)
                                                                       tCCACGCCGTGGTCGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTAATTTACGCAATTGCATCggg                                                                              <  3:663505/149‑1 (MQ=255)
                                                                        ccACGCCGTGGTCGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTAATTTACGCAATTGCATCGGGt                                                                             >  1:669949/1‑149 (MQ=255)
                                                                         cACGCCGTGGTCGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTAATTTACGCAATTGCATCGGGTc                                                                            >  1:684385/1‑149 (MQ=255)
                                                                         cACGCCGTGGTCGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTAATTTACGCAATTGCATCGGGTc                                                                            <  2:523882/149‑1 (MQ=255)
                                                                              cGTGGTCGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTAATTTACGCAATTGCATCGGGTCCctct                                                                       >  1:330529/1‑149 (MQ=255)
                                                                                  gTCGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTAATTTAc                                                                                            <  1:44877/124‑1 (MQ=255)
                                                                                  gTCGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTAATTTAc                                                                                            >  2:44877/1‑124 (MQ=255)
                                                                                   tCGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTAATTTGCGCAATTGCATCGGGTCCCTCTCAGAt                                                                  <  4:713200/149‑1 (MQ=255)
                                                                                      gCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTAATTTACGCAATTGCATCGGGTCCCTCTCAGATGGa                                                               <  4:237660/149‑1 (MQ=255)
                                                                                         tAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTAATTTACGCAATTGCATCGGGTCCCTCTCAGATGGACaa                                                            <  2:532714/149‑1 (MQ=255)
                                                                                          aaaCTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTAATGTACGCAATTGCATCGGGTCCCTCTCAGACGGACaaa                                                           <  2:420329/149‑1 (MQ=255)
                                                                                           aaCTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAgg                                                                                                                              <  3:263984/81‑1 (MQ=255)
                                                                                           aaCTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAgg                                                                                                                              >  4:263984/1‑81 (MQ=255)
                                                                                            aCTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAAc                                                                                                                         <  2:104918/85‑1 (MQ=255)
                                                                                            aCTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAAc                                                                                                                         >  1:104918/1‑85 (MQ=255)
                                                                                             cTGGTTTTCGACGGTTGGCAAAGGGGACATCTTTACCACACCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTAACTTACGCAATTGCATCGGGTCCCTCTCAGATGGCCAAatta                                                        <  2:636814/149‑1 (MQ=255)
                                                                                              tGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTAATTTACGCAATTGCATCGGGTcc                                                                           <  2:458671/129‑1 (MQ=255)
                                                                                              tGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTAATTTACGCAATTGCATCGGGTcc                                                                           >  1:458671/1‑129 (MQ=255)
                                                                                              tGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTAATTTACGCAATTGCATCGGGTCCCTCTCAGATGGACAAattat                                                       <  1:93134/149‑1 (MQ=255)
                                                                                                gTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTAATTTACGCAATTGCATCGGGTCCCTCTCAGATGGACAAAttatta                                                     >  4:621304/1‑149 (MQ=255)
                                                                                                 ttttCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTAATTTACGCAATTGCATCGGGTCCCTCTCAGATGGACAAATTATTaa                                                    >  2:639744/1‑149 (MQ=255)
                                                                                                      aaCGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTAATTTACGCAATTGCATCGGGTCCCTCTCAGATGGACAAATTATTAATTTTc                                               >  2:573670/1‑149 (MQ=255)
                                                                                                      aaCGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTAATTTACGCAATTGCATCGGGTCCCTCTCAGATGGACAAATTATTAATTTTc                                               >  3:721621/1‑149 (MQ=255)
                                                                                                         ggTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTAATTTACGCAATTGCATCGGGTCCCTCTCAGATGGACAAATTATTAATTTTCATc                                            >  3:703854/1‑149 (MQ=255)
                                                                                                           ttGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTAATTTACGCAATTGCATCGGGTCCCTCTCAGATGGACAAATTATTAATTTTCATCTg                                          >  1:526763/1‑149 (MQ=255)
                                                                                                            tGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTAATTTACGCAATTGCATCGGGTCCCTCTCAGATGGACAAATTATTAATTTTCATCTGc                                         <  4:161656/149‑1 (MQ=255)
                                                                                                                aaTGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTAATTTACGCAATTGCATCGGGTCCCTCTCAGATGGACAAATTATTAATTTTCATCTGcatca                                     <  1:119456/149‑1 (MQ=255)
                                                                                                                   ggggACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTAATTTACGCAATTGCACCGGGTCCCTCTCAGATGGACAAATTATTAATTTTCATCTGCATCAttt                                  >  1:640188/1‑149 (MQ=255)
                                                                                                                     ggACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTAATTTACGCAATTGCATCGGGTCCCTCTCAGATGGACAAATTATTAATTTTCATCTGCATCATTTCg                                >  2:577435/1‑149 (MQ=255)
                                                                                                                      gACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTAATTTACGCAATTGCATCGGGTCCCTCTCAGATGGACAAATTATTAATTTTCATCTGCATCATTTCGc                               <  1:760032/149‑1 (MQ=255)
                                                                                                                        cATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTAGCAGCAGTAATTTACGCAATTGCATCGGGTCCCTCTCAGATGGACAAATTATTAATTTTCACCTGCATCATTTCGCTc                             >  1:597517/1‑149 (MQ=255)
                                                                                                                           cTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTAATTTACGCAATTGCATCGGGTCCCTCTCAGATGGACAAATTATTAATTTTCATCTGCATCATTTCGCTCAGt                          <  1:36481/149‑1 (MQ=255)
                                                                                                                            tttACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTAATTTACGCAATTGCATCGGGTCCCTCTCAGATGGACAAATTATTAATTTTCATCTGCATCATTTCGCTCAGTg                         >  4:707063/1‑149 (MQ=255)
                                                                                                                                 cacaTCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTAATTTACGCAATTGCATCGGGTCCCTCTCAGATGGACAAATTATTAATTTTCATCTGCATCATTTCGCTCAGTGCCGGa                    <  1:149509/149‑1 (MQ=255)
                                                                                                                                   caTCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTAATTTACGCAATTGCATCGGGTCCCTCTCAGATGGACAAATTATTAATTTTCATCTGCATCATTTCGCTCAGTGCCGGAGt                  >  4:272310/1‑149 (MQ=255)
                                                                                                                                    aTCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAGAGGGTAAGGGAAAGCAGTCCTGGCAGCAGTAATTTACGCAATTGCATCGGCTCCCTCTCAGATGGACAAATGATTAATTGTCATCTGCATCATTTCGCTCAGTGCCGGAGTa                 >  3:472365/1‑149 (MQ=255)
                                                                                                                                       ttcttcGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTAATTTACGCAATTGCATCGGGt                                                                             >  1:695545/1‑86 (MQ=255)
                                                                                                                                       ttcttcGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTAATTTACGCAATTGCATCGGGt                                                                             >  1:161381/1‑86 (MQ=255)
                                                                                                                                       ttcttcGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTAATTTACGCAATTGCATCGGGt                                                                             <  2:695545/86‑1 (MQ=255)
                                                                                                                                       ttcttcGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTAATTTACGCAATTGCATCGGGt                                                                             <  2:161381/86‑1 (MQ=255)
                                                                                                                                       ttcttcGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTAATTTACGCAATTGCATCGGGTCCCTCTCAGATGGACAAATTATTAAtt                                                  >  1:656854/1‑113 (MQ=255)
                                                                                                                                       ttcttcGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTAATTTACGCAATTGCATCGGGTCCCTCTCAGATGGACAAATTATTAATTTTCATCTGCATCATTTCGCTCAGTGCCGGAGTAACa              >  4:720408/1‑149 (MQ=255)
                                                                                                                                       ttcttcGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTAATTTACGCAATTGCATCGGGTCCCTCTCAGATGGACAAATTATTAATTTTCATCTGCATCATTTCGCTCAGTGCCGGAGTAACa              >  3:501947/1‑149 (MQ=255)
                                                                                                                                       ttcttcGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTAACTTACGCAATTGCATCGGGTCCCTCTCAGATGGACAAATTATTAAtg                                                  <  2:656854/113‑2 (MQ=255)
                                                                                                                                        tgttcGCTGTTTAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTAATTTACGCAATTGCATCGGGTCCCTCTCAGATGGACAAATTATTAATTTTCATCTGCATCATTTCGCTCAGTGCCGGAGTAACAt             >  1:42671/3‑149 (MQ=255)
                                                                                                                                        tgttcGCTGTTTAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTAATTTACGCAATTGCATCGGGTCCCTCTCAGATGGACAAATTATTAATTTTCATCTGCATCATTTCGCTCAGTGCCGGAGTAACAt             >  2:139046/3‑149 (MQ=255)
                                                                                                                                        tcttcGCTGTTAAACAGCGAACAGTCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTAATTTACGCAATTGCATCGGGTCCCTCTCAGATGGACAAATTATTAATTTTCATCTGCATCATTTCGCTCAGTGTCGGAGTAACAt             >  3:657826/1‑149 (MQ=255)
                                                                                                                                          ttcGCTGTTTAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTAATTTACGCAATTGCATCGGGTCCCTCTCAGATGGACAAATTATTAATTTTCATCTGCATCATTTCGCTCAGTGCCGGAGTAACATCg           >  1:328987/1‑149 (MQ=255)
                                                                                                                                           tcGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTAATTTACGCAATTGCATCGGGTCCCTCTCAGATGGACAAATTATTAATTTTCATCTGCATCATTTCGCTCAGTGCCGGAGTAACATCGc          <  2:256025/149‑1 (MQ=255)
                                                                                                                                             gCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTAATTTACGCAATTGCATCGGGTCCCTCTCAGATGGACAAATTATTAATTTTCATCTGCATCATTTCGCTCAGTGCCGGAGTAACATCGCtt        >  1:93861/1‑149 (MQ=255)
                                                                                                                                               tGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTAATTTACGCAATTGCATCGGGTCCCTCTCAGATGGACAAATTATTAATTTTCATCTGCATCATTTCGCTCAGTGCCGGAGTAACATCGCtttt      >  4:41197/1‑149 (MQ=255)
                                                                                                                                               tGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTAATTTACGCAATTGCATCGGGTCCCTCTCAGATGGACAAATTATTAATTTTCATCTGCATCATTTCGCTCAGTGCCGGAGTAACATCGCtttt      >  2:176041/1‑149 (MQ=255)
                                                                                                                                               tGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTAATTCACGCAATTGCATCGGGTCCCTCTCAGATGGACAAATTATTAATTTTCATCTGCATCATTTCGCTCAGTGCCGGAGTAACATCGCtttt      <  3:525219/149‑1 (MQ=255)
                                                                                                                                                gTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTAATTTACGCAATTGCATCGGGTCCCTCTCAGATGGACAAATTATTAATTTTCATCTGCATCATTTCGCTCAGTGCCGGAGTAACATCGCTTTTc     <  2:545956/149‑1 (MQ=255)
                                                                                                                                                  tAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTAATTTACGCAATTGCATCGGGTCCCTCTCATATGGACAAATTATTAATTTTCATCTGCATCATTTCGCTCAGTGCCGGAGTAACATCGCTTTTCAc   <  4:437596/149‑1 (MQ=255)
                                                                                                                                                   aaaCAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTAATTTACGCAATTGCATCGGGTCCCTCTCAGATGGACAAATTATTAATTTTCATCTGCATCATTTCGCTCAGTGCCGGAGTAACATCGCTTTTCACg  >  1:488637/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
GACAACGTTGTCCGCCAGTGCCGGATGAAGATTGGAATAGAAGTTGCCAATGCCGCTACCAACGGAAGTGCTCCACGCCGTGGTCGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTAATTTACGCAATTGCATCGGGTCCCTCTCAGATGGACAAATTATTAATTTTCATCTGCATCATTTCGCTCAGTGCCGGAGTAACATCGCTTTTCACG  >  NC_000913/2638437‑2638732

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