Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,249,431 T→A 5.9% intergenic (+120/‑123) lamB → / → malM maltose outer membrane porin (maltoporin)/maltose regulon periplasmic protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,249,4310TA 5.9% 210.2 / 2.9 77intergenic (+120/‑123)lamB/malMmaltose outer membrane porin (maltoporin)/maltose regulon periplasmic protein
Reads supporting (aligned to +/- strand):  ref base T (33/39);  new base A (2/3);  total (35/42)
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.24e-01

TCGGTGCCCAGATGGAAATCTGGTGGTAATAGCAAAACCTGGGCCGGATAAGGCGTTTACGCCGCATTCGGCAACCAACGCCTGATGCGACGCTTGCGCGTCTTATCAGGCCTACAACGGCTGTCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATAAAAACAGGTTGTCATTATCTGAAAGGGGCGAAAGCCCCTCTGATTATCGGGTTTAGCGCGCTATTGCCTGGCTACCGCTGAGCTCCAGATTTTGAGGTGAAAACAATGAAAATGAATAAAAGTCTCATCG  >  NC_000913/4249283‑4249578
                                                                                                                                                    |                                                                                                                                                   
tCGGTGCCCAGATGGAAATCTGGTGGTAATAGCAAAACCTGGGCCGGATAAGGCGTTTACGCCGCATTCGGCAACCAACGCCTGATGCGACGCTTGCGCGTCTTATCAGGCCTACAACGGCTGTCAAATGTAGGCCGGATAAGGCGttt                                                                                                                                                     >  1:277749/1‑149 (MQ=255)
  ggTGCCCAGAGGGAAATCTGGAGGTAATAGCAAAACCTGGGCCGGATAAGGCGTTTACGCCGCATTCGGCAACCAACGCCTGATGCGACGCTTGCTCGTCTTATCAGGCCTACAACGGCTGTCAAATGTAGGCCGGATAAGGCGTTTAc                                                                                                                                                   <  2:367692/149‑1 (MQ=255)
    tGCCCAGATGGAAATCTGGTGGTAATAGCAAAACCTGGGCCGGATAAGGCGTTTACGCCGCATTCGGCAACCAACGCCTGATGCGACGCTTGCGCGTCTTATCAGGCCTACAACGGCTGTCAAATGTAGGCCGGATAAGGCGTTTAcgc                                                                                                                                                 <  2:1113391/149‑1 (MQ=255)
     gCCCAGATGGAAATCTGGTGGTAATAGCAAAACCTGGGCCGGATAAGGCGTTTACGCCGCATTCGGCAACCAACGCCTGATGCGACGCTTGCGCGTCTTATCAGGCCTACAACGGCTGTCAAATGTAGGCCGGATAAGGCGTTTAcgcc                                                                                                                                                <  2:10872/149‑1 (MQ=255)
          gATGGAAATCTGGTGGTAATAGCAACACCTGGGCCGGATAAGGCGTTTACGCCGCATTCGGCAACCAACGCCTGCTGCGACGCTTGCGCGTCTTATCAGGCCTACAACGGCTGTCAAATGTAGGCCGGATGCGGCGTAAAcgccttatc                                                                                                                                           <  1:129540/149‑6 (MQ=25)
          gATGGAAATCTGGTGGTAATAGCAAAACCTGGGCCGGATAAGGCGTTTACGCCGCATTCGGCAACCAACGCCTGATGCGACGCTTGCGCGTCTTATCAGGCCTACAACGGCTGTCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATc                                                                                                                                           >  2:476114/1‑149 (MQ=255)
          gATGGAAATCTGGTGGTAATAGCAAAACCTGGGCCGGATAAGGCGTTTACGCCGCATTCGGCAACCAACGCCTGATGCGACGCTTGCGCGTCTTATCAGGCCTACAACGGCTGTCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATc                                                                                                                                           >  1:736497/1‑149 (MQ=255)
           aTGGAAATCTGGTGGTAATAGCAAAACCTGGGCCGGATAAGGCGTTTACGCCGCATTCGGCAACCAACGCCTGATGCGACGCTTGCGCGTCTTATCAGGCCTACAACGGCTGTCAAATGTAGGCCGGATGCGGCGTAAACGCCTTATcc                                                                                                                                          <  2:1050242/149‑1 (MQ=25)
            tGGAAATCTGGTGGTAAAAGCAAACCCTGGGCCGGATAAGGCGTTTACGCCGCATTCGGCAACCAACGCCTGATGCGACGCTTGCGCGTCTTATCAGGCCTACAACGGCTGTCAAATGTAGGCCGGATGCGGCGTAAACGCCTTATCCg                                                                                                                                         <  1:237743/149‑1 (MQ=25)
             ggAAATCTGGTGGTAATCGCAAAACCTGGGCCGGATAAGGCGTTTACGCCGCATTCGGCAACCAACGCCTGATGCGACGCTTGCGCGTCTTATCAGGCCTACAACGGCTGACAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCgg                                                                                                                                        <  1:828859/149‑1 (MQ=255)
             ggAAATCTGGTGGTAATAGCAAAACCTGGGCCGGATAAGGCGTTTACGCCGCATTCGGCAACCAACGCCTGATGCGACGCTTGCGCGTCTTATCAGGCCTACAACGGCTGTCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCgg                                                                                                                                        <  2:231763/149‑1 (MQ=255)
               aaaTCTGGTGGTAATAGCAAACCCTGGGCCGGATAAGGCGTTTACGCCGCATTCGGCAACCAACGCCTGATGCGACGCTTGCGCGTCTTATCAGGCCTACAACGGCTGTCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCa                                                                                                                                      >  1:67558/1‑149 (MQ=255)
                    tggtggTAATAGCAAAACCTGGGCCGGATAAGGCGTTTACGCCGCATTCGGCAACCAACGCCTGATGCGACGCTTGCGCGTCTTATCAGGCCTACAACGGCTGTCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATaaaa                                                                                                                                 >  1:677032/1‑149 (MQ=255)
                     ggtggtAATAGCAAAACCTGGGCCGGATCAGGCGGTTACGCCGCATTCGGCAACCAACGCCTGATGCGACGCTTGCGCGTCTTATCAGGCCTACAACGGCTGTCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATaaaaa                                                                                                                                <  1:668380/149‑1 (MQ=255)
                        ggtAATAGCAAAACCTGGGCCGGATAAGGCGTTTACGCCGCATTCGGCAACCAACGCCTGATGCGACGCTTGCGCGTCTTATCAGGCCTACAACGGCTGTCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATAAAAACAg                                                                                                                             <  1:38637/149‑1 (MQ=255)
                         gtAATAGCAAAACCTGGGCCGGATAAGGCGTTTACGCCGCATTCGGCAACCAACGCCTGATGCGACGCTTGCGCGTCTTATCAGGCCTACAACGGCTGTCAAATGTAGGCCGGATAAGGTGTTTACGCCGCATCCGGCATAAAAACAgg                                                                                                                            >  2:494685/1‑149 (MQ=255)
                         gtAAGAGCAAAACCCGGCCCCGATCAGGCGTTTACGCCGCATTCGGCAACCAACGCCTGATGCGAAGCTTGCGCGTCTTATCAGGCCTACAACGGCTGTCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATAAAAACAgg                                                                                                                            <  2:535083/149‑1 (MQ=21)
                              aGCAAAACCTGGGCCGGATAAGGCGTTTACGCCGCATTCGGCAACCAACGCCTGATGCGACGCTTGCGCGTCTTATCAGGCCTACAACGGCTGTCAAATGTAGGCCGGAGAAGGCGTTTACGCCGCATCCGGCATATAAACAGGTTGTc                                                                                                                       >  1:352931/1‑149 (MQ=255)
                               gCAAAACCTGGGCCGGATAAGGCGTTTACGCCGCATTCGGCAACCAACGCCTGATGCGACGCTTGCGCGTCTTATCAGGCCTACAACGGCTGTCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATAAAAACAGGTTGTCa                                                                                                                      >  1:67118/1‑149 (MQ=255)
                                cAAAACCTGGGCCGGATAAGGCGTTTACGCCGCATTCGGCAACCAACGCCTGATGCGACGCTTGCGCGTCTTATCAGGCCTACAACGGCTGTCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATAAAAACAGGTTGTCAt                                                                                                                     >  2:491085/1‑149 (MQ=255)
                                   aaCCTGGGCCGGATAAGGCGTTTACGCCGCATTCGGCAACCAACGCCTGATGCGACGCTTGCGCGTCTTATCAGGCCTACAACGGCTGTCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATAAAAACAGGTTGTCATTAt                                                                                                                  <  1:704358/149‑1 (MQ=255)
                                    aCCTGGGCCGGAGAAGGAGTTTACGCCGCATTCGGCAAACCTCGCCTGATGCGACGCTTGCGCGTCTTATCAGGCCTACAACGGCTGTCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATAAAAACAGGTTGTCATTATc                                                                                                                 <  1:1140851/149‑1 (MQ=25)
                                         ggCCGGATAAGGCGTTTACGCCGCATTCGGCAACCAACGCCTGATGCGACGCTTGCGCGTCTTATCAGGCCTACAACGGCTGTCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATAAAAACAGGTTGTCATTATCTGaaa                                                                                                            <  2:1116884/149‑1 (MQ=255)
                                          gCCGTATAAGTCGTTTACGCCGCATTCGGCAACCAACGCCTGATGCGTCGCTTGCGTGTCTTATCAGGCCTACAACGGCTGTCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATAAATACAGGTTGTCATTATCTGAAAg                                                                                                           <  2:1087219/149‑1 (MQ=25)
                                              gATAAGGCGTTTACGCCGCATTCGGCAACCAACGCCTGATGCGACGCTTGCGCGTCTTATCAGGCCTACAACGGCTGTCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATAAAAACAGGTTGTCATTATCTGAAAGGGGc                                                                                                       <  1:582444/149‑1 (MQ=255)
                                              gATAAGGCGTTTACGCCGCATTCGGCAACCAACGCCTGATGCGACGCTTGCGCGTCTTATCAGGCCTACAACGGCTGTCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATAAAAACAGGTTGTCATTATCTGAAAGGGGc                                                                                                       <  1:764188/149‑1 (MQ=255)
                                                 aaGGCGTTTACGCCGCATTCGGCAACCAACGCCTGATGCGACGCTTGCGCGTCTTATCAGGCCTACAACGGCTGTCAAATGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCATAAAAACAGGTTGTCATTATCTGAAAGGTGCGaa                                                                                                    >  2:717211/1‑149 (MQ=255)
                                                   ggCGTTTACGCCGCATTCGGCAACCAACGCCTGATGCGACGCTTGCGCGTCTTATCAGGCCTACAACGGCTGTCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATAAAAACAGGTTGTCATTATCTGAAAGGGGCGAAAg                                                                                                  <  2:526014/149‑1 (MQ=255)
                                                   ggCGTTTACGCCGCATTCGGCAACCAACGCCTGATGCGACGCTTGCGCGTCTTATCAGGCCTACAACGGCTGTCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATAAAAACAGGTTGTCATTATCTGAAAGGGGCGAAAg                                                                                                  <  2:808110/149‑1 (MQ=255)
                                                    gCGTTTACGCCGCATTCGGCAACCAACGCCTGATGCGACGCTTGCGCGTCTTATCAGGCCTACAACGGCTGTCAACTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATAAAAACAGGTTGTCATTATCTGAAAGGGGCGAAAGc                                                                                                 <  1:385176/149‑1 (MQ=255)
                                                      gTTTACGCCGCATTCGGCAACCAACGCCTGATGCGACGCTTGCGCGTCTTATCAGGCCTACAACGGCTGTCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATAAAAACAGGTTGTCATTATCTGAAAGGGGCGAAAGccc                                                                                               >  2:892402/1‑149 (MQ=255)
                                                        ttACGCCGCATTCGGCAACCAACGCCTGATGCGACGCTTGCGCGTCTTATCAGGCCTACAACGGCTGTCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATAAAAACAGGTTGTCATTATCTGAAAGGGGCGAAAGCCCct                                                                                             >  2:354188/1‑149 (MQ=255)
                                                          aCGCCGCATTCGGCAACCAACGCCTGATGCGACGCGTGCGCGTCTTATCAGGCCTACAACGGCTGTCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATAAAAACAGGTTGTCATTATCTGAAAGGGGCGAAAGCCCctct                                                                                           <  1:130613/149‑1 (MQ=255)
                                                            gccgcATTCGTCAACCAACGCCTGATGCGACGCTTGCGCGTCTTATCAGGCCTACAACGGCTGTCAAATGAAGGCCGGAAAAGGCGTTTACGCCGCATCCGGCATAAAAACAGGTTGTCATTATCAGAAAGGGGCGAACGCCCCTCTGa                                                                                         <  2:1005722/149‑1 (MQ=37)
                                                                  ttCGGCAACCAACGCCTCATGCGACGCTTGCGCGTCTTATGAGGCCTACAACGGCTGTCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATAAAAACAGGTTGTCATTATCTGAAAGGGGCGAAAGCGCCTCTGATTATCg                                                                                   <  2:795306/149‑1 (MQ=255)
                                                                   tCGGCAACCAACGCCTGATGCGACGCTTGCGCGTCTTATCAGGCCTACATCGGCTGTCAAATGTAGTCCGGATAAGGCGTTTACGCCGCATCCGGCATAAAAACAGGTTGTCATTATCTGAAAGGGGCGAAAGCCCCTCTGATTATCgg                                                                                  <  2:768843/149‑1 (MQ=255)
                                                                    cGGCAACCAACGCCTGATGCGACGCTTGCGCGTCTTATCAGGCCTACAACGGCTGTCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATAAAAACAGGTTGTCATTATCTGAAAGGGGCGAAAGCCCCTCTGATTATCggg                                                                                 >  1:1165331/1‑149 (MQ=255)
                                                                      gCACCCAACGCCTGATGCGACGCTTGCGCGTCTTATCAGGCCTACAACGGCTGTCAAATGTAGGCCGGATAAGGCGTTTACGCCGCCTCCGGCATAAAAACAGGTTGTCATTATCTGAAAGGGGCGAAAGCCCCGCTGATTATCGGGtt                                                                               >  1:1040853/1‑149 (MQ=255)
                                                                       caaccaacGCCTGATGCGACGCTTGCGCGTCTTATCAGGCCTACAACGGCTGTCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATAAAAACAGGTTGTCATTATCTGAAAGGGGCGAAAGCCCCTCTGATTATCGGGttt                                                                              >  2:719842/1‑149 (MQ=255)
                                                                                 cTGATGCGACGCTTGCGCGTCTTATCAGGCCTACAACGGCTGTCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATAAAAACAGGTTGTCATTATCTGAAAGGGGCGAAAGCCCCTCTGATTATCGGGTTTAGCGCGCTAt                                                                    >  1:987822/1‑149 (MQ=255)
                                                                                  tGATGCGACGCTTGCGCGTCTTATCAGGCCTACAACGGCTGTCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATAAAAACAGGTTGTCATTATCTGAAAGGGGCGAAAGCCCCTCTGATTATCGGGTTTAGCGCGCTAtt                                                                   <  2:364744/149‑1 (MQ=255)
                                                                                     tGCGACGCTTGCGCGTCTTATCAGGCCTACAACGGCTGTCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATAAAAACAGGTTGTCATTATCTGAAAGGGGCGAAAGCCCCTCTGATTATCGGGTTTAGCGCGCTATTGcc                                                                >  2:64763/1‑149 (MQ=255)
                                                                                      gctcagCTTGCGCGTCTTATCAGGCCTACAACGGCTGTCAAATGTATGCCGTATAAGGCGTTTACGCCGCATCGGGCATAAAAACAGGTTTTCATTATCTGAAAGGGGCGAATGCCCCTCTGATTATCGGTTTTAGCGCGCTATTGCCt                                                               <  1:883402/144‑1 (MQ=38)
                                                                                      gagACGCTTGCGCGTCTTTTCAGGCCTACAACGGCTGTCAAATGTAGGCCGGATAAGGCGTTTACGGCGCATCCGGCATAAAAACAGGTTGTCATTATCTGAAAGGGGCGAAAGCCCCTCTGATTATCGGGTTTAGCGCGCTATTGCCt                                                               <  1:878882/147‑1 (MQ=255)
                                                                                      gactcGCTTGCGCGTCTTATCAGCCCTACAACTGCTGTAAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATAAAAACAGGTTGTCATTATCTGAAAGGGGCGAAAGCCCCTCTGATTATCGGGTTTAGCGCGCTATTGCCt                                                               <  1:1181477/145‑1 (MQ=255)
                                                                                      gCGACGCTTGCGCGTCTTATCAGGCCTACAACGGCTGTCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATAAAAACAGGTTGTCATTATCTGAAAGGGGCGAAAGCCCCTCTGATTATCGGGTTTAGCGCGCTATTGCCt                                                               <  2:53378/149‑1 (MQ=255)
                                                                                       cGACGCATGCGCGTCGTATCAGGCCTACAACGGCTGTCAAAGGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATAAAAACAGGTTGTCATTATCTGAAAGGGGCGAAAGCCCCTCTGATTATCGGGTTTAGCGCGCTATTGCCTg                                                              <  1:1061141/149‑1 (MQ=255)
                                                                                         ccGCTTGCTCGACTTATCAAGCCTACAACGGCTGTCAAATGTAGCCCGGATTAGGCGTTTACGCCGCATCAGGGATAATAACAGGATGACATTATATGAAAGTGGCGAAAGACCCTCTGATTGTCGTGTTTAGCGCGCTATTGCCTGGc                                                            >  2:780836/2‑149 (MQ=255)
                                                                                                 gcgTCTTATCCGGCCTACAACGGCTGTCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATAAAAACAGGTTGTCATTATCTGAAAGGGGCGAAAGCCCCTCTTCTTATCGGGTTTAGCGCGCTATTGCCTGGCTACCGCTg                                                    >  1:92588/1‑149 (MQ=255)
                                                                                                 gcgTCTTATCAGGCCTACTACGGCTGTCAATTGTAGGCCGGATAATGCGTTTACGCCGCATCCGGCATAAAAACAGGTTGTCATTATCTGAAAGGGGCGAATGCCCCTCTGATTATCGGGTTTAGCGCGCTATTGCCTGGCTACCGCTg                                                    <  2:735950/149‑1 (MQ=255)
                                                                                                 gcgTCTTATCAGGCCTACAACGGCTGTCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATAAAAACAGGTTGTCATTATCTGAAAGGGGCGAAAGCCCCTCTGATTATCGGGTTTAGCGCGCTATTGCCTGGCTACCGCTg                                                    >  1:431419/1‑149 (MQ=255)
                                                                                                   gTCTTATCAGGCCTACAACGGCTGTCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATAAAAACAGGTTGTCATTATCTGAAAGGGGCGAAAGCCCCTCTGATTATCGGGTTTAGCGCGCTATTGCCTGGCTACCGCTGAg                                                  <  2:67558/149‑1 (MQ=255)
                                                                                                     cTTATCAGGCCTACAACGGCTGTCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATAAAAACAGGTTGTCATTATCTGAAAGGGGCGAAAGCCCCTCTGATTATCGGGTTTAGCGCGCTATTGCCTGGCTACCGCTGAGCt                                                >  2:1020546/1‑149 (MQ=255)
                                                                                                      ttATCAGGCCTACAACGGCTGTCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATAAAAACAGGTTGTCATTATCTGAAAGGGGCGAAAGCCCCTCTGATTATCGGGTTTAGCGCGCTATTGCCTGGCTACCGCTGAGCTc                                               >  2:292887/1‑149 (MQ=255)
                                                                                                        aTCAGGCCTACAACGGCTGTCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATAAAAACAGGTTGTCATTATCTGAAAGGGGCGAAAGCCCCTCTGATTTTCTGGTTTAGCGCGCTATTGCCTGGCTACCGCTGAGCTcaa                                             >  1:201617/1‑147 (MQ=255)
                                                                                                        aTCAGGCCTACAACGGCTGTCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATAAAAACAGGTTGTCATTATCTGAAAGGGGCGAAAGCCCCTCTGATTATCGGGTTTAGCGCGCTATTGCCTGGCTACCGCTGAGCTCCa                                             >  2:423140/1‑149 (MQ=255)
                                                                                                         tCAGGCCTACAACGGCTGTCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATAAAAACAGGTTGTCATTATCTGAAAGGGGCGAAAGCCCCTCTGATTATCGGGTTTAGCGCGCTATTGCCTGGCTACCGCTGCGCTCCAg                                            >  2:1131339/1‑149 (MQ=255)
                                                                                                          cAGGCCTACAACGGCTGTCAAATGTAGGACGGATAAGGCGTTTACGCCGCATCCGGCATAAAAACAGGTTGTCATTATCTGAAAGGGGCGAAAGCCCCTCTGATTATCGGGTTTAGCGCGCTATTGCCTGGCTACCGCTGAGCTCCAGa                                           <  2:204275/149‑1 (MQ=255)
                                                                                                           aGGCCTACAACGGCTGTCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATAAAAACAGGTTGTCATTATCTGAAAGGGGCGAAAGCCCCTCTGATTATCGGGTTTAGCGCGCTATTGCCTGGCTACCGCTGAGCTCCAGAt                                          >  2:116070/1‑149 (MQ=255)
                                                                                                                tACAACGGCTGTCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATAAAAACAGGTTGTCATTATCTGAAAGGGGCGAAAGCCCCTCTGATTATCGGGTTTAGCGCGCTATTGCCTGGCTACCGCTGAGCTCCAGATTTTGa                                     <  2:736497/149‑1 (MQ=255)
                                                                                                                 aCAACGGCTGTCAACTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATAAAAACAGGTTGTCATTATCTGAAAGGGGCGAAAGCCCCTCTGATTATCGGGTTTAGCGCGCTATTGCCTGGCTACCGCTGAGCTCCAGATTTTGAg                                    <  2:787734/149‑1 (MQ=255)
                                                                                                                  cAACGGCTGTCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATAAAAACAGGTTGTCATTATCTGAAAGGGGCGAAAGCCCCTCTGATTATCGGGTTTAGCGCGCTATTGCCTGGCTACCGCTGAGCTCCAGATTTTGAgg                                   >  1:636707/1‑149 (MQ=255)
                                                                                                                   aaCGGCTGTCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATAAAAACAGGTTGTCATTATCTGAAAGGGGCGAAAGCCCCTCTGATTATCGGGTTTAGCGCGCTATTGCCTGGCTACCGCTGAGCTCCAGATTTTGAGGt                                  >  1:1111735/1‑149 (MQ=255)
                                                                                                                   aaCGGCTGTCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATAAAAACAGGTTGTCATTATCTGAAAGGGGCGAAAGCCCCTCTGATTATCGGGTTTAGCGCGCTATTGCCTGGCTACCGCTGAGCTCCAGATTTTGAGGt                                  <  1:446185/149‑1 (MQ=255)
                                                                                                                      ggCTGTCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATAAAAACAGGTTGTCATTATCTGAAAGGGGCGAAAGCCCCTCTGATTATCGGGTTTAGCGCGCTATTGCCTGGCTACCGCTGAGCTCCAGATTTTGAGGTGaa                               <  1:482107/149‑1 (MQ=255)
                                                                                                                           tCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATAAAAACAGGTTGTCATTATCTGAAAGGGGCGAAAGCCCCTCTGATTATCGGGTTTAGCGCGCTATTGCCTGGCTACCGCTGAGCTCCAGATTTTGAGGTGAAAACaa                          <  1:476114/149‑1 (MQ=255)
                                                                                                                           gcAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATAAAAACAGGTTGTCATTATCTGAAAGGGGCGAAAGCCCCTCTGATTATCGGGTTTAGCGCGCTATTGCCTGGCTACCGCTGAGCTCCAGATTTTGAGGTGAAAACaa                          <  2:780761/148‑1 (MQ=255)
                                                                                                                                tGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATAAAAACAGGTTGTCATTATCTGAAAGGGGCGAAAGCCCCTCTGATTATCGGGTTTAGCGCGCTATTGCCTGGCTACCGCTGAGCTCCAGATTTTGAGGTGAAAACaatgaaa                     <  2:704870/149‑1 (MQ=255)
                                                                                                                                 gTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATAAAAACAGGTTGTCATTATCTGAAAGGGGCGAAAGCCCCTCTGATTATCGGGTTTAGCGCGCTATTGCCTGGCTACCGCTGAGCTCCAGATTTTGAGGTGAAAACaatgaaaa                    >  1:624526/1‑149 (MQ=255)
                                                                                                                                  tatgCCGGATGCGGCGTAAACGCCTTATCCGGCATAAAAACAGGTTGTCATTATCTGAAAGGGGCGAAAGCCCCTCTGATTATCGGGTTTAGCGCGCTATTGCCTGGCTACCGCTGAGCTCCAGATTTTGAGGTGAAAACaatgaaaat                   >  1:149339/4‑149 (MQ=255)
                                                                                                                                  tAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATAAAACCAGGGTGTCATTATCTGACAGGGGCGAAAGCCCCTCTGATTATCGGGTTTAGCGCGCTATTGCCTGGCTACCGCTGAGCTCCAGATTTTGAGGTGAAAACaatgaaaat                   <  1:811572/149‑1 (MQ=255)
                                                                                                                                  tAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATAAAAACAGGTTGTCATTATCTGAAAGGGGCGAAAGCCCCTCTGATTATCGGGTTTAGCGCGCTATTGCCTGGCTACCGCTGAGCTCCAGATTTTGAGGTGAAAACaatgaaaat                   >  2:835307/1‑149 (MQ=255)
                                                                                                                                  tAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATAAAAACAGGTTGTCATTATCTGAAAGGGGCGAAAGCCCCTCTGATTATCGGGTTTAGCGCGCTATTGCCTGGCTACCGCTGAGCTCCAGATTTTGAGGTGAAAACaatgaaaat                   >  2:1184283/1‑149 (MQ=255)
                                                                                                                                    ggCCGGATAAGGCGTTTACGCCGCATCCGGCATAAAAACAGGTTGTCATTATCTGAAAGGGGCGAAAGCCCCTCTGATTATCGGGTTTAGCGCGCTATTGCCTGGCTACCGCTGAGCTCCAGATTTTGAGGTGAAAACaatgaaaatga                 <  1:719842/149‑1 (MQ=255)
                                                                                                                                              ggCGTTTACGCCGCATCCGGCATAAAAACAGGTTGTCATTATCTGAAAGGGGCGAAAGCCCCTCTGATTATCGGGTTTAGCGCGCTATTGCCTGGCTACCGCTGAGCTCCAGATTTTGAGGTGAAAACAATGAAAATGATTAAAAGtct       >  1:516210/1‑149 (MQ=255)
                                                                                                                                              ggCGTTTACGCCGCATCCGGCATAAAAACAGGTTGTCATTATCTGAAAGGGGCGAAAGCCCCTCTGATTATCGGGTTTAGCGCGCTATTGCCTGGCTACCGCTGAGCTCCAGATTTTGAGGTGAAAACAATGAAAATGAATAAAAGtct       >  2:445679/1‑149 (MQ=255)
                                                                                                                                               gCGTTTACGCCGCATCCGGCATAAAAACAGGTTGTCATTATCTGAAAGGGGCGAAAGCCCCTCTGATTATCGGGTTTAGCGCGCTATTGCCTGGCTACCGCTGAGCTCCAGATTTTGAGGTGAAACCAATGAAAATGAATAAAAGtctc      >  2:69924/1‑149 (MQ=255)
                                                                                                                                                cGTTTACGCCGCATCCGGCATAAAAACAGGTTGTCATTATCTGAAAGGGGCGAAAGCCCCTCTGATTATCGGGTTTAGCGCGCTATTGCCTGGCTACCGCTGAGCTCCAGATTTTGAGGTGAAAACAATGAAAATGAATAAAAGTCTCa     >  1:507845/1‑149 (MQ=255)
                                                                                                                                                   ttACGCCGCATCCGGCATAAAAACAGGTTGTCATTATCTGAAAGGGGCGAAAGCCCCTCTGATTATCGGGTTTAGCGCGCTATTGCCTGGCTACCGCTGAGCTCCAGATTTTGAGGTGAAAACAATGAAAATGAATAAAAGTCTCATCg  <  2:277749/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
TCGGTGCCCAGATGGAAATCTGGTGGTAATAGCAAAACCTGGGCCGGATAAGGCGTTTACGCCGCATTCGGCAACCAACGCCTGATGCGACGCTTGCGCGTCTTATCAGGCCTACAACGGCTGTCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATAAAAACAGGTTGTCATTATCTGAAAGGGGCGAAAGCCCCTCTGATTATCGGGTTTAGCGCGCTATTGCCTGGCTACCGCTGAGCTCCAGATTTTGAGGTGAAAACAATGAAAATGAATAAAAGTCTCATCG  >  NC_000913/4249283‑4249578

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑