Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,952,594 C→T 6.0% A110V (GCG→GTG)  yecN → MAPEG family inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,952,5940CT 6.0% 262.2 / 4.5 84A110V (GCG→GTG) yecNMAPEG family inner membrane protein
Reads supporting (aligned to +/- strand):  ref base C (37/41);  new base T (2/3);  total (40/44)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.52e-01

GCTGTTTATGGAAATGAATGGCGCAGAAACCTGGATGGTGCATATTTGCGGCATCGTTTTGCTTGCTGGTCGTCTGATGCATTATTACGGTTTTCATCACCGTCTGTTCCGCTGGCGACGTTCCGGCATGAGCGCCACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTTGGTATATGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAGCGCACAATACGCCACTTTCTTTTTCCCGGATTTTTACGTTATGTCTCACCGCGACACGCTATTTTCTGCCCCTATCGC  >  NC_000913/1952448‑1952739
                                                                                                                                                  |                                                                                                                                                 
gCTGTTTATGGAAATGAATGGCGCAGAAACCTGGATGGTGCATATTTGCGGCATCGTTTTGCTTGCTGGTCGTCTGATGCATTATTACGGTTTTCATCACCGTCTGTTCCGCTGGCGACGTTCCGGCATGAGCGCCACCTGGTGTgcgc                                                                                                                                                 >  7:181563/1‑149 (MQ=255)
   gTTTATGGAAATGAATGGCGCAGAAACCTGGATGGTGCATATTTGCGGCATCGTTTTGCTTGCTGGTCGTCTGATGCATTATTACGGTTTTCATCACCGTCTGTTCCGCTGGCGACGTTCCGGCATGAGCGCCACCTGGTGTGCGCTGt                                                                                                                                              >  8:186978/1‑149 (MQ=255)
    tttATGGAAATGAATGGCGCAGAAACCTGGATGGTGCATATTTGCGGCATCGTTTTGCTTGCTGGTCGTCTGATGCATTATTACGGTTTTCATCACCGTCTGTTCCGCTGGCGACGTTCCGGCATGAGCGCCACCTGGTGTGCGCTGtt                                                                                                                                             <  8:142559/149‑1 (MQ=255)
       aTGGAAATGAATGGCGCAGAAACCTGGATGGTGCATATTTGCGGCATCGTTTTGCTTGCTGGTCGTCTGATGCATTATTACGGTTTTCATCACCGTCTGTTCCGCTGGCGACGTTCCGGCATGAGCGCCACCTGGTGTGCGCTGTTGCt                                                                                                                                          >  3:164504/1‑149 (MQ=255)
       aTGGAAATGAATGGCGCAGAAACCTGGATGGTGCATATTTGCGGCATCGTTTTGCTTGCTGGTAGTCTGATGCATTATTACGGTTTTCATCACCGTCTGTTCCGCTGGCGACGTTCCGGCATGAGCGCAACCTGTTGTGAGATGTTGCt                                                                                                                                          >  4:223994/1‑149 (MQ=255)
           aaaTGAATGGCGCAGAAACCTGGATGGTGCATATTTGCGGCATCGTTTTGCTTGCTGGTCGTCTGATGCATTATTACGGTTTTCATCACCGTCTGTTCCGCTGGCGACGTTCCGGCATGAGCGCCACCTGGTGTGTGCTGTTGCTGATg                                                                                                                                      <  1:64193/149‑1 (MQ=255)
                aatgGCGCAGAAACCTGGATGGTGCATATTTGCGGCATCGTTTTGCTTGCTGGTCGTCTGATGCATTATTACGGTTTTCATCACCGTCTGTTCCGCTGGCGACGTTCCGGCATGAGCGCCCCCTGGTGTGCGCTGTTGCTGATGGTGCt                                                                                                                                 >  6:116584/1‑149 (MQ=255)
                  tgGCGCAGAAACCTGGATGGTGCATATTTGCGGCATCGTTTTGCTTGCTGGTCGTCTGATGCATTATTACGGTTTTCATCACCGTCTGTTCCGCTGGCGACGTTCCGGCATGAGCGCCACCTGGTGTGCGCTGTTGCTGATGGTGCTgg                                                                                                                               <  3:350258/149‑1 (MQ=255)
                  tgGCGCAGAAACCTGGATGGTGCATATTTGCGGCATCGTTTTGCTTGCTGGTCGTCTGATGCATTATTACGGTTTTCATCACCGTCTGTTCCGCTGGCGACGTTCCGGCATGAGCGCCACCTGGTGTGCGCTGTTGCTGATGGTGCTgg                                                                                                                               >  8:165448/1‑149 (MQ=255)
                  tgGCGCAGAAACCTGGATGGTGCATATTTGCGGCATCGTTTTGCTTGCTGGTCGTCTGATGCATTATTACGGTTTTCATCACCGTCTGTTCCGCTGGCGACGTTCCGGCATGAGCGCCACCTGGTGTGCGCTGTTGCTGATGGTGCTgg                                                                                                                               <  3:156803/149‑1 (MQ=255)
                    gcgcAGAAACCTGGATGGTGCATATTTGCGGCATCGTTTTGCTTGCTGGTCGTCTGATGCATTATTACGGTTTTCATCACCGTCTGTTCCGCTGGCGACGTTCCGGCATGAGCGCCACCTGGTGTGCGCTGTTGCTGATGGTGCTGGc                                                                                                                              <  2:126787/148‑1 (MQ=255)
                             ccTGGATGGTGCATATTTGCGGCATCGTTTTGCTTGCTGGTCGTCTGATGCATTATTACGGTTTTCATCACCGTCTGTTCCGCTGGCGACGTTCCGGCATGAGCGCCACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTTgg                                                                                                                    <  2:8890/149‑1 (MQ=255)
                              cTGGATGGTGCATATTTGCGGCATCGTTTTGCTTGCTGGTCGTCTGATGCATTATTACGGTTTTCATCACCGTCTGTTCCGCTGGCGACGTTCCGGCATGAGCGCCACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTTgg                                                                                                                    >  4:180782/1‑148 (MQ=255)
                              cTGGATGGTGCATATTTGCGGCATCGTTTTGCTTGCTGGTCGTCTGATGCATTATTACGGTTTTCATCACCGTCTGTTCCGCTGGCGACGTTCCGGCATGAGCGCCACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTTGGt                                                                                                                   <  3:317780/149‑1 (MQ=255)
                              cTGGATGGTGCATATTTGCGGCATCGTTTTGCTTGCTGGTCGTCTGATGCATTATTACGGTTTTCATCACCGTCTGTTCCGCTGGCGACGTTCCGGCATGAGCGCCACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTTGGt                                                                                                                   <  3:180782/149‑1 (MQ=255)
                              cTGGATGGTGCATATTTGCGGCATCGTTTTGCTTGCTGGTCGTCTGATGCATTATTACGGTTTTCATCACCGTCTGTTCCGCTGGCGACGTTCCGGCATGAGCGCCACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTTGGt                                                                                                                   <  8:181563/149‑1 (MQ=255)
                               tGGATGGTGCATATTTGCGGCATCGTTTTGCTTGCTGGTCGTCTGATTCATTATTACGGTTTTCATAACACTCTGTTACGCTGGCGACGTTCCGGCATGAGCGCCACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTTGGta                                                                                                                  >  3:40824/1‑149 (MQ=255)
                                ggATGGTGCATATTTGCGGCATCGTTTTGCTTGCTGGTCGTCTGATGCATTATTACGGTTTTCATCACCGTCTGTTCCGCTGGCGACGTTCCGGCATGAGCGCCACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTTGGta                                                                                                                  >  2:279268/1‑148 (MQ=255)
                                 gATGGTGCATATTTGCGGCATCGTTTTGCTTGCTGGTCGTCTGATGCATTATTACGGTTTTCATCACCGTCTGTTCCGCTGGCGACGTTCCGGCATGAGCGCCACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTTGGta                                                                                                                  <  4:275514/147‑1 (MQ=255)
                                 gATGGTGCATATTTGCGGCAGCGTTTTGCTTGCTGGTCGTCTGATGCATTATTACGGTTTTCATCACCGTCTGTTCCGCTGGCGACGTTCCGGCATGAGCGCCACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTTGGtata                                                                                                                >  7:382534/1‑149 (MQ=255)
                                     gTGCATATTTGCGGCATCGTTTTGCTTGCTGGTCGTCTGAGGCATTATTACGGTTTTCATCACCGTCTGTTCCGCTGGCGACGTTCCGGCATGAGCGCCACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTTGGTATATGcc                                                                                                            >  1:6128/1‑149 (MQ=255)
                                      tGCATATTTGCGGCATCGTTTTGCTTGCTGGTCGTCTGATGCATTATTACGGTTTTCATCACCGTCTGTTCCGCTGGCGACGTTCCGGCATGAGCGCCACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTTGGTATATGccc                                                                                                           >  4:187528/1‑149 (MQ=255)
                                       gCATTTTTGCGGCATCGTTTTGCTTGCTGGTCGTCTGATGCATTATTACGGTTTTCATCACCGTCTGTTCCGCTGGCGACGTTCCGGCATGAGCGCCACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTTgg                                                                                                                    <  4:222592/139‑1 (MQ=255)
                                       gCATATTTGCGGCATCGTTTTGCTTGCTGGTCGTCTGATGCATTATTACGGTTTTCATCACCGTCTGTTCCGCTGGCGACGTTCCGGCATGAGCGCCACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTTgg                                                                                                                    >  3:222592/1‑139 (MQ=255)
                                        cATATTTGCGGCATCGTTTTGCTTGCTGGTCGTCTGATGCATTATTACGGTTTTCATCACCGTCTGTTCCGCTGGCGACGTTCCGGCATGAGCGCCACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTTGGTATATGCCCt                                                                                                          <  4:240935/148‑1 (MQ=255)
                                         atatTTGCGGCATCGTTTTGCTTGCTGGTCGTCTGATGCATTATTACGGTTTTCATCACCGTCTGTTCCGCTGGCGACGTTCCGGCATGAGCGCCACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTTGGTATATGc                                                                                                             >  3:154467/1‑144 (MQ=255)
                                         atatTTGCGGCATCGTTTTGCTTGCTGGTCGTCTGATGCATTATTACGGTTTTCATCACCGTCTGTTCCGCTGGCGACGTTCCGGCATGAGCGCCACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTTGGTATATGc                                                                                                             <  4:154467/144‑1 (MQ=255)
                                         atatTTGCGGCATCGTTTTGCTTGCTGGTCGTCTGATGCATTATTACGGTTTTCATCACCGTCTGTTCCGCTGGCGACGTTCCGGCATGAGCGCCACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTTGGTATATGCCCTgg                                                                                                        >  3:216295/1‑149 (MQ=255)
                                         atatTTGCGGCATCGTTTTGCTTGCTGGTCGTCTGATGCATTATTACGGTTTTCATCACCGTCTGTTCCGCTGGCGACGTTCCGGCATGAGCGCCACCTGGTGTGCGCTGTTGCTGATGGCGCTGGCGAATCTTTGGTATATGTCCTgg                                                                                                        <  7:346801/149‑1 (MQ=255)
                                            tttGCGGCATCGTTTTGCTTGCTGGTCGTCTGATGCATTATTACGGTTTTCATCACCGTCTGTTCCGCTGGCGACGTTCCGGCATGAGCGCCACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTTGGTATATGCCCTGGGa                                                                                                      <  4:284457/148‑1 (MQ=255)
                                                cGGCATCGTTTTGCTTGCTGGTCGTCTGATGCATTATTACGGTTTTCATCACCGTCTGTTCCGCTGGCGACGTTCCGGCATGAGCGCCACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTTGGTATATGCCCTGGGAGtgg                                                                                                  <  3:171997/148‑3 (MQ=255)
                                                cGGCATCGTTTTGCTTGCTGGTCGTCTGATGCATTATTACGGTTTTCATCACCGTCTGTTCCGCTGGCGACGTTCCGGCATGAGCGCCACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTTGGTATATGCCCTGGGAGTTg                                                                                                  <  7:49248/148‑1 (MQ=255)
                                                 ggCATCGTTTTGCTTGCTGGTCGTCTGATGCATTATTACGGTTTTCATCACCGTCTGTTCCGCTGGCGACGTTCCGGCATGAGCGCCACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTTGGTATATGCCCTGGGAGTTGGt                                                                                                >  2:65863/1‑149 (MQ=255)
                                                   cATCGTTTTGCTTGCTGGTCGTCTGATGCATTATTACGGTTTTCATCACCGTCTGTTCCGCTGGCGACGTTCCGGCATGAGCGCCACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTTGGTATATGCCCTGGGAGTTGGttt                                                                                              <  8:233414/149‑1 (MQ=255)
                                                       gTTTTGCTTGCTGGTCGTCTGATGCATTATTACGGTTTTCATCACCGTCTGTTCCGCTGGCGACGTTCCTGCATGAGCGCCACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTTGGTATATGCCCTGGGAGTTGGTTTtct                                                                                           <  7:285537/148‑1 (MQ=255)
                                                        ttttGCTTGCTGGTCGTCTGATGCATTATTACGGTTTTCATCACCGTCTGTTCCGCTGGCGACGTTCCGGCATGAGCGCCACCTGGTGTGCGCTGTTGCTGATGGTGCTTGCGAATCTTTGGTATATGCCCTGGGAGTTGGTTTTCTcc                                                                                         >  3:257946/1‑149 (MQ=255)
                                                        ttttGCTTGCTGGTCGTCTGATGCATTATTACGGTTTTCATCACCGTCTGTTCCGCTGGCGACGTTCCGGCATGAGCGCCACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTTGGTATATGCCCTGGGAGTTGGTTTTCTcc                                                                                         >  2:22933/1‑149 (MQ=255)
                                                        ttttGCTTGCTGGTCGTCTGATGCATTATTACGGTTTTCATCACCGTCTGTTCCGCTGGCGACGTTCCGGCATGAGCGCCACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTTGGTATATGCCCTGGGAGTTGGTTTTCTcc                                                                                         >  7:95358/1‑149 (MQ=255)
                                                           tgcttgctGGTCGTCTGATGCATTATTACGGTTTTCATCACCGTCTGTTCCGCTGGCGACGTTCCGGCATGAGCGCCACCTGGTGTGCGCTGTTGCTGATAGTGCTGGCGAATCTTTGGTATATGCCCTGGGAGTTGGTTTTCTCCCt                                                                                       <  1:207628/148‑1 (MQ=255)
                                                             cttgctGGTCGTCTGATGCATTATTACGGTTTTCATCACCGTCTGTTCCGCTGGCGACGTTCCGGCATGAGCGCCACCTGGTGTGCGCTGTTGCTGATGGTGCGGGCGAATCTTTGGTATATGCCCTGGGAGTTGGTTTTCTCCCTGc                                                                                     <  3:276824/148‑1 (MQ=255)
                                                                  tGGTCGTCTGATGCATTATTACGGTTTTCATCACCGTGTGTTCCGCTGGCGACGTTCCGGCATGAGCGCCACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTTGGGATATGCCCGGGGAGTTGGTTTTCTCCCTGCGTTAgc                                                                               >  4:302308/1‑149 (MQ=255)
                                                                    gtcgtcTGATGCATTATTACGGTTTTCATCACCGTCTGTTCCGCTGGCGACGTTCCGGCATGAGCGCCACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTTGGTATATGCCCTGGGAGTTGGTATTCTCCCTGCGTTAgcg                                                                              >  8:108132/1‑148 (MQ=255)
                                                                          tGATGCATTATTACGGTTTTCATCACCGTCTGTTCCGCTGGCGACGTTCCGGCATGAGCGCCACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTTGGTATATGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAGCGCACAATa                                                                       >  2:288667/1‑149 (MQ=255)
                                                                          tGATGCATTATTACGGTTTTCATCACCGTCTGTTCCGCTGGCGACGTTCCGGCATGAGCGCCACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTTGGTATATGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAGCGCACAATa                                                                       >  7:179392/1‑149 (MQ=255)
                                                                            aTGCATTATTACGGTTTTCATCACCGTCTGTTCCGCTGGCGACGTTCCGGCATGAGCGCCACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTTGGTATATGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAGCGCACAATAc                                                                      <  4:269904/148‑1 (MQ=255)
                                                                            aTGCATTATTACGGTTTTCATCACCGTCTGTTCCGCTGGCGACGTTCCGGCATGAGCGCCACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTTGGTATATGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAGCGCACAATACg                                                                     >  5:304068/1‑149 (MQ=255)
                                                                               cATTATTACGGTTTTCATCACGGTCTGGTCCGCTGGCGTCGTTCCTGCATGAGCGCCACCTGGTGTGCGCTGTTGCTGATGGTGCTGGAGAATCTTTGGTATATGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAGCGCGCAATACGCCa                                                                  <  4:40824/149‑1 (MQ=255)
                                                                               cATTATTACGGTTTTCATCACCGTCTGTTCCGCTGGCGACGTTCCGGCATGAGCGCCACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTTGGTATATGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAGCGCACAATACGCCa                                                                  <  7:165448/149‑1 (MQ=255)
                                                                                attattACGGTTTTCATAACCGTCTGTTCCGCTGGCGACGTTCCGGCATGAGCGCCACCTGGTGTGCGCTGTTGCTGATGGTGCGGGCGAATCTTGGGTATATGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAGCGCACAATACGCCAc                                                                 <  6:155170/149‑1 (MQ=255)
                                                                                    ttaCGGTTTTCATCACCGTCTGTTCCGCTGGCGACGTTCCGGCATGAGCGCCACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTTGGTATATGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAGCGCACAATACGCCActttc                                                             >  1:164920/1‑149 (MQ=255)
                                                                                       cGGTTTTCATCACCGTCTGTTCCGCTGGCGACGTTCCGGCATGAGCGCCACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTTGGTATATGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAGCGCACAATACGCCActttctt                                                           <  8:179392/148‑1 (MQ=255)
                                                                                        ggTGTTCATCACCGTCTGTTCCGCTGGCGACGTTCCGGCATGAGCGCCTCCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTTGGTATATGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAGCGCACAATACGCCACTTTCtttt                                                         <  3:223994/149‑1 (MQ=255)
                                                                                         gTTTTCATCACCGTCTGTTCCGCTGGCGACGTTCCGGCATGAGCGCCACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTTGGTATATGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAGCGCACAATACGCCACTTTCttttt                                                        >  5:61503/1‑149 (MQ=255)
                                                                                          ttttCATCACCGTCTGTTCCGCTGGCGACGTTCCGGCATGAGCGCCACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTTGGTATATGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAGCGCACAATACGCCACTTTCTTTTTc                                                       <  1:288667/149‑1 (MQ=255)
                                                                                              catcaCCGTCTGTTCCGCTGGCGACGTTCCGGCATGAGCGCCACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTTGGTATATGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAGCGCACAATACGCCACTTTCTTTTTCCCg                                                    >  6:58589/1‑148 (MQ=255)
                                                                                               atcaCCGTCTGTTCCGCTGGCGACGTTCCGGCATGAGCGCCACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTTGGTATATGCCCTGGGAGTTGGTTTtct                                                                                           >  1:329645/1‑108 (MQ=255)
                                                                                               atcaCCGTCTGTTCCGCTGGCGACGTTCCGGCATGAGCGCCACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTTGGTATATGCCCTGGGAGTTGGTTTtct                                                                                           <  2:329645/108‑1 (MQ=255)
                                                                                               atcaCCGTCTGTTCCGCTGGCGACGTTCCGGCATGAGCGCCACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTTGGTATATGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAGCGCACAATACGCCACTTTCTTTTTCCCGGa                                                  >  3:52129/1‑149 (MQ=255)
                                                                                               atcaCCGTCTGTTCCGCTGGCGACGTTCCGGCATGAGCGCCACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTTGGTATATGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAGCGCACAATACGCCACTTTCTTTTTCCCGGa                                                  >  7:201867/1‑149 (MQ=255)
                                                                                                    cGTCTGTTCCGCTGGCGACGTTCCGGCATGAGCGCCACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTTGGTATATGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAGCGCACAATACGCCACTTTCTTTTTCCCGGAtttt                                              >  2:6721/1‑148 (MQ=255)
                                                                                                          ttCCGCTGGCGACGTTCCGGCATGAGCGCCACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTTGGTATATGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAGCGCACAATACGCCACTTTCTTTTTCCCGGATTTTTACGtt                                        <  1:65863/148‑1 (MQ=255)
                                                                                                            ccGCTGGCGACGTTCCGGCATGAGCGCCACCTGGTGTGCGCTGTTTCTGATGGTGCTGGCGAATCTTTGGTATATGCCCTGGGTGTGGGTTTTCTCCCTGCGTTAGCGCACAATACGCCACTTTCTTTTTCCCGGATTTTTACGTTAt                                      <  1:279268/148‑1 (MQ=255)
                                                                                                             cGCTGGCGACGTTCCGGCATGAGCGCCACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTTGGTATATGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAGCGCACAATACGCCACTTTCTTTTTCCCGGATTTTTACGTTATGt                                    <  6:304068/149‑1 (MQ=255)
                                                                                                             cGCTGGCGACGTTCCGGCATGAGCGCCACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTTGGTATATGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAGCGCACAATACGCCACTTTCTTTTTCCCGGATTTTTACGTTATGt                                    >  1:386211/1‑149 (MQ=255)
                                                                                                                tGGCGACGTTCCGGCATGAGCGCCACCTGGTGTGTGCTGTTGCTGATGGTGCTGGCGAATCTTTGGTATATGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAGCGCACAATACGCCAc                                                                 <  4:132694/117‑1 (MQ=255)
                                                                                                                tGGCGACGTTCCGGCATGAGCGCCACCTGGTGTGTGCTGTTGCTGATGGTGCTGGCGAATCTTTGGTATATGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAGCGCACAATACGCCAc                                                                 <  4:127359/117‑1 (MQ=255)
                                                                                                                tGGCGACGTTCCGGCATGAGCGCCACCTGGTGTGTGCTGTTGCTGATGGTGCTGGCGAATCTTTGGTATATGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAGCGCACAATACGCCAc                                                                 >  3:127359/1‑117 (MQ=255)
                                                                                                                tGGCGACGTTCCGGCATGAGCGCCACCTGGTGTGTGCTGTTGCTGATGGTGCTGGCGAATCTTTGGTATATGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAGCGCACAATACGCCAc                                                                 >  3:132694/1‑117 (MQ=255)
                                                                                                                   cGACGTTCCGGCATGAGCGCCACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTTGGTATATGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAGCGCACAATACGCCACTTTCTTTTTCCCGGATTTTTACGTTATGTCTCAcc                              <  7:186978/149‑1 (MQ=255)
                                                                                                                       gTTCCGGCATGAGCGCCACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTTGGTATATGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAGCGCACAATACGCCACTTTCTTTTTCCCGGATTTTTACGTTATGTCTCACCGCGa                          >  1:300882/1‑149 (MQ=255)
                                                                                                                           cGGCATGAGCGCCACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTTGGTATATGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAGCGCACAATACGCCACTTTCTTTTTCCCGGATTTTTACGTTATGTCTCACCGCGacac                       <  3:187528/148‑1 (MQ=255)
                                                                                                                            ggCATGAGCGCCACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTTGGTATATGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAGCGCACAATACGCCACTTTCTTTTTCCCGGATTTTTACGTTATGTCTCACCGCGACACGc                     >  1:61678/1‑149 (MQ=255)
                                                                                                                             gCATGAGCGCCACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTTGGTATATGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAGCGCACAATACGCCACTTTCTTTTTCCCGGATTTTTACGTTATGTCTCACCGCGACACGc                     >  1:240063/1‑148 (MQ=255)
                                                                                                                              cATGAGCGCCACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTTGGTATATGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAGCGCACAATACGCCACTTTCTTTTTCCCGGATTTTTACGTTATGTCTCACCGCGACACGCTa                   <  4:164504/149‑1 (MQ=255)
                                                                                                                              cATGAGCGCCACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTTGGTATATGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAGCGCACAATACGCCACTTTCTTTTTCCCGGATTTTTACGTTATGTCTCACCGCGACACGCTa                   <  1:22933/149‑1 (MQ=255)
                                                                                                                                  aGCGCCACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTTGGTATATGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAGCGCACAATACGCCACTTTCTTTTTCCCGGATTTTTACGTTATGTCTCACCGCGACACGCTAtttt               <  7:137674/149‑1 (MQ=255)
                                                                                                                                    cgcCACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTTGGTATATGCCCTGGTAGTTGGTTTTCTCCCTGCGTAAGCGCACAATACGCCACTTTCTTTTTCCCGGATTTTTACGTTATGTCTCACCGCGACACGCTATTTTCt             <  8:382534/149‑1 (MQ=255)
                                                                                                                                    cgcCACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTTGGTATATGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAGCGCACAATACGCCACTTTCTTTTTCCCGGATTTTTACGTTATGTCTCACCGCGACACGCTATTTTCt             >  3:306896/1‑149 (MQ=255)
                                                                                                                                      ccACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTTGGTATATGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAGCGCACAATACGCCACTTTCTTTTTCCCGGATTTTTACGTTATGTCTCACCGCGACACGCTATTTTCTGc           <  5:58589/149‑1 (MQ=255)
                                                                                                                                       cACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTTGGTATATGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAGCGCACAATACGCCACTTTCTTTTTCCCGGATTTTTACGTTATGTCTCACCGCGACACGCTATTTTCTGcc          >  3:2729/1‑149 (MQ=255)
                                                                                                                                        aCCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTTGGTATATGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAGCGCACAATACGCCACTTTCTTTTTCCCGGATTTTTACGTTATGTCTCACCGCGACACGCTATTTTCTGccc         <  2:164920/149‑1 (MQ=255)
                                                                                                                                        aCCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTTGGTATATGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAGCGCACAATACGCCACTTTCTTTTTCCCGGATTTTTACGTTATGTCTCACCGCGACACGCTATTTTCTGccc         >  1:37032/1‑149 (MQ=255)
                                                                                                                                        aCCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTTGGTATATGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAGCGCACAATACGCCACTTTCTTTTTCCCGGATTTTTACGTTATGTCTCACCGCGACACGCTATTTTCTGccc         >  3:161746/1‑149 (MQ=255)
                                                                                                                                        aCCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTTGGTATATGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAGCGCACAATACGCCACTTTCTTTTTCCCGGATTTTTACGTTATGTCTCACCGCGACACGCTATTTTCTGccc         <  4:2729/149‑1 (MQ=255)
                                                                                                                                               gtgCGCTGTTGCTGATGGTGCTGGCGAATCTTTGGTATATGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAGCGCACAATACGCCACTTTCTTTTTCCCGGATTTTTACGTTATGTCTCACCGCGACACGCTATTTTCTGCCCCTATCGc  <  7:354185/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
GCTGTTTATGGAAATGAATGGCGCAGAAACCTGGATGGTGCATATTTGCGGCATCGTTTTGCTTGCTGGTCGTCTGATGCATTATTACGGTTTTCATCACCGTCTGTTCCGCTGGCGACGTTCCGGCATGAGCGCCACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTTGGTATATGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAGCGCACAATACGCCACTTTCTTTTTCCCGGATTTTTACGTTATGTCTCACCGCGACACGCTATTTTCTGCCCCTATCGC  >  NC_000913/1952448‑1952739

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