Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,355,187 A→G 5.0% V95A (GTG→GCG)  sapB ← antimicrobial peptide transport ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,355,1870AG 5.0% 730.1 / 14.4 224V95A (GTG→GCG) sapBantimicrobial peptide transport ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base A (100/112);  new base G (4/8);  total (104/120)
Fisher's exact test for biased strand distribution p-value = 3.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.34e-01

Only 100 of 228 total aligned reads displayed.

TGAAACGGGCAACCAACCTAGCGTGAGTGAACAAAACAGCGTCAACAGAAGCGCCAGCCAGAACACCGGGATTGAAAAGCCCAGTAAGGCGATGGCGTTGATCAAGTTATCCTGCCACTTATGGCGCGTAATCCCGGCAATCATCCCCACGGGGATGCCAACAATCAGCGCAAAACCGAAGGCGAGGATGCACAGTTCCATCGTCGCCGGGAAGACCTCTTTTAACTGTTCGGCAATCGGTTGACCGTTAATGCTGGAAACACCAAAATCCCAGTGGATCAGGCCGTTAAACCAGAA  >  NC_000913/1355039‑1355335
                                                                                                                                                    |                                                                                                                                                    
       ggCAACCAACCTAGCGTGAGTGAACAAAACAGCGTCAACAGAAGCGCCAGCCAGAACACCGGGATTGAAAAGCCCAGTAAGGCGATGGCGTTGATCAAGTTATCCTGCCACTTATGGCGCGTAATCCCGGCAATCATCCCCACGGGGAt                                                                                                                                               >  4:815934/1‑149 (MQ=255)
       ggCAACCAACCTAGCGTGAGTGAACAAAACAGCGTCAACAGAAGCGCCAGCCAGAACACCGGGATTGAAAAGCCCAGTAAGGCGATGGCGTTGATCAAGTTATCCTGCCACTTATGGCGCGTAATCCCGGCAATCATCCCCACGGGGAt                                                                                                                                               <  1:397295/149‑1 (MQ=255)
        gCAACCAACCTAGCGTGAGTGAACAAAACAGCGTCAACAGAAGCGCCAGCCAGAACACCGGGATTGAAAAGCCCAGTAAGGCGATGGCGTTGATCAAGTTATCCTGCCACTTATGGCGCGTAATCCCGGCAATCATCCCCACGGGGATg                                                                                                                                              >  1:1229674/1‑149 (MQ=255)
        gCAACCAACCTAGCGTGAGTGAACAAAACAGCGTCAACAGAAGCGCCAGCCAGAACACCGGGATTGAAAAGCCCAGTAAGGCGATGGCGTTGATCAAGTTATCCTGCCACTTATGGCGCGTAATCCCGGCAATCATCCCCACGGGGATg                                                                                                                                              <  2:384958/149‑1 (MQ=255)
        gCAACCAACCTAGCGTGAGTGAACAAAACAGCGTCAACAGAAGCGCCAGCCAGAACACCGGGATTGAAAAGCCCAGTAAGGCGATGGCGTTGATCAAGTTATCCTGCCACTTATGGCGCGTAATCCCGGCAATCATCCCCACGGGGATg                                                                                                                                              >  3:766044/1‑149 (MQ=255)
         caaccaacCTAGCGTGAGTGAACAAAACAGCGTCAACAGAAGCGCCAGCCAGAACACCGGGATTGAAAAGCCCAGTAAGGCGATGGCGTTGATCAAGTTATCCTGCCACTTATGGCGCGTAATCCCGGCAATCATCCCCGTGGGGATga                                                                                                                                             <  2:899327/149‑2 (MQ=255)
          aaccaaccTAGCGTGAGTGAACAAAACAGCGTCAACAGAAGCGCCAGCCAGAACACCGGGATTGAAAAGCCCAGTAAGGCGATGGCGTTGATCAAGTTATCCTGCCACTTATGGCGCGTAATCCCGGCAATCATCCCCGTGGGGATgat                                                                                                                                            <  1:1270709/149‑3 (MQ=255)
          aaccaaccTAGCGTGAGTGAACAAAACAGCGTCAACAGAAGCGCCAGCCAGAACACCGGGATTGAAAAGCCCAGTAAGGCGATGGCGTTGATCAAGTTATCCTGCCACTTATGGCGCGTAATCCCGGCAATCATCCCCACGGGGATGcc                                                                                                                                            <  4:1058155/149‑1 (MQ=255)
           accaaccTAGCGTGAGTGAACAAAACAGCGTCAACAGAAGCGCCAGCCAGAACACCGGGATTGAAAAGCCCAGTAAGGCGATGGCGTTGATCAAGTTATCCTGCCACTTATGGCGCGTAATCCCGGCAATCATCCCCGTGGGGATgatt                                                                                                                                           <  1:1321822/149‑4 (MQ=255)
           accaaccTAGCGTGAGTGAACAAAACAGCGTCAACAGAAGCGCCAGCCAGAACACCGGGATTGAAAAGCCCAGTAAGGCGATGGCGTTGATCAAGTTATCCTGCCACTTATGGCGCGTAATCCCGGCAATCATCCCCACGGGGATGCca                                                                                                                                           >  4:361771/1‑149 (MQ=255)
           accaaccTAGCGTGAGTGAACAAAACAGCGTCAACAGAAGCGCCAGCCAGAACACCGGGATTGAAAAGCCCAGTAAGGCGATGGCGTTGATCAAGTTATCCTGCCACTTATGGCGCGTAATCCCGGCAATCATCCCCACGGGGATGCca                                                                                                                                           <  2:980820/149‑1 (MQ=255)
               accTAGCGTGAGTGAACAAAACAGCGTCAACAGAAGCGCCAGCCAGAACACCGGGATTGAAAAGCCCAGTAAGGCGATGGCGTTGATCAAGTTATCCTGCCACTTATGGCGCGTAATCCCGGCAATCATCCCCACGGGGATGCcaacaa                                                                                                                                       >  4:667997/1‑149 (MQ=255)
               accTAGCGTGAGTGAACAAAACAGCGTCAACAGAAGCGCCAGCCAGAACACCGGGATTGAAAAGCCCAGTAAGGCGATGGCGTTGATCAAGTTATCCTGCCACTTATGGCGCGTAATCCCGGCAATCATCCCCACGGGGATGCcaacaa                                                                                                                                       <  3:461853/149‑1 (MQ=255)
                ccTAGCGTGAGTGAACAAAACAGCGTCAACAGAAGCGCCAGCCAGAACACCGGGATTGAAAAGCCCAGTAAGGCGATGGCGTTGATCAAGTTATCCTGCCACTTATGGCGCGTAATCCCGGCAATCATCCCCACGGGGATGCCAACAAt                                                                                                                                      <  4:1417357/149‑1 (MQ=255)
                 cTAGCGTGAGTGAACAAAACAGCGTCAACAGAAGCGCCAGCCAGAACACCGGGATTGAAAAGCCCAGTAAGGCGATGGCGTTGATCAAGTTATCCTGCCACTTATGGCGCGTAATCCCGGCAATCATCCCCACGGGGATGCCAACAATc                                                                                                                                     <  2:735821/149‑1 (MQ=255)
                  tAGCGTGAGTGAACAAAACAGCGTCAACAGAAGCGCCAGCCAGAACACCGGGATTGAAAAGCCCAGTAAGGCGATGGCGTTGATCAAGTTATCCTGCCACTTATGGCGCGTAATCCCGGCAATCATCCCCACGGGGATGCCAACAATCa                                                                                                                                    <  3:1543059/149‑1 (MQ=255)
                           tgaACAAAACAGCGTCAACAGAAGCGCCAGCCAGAACACCGGGATTGAAAAGCCCAGTAAGGCGATGGCGTTGATCAAGTTATCCTGCCACTTATGGCGCGTAATCCCGGCAATCATCCCCACGGGGATGCCAACAATCAGCGCAAAAc                                                                                                                           >  1:943499/1‑149 (MQ=255)
                             aaCAAAACAGCGTCAACAGAAGCGCCAGCCAGAACACCGGGATTGAAAAGCCCAGTAAGGCGATGGCGTTGATCAAGTTATCCTGCCACTTATGGCGCGTAATCCCGGCAATCATCCCCACGGGGATGCCAACAATCAGCGCAAAACCg                                                                                                                         >  3:898398/1‑149 (MQ=255)
                              aCAAAACAGCGTCAACAGAAGCGCCAGCCAGAACACCGGGATTGAAAAGCCCAGTAAGGCGATGGCGTTGATCAAGTTATCCTGCCACTTATGGCGCGTAATCCCGGCAATCATCCCCACgg                                                                                                                                                   <  3:1658516/122‑1 (MQ=255)
                              aCAAAACAGCGTCAACAGAAGCGCCAGCCAGAACACCGGGATTGAAAAGCCCAGTAAGGCGATGGCGTTGATCAAGTTATCCTGCCACTTATGGCGCGTAATCCCGGCAATCATCCCCACGGGGATGCCAACAATCAGCGCAAAACCGa                                                                                                                        >  3:1482195/1‑149 (MQ=255)
                              aCAAAACAGCGTCAACAGAAGCGCCAGCCAGAACACCGGGATTGAAAAGCCCAGTAAGGCGATGGCGTTGATCAAGTTATCCTGCCACTTATGGCGCGTAATCCCGGCAATCATCCCCACGGGGATGCCAACAATCAGCGCAAAACCGa                                                                                                                        >  4:899993/1‑149 (MQ=255)
                               cAAAACAGCGTCAACAGAAGCGCCAGCCAGAACACCGGGATTGAAAAGCCCAGTAAGGCGATGGCGTTGATCAAGTTATCCTGCCACTTATGGCGCGTAATCCCGGCAATCATCCCCACGGGGATGCCAACAATCAGCGCAAAACCGaa                                                                                                                       <  3:596159/149‑1 (MQ=255)
                               cAAAACAGCGTCAACAGAAGCGCCAGCCAGAACACCGGGATTGAAAAGCCCAGTAAGGCGATGGCGTTGATCAAGTTATCCTGCCACTTATGGCGCGTAATCCCGGCAATCATCCCCACGGGGATGCCAACAATCAGCGCAAAACCGaa                                                                                                                       <  3:668129/149‑1 (MQ=255)
                               cAAAACAGCGTCAACAGAAGCGCCAGCCAGAACACCGGGATTGAAAAGCCCAGTAAGGCGATGGCGTTGATCAAGTTATCCTGCCACTTATGGCGCGTAATCCCGGCAATCATCCCCACGGGGATGCCAACAATCAGCGCAAAACCGaa                                                                                                                       <  3:84511/149‑1 (MQ=255)
                                aaaaCAGCGTCAACAGAAGCGCCAGCCAGAACACCGGGATTGAAAAGCCCAGTAAGGCGATGGCGTTGATCAAGTTATCCTGCCACTTATGGCGCGTAATCCCGGCAATCATCCCCACGGGGATGCCAACAATCAGCGCAAAACCGAAg                                                                                                                      >  3:1600459/1‑149 (MQ=255)
                                aaaaCAGCGTCAACAGAAGCGCCAGCCAGAACACCGGGATTGAAAAGCCCAGTAAGGCGATGGCGTTGATCAAGTTATCCTGCCACTTATGGCGCGTAATCCCGGCAATCATCCCCACGGGGATGCCAACAATCAGCGCAAAACCGAAg                                                                                                                      >  1:1348910/1‑149 (MQ=255)
                                 aaaCAGCGTCAACAGAAGCGCCAGCCAGAACACCGGGATTGAAAAGCCCAGTAAGGCGATGGCGTTGATCAAGTTATCCTGCCACTTATGGCGCGTAATCCCGGCAATCATCCCCACGGGGATGCCAACAATCAGCGCAAAACCGAAgg                                                                                                                     <  3:199809/149‑1 (MQ=255)
                                    cAGCGTCAACAGAAGCGCCAGCCAGAACACCGGGATTGAAAAGCCCAGTAAGGCGATGGCGTTGATCAAGTTATCCTGCCACTTATGGCGCGTAATCCCGGCAATCATCCCCACGGGGATGCCAACAATCAGCGCAAAACCGAAGGCGa                                                                                                                  <  2:1437619/149‑1 (MQ=255)
                                      gCGTCAACAGAAGCGCCAGCCAGAACACCGGGATTGAAAAGCCCAGTAAGGCGATGGCGTTGATCAAGTTATCCTGCCACTTATGGCGCGTAATCCCGGCAATCATCCCCACGGGGATGCCAACAATCAGCGCAAAACCGAAGGCGAgg                                                                                                                <  4:295316/149‑1 (MQ=255)
                                        gTCAACAGAAGCGCCAGCCAGAACACCGGGATTGAAAAGCCCAGTAAGGCGATGGCGTTGATCAAGTTATCCTGCCACTTATGGCGCGTAATCCCGGCAATCATCCCCACGGGGATGCCAACAATCAGCGCAAAACCGAAGGCGAGGAt                                                                                                              <  3:1178202/149‑1 (MQ=255)
                                           aaCAGAAGCGCCAGCCAGAACACCGGGATTGAAAAGCCCAGTAAGGCGATGGCGTTGATCAAGTTATCCTGCCACTTATGGCGCGTAATCCCGGCAATCATCCCCACGGGGATGCCAACAATCAGCGCAAAACCGAAGGCGAGGATGca                                                                                                           <  1:1159174/149‑1 (MQ=255)
                                             cAGAAGCGCCAGCCAGAACACCGGGATTGAAAAGCCCAGTAAGGCGATGGCGTTGATCAAGTTATCCTGCCACTTATGGCGCGTAATCCCGGCAATCATCCCCACGGGGATGCCAACAATCAGCGCAAAACCGGAGGCGAGGATGcaca                                                                                                         <  4:659127/149‑1 (MQ=255)
                                              aGAAGCGCCAGCCAGAACACCGGGATTGAAAAGCCCAGTAAGGCGATGGCGTTGATCAAGTTATCCTGCCACTTATGGCGCGTAATCCCGGCAATCATCCCCACGGGGATGCCAACAATCAGCGCAAAACCGAAGGCGAGGATGCACAg                                                                                                        >  2:1196694/1‑149 (MQ=255)
                                                 aGCGCCAGCCAGAACACCGGGATTGAAAAGCCCAGTAAGGCGATGGCGTTGATCAAGTTATCCTGCCACTTATGGCGCGTAATCCCGGCAATCATCCCCACGGGGATGCCAACAATCAGCGCAAAACCGAAGGCGAGGATGCACAGTTc                                                                                                     <  2:1327873/149‑1 (MQ=255)
                                                  gcgcCAGCCAGAACAGCGGGATTGAAAAGCCCAGTAAGGCGATGGCGTTGATCAAGTTATCCTGCCACTTATGGCGCGTAATCCCGGCAATCATCCCCGCGGGGATGCCAACAATCAGCGCAAAACCGAAGGCGAGGATGCACAGTTcc                                                                                                    <  1:1008654/149‑1 (MQ=255)
                                                  gcgcCAGCCAGAACACCGGGATTGAAAAGCCCAGTAAGGCGATGGCGTTGATCAAGTTATCCTGCCACTTATGGCGCGTAATCCCGGCAATCATCCCCACGGGGATGCCAACAATCAGCGCAAAACCGAAGGCGAGGATGCACAGTTcc                                                                                                    <  3:333974/149‑1 (MQ=255)
                                                  gcgcCAGCCAGAACACCGGGATTGAAAAGCCCAGTAAGGCGATGGCGTTGATCAAGTTATCCTGCCACTTATGGCGCGTAATCCCGGCAATCATCCCCACGGGGATGCCAACAATCAGCGCAAAACCGAAGGCGAGGATGCACAGTTcc                                                                                                    <  1:1697495/149‑1 (MQ=255)
                                                    gccagccaGAACACCGGGATTGAAAAGCCCAGTAAGGCGATGGCGTTGATCAAGTTATCCTGCCACTTATGGCGCGTAATCCCGGCAATCATCCCCACGGGGATGCCAACAATCAGCGCAAAACCGAAGGCGAGGATGCACAGTTCCAt                                                                                                  >  2:31815/1‑149 (MQ=255)
                                                     ccagccagAACACCGGGATTGAAAAGCCCAGTAAGGCGATGGCGTTGATCAAGTTATCCTGCCACTTATGGCGCGTAATCCCGGCAATCATCCCCACGGGGATGCCAACAATCAGCGCAAAACCGAAAGCGAGGATGCACAGTTCCAtc                                                                                                 <  2:79159/149‑1 (MQ=255)
                                                       agccagAACACCGGGATTGAAAAGCCCAGTAAGGCGATGGCGTTGATCAAGTTATCCTGCCACTTATGGCGCGTAATCCCGGCAATCATCCCCACGGGGATGCCAACAATCAGCGCAAAACCGAAGGCGAGGATGCACAGTTCCAtcgt                                                                                               >  3:1413641/1‑149 (MQ=255)
                                                        gccagAACACCGGGATTGAAAAGCCCAGTAAGGCGATGGCGTTGATCAAGTTATCCTGCCACTTATGGCGCGTAATCCCGGCAATCATCCCCACGGGGATGCCAACAATCAGCGCAAAACCGAAGGCGAGGATGCACAGTTCCAtcgtc                                                                                              >  4:60305/1‑149 (MQ=255)
                                                                acCGGGATTGAAAAGCCCAGTAAGGCGATGGCGTTGATCAAGTTATCCTGCCACTTATGGCGCGTAATCCCGGCAATCATCCCCACGGGGATGCCAACAATCAGCGCAAAACCGAAGGCGAGGATGCACAGTTCCATCGTCGCCGGGaa                                                                                      >  2:1134833/1‑149 (MQ=255)
                                                                       ttGAAAAGCCCAGTAAGGCGATGGCGTTGATCAAGTTATCCTGCCACTTATGGCGCGTAATCCCGGCAATCATCCCCACGGGGATGCCAACAATCAGCGCAAAACCGAAGGCGAGGATGCACAGTTCCATCGTCGCCGGGAAGACctct                                                                               >  3:477268/1‑149 (MQ=255)
                                                                        tGAAAAGCCCAGTAAGGCGATGGCGTTGATCAAGTTATCCTGCCACTTATGGCGCGTAATCCCGGCAATCATCCCCACGGGGATGCCAACAATCAGCGCAAAACCGAAGGCGAGGATGCACAGTTCCATCGTCGCCGGGAAGACCTCtt                                                                              <  3:1716391/149‑1 (MQ=255)
                                                                          aaaaGCCCAGTAAGGCGATGGCGTTGATCAAGTTATCCTGCCACTTATGGCGCGTAATCCCGGCAATCATCCCCACGGGGATGCCAACAATCAGCGCAAAACCGAAGGCGAGGATGCACAGTTCCATCGTCGCCGGGAAGACCTCtttt                                                                            >  1:1419872/1‑149 (MQ=255)
                                                                          aaaaGCCCAGTAAGGCGATGGCGTTGATCAAGTTATCCTGCCACTTATGGCGCGTAATCCCGGCAATCATCCCCACGGGGATGCCAACAATCAGCGCAAAACCGAAGGCGAGGATGCACAGTTCCATCGTCGCCGGGAAGACCTCtttt                                                                            >  4:550174/1‑149 (MQ=255)
                                                                            aaGCCCAGTAAGGCGATGGCGTTGATCAAGTTATCCTGCCACTTATGGCGCGTAATCCCGGCAATCATCCCCACGGGGATGCCAACAATCAGCGCAAAACCGAAGGCGAGGATGCACAGTTCCATCGTCGCCGGGAAGACCTCTTTTaa                                                                          <  1:245353/149‑1 (MQ=255)
                                                                             aGCCCAGTAAGGCGATGGCGTTGATCAAGTTATCCTGCCACTTATGGCGCGTAATCCCGGCAATCATCCCCACGGGGATGCCAACAATCAGCGCAAAACCGAAGGCGAGGATGCACAGTTCCATCGTCGCCGGGAAGACCTCTTTTAAc                                                                         <  3:550174/149‑1 (MQ=255)
                                                                              gCCCAGTAAGGCGATGGCGTTGATCAAGTTATCCTGCCACTTATGGCGCGTAATCCCGGCAATCATCCCCACGGGGATGCCAACAATCAGCGCAAAACCGAAGGCGAGGATGCACAGTTCCATCGTCGCCGGGAAGACCTCTTTTAACt                                                                        >  4:1716484/1‑149 (MQ=255)
                                                                                ccAGTAAGGCGATGGCGTTGATCAAGTTATCCTGCCACTTATGGCGCGTAATCCCGGCAATCATCCCCACGGGGATGCCAACAATCAGCGCAAAACCGAAGGCGAGGATGCACAGTTCCATCGTCGCCGGGAAGACCTCTTTTAACTGt                                                                      >  1:112192/1‑149 (MQ=255)
                                                                                 cAGTAAGGCGATGGCGTTGATCAAGTTATCCTGCCACTTATGGCGCGTAATCCCGGCAATCATCCCCACGGGGATGCCAACAATCAGCGCAAAACCGAAGGCGAGGATGCACAGTTCCATCGTCGCCGGGAAGACCTCTTTTAACTGtt                                                                     <  1:848527/149‑1 (MQ=255)
                                                                                  aGTAAGGCGATGGCGTTGATCAAGTTATCCTGCCACTTATGGCGCGTAATCCCGGCAATCATCCCCACGGGGATGCCAACAATCAGCGCAAAACCGAAGGCGAGGATGCACAGTTCCATCGTCGCCGGGAAGACCTCTTTTAACTGTTc                                                                    <  4:456132/149‑1 (MQ=255)
                                                                                    tAAGGCGATGGCGTTGATCAAGTTATCCTGCCACTTATGGCGCGTAATCCCGGCAATCATCCCCACGGGGATGCCAACAATCAGCGCAAAACCGAAGGCGAGGATGCACAGTTCCATCGTCGCCGGGAAGACCTCTTTTAACTGTTCgg                                                                  <  3:60305/149‑1 (MQ=255)
                                                                                      aGGCGATGGCGTTGATCAAGTTATCCTGCCACTTATGGCGCGTAATCCCGGCAATCATCCCCACGGGGATGCCAACAATCAGCGCAAAACCGAAGGCGAGGATGCACAGTTCCATCGTCGCCGGGAAGACCTCTTTTAACTGTTCGGCa                                                                >  2:673951/1‑149 (MQ=255)
                                                                                         cGATGGCGTTGATCAAGTTATCCTGCCACTTATGGCGCGTAATCCCGGCAATCATCCCCACGGGGATGCCAACAATCAGCGCAAAACCGAAGGCGAGGATGCACAGTTCCATCGTCGCCGGGAAGACCTCTTTTAACTGTTCGGCAATc                                                             <  3:1394555/149‑1 (MQ=255)
                                                                                             ggCGTTGATCAAGTTATCCTGCCACATATGGCGCGTAATCCCGGCCATCTTCCCCACGGGCATGCCAACCACCAGCGCAAAACCGTAGGCGAGGATGCCCCATTCCATCGTCGCGGGGCACACCTCTTTTAACTATTCGGCAATCGGtt                                                         <  4:404377/149‑1 (MQ=255)
                                                                                                gTTGATCAAGTTATCCTGCCACTTATGGCGCGTAATCCCGGCAATCATCCCCACGGGGATGCCAACAATCAGCGCAAAACCGAAGGCGAGGATGCACAGTTCCATCGTCGCCGGGAAGAc                                                                                   >  4:1668688/1‑120 (MQ=255)
                                                                                                gTTGATCAAGTTATCCTGCCACTTATGGCGCGTAATCCCGGCAATCATCCCCACGGGGATGCCAACAATCAGCGCAAAACCGAAGGCGAGGATGCACAGTTCCATCGTCGCCGGGAAGACCTCTTTTAACTGTTCGGCAATCGGTTGAc                                                      <  1:539715/149‑1 (MQ=255)
                                                                                                 ttGATCAAGTTATCCTGCCACTTATGGCGCGTAATCCCGGCAATCATCCCCACGGGGATGCCAACAATCAGCGCAAAACCGAAGGCGAGGATGCACAGTTCCATAGTCGCCGGGAAGACCTCTTTTAACTGTTCGGCAATCGGTTGAcc                                                     <  4:1061788/149‑1 (MQ=255)
                                                                                                   gATCAAGTTATCCTGCCACTTATGGCGCGTAATCCCGGCAATCATCCCCACGGGGATGCCAACAATCAGCGCAAAACCGAAGGCGAGGATGCACAGTTCCATCGTCGCCGGGAAGACCTCTTTTAACTGTTCGGCAATCGGTTGACCGt                                                   >  1:462193/1‑149 (MQ=255)
                                                                                                   gATCAAGTTATCCTGCCACTTATGGCGCGTAATCCCGGCAATCATCCCCACGGGGATGCCAACAATCAGCGCAAAACCGAAGGCGAGGATGCACAGTTCCATCGTCGCCGGGAAGACCTCTTTTAACTGTTCGGCAATCGGTTGACCGt                                                   >  3:392228/1‑149 (MQ=255)
                                                                                                   gATCAAGTTATCCTGCCACTTATGGCGCGTAATCCCGGCAATCATCCCCACGGGGATGCCAACAATCAGCGCAAAACCGAAGGCGAGGATGCACAGTTCCATCGTCGCCGGGAAGACCTCTTTTAACTGTTCGGCAATCGGTTGACCGt                                                   >  2:664815/1‑149 (MQ=255)
                                                                                                     tCAAGTTATCCTGCCACTTATGGCGCGTAATCCCGGCAATCATCCCCACGGGGATGCCAACAATCAGCGCAAAACCGAAGGCGAGGATGCACAGTTCCATCGTCGCCGGGAAGAc                                                                                   <  4:1477656/115‑1 (MQ=255)
                                                                                                      cAAGTTATCCTGCCACTTATGGCGCGTAATCCCGGCAATCATCCCCACGGGGATGCCAACAATCAGCGCAAAACCGAAGGCGAGGATGCACAGTTCCATCGTCGCCGGGAAGACCTCTTTTAACTGTTCGGCAATCGGTTGACCGTTaa                                                <  1:676352/149‑1 (MQ=255)
                                                                                                       aaGTTATCCTGCCACTTATGGCGCGTAATCCCGGCAATCATCCCCACGGGGATGCCAACAATCAGCGCAAAACCGAAGGCGAGGATGCACAGTTCCATCGTCGCCGGGAAGACCTCTTTTAACTGTTCGGCAATCGGTTGACCGTTAAt                                               >  2:124193/1‑149 (MQ=255)
                                                                                                        aGTTATCCTGCCACTTATGGCGCGTAATCCCGGCAATCATCCCCACGGGGATGCCAACAATCAGCGCAAAACCGAAGGCGAGGATGCACAGTTCCATCGTCGCCGGGAAGACCTCTTTTAACTGTTCGGCAATCGGTTGACCGTTAATg                                              >  3:664412/1‑149 (MQ=255)
                                                                                                         gTTATCCTGCCACTTATGGCGCGTAATCCCGGCAATCATCCCCACGGGGATGCCAACAATCAGCGCAAAACCGAAGGCGAGGATGCACAGTTCCATCGTCGCCGGGAAGACCTCtt                                                                              >  4:776011/1‑116 (MQ=255)
                                                                                                         gTTATCCTGCCACTTATGGCGCGTAATCCCGGCAATCATCCCCACGGGGATGCCAACAATCAGCGCAAAACCGAAGGCGAGGATGCACAGTTCCATCGTCGCCGGGAAGACCTCTTTTAACTGTTCg                                                                   <  3:939543/127‑1 (MQ=255)
                                                                                                          ttATCCTGCCACTTATGGCGCGTAATCCCGGCAATCATCCCCACGGGGATGCCAACAATCAGCGCAAAACCGAAGGCGAGGATGCACAGTTCCATCGTCGCCGGGAAGACCTCTTTTAACTGTTCGGCAATCGGTTGACCGTTAATGCt                                            <  2:464830/149‑1 (MQ=255)
                                                                                                           tATCCTGCCACTTATGGCGCGTAATCCCGGCAATCATCCCCACGGGGATGCCAACAATCAGCGCAAAACCGAAGGCGAGGATGCACAGTTCCATCGTCGCCGGGAAGACCTCTTTTAACTGTTCGGCAATCGGTTGACCGTTAATGCTg                                           <  3:1269438/149‑1 (MQ=255)
                                                                                                            aTCCTGCCACTTATGGCGCGTAATCCCGGCAATCATCCCCACGGGGATGCCAACAATCAGCGCAAAACCGAAGGCGAGGATGCACAGTTCTATCGTCGCCGGGAAGACCTCTTTTAACTGTTCGGCAATCGGTTGACCGTTAATGCTgg                                          <  3:815934/149‑1 (MQ=255)
                                                                                                            aTCCTGCCACTTATGGCGCGTAATCCCGGCAATCATCCCCACGGGGATGCCAACAATCAGCGCAAAACCGAAGGCGAGGATGCACAGTTCCATCGTCGCCGGGAAGACCTCTTTTAACTGTTCGGCTATCGGTTGACCGTTAATGCTgg                                          <  1:1066803/149‑1 (MQ=255)
                                                                                                            aTCCTGCCACTTATGGCGCGTAATCCCGGCAATCATCCCCACGGGGATGCCAACAATCAGCGCAAAACCGAAGGCGAGGATGCACAGTTCCATCGTCGCCGGGAAGACCTCTTTTAACTGTTCGGCAATCGGTTGACCGTTAATGCTgg                                          <  3:1081033/149‑1 (MQ=255)
                                                                                                                tGCCACTTATGGCGCGTAATCCCGGCAATCATCCCCACGGGGATGCCAACAATCAGCGCAAAACCGAAGGCGAGGATGCACAGTTCCATCGTCGCCGGGAAGACCTCTTTTAACTGTTCGGCAATCGGTTGACCGTTAATGCTGGAAac                                      >  1:202917/1‑149 (MQ=255)
                                                                                                                    aCTTATGGCGCGTAATCCCGGCAATCATCCCCACGGGGATGCCAACAATCAGCGCAAAACCGAAGGCGAGGATGCACAGTTCCATCGTCGCCGGGAAGACCTCTTTTAACTGTTCGGCAATCGGTTGACCGTTAATGCTGGAAaca                                     <  4:1543722/146‑1 (MQ=255)
                                                                                                                      ttATGGCGCGTAATCCCGGCAATCATCCCCACGGGGATGCCAACAATCAGCGCAAAACCGAAGGCGAGGATGCACAGTTCCATCGTCGCCGGGAAGACCTCTTTTAACTGTTCGGCAATCGGTTGACCGTTAATGCTGGAAACACCaaa                                >  1:1745530/1‑149 (MQ=255)
                                                                                                                      ttATGGCGCGTAATCCCGGCAATCATCCCCACGGGGATGCCAACAATCAGCGCAAAACCGAAGGCGAGGATGCACAGTTCCATCGTCGCCGGGAAGACCTCTTTTAACTGTTCGGCAATCGGTTGACCGTTAATGCTGGAAACACCaaa                                >  2:1370514/1‑149 (MQ=255)
                                                                                                                          ggCGCGTAATCCCGGCAATCATCCCCACGGGGATGCCAACAATCAGCGCAAAACCGAAGGCGAGGATGCACAGTTCCATCGTCGCCGGGAAGACCTCTTTTAACTGTTCGGCAATCGGTTGACCGTTAATGCTGGAAACACCAAAATcc                            <  2:604256/149‑1 (MQ=255)
                                                                                                                            cgcgTAATCCCGGCAATCATCCCCACGGGGATGCCAACAATCAGCGCAAAACCGAAGGCGAGGATGCACAGTTCCATCGTCGCCGGGAAGACCTCTTTTAACTGTTCGGCAATCGGTTGACCGTTAATGCTGGAAACACCAAAATCCCa                          >  1:180619/1‑149 (MQ=255)
                                                                                                                            cgcgTAATCCCGGCAATCATCCCCACGGGGATGCCAACAATCAGCGCAAAACCGAAGGCGAGGATGCACAGTTCCATCGTCGCCGGGAAGACCTCTTTTAACTGTTCGGCAATCGGTTGACCGTTAATGCTGGAAACACCAAAATCCCa                          <  4:38524/149‑1 (MQ=255)
                                                                                                                            cgcgTAATCCCGGCAATCATCCCCACGGGGATGCCAACAATCAGCGCAAAACCGAAGGCGAGGATGCACAGTTCCATCGTCGCCGGGAAGACCTCTTTTAACTGTTCGGCAATCGGTTGACCGTTAATGCTGGAAACACCAAAATCCCa                          <  2:180619/149‑1 (MQ=255)
                                                                                                                            cgcgTAATCCCGGCAATCATCCCCACGGGGATGCCAACAATCAGCGCAAAACCGAAGGCGAGGATGCACAGTTCCATCGTCGCCGGGAAGACCTCTTTTAACTGTTCGGCAATCGGTTGACCGTTAATGCTGGAAACACCAAAATCCCa                          <  2:1625668/149‑1 (MQ=255)
                                                                                                                             gcgTAATCCCGGCAATCATCCCCACGGGGATGCCAACAATCAGCGCAAAACCGAAGGCGAGGATGCACAGTTCCATCGTCGCCGGGAAGACCTCTTTTAACTGTTCGGCAATCGGTTGACCGTTAATGCTGGAAACACCAAAATCCCAg                         >  4:369716/1‑149 (MQ=255)
                                                                                                                               gTAATCCCGGCAATCATCCCCACGGGGATGCCAACAATCAGCGCAAAACCGAAGGCGAGGATGCACAGTTCCATCGTCGCCGGGAAGACCTCTTTTAACTGTTCg                                                                   <  4:1212123/105‑1 (MQ=255)
                                                                                                                               gTAATCCCGGCAATCATCCCCACGGGGATGCCAACAATCAGCGCAAAACCGAAGGCGAGGATGCACAGTTCCATCGTCGCCGGGAAGACCTCTTTTAACTGTTCg                                                                   >  3:1212123/1‑105 (MQ=255)
                                                                                                                                  aTCCCGGCAATCATCCCCACGGGGATGCCAACAATCAGCGCAAAACCGAAGGCGAGGATGCACAGTTCCATCGTCGCCGGGAAGACCTCTTTTAACTGTTCGGCAATCGGTTGACCGTTAATGCTGGAAACACCAAAATCCCAGTGGAt                    >  3:169583/1‑149 (MQ=255)
                                                                                                                                   tCCCGGCAATCATCCCCACGGGGATGCCAACAATCAGCGCAAAACCGAAGGCGAGGATGCACAGTTCCATCGTCGCCGGGAAGACCTCTTTTAACTGTTCGGCAATCGGTTGACCGTTAATGCTGGAAACACCAAAATCCCa                          >  2:1399414/1‑142 (MQ=255)
                                                                                                                                    cccGGCAATCATCCCCACGGGGATGCCAACAATCAGCGCAAAACCGAAGGCGAGGATGCACAGTTCCATCGTCGCCGGGAAGACCTCTTTTAACTGTTCGGCAATCGGTTGACCGTTAATGCTGGAAACACCAAAATCCCAGTGGATCa                  <  4:1699598/149‑1 (MQ=255)
                                                                                                                                     ccGGCAATCATCCCCACGGGGATGCCAACAATCAGCGCAAAACCGAAGGCGAGGATGCACAGTTCCAt                                                                                                  <  3:393228/68‑1 (MQ=255)
                                                                                                                                     ccGGCAATCATCCCCACGGGGATGCCAACAATCAGCGCAAAACCGAAGGCGAGGATGCACAGTTCCATCGTCGCCGGGAAGACCTCTTTTAACTGTTCGGCAATCGGTTGACCGTTAATGCTGGAAACACCAAAATCCCAGTGGATCAg                 >  3:1041192/1‑149 (MQ=255)
                                                                                                                                      cGGCAATCATCCCCACGGGGATGCCAACAATCAGCGCAAAACCGAAGGCGAGGATGCACAGTTCCATCGTCGCCGGGAAGACCTCTTTTAACTGTTCGGCAATCGGtt                                                         <  4:1112656/108‑1 (MQ=255)
                                                                                                                                        gCAATCATCCCCACGGGGATGCCAACAATCAGCGCAAAACCGAAGGCGAGGATGCACAGTTCCATCGTCGCCGGGAAGACCTCTTTTAACTGTTCGGCAATCGGTTGACCGTTAATGCTGGAAACACCAAAATCCCAGTGGATCAGGcc              <  2:197220/149‑1 (MQ=255)
                                                                                                                                          tggcatcCCCGTGGGGATGCCAACAATCAGCGCAAAACCGAAGGCGAGGATGCACAGTTCCATCGTCGCCGGGAAGACCTCTTTTAACTGTTCGGCAATCGGTTGACCGTTAATGCTGGAAACACCAAAATCCCAGTGGATCAGGCCGt            >  1:669448/4‑149 (MQ=255)
                                                                                                                                             catcCCCACGGGGATGCCAACAATCAGCGCAAAACCGAAGGCGAGGATGCACAGTTCCATCGTCGCCGGGAAGACCTCTTTTAACTGTTCGGCAATCGGTTGACCGTTAATGCTGGAAACACCAAAATCCCAGTGGATCAGGCCGTTaa         <  4:700409/149‑1 (MQ=255)
                                                                                                                                                  ccACGGGGATGCCAACAATCAGCGCAAAACCGAAGGCGAGGATGCACAGTTCCATCGTCGCCGGGAAGACctc                                                                                >  3:1630015/1‑73 (MQ=255)
                                                                                                                                                  ccACGGGGATGCCAACAATCAGCGCAAAACCGAAGGCGAGGATGCACAGTTCCATCGTCGCCGGGAAGACctc                                                                                <  4:1630015/73‑1 (MQ=255)
                                                                                                                                                  ccACGGGGATGCCAACAATCAGCGCAAAACCGAAGGCGAGGATGCACAGTTCCATCGTCGCCGGGAAGACCTCTTTTAACTGTTCGGCAATCGGTTGACCGTTAATGCTGGAAACACCAAAATCCCAGTGGATCAGGCCGTTAAACCAg    >  2:279888/1‑149 (MQ=255)
                                                                                                                                                  ccACGGGGATGCCAACAATCAGCGCAAAACCGAAGGCGAGGATGCACAGTTCCATCGTCGCCGGGAAGACCTCTTTTAACTGTTCGGCAATCGGTTGACCGTTAATGCTGGAAACACCAAAATCCCAGTGGATCAGGCCGTTAAACCAg    >  2:149337/1‑149 (MQ=255)
                                                                                                                                                   cACGGGGATGCCAACAATCAGCGCAAAACCGAAGGCGAGGATGCACAGTTCCATCGTCGCCGGGAAGACCTCTTTAACTGTTTGGCAATCGGTTGAACGGTAAAGCTGGAAAAACCAAAATCCCAGTGGATCAGGCCGGTAAACCAGaa  >  1:1657288/1‑149 (MQ=255)
                                                                                                                                                    aCGGGGATGCCAACAATCAGCGCAAAACCGAAGGCGAGGATGCACAGTTCCATCGTCGCCGGGAAGACCTCTTTTAACTGTTCGGCAATCGGTTGACCGTTAATGCTGGAAACACCAAAATCCCAGTGGATCAGGCCGTTAAACCAGaa  <  2:693636/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
TGAAACGGGCAACCAACCTAGCGTGAGTGAACAAAACAGCGTCAACAGAAGCGCCAGCCAGAACACCGGGATTGAAAAGCCCAGTAAGGCGATGGCGTTGATCAAGTTATCCTGCCACTTATGGCGCGTAATCCCGGCAATCATCCCCACGGGGATGCCAACAATCAGCGCAAAACCGAAGGCGAGGATGCACAGTTCCATCGTCGCCGGGAAGACCTCTTTTAACTGTTCGGCAATCGGTTGACCGTTAATGCTGGAAACACCAAAATCCCAGTGGATCAGGCCGTTAAACCAGAA  >  NC_000913/1355039‑1355335

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