Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,638,584 A→T 5.1% 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.1% 423.9 / 8.9 118F23L (TTT→TTA) bamBBamABCDE complex OM biogenesis lipoprotein
Reads supporting (aligned to +/- strand):  ref base A (49/63);  new base T (3/3);  total (52/66)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 119 total aligned reads displayed.

AGACAACGTTGTCCGCCAGTGCCGGATGAAGATTGGAATAGAAGTTGCCAATGCCGCTACCAACGGAAGTGCTCCACGCCGTGGTCGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTAATTTACGCAATTGCATCGGGTCCCTCTCAGATGGACAAATTATTAATTTTCATCTGCATCATTTCGCTCAGTGCCGGAGTAACATCGCTTTTCA  >  NC_000913/2638436‑2638730
                                                                                                                                                    |                                                                                                                                                  
aGACAACGTTGTCCGCCAGTGCCGGATGAAGATTGGAATAGAAGTTGCCAATGCCGCTACCAACGGAAGTGCTCCACGCCGTGGTCGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTa                                                                                                                                                    >  4:96088/1‑149 (MQ=255)
    aaCGTTGTCCGCCAGTGCCGGATGAAGATTGGAATAGAAGTTGCCAATGCCGCTACCAACGGAAGTGCTCCACGCCGTGGTCGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACa                                                                                                                                                >  3:14829/1‑149 (MQ=255)
    aaCGTTGTCCGCCAGTGCCGGATGAAGATTGGAATAGAAGTTGCCAATGCCGCTACCAACGGAAGTGCTCCACGCCGTGGTCGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACa                                                                                                                                                >  4:902638/1‑149 (MQ=255)
    aaCGTTGTCCGCCAGTGCCGGATGAAGATTGGAATAGAAGTTGCCAATGCCGCTACCAACGGAAGTGCTCCACGCCGTGGTCGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACa                                                                                                                                                >  3:516671/1‑149 (MQ=255)
    aaCGTTGTCCGCCAGTGCCGGATGAAGATTGGAATAGAAGTTGCCAATGCCGCTACCAACGGAAGTGCTCCACGCCGTGGTCGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACa                                                                                                                                                <  4:516671/149‑1 (MQ=255)
     aCGTTGTCCGCCAGTGCCGGATGAAGATTGGAATAGAAGTTGCCAATGCCGCTACCAACGGAAGTGCTCCACGCCGTGGTCGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAg                                                                                                                                               >  1:104877/1‑149 (MQ=255)
      cGTTGTCCGCCAGTGCCGGATGAAGATTGGAATAGAAGTTGCCAATGCCGCTACCAACGGAAGTGCTCCACGCCGTGGTCGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGc                                                                                                                                              <  4:530864/149‑1 (MQ=255)
        ttGTCCGCCAGTGCCGGATGAAGATTGGAATAGAAGTTGCCAATGCCGCTACCAACGGAAGTGCTCCACGCCGTGGTCGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGa                                                                                                                                            <  3:1136965/149‑1 (MQ=255)
         tGTCCGCCAGTGCCGGATGAAGATTGGAATAGAAGTTGCCAATGCCGCTACCAACGGAAGTGCTCCACGCCGTGGTCGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGaa                                                                                                                                           >  1:315722/1‑149 (MQ=255)
          gTCCGCCAGTGCCGGATGAAGATTGGAATAGAAGTTGCCAATGCCGCTACCAACGGAAGTGCTCCACGCCGTGGTCGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTTAACAGCGaag                                                                                                                                          <  2:87349/149‑2 (MQ=255)
          gTCCGCCAGTGCCGGATGAAGATTGGAATAGAAGTTGCCAATGCCGCTACCAACGGAAGTGCTCCACGCCGTGGTCGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAAc                                                                                                                                          <  3:385011/149‑1 (MQ=255)
           tCCGCCAGTGCCGGATGAAGATTGGAATAGAAGTTGCCAATGCCGCTACCAACGGAAGTGCTCCACGCCGTGGTCGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACa                                                                                                                                         <  2:406875/149‑1 (MQ=255)
            ccGCCAGTGCCGGATGAAGATTGGAATAGAAGTTGCCAATGCCGCTACCAACGGAAGTGCTCCACGCCGTGGTCGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTTAACAGCGaagaa                                                                                                                                        <  3:842534/149‑4 (MQ=255)
              gCCAGTGCCGGATGAAGATTGGAATAGAAGTTGCCAATGCCGCTACCAACGGAAGTGCTCCACGCCGTGGTCGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGcc                                                                                                                                      <  1:336664/149‑1 (MQ=255)
                      cGGATGAAGATTGGAATAGAAGTTGCCAATGCCGCTACCAACGGAAGTGCTCCACGCCGTGGTCGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTaaaa                                                                                                                              <  1:910351/149‑1 (MQ=255)
                        gATGAAGATTGGAATAGAAGTTGCCAATGCCGCTACCAACGGAAGTGCTCCACGCCGTGGTCGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAgg                                                                                                                            >  3:1059397/1‑149 (MQ=255)
                         aTGAAGATTGGAATAGAAGTTGCCAATGCCGCTACCAACGGAAGTGCTCCACGCCGTGGTCGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAggg                                                                                                                           <  3:96088/149‑1 (MQ=255)
                         aTGAAGATTGGAATAGAAGTTGCCAATGCCGCTACCAACGGAAGTGCTCCACGCCGTGGTCGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAggg                                                                                                                           >  4:742345/1‑149 (MQ=255)
                          tGAAGATTGGAATAGAAGTTGCCAATGCCGCTACCAACGGAAGTGCTCCGCGCCGTGGTCGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGt                                                                                                                          >  3:594588/1‑149 (MQ=255)
                          tGAAGATTGGAATAGAAGTTGCCAATGCCGCTACCAACGGAAGTGCTCCACGCCGTGGTCGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGt                                                                                                                          <  3:827705/149‑1 (MQ=255)
                           gAAGATTGGAATAGAAGTTGCCAATGCCGCTACCAACGGAAGTGCTCCACGCCGTGGTCGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTa                                                                                                                         >  4:512489/1‑149 (MQ=255)
                             aGATTGGAATAGAAGTTGCCAATGCCGCTACCAACGGAAGTGCTCCACGCCGTGGTCGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAAc                                                                                                                       >  3:562261/1‑149 (MQ=255)
                              gATTGGAATAGAAGTTGCCAATGCCGCTACCAACGGAAGTGCTCCACGCCGTGGTCGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACg                                                                                                                      >  2:1150789/1‑149 (MQ=255)
                               aTTGGAATAGAAGTTGCCAATGCCGCTACCAACGGAAGTGCTCCACGCCGTGGTCGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACgg                                                                                                                     >  1:332272/1‑149 (MQ=255)
                               aTTGGAATAGAAGTTGCCAATGCCGCTACCAACGGAAGTGCTCCACGCCGTGGTCGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACgg                                                                                                                     >  4:893104/1‑149 (MQ=255)
                                  ggAATAGAAGTTGCCAATGCCGCTACCAACGGAAGTGCTCCACGCCGTGGTCGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGaaa                                                                                                                  <  4:349295/149‑1 (MQ=255)
                                   gAATAGAAGTTGCCAATGCCGCTACCAACGGAAGTGCTCCACGCCGTGGTCGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAg                                                                                                                 <  3:193706/149‑1 (MQ=255)
                                      tAGAAGTTGCCAATGCCGCTACCAACGGAAGTGCTCCACGCCGTGGTCGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAg                                                                                                              >  1:238797/1‑149 (MQ=255)
                                        gAAGTTGCCAATGCCGCTACCAACGGAAGTGCTCCACGCCGTGGTCGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTc                                                                                                            <  3:311469/149‑1 (MQ=255)
                                         aaGTTGCCAATGCCGCTACCAACGGAAGTGCTCCACGCCGTGGTCGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGaaa                                                                                                                  >  4:739950/1‑142 (MQ=255)
                                         aaGTTGCCAATGCCGCTACCAACGGAAGTGCTCCACGCCGTGGTCGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTcc                                                                                                           >  3:476128/1‑149 (MQ=255)
                                         aaGTTGCCAATGCCGCTACCAACGGAAGTGCTCCACGCCGTGGTCGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAACGGTAACGGaaa                                                                                                                  <  3:739950/142‑1 (MQ=255)
                                          aGTTGCCACTGCCGCTAACAACGGAAGTGCTCCACGCCGTGGTCGGCGTAACTTGGTTTTTAACGGTTGGCAATGGGGAAATCTTTACCCAATCTAATTCGCTGTTAAACAGAGACCAGCCGCTTACAAGGGTAACGGAACGCAGTcat                                                                                                          <  2:326399/149‑3 (MQ=255)
                                            ttGCCAATGCCGCTACCAACGGAAGTGCTCCACGCCGTGGTCGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTgg                                                                                                        >  3:526954/1‑149 (MQ=255)
                                            ttGCCAATGCCGCTACCAACGGAAGTGCTCCACGCCGTGGTCGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTgg                                                                                                        <  3:893104/149‑1 (MQ=255)
                                              gCCAATGCCGCTACCAACGGAAGTGCTCCACGCCGTGGTCGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTTAACAGCGaagaa                                                                                                                                        <  2:805573/115‑4 (MQ=255)
                                              gCCAATGCCGCTACCAACGGAAGTGCTCCACGCCGTGGTCGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTTAACAGCGaagaa                                                                                                                                        >  1:805573/1‑112 (MQ=255)
                                               ccAATGCCGCTACCAACGGAAGTGCTCCACGCCGTGGTCGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGgcag                                                                                                     <  2:496318/149‑1 (MQ=255)
                                                cAATGCCGCTACCAACGGAAGTGCTCCACGCCGTGGTCGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGgcagc                                                                                                    >  1:1037432/1‑149 (MQ=255)
                                                     ccGCTACCAACGGAAGTGCTCCACGCCGTGGTCGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTaa                                                                                               >  2:72189/1‑149 (MQ=255)
                                                       gCTACCAACGGAAGTGCTCCACGCCGTGGTCGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTg                                                                                                         >  2:1083194/1‑137 (MQ=255)
                                                          aCCAACGGAAGTGCTCCACGCCGTGGTCGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTAATTTAc                                                                                          >  3:1077278/1‑149 (MQ=255)
                                                          aCCAACGGAAGTGCTCCACGCCGTGGTCGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTAATTTAc                                                                                          <  1:888850/149‑1 (MQ=255)
                                                           ccAACGGAAGTGCTCCACGCCGTGGTCGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTAATTTACg                                                                                         <  4:594588/149‑1 (MQ=255)
                                                             aaCGGAAGTGCTCCACGCCGTGGTCGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAgg                                                                                                                            <  1:1024693/112‑1 (MQ=255)
                                                             aaCGGAAGTGCTCCACGCCGTGGTCGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAgg                                                                                                                            >  2:1024693/1‑112 (MQ=255)
                                                               cGGAAGTGCTCCACGCCGTGGTCGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTAATTTACGCAAt                                                                                     <  4:513516/149‑1 (MQ=255)
                                                                   aGTGCTCCACGCCGTGGTCGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTAATTTACGCAATTGCa                                                                                 <  4:127604/149‑1 (MQ=255)
                                                                   aGTGCTCCACGCCGTGGTCGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTAATTTACGCAATTGCa                                                                                 >  3:945174/1‑149 (MQ=255)
                                                                         ccACGCCGTGGTCGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTAATTTACGCAATTGCATCGGGt                                                                           >  1:1066280/1‑149 (MQ=255)
                                                                           aCGCCGTGGTCGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTAATTTACGCAATTGCATCGGGTcc                                                                         <  2:238797/149‑1 (MQ=255)
                                                                           aCGCCGTGGTCGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTAATTTACGCAATTGCATAGGGTcc                                                                         <  4:1048342/149‑1 (MQ=255)
                                                                             gCCGTGGTCGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTAATTTACGCAATTGCATCGGGTCCct                                                                       >  3:683296/1‑149 (MQ=255)
                                                                                gTGGTCGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTAATTTACGCAATTGCATCGGGTCCCtctc                                                                    <  4:476128/149‑1 (MQ=255)
                                                                                    tCGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTAATTTACGCAATTGCATCGGGTCCCTCTCAGAt                                                                <  1:72189/149‑1 (MQ=255)
                                                                                    tCGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTAATTTACGCAATTGCATCGGGTCCCTCTCAGAt                                                                <  4:1059397/149‑1 (MQ=255)
                                                                                       gCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTAATTTACGCAATTGCATCGGGTCCCTCTCAGATGGa                                                             <  3:1145693/149‑1 (MQ=255)
                                                                                           aaaCTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTAATTTACGCAATTGCATCGGGTCCCTCTCAGATGGACaaa                                                         <  3:742345/149‑1 (MQ=255)
                                                                                              cTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTAATTTACGCAATTGCATCGGGTCCctct                                                                     <  3:1066468/134‑1 (MQ=255)
                                                                                              cTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTAATTTACGCAATTGCATCGGGTCCctct                                                                     >  4:1066468/1‑134 (MQ=255)
                                                                                              cTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTAATTTACGCAATTGCATCGGGTCCCTCTCAGATGGACAAatta                                                      <  1:154775/149‑1 (MQ=255)
                                                                                              cTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTAATTTACGCAATTGCATCGGGTCCCTCTCAGATGGACAAatta                                                      <  3:902638/149‑1 (MQ=255)
                                                                                               tGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTAATTTACGCAATTGCATCGGGTCCCTCTCAGATGGACAAattct                                                     <  1:268590/149‑3 (MQ=255)
                                                                                                ggTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTAATTTa                                                                                           <  4:307128/110‑1 (MQ=255)
                                                                                                ggTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTAATTTa                                                                                           >  3:307128/1‑110 (MQ=255)
                                                                                                  ttttCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGgc                                                                                                       >  4:862940/1‑96 (MQ=255)
                                                                                                    ttCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTAATTTACGCAATTGCATCGGGTCCCTCTCAGATGGACAAATTATTAAtt                                                <  3:86986/149‑1 (MQ=255)
                                                                                                       aaCGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTAATTTACGCAATTGCATCGGGTCCCTCTCAGATGGACAAATTATTAATTTTc                                             <  4:564861/149‑1 (MQ=255)
                                                                                                       aaCGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTAATTTACGCAATTGCATCGGGTCCCTCTCAGATGGACAAATTATTAATTTTc                                             >  1:328113/1‑149 (MQ=255)
                                                                                                       aaCGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTAATTTACGCAATTGCATCGGGTCCCTCTCAGATGGACAAATTATTAATTTTc                                             <  3:875193/149‑1 (MQ=255)
                                                                                                        aCGGTTGGAAATGGGGACATCTTTACCAAATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTAATTTACGCAATTGCATCGGGTCCCTCTCAGATGGACAAATTATTAATTTTCa                                            <  3:409545/149‑1 (MQ=255)
                                                                                                             tGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTAATTTACGCAATTGCATCGGGTCCCTCTCAGATGGACAAATTATTAATTTTCATCTGc                                       <  4:526954/149‑1 (MQ=255)
                                                                                                               gCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTAATTTACGCAATTGCATCGGGTCCCTCTCAGATGGACAAATTATTAATTTTCATCTGcat                                     >  3:417350/1‑149 (MQ=255)
                                                                                                                    ggggACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTAATTTACGCAATTGCATCGGGTCCCTCTCAGATGGACAAATTATTAATTTTCATCTGCATCAttt                                <  2:1128521/149‑1 (MQ=255)
                                                                                                                     gggACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTAATTTACGCAATTGCATCGGGTCCCTCTCAGATGGACa                                                           <  2:667921/121‑1 (MQ=255)
                                                                                                                      ggACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTAATTTACGCAAtt                                                                                    >  1:1007157/1‑95 (MQ=255)
                                                                                                                      ggACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTAATTTACGCAAtt                                                                                    <  2:1007157/95‑1 (MQ=255)
                                                                                                                       gACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGTAGTCCTGGCAGCAGTAATTTACGCAATTGCATCGGGTCCCTCTCAGATGGACAAATTATTAATTTTCATCTGCATCATTTCGc                             >  4:427787/1‑149 (MQ=255)
                                                                                                                       gACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTAATTTACGCAATTGCATCGGGTCCCTCTCAGATGGACAAATTATTAATTTTCATCTGCATCATTTCGc                             >  2:48343/1‑149 (MQ=255)
                                                                                                                        aCATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTAATTTACGCAATTGCATCGGGTCCCTCTCAGATGGACAAATTATTAATTTTCATCTGCATCATTTCGCt                            >  3:998994/1‑149 (MQ=255)
                                                                                                                          aTCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTAATTTACGCAATTGCATCGGGTCCCTCTCAGATGGACAAATTATTAATTTTCATCTGCATCATTTCg                              >  1:107561/1‑145 (MQ=255)
                                                                                                                          aTCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTAATTTACGCAATTGCATCGGGTCCCTCTCAGATGGACAAATTATTAATTTTCATCTGCATCATTTCg                              <  2:107561/145‑1 (MQ=255)
                                                                                                                                aCCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTAATTTACGCAATTGCATCGGGTCCCTCTCAGATGGACAAat                                                        <  4:952326/113‑1 (MQ=255)
                                                                                                                                aCCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTAATTTACGCAATTGCATCGGGTCCCTCTCAGATGGACAAATTATTAATTTTCATCTGCATCATTTCGCTCAGTGCCg                    <  2:328113/149‑1 (MQ=255)
                                                                                                                                aCCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTAATTTACGCAATTGCATCGGGTCCCTCTCAGATGGACAAATTATTAATTTTCATCTGCATCATTTCGCTCAGTGCCg                    <  4:562261/149‑1 (MQ=255)
                                                                                                                                 ccACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTAATTTACGCAATTGCATCGGGTCCCTCTCAGATGGACAAATTATTAATTTTCATCTGCATCATTTCGCTCAGTGCCgg                   >  4:732902/1‑149 (MQ=255)
                                                                                                                                  cacaTCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTAATTTACGCAATTGCATCGGGTCCCTCTCAGATGGACAAATTATTAATTTTCATCTGCATCATTTCGCTCAGTGCCGGa                  <  3:732902/149‑1 (MQ=255)
                                                                                                                                  cacaTCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTAATTTACGCAATTGCATCGGGTCCCTCTCAGATGGACAAATTATTAATTTTCATCTGCATCATTTCGCTCAGTGCCGGa                  <  4:149741/149‑1 (MQ=255)
                                                                                                                                   acaTCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTAATTTACGCAATTGCATCGGGTCCCTCTCAGATGGACAAATTATTAATTTTCATCTGCATCATTTCGCTCAGTGCCGGAg                 >  4:352227/1‑149 (MQ=255)
                                                                                                                                   acaTCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTAATTTACGCAATTGCATCGGGTCCCTCTCAGATGGACAAATTATTAATTTTCATCTGCATCATTTCGCTCAGTGCCGGAg                 >  2:741399/1‑149 (MQ=255)
                                                                                                                                    caTCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTAATTTACGCAATTGCATCGGGTCCCTCTCAGATGGACAAATTATTAATTTTCATCTGCATCATTTCGCTCa                          <  1:664054/139‑1 (MQ=255)
                                                                                                                                      tcttctTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTAATTTACGCAATTGCATCGGGTCCCTCTCAGATGGACAAATTATTAATTTTCATCTGCATCATTTCGCTCAGTGCCGGAGTaa              >  1:235738/1‑149 (MQ=255)
                                                                                                                                        ttcttcGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTAATTTACGCAATTGCATCGGGTCCCTCTCAGATGGACAAATTATTAATTTTCATCTGCATCATTTCGCTCAGTGCCGGAGTAACa            >  3:191483/1‑149 (MQ=255)
                                                                                                                                         tgttcGCTGTTTAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTAATTTACGCAATTGCATCGGGTCCCTCTCAGATGGACAAATTATTAATTTTCATCTGCATCATTTCGCTCAGTGCCGGAGTAACAt           >  1:453512/3‑149 (MQ=255)
                                                                                                                                         tgttcGCTGTTTAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTAATTTACGCAATTGCATCGGGTCCCTCTCAGATGGACAAATTATTAATTTTCATCTGCATCATTTCGCTCAGTGCCGGAGTAACAt           >  3:300617/3‑149 (MQ=255)
                                                                                                                                          cttcGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGTAGTCCTGGCAGCAGTAATTTACGCAATTGCATCGGGTCCCTCTCAGATGGACAAATTATTAATTTTCATCTGCATCATTTCGCTCAGTGCCGGAGTAACATc          <  3:427787/149‑1 (MQ=255)
                                                                                                                                            tcGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTAATTTACGCAATTGCATCGGGTCCCTCTCAGATGGACAAATTATTAATTTTCATCTGCATCATTTCGCTCAGTGCCGGa                  <  3:519599/139‑1 (MQ=255)
                                                                                                                                            tcGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTAATTTACGCAATTGCATCGGGTCCCTCTCAGATGGACAAATTATTAATTTTCATCTGCATCATTTCGCTCAGTGCCGGa                  >  4:519599/1‑139 (MQ=255)
                                                                                                                                                 gTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTAATTTACGCAATTGCATCGGGTCCCTCTCAGATGGACAAATTATTAATTTTCATCTGCATCATTTCGCTCAGTGCCGGAGTAACATCGCTTTTc   <  4:945174/149‑1 (MQ=255)
                                                                                                                                                  ttAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTAATTTACGCAATTGCATCGGGTCCCTCTCAGATGGACAAATTATTAATTTTCATCTGCATCATTTCGCTCAGTGCCGGAGTAACATCGCTTTTCa  <  4:683296/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
AGACAACGTTGTCCGCCAGTGCCGGATGAAGATTGGAATAGAAGTTGCCAATGCCGCTACCAACGGAAGTGCTCCACGCCGTGGTCGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTAATTTACGCAATTGCATCGGGTCCCTCTCAGATGGACAAATTATTAATTTTCATCTGCATCATTTCGCTCAGTGCCGGAGTAACATCGCTTTTCA  >  NC_000913/2638436‑2638730

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