Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,101,655 T→A 9.4% intergenic (‑32/+152) yggN ← / ← yggL DUF2884 family putative periplasmic protein/DUF469 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,101,6550TA 9.4% 243.1 / 12.2 74intergenic (‑32/+152)yggN/yggLDUF2884 family putative periplasmic protein/DUF469 family protein
Reads supporting (aligned to +/- strand):  ref base T (37/30);  new base A (2/5);  total (39/35)
Fisher's exact test for biased strand distribution p-value = 2.45e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.49e-01

ACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGTAAACACCTTATCCCGCCTACTCGATCGTGCATCTTTCCAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCATTTGTCATCGCCCGT  >  NC_000913/3101508‑3101798
                                                                                                                                                   |                                                                                                                                               
aCGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCAtt                                                                                                                                                >  1:523990/1‑149 (MQ=255)
     ttGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTg                                                                                                                                           <  1:360650/149‑1 (MQ=255)
       gCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGc                                                                                                                                         <  3:347493/149‑1 (MQ=255)
       gCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCAAATGCTGGc                                                                                                                                         <  4:137410/149‑1 (MQ=255)
       gCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCAAATGCTGGc                                                                                                                                         <  4:212913/149‑1 (MQ=255)
        cGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGcc                                                                                                                                        >  3:496456/1‑149 (MQ=255)
          ggCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCAAATGCTGGcctg                                                                                                                                      <  1:396007/149‑3 (MQ=255)
          ggCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCAAATGCTGGcctg                                                                                                                                      <  2:738264/149‑3 (MQ=255)
             tGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGa                                                                                                                                   >  3:60500/1‑149 (MQ=255)
                 aaTCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGa                                                                                                                                   >  3:316809/1‑145 (MQ=255)
                 aaTCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGa                                                                                                                                   <  4:316809/145‑1 (MQ=255)
                  aTCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGAt                                                                                                                              <  1:70657/149‑1 (MQ=255)
                   tCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATa                                                                                                                             <  3:65334/149‑1 (MQ=255)
                       agcGTCACGCCGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACaa                                                                                                                         >  3:534412/3‑149 (MQ=255)
                        tCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACaaa                                                                                                                        <  4:648811/149‑1 (MQ=255)
                        tCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACaaa                                                                                                                        <  3:496359/149‑1 (MQ=255)
                         cGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAg                                                                                                                       <  3:467323/149‑1 (MQ=255)
                           tCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTa                                                                                                                     >  4:504589/1‑149 (MQ=255)
                             aCGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTAcg                                                                                                                   >  4:368006/1‑149 (MQ=255)
                                  gcgTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGttt                                                                                                              <  3:368006/149‑1 (MQ=255)
                                    gTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTc                                                                                                            >  1:291330/1‑149 (MQ=255)
                                    gTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGACCAGCATTGGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTc                                                                                                            <  3:438963/149‑1 (MQ=255)
                                     tGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTcc                                                                                                           <  4:322/149‑1 (MQ=255)
                                                tGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATgc                                                                                                >  4:575921/1‑149 (MQ=255)
                                                     gTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGt                                                                                           <  4:253934/149‑1 (MQ=255)
                                                      tCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGTa                                                                                          >  2:608082/1‑149 (MQ=255)
                                                         gCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGTAAac                                                                                       <  3:257898/149‑1 (MQ=255)
                                                               gCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGTAAACACCTTa                                                                                 <  3:575921/149‑1 (MQ=255)
                                                                 ggTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGTAAACACCTTATc                                                                               >  1:230705/1‑149 (MQ=255)
                                                                   tCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGTAAACACCTTATccc                                                                             >  1:722012/1‑149 (MQ=255)
                                                                      ttGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGTAAACACCTTATCCCGcc                                                                          >  3:416611/1‑149 (MQ=255)
                                                                          cGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCGTCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGTAAACACCTTATCCCGCCTACt                                                                      <  3:493723/149‑1 (MQ=255)
                                                                           gTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGTAAACACCTTAACCCGCCTACTc                                                                     <  3:248964/149‑1 (MQ=255)
                                                                             cACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGTAAACACCTTATCCCGCCTACTCGa                                                                   >  1:600224/1‑149 (MQ=255)
                                                                               cTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGTAAACACCTTATCCCGCCTACTCGATc                                                                 <  1:405113/149‑1 (MQ=255)
                                                                                tGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGTAAACACCTTATCCCGCCTACTCGATCg                                                                >  3:273522/1‑149 (MQ=255)
                                                                                 gAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGTAAACACCTTATCCCGCCTACTCGATCGt                                                               <  4:749143/149‑1 (MQ=255)
                                                                                 gAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGTAAACACCTTATCCCGCCTACTCGATCGt                                                               <  4:208959/149‑1 (MQ=255)
                                                                                 caaaaTGCTGCTGGCGGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGTAAACACCTTATCCCGCCTACTCGATCGt                                                               >  1:315929/2‑149 (MQ=255)
                                                                                  aaaaTGCTGCTGGCGGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGTAAACACCTTATCCCGCCTACTCGATCGTg                                                              >  4:28621/1‑149 (MQ=255)
                                                                                   aaGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGTAAACACCTTATCCCGCCTACTCGATCGTGc                                                             <  1:2251/149‑1 (MQ=255)
                                                                                   aaGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGTAAACACCTTATCCCGCCTACTCGATCGTGc                                                             <  3:735667/149‑1 (MQ=255)
                                                                                    aGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGTAAACACCTTATCCCGCCTACTCGATCGTGCa                                                            <  4:402945/149‑1 (MQ=255)
                                                                                       gctgcCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGTAAACACCTTATCCCGCCTACTCGATCGTGCATCt                                                         <  4:534412/149‑1 (MQ=255)
                                                                                         tgcCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTACATGCCGGATGCGCCGTAAACACCTTATCCCGCCTACTCGATCGTGCATCttt                                                       >  4:10939/1‑149 (MQ=255)
                                                                                          gccgccAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGTAAACACCTTATCCCGCCTACTCGATCGTGCATCTTTc                                                      >  2:82105/1‑149 (MQ=255)
                                                                                                   agcaTTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGTAAACACCTTATCCCGCCTACTCGATCGTGCATCTTTCCAGGCCTGa                                             <  3:336581/149‑1 (MQ=255)
                                                                                                    gcaTTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGTAAACACCTTATCCCGCCTACTCGATCGTGCATCTTTCCAGGCCTGAt                                            >  1:618000/1‑149 (MQ=255)
                                                                                                    gcaTTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGTAAACACCTTATCCCGCCTACTCGATCGTGCATCTTTCCAGGCCTGAt                                            <  4:103137/149‑1 (MQ=255)
                                                                                                             gcATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGTAAACACCTTATCCCGCCTACTCGATCGTGCATCTTTCCAGGCCTGATAAGCGTAgc                                   >  3:230696/1‑149 (MQ=255)
                                                                                                             gcATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGTAAACACCTTATCCCGCCTACTCGATCGTGCATCTTTCCAGGCCTGATAAGCGTAgc                                   >  1:300724/1‑149 (MQ=255)
                                                                                                                  atAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGTAAACACCTTATCCCGCCTACTCGATCGTGCATCTTTCCAGGCCTGATAAGCGTAGCGCATc                              >  1:200095/1‑149 (MQ=255)
                                                                                                                  atAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGTAAACACCTTATCCCGCCTACTCGATCGTGCATCTTTCCAGGCCTGATAAGCGTAGCGCATc                              <  1:579439/149‑1 (MQ=255)
                                                                                                                  atAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGTAAACACCTTATCCCGCCTACTCGATCGTGCATCTTTCCAGGCCTGATAAGCGTAGCGCATc                              >  2:579439/1‑149 (MQ=255)
                                                                                                                   tAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGTAAACACCTTATCCCGCCTACTCGATCGTGCATCTTTCCAGGCCTGATAAGCGTAGCGCATCa                             >  2:746528/1‑149 (MQ=255)
                                                                                                                   aaGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCTGATGCGGCGTAAACTCCTTATCCCGCCTACTCGATCGTGCATCTTTCCAGGCCTGATAAGCGTAGCCCATCa                             >  4:491093/2‑149 (MQ=255)
                                                                                                                      tCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGTAAACACCTTATCCCGCCTACTCGATCGTGCATCTTTCCAGGCCTGATAAGCGTAg                                    >  4:12556/1‑139 (MQ=255)
                                                                                                                      tCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGTAAACACCTTATCCCGCCTACTCGATCGTGCATCTTTCCAGGCCTGATAAGCGTAg                                    <  3:12556/139‑1 (MQ=255)
                                                                                                                        ttCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGTAAACACCTTAt                                                                                <  4:22939/93‑1 (MQ=255)
                                                                                                                        ttCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGTAAACACCTTAt                                                                                >  3:22939/1‑93 (MQ=255)
                                                                                                                        ttCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGTAAACACCTTATCCCGCCTACTCGATCGTGCATCTTTCCAGGCCTGATAAGCGTAGCGCATCAGGCAg                        >  3:175348/1‑149 (MQ=255)
                                                                                                                           ccTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGTAAACACCTTATCCCGCCTACTCGATCGTGCATCTTTCCAGGCCTGATAAGCGTAGCGCATCAGGCAGttt                     >  4:208343/1‑149 (MQ=255)
                                                                                                                            cTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTCTTCCATGCCGGATGCGGCGTAAACACCTTATCCCGCCTACTCGATCGTGCATCTTTCCAGGCCTGATAAGCGTAGCGCATCAGGCAGtttc                    >  3:68189/1‑148 (MQ=255)
                                                                                                                             tCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGTAAACACCTTATCCCGCCTACTCGATCGTGCATCTTTCCAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTg                   >  3:153516/1‑149 (MQ=255)
                                                                                                                                 gAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGTAAACACCTTATCCCGCCTACTCGATCGTGCATCTTTCCAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCAtt               >  3:199663/1‑149 (MQ=255)
                                                                                                                                 gAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGTAAACACCTTATCCCGCCTACTCGATCGTGCATCTTTCCAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCAtt               >  3:302640/1‑149 (MQ=255)
                                                                                                                                  aaaaaaTCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGTAAACACCTTATCCCGCCTACTCGATCGTGCATCTTTCCAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCAttt              <  3:338200/149‑1 (MQ=255)
                                                                                                                                     atctggggCCAGCAAATGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGTAAACACCTTATCCCGCCTACTCGATCGTGCATCTTTCCAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCATTTGTc           >  2:161943/7‑149 (MQ=255)
                                                                                                                                        tggggCCAGCAAATGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGTAAACACCTTATCCCGCCTACTCGATCGTGCATCTTTCCAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCATTTGt            >  4:37748/4‑145 (MQ=255)
                                                                                                                                        tggggCCAGCAAATGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGTAAACACCTTATCCCGCCTACTCGATCGTGCATCTTTCCAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCATTTGt            <  3:37748/142‑1 (MQ=255)
                                                                                                                                         cAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGTAAACACCTTATCCCGCCTACTCGATCGTGCATCTTTCCAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCATTTGTCATCg       >  4:626514/1‑149 (MQ=255)
                                                                                                                                              cAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGTAAACACCTTATCCCGCCTACTCGATCGTGCATCTTTCCAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCATTTGTCATCGCCCGt  >  2:161357/1‑149 (MQ=255)
                                                                                                                                              cAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGTAAACACCTTATCCCGCCTACTCGATCGTGCATCTTTCCAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCATTTGTCATCGCCCGt  <  2:523990/149‑1 (MQ=255)
                                                                                                                                              cAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGTAAACACCTTATCCCGCCTACTCGATCGTGCATCTTTCCAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCATTTGTCATCGCCCGt  >  4:607261/1‑149 (MQ=255)
                                                                                                                                              cAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGTAAACACCTTATCCCGCCTACTCGATCGTGCATCTTTCCAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCATTTGTCATCGCCCGt  >  1:224717/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                               
ACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGTAAACACCTTATCCCGCCTACTCGATCGTGCATCTTTCCAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCATTTGTCATCGCCCGT  >  NC_000913/3101508‑3101798

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