Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,637,444 A→G 6.3% intergenic (‑88/+30) der ← / ← bamB GTPase; multicopy suppressor of ftsJ/BamABCDE complex OM biogenesis lipoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,637,4440AG 6.3% 241.6 / 3.8 76intergenic (‑88/+30)der/bamBGTPase; multicopy suppressor of ftsJ/BamABCDE complex OM biogenesis lipoprotein
Reads supporting (aligned to +/- strand):  ref base A (31/40);  new base G (3/2);  total (34/42)
Fisher's exact test for biased strand distribution p-value = 6.51e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

GGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGACCACAGGTACCATGTTTAAAGCCTCATTTTTATAAATCATCAGACAACGCACGCTATATTCGCGTCGTTGTTAAAAACAGGAAAACGGCCCCTGTCCAGGAGCCGTTTTCAAAGTGAACGACAGAGACGATTAACGTGTAATAGAGTACACGGTTCCGTCTTTTGCCTGGATCAGCAGTTTGCCGTCAGCGGCAACCGGTTCAGTCTGGAAACCGGAACTATCAACTTTTTGCTGGGCAACGAAACGA  >  NC_000913/2637296‑2637591
                                                                                                                                                    |                                                                                                                                                   
ggTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGACCACAGGTACCATGTTTAAAGCCTCATTTTTATAAATCATCAGACAACGCACGCTATATTCGCGTCGTTGTTAAAAACAGGAAAACGGCCCCTGTCCAGGa                                                                                                                                                     <  2:146753/149‑1 (MQ=255)
     aaTAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGACCACAGGTACCATGTTTAAAGCCTCATTTTTATAAATCATCAGACAACGCACGCTATATTCGCGTCGTTGTTAAAAACAGGAAAACGGCCCCTGTCCAGGAGCCGt                                                                                                                                                <  3:20843/149‑1 (MQ=255)
     aaTAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGACCACAGGTACCATGTTTAAAGCCTCATTTTTATAAATCATCAGACAACGCACGCTATATTCGCGTCGTTGTTAAAAACAGGAAAACGGCCCCTGTCCAGGAGCCGt                                                                                                                                                <  3:74233/149‑1 (MQ=255)
      aTAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGACCACAGGTACCATGTTTAAAGCCTCATTTTTATAAATCATCAGACAACGCACGCTATATTCGCGTCGTTGTTAAAAACAGGAAAACGGCCCCTGTCCAGGAGCCGtt                                                                                                                                               <  4:408681/149‑1 (MQ=255)
      aTAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGACCACAGGTACCATGTTTAAAGCCTCATTTTTATAAATCATCAGACAACGCACGCTATATTCGCGTCGTTGTTAAAAACAGGAAAACGGCCCCTGTCCAGGAGCCGtt                                                                                                                                               >  1:44865/1‑149 (MQ=255)
      aTAACGTGGATTTTCCCACGTTAGGGCGCCCGACAAGCGCGACCACAGGTACCGTGTTTAAAGCCTCATTTTTATAAATCATCAGACAACGCACGCTATATTCGCGTCGTTGTTAAAAACAGGAAAACGGCTCCTGGACAGGGGCCGtt                                                                                                                                               <  1:460935/149‑1 (MQ=255)
        aaCGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGACCACAGGTACCATGTTTAAAGCCTCATTTTTATAAATCATCAGACAACGCACGCTATATTCGCGTCGTTGTTAAAAACAGGAAAACGGCTCCTGGACAGGGGCCGtttt                                                                                                                                             <  2:182409/149‑1 (MQ=255)
        aaCGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGACCACAGGTACCATGTTTAAAGCCTCATTTTTATAAATCATCAGACAACGCACGCTATATTCGCGTCGTTGTTAAAAACAGGAAAACGGCCCCTGTCCAGGAGCCGtttt                                                                                                                                             >  4:194267/1‑149 (MQ=255)
         aCGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGACCACAGGTACCATGTTTAAAGCCTCATTTTTATAAATCATCAGACAACGCACGCTATATTCGCGTCGTTGTTAAAAACAGGAAAACGGCCCCTGTCCAGGAGCCGTTTTc                                                                                                                                            >  4:265456/1‑149 (MQ=255)
           gTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGACCACAGGTACCATGTTTAAAGCCTCATTTTTATAAATCATCAGACAACGCACGCTATATTCGCGTCGTTGTTAAAAACAGGAAAACGGCCCCTGTCCAGGAGCCGTTTTCaa                                                                                                                                          >  4:520756/1‑149 (MQ=255)
            tGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGACCACAGGTACCATGTTTAAAGCCTCATTTTTATAAATCATCAGACAACGCACGCTATATTCGCGTCGTTGTTAAAAACAGGAAAACGGCCCCTGTCCAGGAGCCGTTTTCaaa                                                                                                                                         >  3:211204/1‑149 (MQ=255)
            tGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGACCACAGGTACCATGTTTAAAGCCTCATTTTTATAAATCATCAGACAACGCACGCTATATTCGCGTCGTTGTTAAAAACAGGAAAACGGCCCCTGTCCAGGAGCCGTTTTCaaa                                                                                                                                         <  2:484840/149‑1 (MQ=255)
            tGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGACCACAGGTACCATGTTTAAAGCCTCATTTTTATAAATCATCAGACAACGCACGCTATATTCGCGTCGTTGTTAAAAACAGGAAAACGGCCCCTGTCCAGGAGCCGTTTTCaaa                                                                                                                                         <  2:330355/149‑1 (MQ=255)
               aTTTTCCTACGTTAGGGCGCCCGACAAGCGCGACCACAGGTACCATGTTTAAAGCCTCATTTTTATAAATCATCAGACAACGCACGCTATATTCGCGTCGTTGTTAAAAACAGGAAAACGGCCCCTGTCCAGGAGCCGTTTTCAAAGTg                                                                                                                                      >  3:171362/1‑149 (MQ=255)
                 tttCCTACGTTAGGGCGCCCGACAAGCGCGACCACAGGTACCATGTTTAAAGCCTCATTTTTATAAATCATCAGACAACGCACGCTATATTCGCGTCGTTGTTAAAAACAGGAAAACGGCCCCTGTCCAGGAGCCGTTTTCAAAGTGaa                                                                                                                                    <  4:470861/149‑1 (MQ=255)
                 tttCCTACGTTAGGGCGCCCGACAAGCGCGACCACAGGTACCATGTTTAAAGCCTCATTTTTATAAATCATCAGACAACGCACGCTATATTCGCGTCGTTGTTAAAAACAGGAAAACGGCCCCTGTCCAGGAGCCGTTTTCAAAGTGaa                                                                                                                                    <  3:444884/149‑1 (MQ=255)
                 tttCCTACGTTAGGGCGCCCGACAAGCGCGACCACAGGTACCATGTTTAAAGCCTCATTTTTATAAATCATCAGACAACGCACGCTATATTCGCGTCGTTGTTAAAAACAGGAAAACGGCCCCTGTCCAGGAGCCGTTCTCAAAGTGaa                                                                                                                                    <  2:196065/149‑1 (MQ=255)
                         gTTAGGGCGCCCGACAAGCGCGACCACAGGTACCATGTTTAAAGCCTCATTTTTATAAATCATCAGACAACGCACGCTATATTCGCGTCGTTGTTAAAAACAGGAAAACGGCCCCTGTCCAGGAGCCGTTTTCAAAGTGAACGACagag                                                                                                                            >  2:222032/1‑149 (MQ=255)
                          ttAGGGCGCCCGACAAGCGCGACCACAGGTACCATGTTTAAAGCCTCATTTTTATAAATCATCAGACAACGCACGCTATATTCGCGTCGTTGTTAAAAACAGGAAAACGGCCCCTGTCCAGGAGCCGTTTTCAAAGTGAACGACAgaga                                                                                                                           <  4:374436/149‑1 (MQ=255)
                          ttAGGGCGCCCGACAAGCGCGACCACAGGTACCATGTTTAAAGCCTCATTTTTATAAATCATCAGACAACGCACGCTATATTCGCGTCGTTGTTAAAAACAGGAAAACGGCCCCTGTCCAGGAGCCGTTTTCAAAGTGAACGACAgaga                                                                                                                           >  4:181086/1‑149 (MQ=255)
                           tAGGGCGCCCGACAAGCGCGACCACAGGTACCATGTTTAAAGCCTCATTTTTATAAATCATCAGACAACGCACGCTATATTCGCGTCGTTGTTAAAAACAGGAAAACGGCCCCTGTCCAGGAGCCGTTTTCAAAGTGAACGACAGAGAc                                                                                                                          <  4:171362/149‑1 (MQ=255)
                             gggCGCCCGACAAGCGCGACCACAGGTACCATGTTTAAAGCCTCATTTTTATAAATCATCAGACAACGCACGCTATATTCGCGTCGTTGTTAAAAACAGGAAAACGGCCCCTGTCCAGGAGCCGTTTTCAAAGTGAACGACAGAGACGa                                                                                                                        >  3:71806/1‑149 (MQ=255)
                              ggCGCCCGACAAGCGCGACCACAGGTACCATGTTTAAAGCCTCATTTTTATAAATCATCAGACAACGCACGCTATATTCGCGTCGTTGTTAAAAACAGGAAAACGGCCCCTGTCCAGGAGCCGTTTTCAAAGTGAACGACag                                                                                                                              <  2:303518/142‑1 (MQ=255)
                              ggCGCCCGACAAGCGCGACCACAGGTACCATGTTTAAAGCCTCATTTTTATAAATCATCAGACAACGCACGCTATATTCGCGTCGTTGTTAAAAACAGGAAAACGGCCCCTGTCCAGGAGCCGTTTTCAAAGTGAACGACag                                                                                                                              >  1:303518/1‑142 (MQ=255)
                                 gcCCGACAAGCGCGACCACAGGTACCATGTTTAAAGCCTCATTTTTATAAATCATCAGACAACGCACGCTATATTCGCGTCGTTGTTAAAAACAGGAAAACGGCCCCTGTCCAGGAGCCGTTTTCAAAGTGAACGACAGAGACGATTaa                                                                                                                    >  4:177368/1‑149 (MQ=255)
                                  cccGACAAGCGCGACCACAGGTACCATGTTTAAAGCCTCATTTTTATAAATCATCAGACAACGCACGCTATATTCGCGTCGTTGTTAAAAACAGGAAAACGGCCCCTGTCCAGGAGCCGTTTTCAAAGTGAACGACAGAGACGATTAAc                                                                                                                   <  3:520756/149‑1 (MQ=255)
                                      aCAAGCGCGACCACAGGTACCATGTTTAAAGCCTCATTTTTATAAATCATCAGACAACGCACGCTATATTCGCGTCGTTGTTAAAAACAGGAAAACGGCCCCTGTCCAGGAGCCGTTTTCAAAGTGAACGACAGAGACGATTAACgtgt                                                                                                               <  2:221459/149‑1 (MQ=255)
                                        aaGCGCGACCACAGGTACCATGTTTAAAGCCTCATTTTTATAAATCATCAGACAACGCACGCTATATTCGCGTCGTTGTTAAAAACAGGAAAACGGCCCCTGTCCAGGAGCCGTTTTCAAAGTGAACGACAGAGACGATTAACGTGTaa                                                                                                             <  2:390989/149‑1 (MQ=255)
                                          gcgcGACCACAGGTACCATGTTTAAAGCCTCATTTTTATAAATCATCAGACAACGCACGCTATATTCGCGTCGTTGTTAAAAACAGGAAAACGGCCCCTGTCCAGGAGCCGTTTTCAAAGTGAACGACAGAGACGATTAACGTGTAATa                                                                                                           <  3:265456/149‑1 (MQ=255)
                                           cgcgACCACAGGTACCATGTTTAAAGCCTCATTTTTATAAATCATCAGACAACGCACGCTATATTCGCGTCGTTGTTAAAAACAGGAAAACGGCCCCTGTCCAGGAGCCGTTTTCAAAGTGAACGACAGAGACGATTAACGTGTAATag                                                                                                          >  3:213017/1‑149 (MQ=255)
                                            gcgACCACAGGTACCATGTTTAAAGCCTCATTTTTATAAATCATCAGACAACGCACGCTATATTCGCGTCGTTGTTAAAAACAGGAAAACGGCCCCTGTCCAGGAGCCGTTTTCAAAGTGAACGACAGAGACGATTAACGTGTAATaga                                                                                                         <  3:448002/149‑1 (MQ=255)
                                              gACCACAGGTACCATGTTTAAAGCCTCATTTTTATAAATCATCAGACAACGCACGCTATATTCGCGTCGTTGTTAAAAACAGGAAAACGGCCCCTGTCCAGGAGCCGTTTTCAAAGTGAACGACAGAGACGATTAACGTGTAATAGAGt                                                                                                       >  2:230442/1‑149 (MQ=255)
                                                         ccATGTTTAAAGCCTCATTTTTATAAATCATCAGACAACGCACGCTATATTCGCGTCGTTGTTAAAAACAGGAAAACGGCCCCTGTCCAGGAGCCGTTTTCAAAGTGAACGACAGAGACGATTAACGTGTAATAGAGTACACGGTTCCg                                                                                            >  1:561871/1‑149 (MQ=255)
                                                           aTGTTTAAAGCCTCATTTTTATAATTCATCAGACAACGCACGCTATATTCGCGTCGTTGTTAAAAACAGGAAAACGGCCCCTGTCCAGGAGCCGTTTTCAAAGTGAACGACAGAGACGATTAACGTGTAATAGAGTACACGggtcggtc                                                                                          >  4:2727/1‑149 (MQ=255)
                                                              tttAAAGCCTCATTTTTATAAATCATCAGACAACGCACGCTATATTCGCGTCGTTGTTAAAAACAGGAAAACGGCCCCTGTCCAGGAGCCGTTTTCAAAGTGAACGACAGAGACGATTAACGTGTAATAGAGTACACGGTTCCGTCttt                                                                                       >  2:536963/1‑149 (MQ=255)
                                                                  aaGCCTCATTTTTATAAATCATCAGACAACGCACGCTATATTCGCGTCGTTGTTAAAAACAGGAAAACGGCCCCTGTCCAGGAGCCGTTTTCAAAGTGAACGACAGAGACGATTAACGTGTAATAGAGTACACGGTTCCGTCTTTTGcc                                                                                   <  3:388528/149‑1 (MQ=255)
                                                                  aaGCCTCATTTTTATAAATCATCAGACAACGCACGCTATATTCGCGTCGTTGTTAAAAACAGGAAAACGGCCCCTGTCCAGGAGCCGTTTTCAAAGTGAACGACAGAGACGATTAACGTGTAATAGAGTACACGGTTCCGTCTTTTGcc                                                                                   <  4:213017/149‑1 (MQ=255)
                                                                  aaGCCTCATTTTTATAAATCATCAGACAACGCACGCTATATTCGCGTCGTTGTTAAAAACAGGAAAACGGCCCCTGTCCAGGAGCCGTTTTCAAAGTGAACGACAGAGACGATTAACGTGTAATAGAGTACACGGTTCCGTCTTTTGcc                                                                                   <  4:456854/149‑1 (MQ=255)
                                                                      cTCATTTTTATAAATCATCAGACAACGCACGCTATATTCGCGTCGTTGTTAAAAACAGGAAAACGGCCCCTGTCCAGGAGCCGTTTTCAAAGTGAACGACAGAGACGATTAACGTGTAATAGAGTACACGGTTCCGTCTTTTGCCTGGa                                                                               >  4:534654/1‑149 (MQ=255)
                                                                       tCATTTTTATATATCATCAGACAACGCACGCTATATTCGCGTCGTTGTTAAAAACAGGAAAACGGCCCCTGTCCTGGAGTCGTTTTCAAAGTGAACGACAGAGACGATTAACGTGTAAAAGAGTACACGGTTCTGTCTTTTGCCTGGAt                                                                              <  4:361876/149‑1 (MQ=255)
                                                                         aTTTTTATAAATCATCAGACAACGCACGCTATATTCGCGTCGTTGTTAAAAACAGGAAAACGGCCCCTGTCCAGGAGCCGGTTTCAAAGTGAACGACAGAGACGATTAACGTGTAATAGAGTACACGGTTCGGTCTTTTGCCTGGATca                                                                            >  2:15948/1‑149 (MQ=255)
                                                                              tataAATCATCAGACAACGCACGCTATATTCGCGTCGTTGTTAAAAACAGGAAAACGGCCCCTGTCCAGGAGCCGTTTTCAAAGTGAACGACAGAGACGATTAACGTGTAATAGAGTACACGGTTCCGTCTTTTGCCTGGATCAGCAGt                                                                       >  1:413266/1‑149 (MQ=255)
                                                                                   atcatcAGACAACGCACGCTATATTCGCGTCGTTGTTAAAAACAGGAAAACGGCCCCTGTCCAGGAGCCGTTTTCAAAGTGAACGACAGAGACGATTAACGTGTAATAGAGTACACGGTTCCGTCTTTTGCCTGGATCAGCAGTTTGcc                                                                  <  4:71806/149‑1 (MQ=255)
                                                                                   atcatcAGACAACGCACGCTATATTCGCGTCGTTGTTAAAAACAGGAAAACGGCCCCTGTCCAGGAGCCGTTTTCAAAGTGAACGACAGAGACGATTAACGTGTAATAGAGTACACGGTTCCGTCTTTTGCCTGGATCAGCAGTTTGcc                                                                  <  3:194267/149‑1 (MQ=255)
                                                                                                  acgcTATATTCGCGTCGTTGTTAAAAACAGGAAAACGGCCCCTGTCCAGGAGCCGTTTTCAAAGTGAACGACAGAGACGATTAACGTGTAATAGAGTACACGGTTCCGTCTTTTGCCTGGATCAGCAGTTTGCCGTCAGCGGCAACCgg                                                   >  2:263881/1‑149 (MQ=255)
                                                                                                    gcTATATTCGCGTCGTTGTTAAAAACGGGAAAACGGCCCCTGTCCAGGAGCCGTTTTCAAAGTGAACGACAGAGCCGATTAACGTGTAATAGAGTACACGGTTCCGTCTTTTGCCTGGATCAGCAGTTTGCCGTCAGCGGCAACCGGtt                                                 >  3:362853/1‑149 (MQ=255)
                                                                                                    gcTATATTCGCGTCGTTGTTAAAAACAGGAAAACGGCCCCTGTCCAGGAGCCGTTTTCAAAGTGAACGACAGAGACGATTAACGTGTAATAGAGTACACGGTTCCGTCTTTTGCCTGGATCAGCAGTTTGCCGTCAGCGGCAACCGGtt                                                 >  4:465040/1‑149 (MQ=255)
                                                                                                    gcTATATTCGCGTCGTTGTTAAAAACAGGAAAACGGCCCCTGTCCAGGAGCCGTTTTCAAAGTGAACGACAGAGACGATTAACGTGTAATAGAGTACACGGTTCCGTCTTTTGCCTGGATCAGCAGTTTGCCGTCAGCGGCAACCGGtt                                                 <  3:181086/149‑1 (MQ=255)
                                                                                                      tataTTCGCGTCGTTGTTAAAAACAGGAAAACGGCCCCTGTCCAGGAGCCGTTTTCAAAGTGAACGACAGAGACGATTAACGTGTAATAGAGTACACGGTTCCGTCTTTTGCCTGGATCAGCAGTTTGCCGTCAGCGGCAACCGGTTCa                                               >  2:320776/1‑149 (MQ=255)
                                                                                                        tatTCGCGTCGTTGTTAAAAACAGGAAAACGGCCCCTGTCCAGGAGCCGTTTTCAAAGTGAACGACAGAGACGATTAACGTGTAATAGAGTACACGGTTCCGTCTTTTGCCTGGATCAGCAGTTTGCCGTCAGCGGCAACCGGTTCAGt                                             <  2:434823/149‑1 (MQ=255)
                                                                                                         atTCGCGTCGTTGTTAAAAACAGGAAAACGGCCCCTGTCCAGGAGCCGTTTTCAAAGTGAACGACAGAGACGATTAACGTGTAATAGAGTacac                                                                                                   >  1:369352/1‑94 (MQ=255)
                                                                                                         atTCGCGTCGTTGTTAAAAACAGGAAAACGGCCCCTGTCCAGGAGCCGTTTTCAAAGTGAACGACAGAGACGATTAACGTGTAATAGAGTacac                                                                                                   <  2:369352/94‑1 (MQ=255)
                                                                                                         atTCGCGTCGTTGTTAAAAACAGGAAAACGGCCCCTGTCCAGGAGCCGTTTTCAAAGTGAACGACAGAGACGATTAACGTGTAATAGAGTACACGGTTCCGTCTTTTGCCTGGATCAGCAGTTTGCCGTCAGCGGCAACCGGTTCAGTc                                            <  1:222032/149‑1 (MQ=255)
                                                                                                          ttCGCGTCGTTGTTAAAAACAGGAAAACGGCCCCTGTCCAGGAGCCGTTTTCAAAGTGAACGACAGAGACGATTAACGTGTAATAGAGTACACGGTTCCGTCTTTTGCCTGGATCAGCAGTTTGCCGTCAGCGGCAACCGGTTCAGTCt                                           <  3:374250/149‑1 (MQ=255)
                                                                                                          ttCGCGTCGTTGTTAAAAACAGGAAAACGGCCCCTGTCCAGGAGCCGTTTTCAAAGTGAACGACAGAGACGATTAACGTGTAATAGAGTACACGGTTCCGTCTTTTGCCTGGATCAGCAGTTTGCCGTCAGCGGCAACCGGTTCAGTCt                                           <  1:230442/149‑1 (MQ=255)
                                                                                                                  gttgttAAAAACAGGAAAACGGCCCCTGTCCAGGAGCCGTTTTCAAAGTGAACGACAGAGACGATTAACGTGTAATAGAGTACACGGTTCCGTCTTTTGCCTGGATCAGCAGTTTGCCGTCAGCGGCAACCGGTTCAGTCTGGAAACCg                                   >  1:154741/1‑149 (MQ=255)
                                                                                                                     gttAAAAACAGGAAAACGGCCCCTGTCCAGGAGCCGTTTTCAAAGTGAACGACAGAGACGATTAACGTGTAATAGAGTACACGGTTCCGTCTTTTGCCTGGATCAGCAGTTTGCCGTCAGCGGCAACCGGTTCAGTCTGGAAACCGGaa                                <  3:177368/149‑1 (MQ=255)
                                                                                                                     gttAAAAACAGGAAAACGGCCCCTGTCCAGGAGCCGTTTTCAAAGTGAACGACAGAGACGATTAACGTGTAATAGAGTACACGGTTCCGTCTTTTGCCTGGATCAGCAGTTTGCCGTCAGCGGCAACCGGTTCAGTCTGGAAACCGGaa                                <  1:452718/149‑1 (MQ=255)
                                                                                                                      ttAAAAACAGGAAAACGGCCCCTGTCCAGGAGCCGTTTTCAAAGTGAACGACAGAGACGATTAACGTGTAATAGAGTACACGGTTCCGTCTTTTGCCTGGATCAGCAGTTTGCCGTCAGCGGCAACCGGTTCAGTCTGGAAACCGGAAc                               >  1:151674/1‑149 (MQ=255)
                                                                                                                         aaaaCAGGAAAACGGCCCCTGTCCAGGAGCCGTTTTCAAAGTGAACGACAGAGACGATTAACGTGTAATAGAGTACACGGTTCCGTCTTTTGCCTGGATCAGCAGTTTGCCGTCAGCGGCAACCGGTTCAGTCTGGAAACCGGAACTAt                            >  3:237520/1‑149 (MQ=255)
                                                                                                                            aCAGGAAAACGGCCCCTGTCCAGGAGCCGTTTTCAAAGTGAACGACAGAGACGATTAACGTGTAATAGAGTACACGGTTCCGTCTTTTGCCTGGATCAGCAGTTTGCCGTCAGCGGCAACCGGTTCAGTCTGGAAACCGGAACTATCaa                         <  3:287318/149‑1 (MQ=255)
                                                                                                                            aCAGGAAAACGGCCCCTGTCCAGGAGCCGTTTTCAAAGTGAACGACAGAGACGATTAACGTGTAATAGAGTACACGGTTCCGTCTTTTGCCTGGATCAGCAGTTTGCCGTCAGCGGCAACCGGTTCAGTCTGGAAACCGGAACTATCaa                         >  2:69859/1‑149 (MQ=255)
                                                                                                                            aCAGGAAAACGGCCCCCGGACAGGAGCCGCTTTCAAAGTGAACGACAGCGTCGAGTAACCTGTAATAGAGTAGACGGTTCCGTCTTTTGCCTGGATCAGCAGTGTGCCGTCAGTGGCAACCGGTTCAGCCTGGAAACCGGGACtctcga                         <  2:99336/149‑3 (MQ=255)
                                                                                                                              aGGAAAACGGCCCCTGTCCAGGAGCCGTTTTCAAAGTGAACGACAGAGACGATTAACGTGTAATAGAGTACACGGTTCCGTCTTTTGCCTg                                                                                 >  2:60746/1‑91 (MQ=255)
                                                                                                                              aGGAAAACGGCCCCTGTCCAGGAGCCGTTTTCAAAGTGAACGACAGAGACGATTAACGTGTAATAGAGTACACGGTTCCGTCTTTTGCCTg                                                                                 <  1:60746/91‑1 (MQ=255)
                                                                                                                                 aaaaCGGCCCCTGTCCAGGAGCCGTTTTCAAAGTGAATGACAGAGACGATTAACGTGTAATa                                                                                                           <  3:7787/62‑1 (MQ=255)
                                                                                                                                 aaaaCGGCCCCTGTCCAGGAGCCGTTTTCAAAGTGAACGACAGAGACGATTAACGTGTAATa                                                                                                           >  4:7787/1‑62 (MQ=255)
                                                                                                                                  aaaCGGCTCCTGGACAGGGGCCGTTTTCAAAGTGAACGACAGAGACGATTAACGTGTAATAGAGTACACGGTTCCGTCTTTTGCCTGGATCAGCAGTTTGCCGTCAGCGGCAACCGGTTCAGTCTGGAAACCGGAACTATCAACttttt                   >  2:391261/1‑149 (MQ=255)
                                                                                                                                  aaaCGGCCCCTGTCCAGGAGCCGTTTTCAAAGTGAACGACAGAGACGATTAACGTGTAATAGAGTACACGGTTCCGTCTTTTGCCTGGATCAGCAGTTTGCCGTCAGCGGCAACCGGTTCAGTCTGGAAACCGGAACTATCAACttttt                   <  1:536963/149‑1 (MQ=255)
                                                                                                                                    aCGGCTCCTGGACAGGGGCCGTTTTCAAAGTGAACGACAGAGACGATTAACGTGTAATAGAGTACACGGTTCCGTCTTTTGCCTGGATCAGCAGTTTGCCGTCAGCGGCAACCGGTTCAGTCTGGAAACCGGAACTATCAACTTTTTGc                 >  2:406362/1‑149 (MQ=255)
                                                                                                                                    aCGGCTCCTGGACAGGGGCCGTTTTCAAAGTGAACGACAGAGACGATTAACGTGTAATAGAGTACACGGTTCCGTCTTTTGCCTGGATCAGCAGTTTGCCGTCAGCGGCAACCGGTTCAGTCTGGAAACCGGAACTATCAACTTTTTGc                 >  4:273000/1‑149 (MQ=255)
                                                                                                                                      ggCCCCTGTCCAGGAGCCGTTTTCAAAGTGAACGACAGAGACGATTAACGTGTAATAGAGTACACGGTTCCGTCTTTTGCCTGGATCAGCAGTTTGCCGTCAGCGGCAACCGGTTCAGTCTGGAAACCGGAACTATCAACTTTTTGCTg               <  4:211204/149‑1 (MQ=255)
                                                                                                                                       gCCCCTGTCCAGGAGCCGTTTTCAAAGTGAACGACAGAGACGATTAACGTGTAATAGAGTACACGGTTCCGTCTTTTGCCTGGATCAGCAGTTTGCCGTCAGCGGCAACCGGTTCAGTCTGGAAACCGGAACTATCAACTTTTTGCTgg              <  2:274092/149‑1 (MQ=255)
                                                                                                                                       gCCCCTGTCCAGGAGCCGTTTTCAAAGTGAACGACAGAGACGATTAACGTGTAATAGAGTACACGGTTCCGTCTTTTGCCTGGATCAGCAGTTTGCCGTCAGCGGCAACCGGTTCAGTCTGGAAACCGGAACTATCAACTTTTTGCTgg              <  1:497669/149‑1 (MQ=255)
                                                                                                                                         cccTGTCCAGGAACCGTTTTCAAAGTGAACGACAGAGACGATTAACGTGTAATAGAGTACACGGTTCCGTCTTTTGCCTGGATCAGCAGTTTGCCGTCAGCGGCAACCGGTTCAGTCTGGAAACCGGAACTATCAACTTTTTGCTGGGc            >  3:277433/1‑149 (MQ=255)
                                                                                                                                                cAGGAGCCGTTTTCAAAGTGAACGACAGAGACGATTAACGTGTAATAGAGTACACGGTTCCGTCTTTTGCCTGGATCAGCAGTTTGCCGTCAGCGGCAACCGGTTCAGTCTGGAAACCGGAACTATCAACTTTTTGCTGGGCaacgaaa     <  2:413266/149‑1 (MQ=255)
                                                                                                                                                   gAGCCGTTTTCAAAGTGAACGACAGAGACGATTAACGTGTAATAGAGTACACGGTTCCGTCTTTTGCCTGGATCAGCAGTTTGCCGTCAGCGGCAACCGGTTCAGTCTGGAAACCGGAACTATCAACTTTTTGCTGGGCaacgaaacga  <  1:320776/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
GGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGACCACAGGTACCATGTTTAAAGCCTCATTTTTATAAATCATCAGACAACGCACGCTATATTCGCGTCGTTGTTAAAAACAGGAAAACGGCCCCTGTCCAGGAGCCGTTTTCAAAGTGAACGACAGAGACGATTAACGTGTAATAGAGTACACGGTTCCGTCTTTTGCCTGGATCAGCAGTTTGCCGTCAGCGGCAACCGGTTCAGTCTGGAAACCGGAACTATCAACTTTTTGCTGGGCAACGAAACGA  >  NC_000913/2637296‑2637591

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