Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,638,584 A→T 5.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 5.3% 391.8 / 7.5 111F23L (TTT→TTAbamBBamABCDE complex OM biogenesis lipoprotein
Reads supporting (aligned to +/- strand):  ref base A (53/52);  new base T (4/2);  total (57/54)
Fisher's exact test for biased strand distribution p-value = 6.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.57e-01

Only 100 of 116 total aligned reads displayed.

GACAACGTTGTCCGCCAGTGCCGGATGAAGATTGGAATAGAAGTTGCCAATGCCGCTACCAACGGAAGTGCTCCACGCCGTGGTCGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTAATTTACGCAATTGCATCGGGTCCCTCTCAGATGGACAAATTATTAATTTTCATCTGCATCATTTCGCTCAGTGCCGGAGTAACATCGCTTTTCACG  >  NC_000913/2638437‑2638732
                                                                                                                                                   |                                                                                                                                                    
gACAACGTTGTCCGCCAGTGCCGGATGAAGATTGGAATAGAAGTTGCCAATGCCGCTACCAACGGAAGTGCTCCACGCCGTGGTCGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTaa                                                                                                                                                     >  1:719443/1‑149 (MQ=255)
 aCAACGTTGTCCGCCAGTGCCGGATGAAGATTGGAATAGAAGTTGCCAATGCCGCTACCAACGGAAGTGCTCCACGCCGTGGTCGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTaaa                                                                                                                                                    <  3:698962/149‑1 (MQ=255)
   aaCGTTGTCCGCCAGTGCCGGATGAAGATTGGAATAGAAGTTGCCAATGCCGCTACCAACGGAAGTGCTCCACGCCGTGGTCGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACa                                                                                                                                                  >  2:348708/1‑149 (MQ=255)
   aaCGTTGTCCGCCAGTGCCGGATGAAGATTGGAATAGAAGTTGCCAATGCCGCTACCAACGGAAGTGCTCCACGCCGTGGTCGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACa                                                                                                                                                  <  1:485088/149‑1 (MQ=255)
      gTTGTCCGCCAGTGCCGGATGAAGATTGGAATAGAAGTTGCCAATGCCGCTACCAACGGAAGTGCTCCACGCCGTGGTCGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCg                                                                                                                                               >  1:28684/1‑149 (MQ=255)
       ttGTCCGCCAGTGCCGGATGAAGATTGGAATAGAAGTTGCCAATGCCGCTACCAACGGAAGTGCTCCACGCCGTGGTCGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGa                                                                                                                                              <  2:28684/149‑1 (MQ=255)
        tGTCCGCCAGTGCCGGATGAAGATTGGAATAGAAGTTGCCAATGCCGCTACCAACGGAAGTGCTCCACGCCGTGGTCGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTTAACAGCGaa                                                                                                                                             <  1:417399/149‑1 (MQ=255)
        tGTCCGCCAGTGCCGGATGAAGATTGGAATAGAAGTTGCCAATGCCGCTACCAACGGAAGTGCTCCACGCCGTGGTCGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGaa                                                                                                                                             <  3:302589/149‑1 (MQ=255)
         gTCCGCCAGTGCCGGATGAAGATTGGAATAGAAGTTGCCAATGCCGCTACCAACGGAAGTGCTCCACGCCGTGGTCGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAAc                                                                                                                                            >  4:801176/1‑149 (MQ=255)
          tCCGCCAGTGCCGGATGAAGATTGGAATAGAAGTTGCCAATGCCGCTACCAACGGAAGTGCTCCACGCCGTGGTCGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACa                                                                                                                                           >  4:247637/1‑149 (MQ=255)
             gCCAGTGCCGGATGAAGATTGGAATAGAAGTTGCCAATGCCGCTACCAACGGAAGTGCTCCACGCCGTGGTCGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGcc                                                                                                                                        <  1:296856/149‑1 (MQ=255)
               cAGTGCCGGATGAAGATTGGAATAGAAGTTGCCAATGCCGCTACCAACGGAAGTGCTCCACGCCGTGGTCGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGc                                                                                                                                      <  2:568612/149‑1 (MQ=255)
                     cGGATGAAGATTGGAATAGAAGTTGCCAATGCCGCTACCAACGGAAGTGCTCCACGCCGTGGTCGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTaaaa                                                                                                                                <  1:110427/149‑1 (MQ=255)
                     cGGATGAAGATTGGAATAGAAGTTGCCAATGCCGCTACCAACGGAAGTGCTCCACGCCGTGGTCGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTaaaa                                                                                                                                <  3:594819/149‑1 (MQ=255)
                       gATGAAGATTGGAATAGAAGTTGCCAATGCCGCTACCAACGGAAGTGCTCCACGCCGTGGTCGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAgg                                                                                                                              <  4:458842/149‑1 (MQ=255)
                       gATGAAGATTGGAATAGAAGTTGCCAATGCCGCTACCAACGGAAGTGCTCCACGCCGTGGTCGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAgg                                                                                                                              <  1:723427/149‑1 (MQ=255)
                       gATGAAGATTGGAATAGAAGTTGCCAATGCCGCTACCAACGGAAGTGCTCCACGCCGTGGTCGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAgg                                                                                                                              <  1:316198/149‑1 (MQ=255)
                       gATGAAGATTGGAATAGAAGTTGCCAATGCCGCTACCAAAGGAAGTGCTCCACGCCGTGGTCGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAgg                                                                                                                              <  3:454777/149‑1 (MQ=255)
                        aTGAAGATTGGAATAGAAGTTGCCAATGCCGCTACCAACGGAAGTGCTCCACGCCGTGGTCGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAggg                                                                                                                             >  4:409860/1‑149 (MQ=255)
                          gAAGATTGGAATAGAAGTTGCCAATGCCGCTACCAACGGAAGTGCTCCACGCCGTGGTCGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTa                                                                                                                           >  2:718295/1‑149 (MQ=255)
                            aGATTGGAATAGAAGTTGCCAATGCCGCTACCAACGGAAGTGCTCCACGCCGTGGTCGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAAc                                                                                                                         >  1:855442/1‑149 (MQ=255)
                             gATTGGAATAGAAGTTGCCAATGCCGCTACCAACGGAAGTGCTCCACGCCGTGGTCGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACg                                                                                                                        >  4:187891/1‑149 (MQ=255)
                              aTTGGAATAGAAGTTGCCAATGCCGCTACCAACGGAAGTGCTCCACGCCGTGGTCGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACgg                                                                                                                       >  2:928551/1‑149 (MQ=255)
                              aTTGGAATAGAAGTTGCCAATGCCGCTACCAACGGAAGTGCTCCACGCCGTGGTCGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACgg                                                                                                                       >  2:816726/1‑149 (MQ=255)
                                 ggAATAGAAGTTGCCAATGCCGCTACCAACGGAAGTGCTCCACGCCGTGGTCGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGaaa                                                                                                                    <  1:274515/149‑1 (MQ=255)
                                  gAATAGAAGTTGCCAATGCCGCTACCAACGGAAGTGCTCCACGCCGTGGTCGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAg                                                                                                                   >  1:495060/1‑149 (MQ=255)
                                  gAATAGAAGTTGCCAATGCCGCTACCAACGGAAGTGCTCCACGCCGTGGTCGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAg                                                                                                                   >  2:227230/1‑149 (MQ=255)
                                     tAGAAGTTGCCAATGCCGCTACCAACGGAAGTGCTCCACGCCGTGGTCGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAg                                                                                                                >  4:801479/1‑149 (MQ=255)
                                          gTTGCCAATGCCGCTACCAACGGAAGTGCTCCACGCCGTGGTCGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTaa                                                                                                                                  <  2:407769/126‑1 (MQ=255)
                                          gTTGCCAATGCCGCTACCAACGGAAGTGCTCCACGCCGTGGTCGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTaa                                                                                                                                  >  1:407769/1‑126 (MQ=255)
                                          gTTGCCAATGCCGCTACCAACGGAAGTGCTCCACGCCGTGGTCGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTg                                                                                                           >  3:213222/1‑149 (MQ=255)
                                           ttGCCAATGCCGCTACCAACGGAAGTGCTCCACGCCGTGGTCGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTgg                                                                                                          >  3:125190/1‑149 (MQ=255)
                                                aaTGCCGCTACCAACGGAAGTGCTCCACGCCGTGGTCGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGTGTAACGGAAAGCAGTCCTGgcagca                                                                                                     >  4:26024/1‑149 (MQ=255)
                                                aaTGCCGCTACCAACGGAAGTGCTCCACGCCGTGGTCGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGgcagca                                                                                                     >  4:694636/1‑149 (MQ=255)
                                                    ccGCTACCAACGGAAGTGCTCCACGCCGTGGTCGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTaa                                                                                                 <  3:187891/149‑1 (MQ=255)
                                                     cGCTACCAACGGAAGTGCTCCACGCCGTGGTCGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTAAt                                                                                                <  2:877038/149‑1 (MQ=255)
                                                     cGCTACCAACGGAAGTGCTCCACGCCGTGGTCGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTAAt                                                                                                <  4:571611/149‑1 (MQ=255)
                                                         aCCAACGGAAGTGCTCCACGCCGTGGTCGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGaaa                                                                                                                    >  3:376728/1‑125 (MQ=255)
                                                         aCCAACGGAAGTGCTCCACGCCGTGGTCGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTAATTTaa                                                                                            <  2:792508/149‑2 (MQ=255)
                                                         aCCAACGGAAGTGCTCCACGCCGTGGTCGGCATAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGaaa                                                                                                                    <  4:376728/125‑1 (MQ=255)
                                                           cAACGGAAGTGCTCCACGCCGTGGTCGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTAATTTACGc                                                                                          >  2:388462/1‑149 (MQ=255)
                                                            aaCGGAAGTGCTCCACGCCGTGGTCGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTAATTTACGCa                                                                                         >  2:866276/1‑149 (MQ=255)
                                                              cGGAAGTGCTCCACGCCGTGGTCGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTAATTTACGCAAt                                                                                       <  4:378694/149‑1 (MQ=255)
                                                               ggAAGTGCTCCACGCCGTGGTCGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTAATTTACGCAAtt                                                                                      <  4:396521/149‑1 (MQ=255)
                                                                  aGTGCTCCACGCCGTGGTCGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTAATTTACGCAATTGCa                                                                                   <  1:348708/149‑1 (MQ=255)
                                                                  aGTGCTCCACGCCGTGGTCGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTAATTTACGCAATTGCa                                                                                   >  3:265091/1‑149 (MQ=255)
                                                                      cTCCACGCCGTGGTCGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTAATTTa                                                                                             <  3:803187/135‑1 (MQ=255)
                                                                      cTCCACGCCGTGGTCGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTAATTTa                                                                                             >  4:803187/1‑135 (MQ=255)
                                                                      cTCCACGCCGTGGTCGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTAATTTACGCAATTGCATCgg                                                                               >  2:486498/1‑149 (MQ=255)
                                                                       tCCACGCCGTGGTCGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTAATTTACGCAATTGCATCggg                                                                              <  1:227230/149‑1 (MQ=255)
                                                                       tCCACGCCGTGGTCGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTAATTTACGCAATTGCATCggg                                                                              >  2:565740/1‑149 (MQ=255)
                                                                       tCCACGCCGTGGTCGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTAATTTACGCAATTGCATCggg                                                                              >  1:628971/1‑149 (MQ=255)
                                                                       tCCACGCCGTGGTCGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTAATTTACGCAATTGCATCggg                                                                              >  3:288646/1‑149 (MQ=255)
                                                                          aCGCCGTGGTCGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGaaa                                                                                                                    <  2:133300/108‑1 (MQ=255)
                                                                          aCGCCGTGGTCGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGaaa                                                                                                                    >  1:133300/1‑108 (MQ=255)
                                                                            gCCGTGGTCGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTAATTTACGCAATTGCATCGGGTCCct                                                                         <  2:818305/149‑1 (MQ=255)
                                                                                tGGTCGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTAATTTACGCAATTGCATCGGGTCCCTCTCa                                                                     <  2:719443/149‑1 (MQ=255)
                                                                                  gTCGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTAATTTACGCAATTGCATCGGGTCCCTCTCAGa                                                                   <  2:394556/149‑1 (MQ=255)
                                                                                   tCGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTAATTTACGCAATTGCATCGGGTCCCTCTCAGAt                                                                  <  1:816726/149‑1 (MQ=255)
                                                                                    cGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTAATTTACGCAATTGCATCGGGTCCCTCTCAGATg                                                                 <  3:575804/149‑1 (MQ=255)
                                                                                    cGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTAATTTACGCAATTGCATCGGGTCCCTCTCAGATg                                                                 <  4:265091/149‑1 (MQ=255)
                                                                                     ggCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTACGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTACTGGCAGCAGTAATTTACGCAATTGCATCGGGTCCCTCTCAGATgg                                                                <  2:694009/149‑1 (MQ=255)
                                                                                       cGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTAATTTACGCAATTGCATCGGGTCCCTCTCAGATGGAc                                                              <  4:287694/149‑1 (MQ=255)
                                                                                           aaCTGGTTTTCAACGGTTGGCAATGGGGATATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTAATTTACGCAATTGCATCGGGTCCCTCTCAGATGGACAAat                                                          <  1:173077/149‑1 (MQ=255)
                                                                                           aaCTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAgg                                                                                                                              >  3:905370/1‑81 (MQ=255)
                                                                                           aaCTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAgg                                                                                                                              <  4:905370/81‑1 (MQ=255)
                                                                                           aaCTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTAATTTACGCAATTGCATCGGGTCCCTCTCAGATGGACAAat                                                          <  1:388462/149‑1 (MQ=255)
                                                                                             cTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTAATTTACGCAATTGCATCGGGTCCCTCTCAGATGGACAAatta                                                        >  3:382622/1‑149 (MQ=255)
                                                                                              tGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTAATTTACGCAATTGCATCGGGTCCCTCTCAGATGGACAAattat                                                       >  4:644471/1‑149 (MQ=255)
                                                                                                 ttttCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTAATTTACGCAATTGCATCGGGTCCCTCTCAGATGGACAAATTATTaa                                                    <  4:125190/149‑1 (MQ=255)
                                                                                                      aaCGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTACACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTAATTTACGCAATTGCATCGGGTCCCTCTCAGATGGACAAATTATTAATTTTc                                               >  1:704036/1‑149 (MQ=255)
                                                                                                      aaCGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTAATTTACGCAATTGCATCGGGTCCCTCTCAGATGGACAAATTATTAATTTTc                                               <  3:73346/149‑1 (MQ=255)
                                                                                                      aaCGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTAATTTACGCAATTGCATCGGGTCCCTCTCAGATGGACAAATTATTAATTTTc                                               <  1:928551/149‑1 (MQ=255)
                                                                                                                      gACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTAATTTACGCAATTGCATCGGGTCCCTCTCAGATGGACAAATTATTAATTTTCATCTGCATCATTTc                                 <  4:133266/147‑1 (MQ=255)
                                                                                                                      gACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTAATTTACGCAATTGCATCGGGTCCCTCTCAGATGGACAAATTATTAATTTTCATCTGCATCATTTc                                 >  3:133266/1‑147 (MQ=255)
                                                                                                                         aTCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTAATTTACGCAATTGCATCGGGTCCCTCTCAGATGGACAAATTATTAATTTTCATCTGCATCATTTCGCTCa                            >  4:254284/1‑149 (MQ=255)
                                                                                                                            tttACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTAATTTACGCAATTGCATCGGGTCCCTCTCAGATGGACAAATTATTAATTTTCATCTGCATCATTTCGCTCAGTg                         >  1:514675/1‑149 (MQ=255)
                                                                                                                                ccACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTAATTTACGCAATTGCATCGGGTCCCTCTCAGATGGACAAATTATTAATTTTCATCTGCATCATTTCGCTCAGTGCCgg                     >  2:109679/1‑149 (MQ=255)
                                                                                                                                 cacaTCTTCTTCGATGTTAAACAGCGAACAGCCCCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTAATTTACGCAATTGCATCGGGTCCCTCTCAGATGGACAAATTATTAATTTTCATCTGCATCATTTCGCTCAGTGCCGGa                    <  3:26024/149‑1 (MQ=255)
                                                                                                                                  acaTCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTAATTTACGCAATTGCATCGGGTCCCTCTCAGATGGACAAATTATTAATTTTCATCTGCATCATTTCGCTCAGTGCCGGAg                   >  4:54408/1‑149 (MQ=255)
                                                                                                                                  acaTCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTAATTTACGCAATTGCATCGGGTCCCTCTCAGATGGACAAATTATTAATTGTCATCTGCATCATTTCGCTCAGTGCCGGAg                   >  3:753196/1‑149 (MQ=255)
                                                                                                                                  acaTCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGAAACGGAAACAAGTCTTGGCAGCAGTAATTTAGCCATTTGCACGGGGCCCTTCTCAGAGGGACAATTTTTAATTTTTCATTTCAATATTTTCGTTCAGGCCggggg                  >  1:149615/1‑146 (MQ=255)
                                                                                                                                   caTCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTAATTTACGCAATTGCATCGGGTCCCTCTCAGATGGACAAATTATTAATTTTCATCTGCATCATTTCGCTCAGTGCCGGAGt                  >  3:303296/1‑149 (MQ=255)
                                                                                                                                    aTCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTAATTTACGCAATTGCATCGGGTCCCTCTCAGATGGACAAATTATTAATTTTCATCTGCATCATTTCGCTCAGTGCCGGAGTa                 >  4:535582/1‑149 (MQ=255)
                                                                                                                                     tcttctTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTAATTTACGCAATTGCATCGGGTCCCTCTCAGATGGACAAATTATTAATTTTCATCTGCATCATTTCGCTCAGTGCCGGAGTaa                <  3:920449/149‑1 (MQ=255)
                                                                                                                                      cttcttCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTAATTTACGCAATTGCATCGGGTCCCTCTCAGATGGACAAATTATTAATTTTCATCTGCATCATTTCGCTCAGTGCCGGAGTAAc               >  1:235610/1‑149 (MQ=255)
                                                                                                                                       ttcttcGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTAATTTACGCAATTGCATCGGGTCCCTCTCAGATGGACAAATTATTAATTTTCATCTGCATCATTTCGCTCAGTGCCGGAGTAACa              >  1:147439/1‑149 (MQ=255)
                                                                                                                                       ctgttcGCTGTTTAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTAATTTACGCAATTGCATCGGGTCCCTCTCAGATGGACAAATTATTAATTTTCATCTGCATCATTTCGCTCAGTGCCGGAGTAACa              >  4:373734/4‑149 (MQ=255)
                                                                                                                                         cttcGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACTGGAAGCAGTCCTGGCAGCAGTAATTTACGCAATTGCATCGGGTCCCTCTCAGATGGACAAATTATTAATTTTCATCTGCATCATTTCGCTCAGTGCCGGAGTAACATc            <  4:803154/149‑1 (MQ=255)
                                                                                                                                         cttcGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTAATTTACGCAATTGCATCGGGTCCCTCTCAGATGGACAAATTATTAATTTTCATCTGCATCATTTCGCTCAGTGCCGGAGTAACATc            <  2:422974/149‑1 (MQ=255)
                                                                                                                                          ttcGCTGTTTAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTAATTTACGCAATTGCATCGGGTCCCTCTCAGATGGACAAATTATTAATTTTCATCTGCATCATTTCGCTCAGTGCCGGAGTAACATCg           >  4:521570/1‑149 (MQ=255)
                                                                                                                                          ttcGCTGTTAAACAGCGAACAGTCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTAATTTACGCAATTGCATCGGGTCCCTCTCAGATGGACAAATTATTAATTTTCATCTGCATCATTTCGCTCAGTGCCGGAGTAACATCg           >  1:517299/1‑149 (MQ=255)
                                                                                                                                          ttcGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTAATTTACGCAATTGCATCGGGTCCCTCTCAGATGGACAAATTATTAATTTTCATCTGCATCATTTCGCTCAGTGCCGGAGTAACATCg           >  3:115746/1‑149 (MQ=255)
                                                                                                                                            cGCTGTTTAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTAATTTACGCAATTGCATCGGGTCCCTCTCAGATGGACAAATTATTAATTTTCATCTGCATCATTTCGCTCAGTGCCGGAGTAACATCGCt         >  1:697396/1‑149 (MQ=255)
                                                                                                                                               tGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTAATTTACGCAATTGCATCGGGTCCCTCTCAGATGGACAAATTATTAATTTTCATCTGCATCATTTCGCTCAGTGCCGGAGTAACATCGCtttt      <  4:765615/149‑1 (MQ=255)
                                                                                                                                                 ttAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTAATTTACGCAATTGCATCGGGTCCCTCTCAGATGGACAAATTATTAATTTTCATCTGCATCATTTCGCTCAGTGCCGGAGTAACATCGCTTTTCa    >  1:722918/1‑149 (MQ=255)
                                                                                                                                                  tAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTAATTTACGCAATTGCATCGGGTCCCTCTCAGATGGACAAATTATTAATTTTCATCTGCATCATTTCGCTCAGTGCCGGAGTAACATCGCTTTTCAc   <  3:801176/149‑1 (MQ=255)
                                                                                                                                                  tAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTAATTTACGCAATTGCATCGGGTCCCTCTCAGATGGACAAATTATTAATTTTCATCTGCATCATTTCGCTCAGTGCCGGAGTAACATCGCTTTTCAc   <  4:332957/149‑1 (MQ=255)
                                                                                                                                                   aaaCAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTAATTTACGCAATTGCATCGGGTCCCTCTCAGATGGACAAATTATTAATTTTCATCTGCATCATTTCGCTCAGTGCCGGAGTAACATCGCTTTTCACg  <  2:628971/149‑1 (MQ=255)
                                                                                                                                                   aaaCAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTAATTTACGCAATTGCATCGGGTCCCTCTCAGATGGACAAATTATTAATTTTCATCTGCATCATTTCGCTCAGTGCCGGAGTAACATCGCTTTTCACg  >  4:361978/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: