Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,538,900 G→A 5.2% T201M (ACG→ATG)  nanC ← N‑acetylnuraminic acid outer membrane channel protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,538,9000GA 5.2% 272.9 / 3.8 77T201M (ACG→ATG) nanCN‑acetylnuraminic acid outer membrane channel protein
Reads supporting (aligned to +/- strand):  ref base G (34/39);  new base A (2/2);  total (36/41)
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.77e-01

TGTGACTGTCTCCTGTCTACTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGTTGTAAACACCCTGACGGTCAAGGTAGTCATATTCTATGTATGGCGTAATATCGGGCGTCATATGGTATTGTAGAACAAATGCATTTTCCGTCGCCCATTTCTTATGGTTTGCATAGCGATAATCGTTCTGTTTGCTGTATAGCGTCGTTTGCCATGCGAAGGTGAAATCACTATTAATATGGTAAGTGACATATC  >  NC_000913/4538766‑4539038
                                                                                                                                      |                                                                                                                                          
tgtgACTGTCTCCTGTCTACTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGTTGTAAACACCCTGACGGTCAAGGTAGTCATATTCTATGTATGGCGTAATATCGGGCGTCata                                                                                                                                       <  2:1280826/140‑1 (MQ=255)
tgtgACTGTCTCCTGTCTACTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGTTGTAAACACCCTGACGGTCAAGGTAGTCATATTCTATGTATGGCGTAATATCGGGCGTCata                                                                                                                                       <  2:1337053/140‑1 (MQ=255)
ggtgACTGTCACCTGTCTACTACAGTTTAAATGACCCACCAATGCGCTAGCTCTGTTGCGAGAAATGATCACTGCCGATGTCATCACCCTGCTGGAGACGGGCGTCAGGTGCTATGTAGTGCGTGATATAGgtcgtccac                                                                                                                                       >  2:1896346/2‑137 (MQ=255)
          tcCTGTCTACTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGTTGTAAACACCCTGACGGTCAAGGTAGTCATATTCTATGTATGGCGTAATATCGGGCGTCATATGGTATTGTa                                                                                                                             >  1:1338444/1‑140 (MQ=255)
                         tttAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGTTGTAAACACCCTGACGGTCAAGGTAGTCATATTCTATGTATGGCGTAATATCGGGCGTCATATGGTATTGTAGAACAAATGCAtttt                                                                                                              >  1:1051595/1‑140 (MQ=255)
                         tttAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGTTGTAAACACCCTGACGGTCAAGGTAGTCATATTCTATGTAAGGCGTAATATCGGGCGTCATATGGTATTGTAGAAAAAATGCAtttt                                                                                                              >  2:1602197/1‑140 (MQ=255)
                         tttAAATGACACACCAATGCGATAACTGGTTTCCGATAAATTATCTCTGCCGTTGTAAACACCCTGACGGTCAAGGTAGTCATCTACAATGTATGGCGACATATCGGGCGTCATATGGTATAGTAGAACAAATGCAattt                                                                                                              >  2:387021/1‑140 (MQ=255)
                           tAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGTTGTAAACACCCTGACGGTCAAGGTAGTCATATTCTATGTATGGCGTAATATCGGGCGTCATATGGTATTGTAGAACAAATGCATTTTcc                                                                                                            >  1:1620137/1‑140 (MQ=255)
                                  cacacCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGTTGTAAACACCCTGACGGTCAAGGTAGTCATATTCTATGTATGGCGTAATATCGGGCGTCATATGGTATTGTAGAACAAATGCATTTTCCGTCGccc                                                                                                     <  1:1005905/140‑1 (MQ=255)
                                   acacCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGTTGTAAACACCCTGTCGGTCAAGGAAGTCTTATTCTATGTTTGGCGTAATATCGGGCGTCATAAGGTATTGAAGAACAAATGCATTATCCTTCGCCCa                                                                                                    >  1:1868010/1‑140 (MQ=255)
                                    cacCAATGCGATAACTGTTTTCCGATAAATTCTCTCTGCCGTTGTAAACACCCTGACGGTCAAGGTAGTCATATTCTATGTATGGCGTAATATCGGGCGTCATATGGTATTGTAGAACAAAGGCATTGTCCGTCGCCCAt                                                                                                   <  2:1634869/140‑1 (MQ=255)
                                         tagCGATATCTGTTTTCCGATAAATTATCTCTGCCGTTGTAAACACCCTGTCGGTCTAGTTAGTCTTATTCTATGTAGGGCGTAATATCGGGCGTCATATGGTATTGTAGAACTAATGCATTTTCCGTCGCCCATTTCtt                                                                                              <  2:911470/138‑1 (MQ=255)
                                         ctGCGATACCTTTTTTCCGATAAATTATCACTGCCGTTTTAAACACCCTGACGGTCAAGGTAGTCATATTCTATGTATGGCGTAATATCGGGCGTCATATGGTATTGTAGAACAAATGCATTTTCCGTCGCCCATTTCtt                                                                                              <  2:923233/139‑1 (MQ=255)
                                         aTGCGAGAACTGTTTTCCGATATAGTATCTCTGCCGTTGTAAAAACCCTGCCGGTCAAGGTAGTCATATTCTATGTATGGCGTAATATCGGGCGTCATATGGTATTGTAGAACAAATGCTTTGTCCGTCGCCCATTTCtt                                                                                              <  2:303330/140‑1 (MQ=255)
                                                 aCTGTTTTCCGATAAATTATCTCTGCCGTTGTAAACACCCTGACGGTCAAGGTAGTCATATTCTATGTATGGCGTAATATCGGGCGTCATATGGTATTGTAGAACAAATGCATTTTCCGTCGCCCATTTCTTATGGTTTg                                                                                      >  2:1385371/1‑140 (MQ=255)
                                                         ccGATAAATTATCTCTGCCGTTGTAAACACCCTGACGGTCAAGGTAGTCATATTCTATGTATGGCGTAATATCGGGCGTCATATGGTATTGTAGAACAAATGCATTTTCCGTCGCCCATTTCTTATGGTTTGCATAGCGa                                                                              <  1:759811/140‑1 (MQ=255)
                                                               aaTTATCTCTGCCGTTGTAAACACCCTGACGGTCAAGGTAGTCATATTCTATGTATGGCGTAATATCGGGCGTCATATGGTATTGTAGAACAAATGCATTTTCCGTCGCCCATTTCTTATGGTTTGCATAGCGATAATCg                                                                        >  2:690244/1‑140 (MQ=255)
                                                               aaTTATCTCTGCCGTTGTAAACACCCTGACGGTCAAGGTAGACATATTCTATGTATGGCGTAATATCGGGCGTCATATGGTATTGTAGAACAAATGCATTTTCCGTCGCCCATTTCTTATGGTTTGCATAGCGATAATCg                                                                        >  2:1089739/1‑140 (MQ=255)
                                                                 ttATCTCTGCCGTTGTAAACACCCTGACGGTCAAGGTAGTCATATTCTATGTATGGCGTAATATCGGGCGTCATATGGTATTGTAGAACAAATGCATTTTCCGTCGCCCATTTCTTATGGTTTGCATAGCGATAATCGtt                                                                      <  1:489209/140‑1 (MQ=255)
                                                                 ttATCTCTGCCGTTGTAAACACCCTGACGGTCAAGGTAGTCATATTCTATGTATGGCGTAATATCGGGCGTCATATGGTATTGTAGAACAAATGCATTTTCCGTCGCCCATTTCTTATGGTTTGCATAGCGATAATCGtt                                                                      <  1:387451/140‑1 (MQ=255)
                                                                      tctGCCGTTGTAAACACCCTGACGGTCAAGGTAGTCATATTCTATGTATGGCGTAATATCGGGCGTCATATGGTATTGTAGAACAAATGCATTGTCCGTCGCCCATTTCTTATGGTTTGCATAGCGATAATCGTTCTGtt                                                                 <  2:1096642/140‑1 (MQ=255)
                                                                       ctGCCGTTGTAATCACCCTGACGTTCAAGGTAGTCATATTCTATGTATGGCGTAATATCGGGCGTCATATGGTTTTGTAGAACAAATGCATTTTCCGTCGCCCATTTCTTATGGTTTGCATAGCGATAATCGTTCTGttt                                                                <  2:1426686/140‑1 (MQ=255)
                                                                       ctGCCGTTGTAAACACCCTGACGGTCAAGGTAGTCATATTCTATGTATGGCGTAATATCGGGCGTCATATGGTATTGTAGAACAAATGCATTTTCCGTCGCCCATTTCTTATGGTTTGCATAGCGATAATCGTTCTGttt                                                                >  2:1509532/1‑140 (MQ=255)
                                                                       ctGCCGTTGTAAACACCCTGACGGTCAAGGTAGTCATATTCTATGTATGGCGTAATATCGGGCGTCATATGGTATTGTAGAACAAATGCATTGTCCGTCGCCCATTTCTTATGGTTTGCATAGCGATAATCGTTCTGttt                                                                <  2:901227/140‑1 (MQ=255)
                                                                        tGCCGTTGTATACACCCAGACGGTCAAGGTAGTCATATTCTATGTATGGCGTAATATCGGGCGTCATATGGTATTGTAGAACAAATGCATTTTCCGTCGCCCATTTCTTATGGTTTGCTTAGCGATAATCGTTCTGTTTg                                                               <  1:895484/140‑1 (MQ=255)
                                                                        tGCCGTTGTAAACACCCTGACGGTCAAGGTAGTCATATTCTATGTATGGCGTAATATCGGGCGTCATATGGTATTGTAGAACAAATGCATTTTCCGTCGCCCATTTCTTATGGTTTGCATAGCGATAATCGTTCTGTTTg                                                               <  2:209703/140‑1 (MQ=255)
                                                                        tGCCGTTGTAAACACCCTGACGGTCAAGGTAGTCATATTCTATGTATGGCGTAATATCGGGCGTCATATGGTATTGTAGAACAAATGCATTTTCCGTCGCCCATTTCTTATGGTTTGCATAGCGATAATCGTTCTGTTTg                                                               <  1:895762/140‑1 (MQ=255)
                                                                        tGCCGTTGTAAACACCCTGACGGTCAAGGTAGTCATATTCTATGTATGGCGTAATATCGGGCGTCATATGGTATTGTAGAACAAATGCATTTTCCGTCGCCCATTTCTTATGGTTTGCATAGCGATAATCGTTCTGTTTg                                                               <  1:1820846/140‑1 (MQ=255)
                                                                         gCCGTTGTAAACACCCTGACGGTCAAGGTAGTCATATTCTATGTATGGCGTAATATCGGGCGTCATATGGTATTGTAGAACAAATGCATTTTCCGTCGCCCATTTCTTATGGTTTGCATAGCGATAATCGTTCTGt                                                                  <  1:678847/136‑1 (MQ=255)
                                                                         gCCGTTGTAAACACCCTGACGGTCAAGGTAGTCATATTCTATGTATGGCGTAATATCGGGCGTCATATGGTATTGTAGAACAAATGCATTTTCCGTCGCCCATTTCTTATGGTTTGCATAGCGATAATCGTTCTGt                                                                  >  2:678847/1‑136 (MQ=255)
                                                                               gTAAACACCCTGACGGTCAAGGTAGTCATATTCTATGTATGGCGTAATATCGGGCGTCATATGGTATTGTAGAACAAATGCATTTTCCGTCGCCCATTTCTTATGGTTTGCATAGCGATAATCGTTCTGTTTGCTGtata                                                        >  2:172500/1‑140 (MQ=255)
                                                                                tAAACACCCTGACGGTCAAGGTAGTCATATTCTATGTATGGCGTAATATCGGGCGTCATATGGTATTGTAGAACAAATGCATTTTCCGTCGCCCATTTCTTATGGTTTGCATAGCGATAATCGTTCTGTTTGCTGTATAg                                                       >  1:780953/1‑140 (MQ=255)
                                                                                      cccTGACGGTCAAGGTAGTCATATTCTATGTATGGCGTAATATCGGGCGTCATATGGTGTTGTAGAACAAATGCATTTTCCGTCGCCCATTTCTTATGGTTTGCATAGCGATAATCGTTCTGTTTGCTGTATAGcgtcgt                                                 <  1:5467/140‑1 (MQ=255)
                                                                                       ccTGACGGTCAAGGTAGTCATATTCTATGTATGGCGTAATATCGGGCGTCATATGGTATTGTAGAACAAATGCATTTTCCGTCGCCCATTTCTTATGGTTTGCATAGCGATAATCGTTCTGTTTGCTGTATAGCGTCGtt                                                >  2:989467/1‑140 (MQ=255)
                                                                                          gACGGTCAAGGTAGTCATATTCTATGTATGGCGTAATATCGGGCGTCATATGGTATTGTAGAACAAATGCATTTTCCGTCGCCCATTTCTTATGGTTTGCATAGCGATAATCGTTCTGTTTGCTGTATAGCGTCGTTTGc                                             >  1:747687/1‑140 (MQ=255)
                                                                                              gTCAAGGTAGTCATATTCTATGTATGGCGTAATATCGGGCGTCATATGGTATTGTAGAACAAATGCATTTTCCGTCGCCCATTTCTTATGGTTTGCATAGCGATAATCGTTCTGTTTGCTGTATAGCGTCGTTTGCCATg                                         <  2:1549419/140‑1 (MQ=255)
                                                                                               tCAAGGTAGTCATATTCTATGTATGGCGTAATATCGGGCGTCATATGGTATTGTAGAACAAATGCATTTTCCGTCGCCCATTTCTTATGGTTTGCATAGCGATAATCGTTCTGTTTGCTGTATAGCGTCGTTTGCCATGc                                        >  2:76912/1‑140 (MQ=255)
                                                                                               tCAAGGTAGTCATATTCTATGTATGGCGTAATATCGGGCGTCATATGGTATTGTAGAACAAATGCATTTTCCGTCGCCCATTTCTTATGGTTTGCATAGCGATAATCGTTCTGTTTGCTGTATAGCGTCGTTTGCCATGc                                        >  2:101575/1‑140 (MQ=255)
                                                                                               tCAAGGTAGTCATATTCTATGTATGGCGTAATATCGGGCGTCATATGGTATTGTAGAACAAATGCATTTTCCGTCGCCCATTTCTTATGGTTTGCATAGCGATAATCGTTCTGTTTGCTGTATAGCGTCGTTTGCCATGc                                        >  2:460435/1‑140 (MQ=255)
                                                                                               tCAAGGTAGTCATATTCTATGTATGGCGTAATATCGGGCGTCATATGGTATTGTAGAACAAATGCATTTTCCGTCGCCCATTTCTTATGGTTTGCATAGCGATAATCGTTCTGTTTGCTGTATAGCGTCGTTTGCCATGc                                        >  2:76883/1‑140 (MQ=255)
                                                                                               tCAAGGTAGTCATATTCTATGTATGGCGTAATATCGGGCGTCATATGGTATTGTAGAACAAATGCATTTTCCGTCGCCCATTTCTTATGGTTTGCATAGCGATAATCGTTCTGTTTGCTGTATAGCGTCGTTTGCCATGc                                        >  2:855670/1‑140 (MQ=255)
                                                                                               tCAAGGTAGTCATATTCTATGTATGGCGTAATATCGGGCGTCATATGGTATTGTAGAACAAATGCATTTTCCGTCGCCCATTTCTTATGGTTTGCATAGCGATAATCGTTCTGTTTGCTGTATAGCGTCGTTAGCCATGc                                        >  2:855639/1‑140 (MQ=255)
                                                                                                    gTAGTCATATTCTATGTATGGCGTAATATCGGGCGTCATATGGTATTGTAGAACAAATGCATTTTCCGTCGCCCATTTCTTATGGTTTGCATAGCGATAATCGTTCTGTTTGCTGTATAGCGTCGTTTGCCATGCGAAgg                                   <  2:1675249/140‑1 (MQ=255)
                                                                                                          atatTCTATGTATGGCGTAATATCGGGCGTCATATGGTATTGTAGAACAAATGCATTTTCCGTCGCCCATTTCTTATGGTTTGCATAGCGATAATCGTTCTGTTTGCTGTATTGCGTCGTTTGCCATGCGAAGGTGAAAt                             <  2:1185629/140‑1 (MQ=255)
                                                                                                          atatTCTATGTATGGCGTAATATCGGGCGTCATATGGTATTGTAGAACAAATGCATTTTCCGTCGCCCATTTCTTATGGTTTGCATAGCGATAATCGTTCTGTTTGCTGTATAGCGTCGTTTGCCATGCGAAGGTGAAAt                             <  1:873495/140‑1 (MQ=255)
                                                                                                          atatTCTATGTATGGCGTAATATCGGGCGTCATATGGTATTGTAGAACAAATGCATTTTCCGTCGCCCATTTCTTATGGTTTGCATAGCGATAATCGTTCTGTTTGCTGTATAGCGTCGTTTGCCATGCGAAGGTGAAAt                             <  1:984372/140‑1 (MQ=255)
                                                                                                          atatTCTATGTATGGCGTAATATCGGGCGTCATATGGTATTGTAGAACAAATGCATTTTCCGTCGCCCATTTCTTATGGTTTGCATAGCGATAATCGTTCTGTTTGCTGTATAGCGTCGTTTGCCATGCGAAGGTGAAAt                             <  2:486381/140‑1 (MQ=255)
                                                                                                             ttCTATGTATGGCGTAATATCGGGCGTCATATGGTATTGTAGAACAAATGCATTTTCCGTCGCCCATTTCTTATGGTTTGCATAGCGATAATCGTTCTGTTTGCTGTATAGCGTCGTTTGCCATGCGAAGGTGAAATCAc                          <  1:1684747/140‑1 (MQ=255)
                                                                                                              tCTATGTATGGCGTAATATCGGGCGTCATATGGTATTGTAGAACAAATGCATTTTCCGTCGCCCATTTCTTATGGTTTGCATAGCGATAATCGTTCTGTTTGCTGTATAGCGTCGTTTGCCATGCGAAGGTGAAATCACt                         >  2:1420278/1‑140 (MQ=255)
                                                                                                               cTATGTATGGCGTAATATCGGGCGTCATATGGTATTGTAGAACAAATGCATTTTCCGTCGCCCATTTCTTATGGTTTGCATAGCGATAATCGTTCTGTTTGCTGTATAGCGTCGTTTGCCATGCGAAGGTGAAATCACTa                        >  2:719485/1‑140 (MQ=255)
                                                                                                               cTATGTATGGCGTAATATCGGGCGTCATATGGTATTGTAGAACAAATGCATTTTCCGTCGCCCATTTCTTATGGTTTGCATAGCGATAATCGTTCTGTTTGCTGTATAGCGTCGTTTGCCATGCGAAGGTGAAATCACTa                        >  2:1690508/1‑140 (MQ=255)
                                                                                                                  tgtatgGCGTAATATCGGGCGTCATATGGTATTGTAGAACAAATGCATTTTCCGTCGCCCATTTCTTATGGTTTGCATAGCGATAATCGTTCTGTTTGCTGTATAGCGTCGTTTGCCATGCGAAGGTGAAATCACTATTa                     >  1:1602716/1‑140 (MQ=255)
                                                                                                                   gtatgGCGTAATATCGGGCGTCATATGGTATTGTAGAACAAATGCATTTTCCGTCGCCCATTTCTTATGGTTTGCATAGCGATAATCGTTCTGTTTGCTGTATAGCGTCGTTTGCCATGCGAAGGTGAAATCACTATTaa                    <  2:780953/140‑1 (MQ=255)
                                                                                                                     atgGCGTAATATCGGGCGTCATATGGTATTGTAGAACAAATGCATTTTCCGTCGCCCATTTCTTATGGTTTGCATAGCGATAATCGTTCTGTTTGCTGTATAGCGTCGTTTGCCATGCGAAGGTGAAATCACTATTAata                  <  1:297249/140‑1 (MQ=255)
                                                                                                                     atgGCGTAATATCGGGCGTCATATGGTATTGTAGAACAAATGCATTTTCCGTCGCCCATTTCTTATGGTTTGCATAGCGATAATCGTTCTGTTTGCTGTATAGCGTCGTTTGCCATGCGAAGGTGAAATCACTATTAata                  <  1:297141/140‑1 (MQ=255)
                                                                                                                      tgGCGTAATATCGGGCATCATATGGTATTGTAGAACAAATGCATTTTCCGTCGCCCATTTCTTATGGTTTGCAt                                                                                   <  1:1038475/74‑1 (MQ=255)
                                                                                                                      tgGCGTAATATCGGGCATCATATGGTATTGTAGAACAAATGCATTTTCCGTCGCCCATTTCTTATGGTTTGCAt                                                                                   >  2:1038588/1‑74 (MQ=255)
                                                                                                                      tgGCGTAATATCGGGCATCATATGGTATTGTAGAACAAATGCATTTTCCGTCGCCCATTTCTTATGGTTTGCAt                                                                                   <  1:1038588/74‑1 (MQ=255)
                                                                                                                      tgGCGTAATATCGGGCATCATATGGTATTGTAGAACAAATGCATTTTCCGTCGCCCATTTCTTATGGTTTGCAt                                                                                   >  2:1038475/1‑74 (MQ=255)
                                                                                                                       ggCGTAATATCGGGCGTCATATGGTATTGTAGAACAAATGCATTTTCCGTCGCCCATTTCTTATGGTTTGCATAGCGATAATCGTTCTGTTTGCTGTATAGCGTCGTTTGCCATGCGAAGGTGAAATCACTATTAATATg                >  1:441042/1‑140 (MQ=255)
                                                                                                                          gTAATATCGGGCGTCATATGGTATTGTAGAACAAATGCATTTTCCGTCGCCCATTTCTTATGGTTTGCATAGCGATAATCGTTCTGTTTGCTGTATAGCGTCGTTTGCCATGCGAAGGTGAAATCACTATTAATATGGTa             <  2:1635374/140‑1 (MQ=255)
                                                                                                                           tAATATCGGGCGTCATATGGTATTGTAGAACAAATGCATTTTCCGTCGCCCATTTCTTATGGTTTGCATAGCGATAATCGTTCTGTTTGCTGTATAGCGTCGTTTGCCATGCGAAGGTGAAATCACTATTAATATGGTaa            <  1:386180/140‑1 (MQ=255)
                                                                                                                           tAATATCGGGCGTCATATGGTATTGTAGAACAAATGCATTTTCCGTCGCCCATTTCTTATGGTTTGCATAGCGATAATCGTTCTGTTTGCTGTATAGCGTCGTTTGCCATGCGAAGGTGAAATCACTATTAATATGGTaa            <  1:1385371/140‑1 (MQ=255)
                                                                                                                              tatCGGGCGTCATATGGTATTGTAGAACAAATGCATTTTCCGTCGCCCATTTCTTATGGTTTGCATAGCGATAATCGTTCTGTTTGCTGTATAGCGTCGTTTGCCATGCGAAGGTGAAATCACTATTAATATGGTAAGTg         >  2:1443232/1‑140 (MQ=255)
                                                                                                                              tatCGGGCGTCATATGGTATTGTAGAACAAATGCATTTTCCGTCGCCCATTTCTTATGGTTTGCATAGCGATAATCGTTCTGTTTGCTGTATAGCGTCGTTTGCCATGCGAAGGTGAAATCACTATTAATATGGTAAGTg         >  2:702616/1‑140 (MQ=255)
                                                                                                                              tatCGGGCGTCATATGGTATTGTAGAACAAATGCATTTTCCGTCGCCCATTTCTTATGGTTTGCATAGCGATAATCGTTCTGTTTGCTGTATAGCGTCGTTTGCCATGCGAAGGTGAAATCACTATTAATATGGTAAGTg         >  2:1140486/1‑140 (MQ=255)
                                                                                                                              tatCGGGCGTCATATGGTATTGTAGAACAAATGCATTTTCCGTCGCCCATTTCTTATGGTTTGCATAGCGATAATCGTTCTGTTTGCTGTATAGCGTCGTTTGCCATGCGAAGGTGAAAACACTATTAATATTGTAAGTg         >  2:1804031/1‑140 (MQ=255)
                                                                                                                              tatCGGGCGTCATATGGTATTGTAGAACAAATGCATTTTCCGTCGCACATTTCTTATGGTTTGCATAGCGATAATCGTTCTGTTTGCTGTATAGCGTCGTTTGCCATGCGAAGGTGAAATCACTATTAATATGGTAAGTg         >  2:1803765/1‑140 (MQ=255)
                                                                                                                                tCGGGCGTCATATGGTATTGTAGAACAAATGCATTTTCCGTCGCCCATTTCTTATGGTTTGCATAGCGTTAATCGTTCTGTTTGCTGTTTAGCGTCGTTTGCCATGCGAAGGTGAATTCACTATTAATATGGTAAGTGAc       <  2:1110162/140‑1 (MQ=255)
                                                                                                                                tCGGGCGTCATATGGTATTGTAGAACAAATGCATTTTCCGTCGCCCATTTCTTATGGTTTGCATAGCGATAATCGTTCTGTTTGCTGTATAGCGTCGTTTGCCATGCGAAGGTGAAATCACTATTAATATGGTAAGTGAc       <  2:295878/140‑1 (MQ=255)
                                                                                                                                tCGGGCGTCATATGGTATTGTAGAACAAATGCATTTTCCGTCGCCCATTTCTTATGGTTTGCATAGCGATAATCGTTCTGTTTGCTGTATAGCGTCGTTTGCCATGCGAAGGTGAAATCACTATTAATATGGTAAGTGAc       <  1:433761/140‑1 (MQ=255)
                                                                                                                                tCGGGCGTCATATGGTATTGTAGAACAAATGCATTTTCCGTCGCCCATTTCTTATGGTTTGCATAGCGATAATCGTTCTGTTTGCTGTATAGCGTCGTTTGCCATGCGAAGGTGAAATCACTATTAATATGGTAAGTGAc       <  1:436587/140‑1 (MQ=255)
                                                                                                                                tCGGGCGTCATATGGTATTGTAGAACAAATGCATTTTCCGTCGCCCATTTCTTATGGTTTGCATAGCGATAATCGTTCTGTTTGCTGTATAGCGTCGTTTGCCATGCGAAGGTGAAATCACTATTAATATGGTAAGTGAc       <  2:1338119/140‑1 (MQ=255)
                                                                                                                                   ggCGTCATATGGTATTGTAGAACAAATGCATTTTCCGTCGCCCATTTCTTATGGTTTGCATAGCGATAATCGTTCTGTTTGCTGTATAGCGTCGTTTGCCATGCGAAGGTGAAATCACTATTAATATGGTAAGTGACata    >  2:1790314/1‑140 (MQ=255)
                                                                                                                                   ggCGTCATATGGTATTGTAGAACAAATGCATTTTCCGTCGCCCATTTCTTATGGTTTGCATAGCGATAATCGTTCTGTTTGCTGTATAGCGTCGTTTGCCATGCGAAGGTGAAATCACTATTAATATGGTAAGTGACata    >  1:944655/1‑140 (MQ=255)
                                                                                                                                   ggCGTCATATGGTATTGTAGAACAAATGCATTTTCCGTCGCCCATTTCTTATGGTTTGCATAGCGATAATCGTTCTGTTTGCTGTATAGCGTCGTTTGCCATGCGAAGGTGAAATCACTATTAATATGGTAAGTGACata    <  1:1373685/140‑1 (MQ=255)
                                                                                                                                   ggCGTCATATGGTATTGTAGAACAAATGCATTTTCCGTCGCCCATTTCTTATGGTTTGCATAGCGATAATCGTTCTGTTTGCTGTATAGCGTCGTTTGCCATGCGAAGGTGAAATCACTATTAATATGGTAAGTGACaaa    >  2:1607040/1‑138 (MQ=255)
                                                                                                                                     cGTCATATGGTATTGTAGAACAAATGCATTTTCCGTCGCCCATTTCTTATGGTTTGCATAGCGATAATCGTTCTGTTTGCTGTATAGCGTCGTTTGCCATGCGAAGGTGAAATCACTATTAATATGGTAAGTGACATatc  <  2:396310/140‑1 (MQ=255)
                                                                                                                                      |                                                                                                                                          
TGTGACTGTCTCCTGTCTACTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGTTGTAAACACCCTGACGGTCAAGGTAGTCATATTCTATGTATGGCGTAATATCGGGCGTCATATGGTATTGTAGAACAAATGCATTTTCCGTCGCCCATTTCTTATGGTTTGCATAGCGATAATCGTTCTGTTTGCTGTATAGCGTCGTTTGCCATGCGAAGGTGAAATCACTATTAATATGGTAAGTGACATATC  >  NC_000913/4538766‑4539038

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: