Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,654,745 C→T 6.5% D133N (GAT→AAT)  rspA ← bifunctional D‑altronate/D‑mannonate dehydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,654,7450CT 6.5% 194.8 / 2.6 61D133N (GAT→AAT) rspAbifunctional D‑altronate/D‑mannonate dehydratase
Reads supporting (aligned to +/- strand):  ref base C (29/28);  new base T (2/2);  total (31/30)
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.00e-01

GCTCTTCCGGCCACTGTCCTTTGGTTGCGGGTTCATAAGCCAGACCTTTACCTTTCGACATGCCGTAGGTGGTTTTCATACCAGGGATTCCGCACTGCACGCGGATGGCTTTGAATCCAAGCTCTTGATGACGGGCATAATCATCCAGAGCTTCATCAATACTGTGACCGGTGGTATGGCAATAAACCATCACCCCTTCACGAGACGCGCCGCCGAGTAACTGGTAAAGCGGCATGTTGGCAGCTTTGGCTTTAATATCCCACAGCGCCATATCAACCGCTGAAATGGCCGA  >  NC_000913/1654601‑1654892
                                                                                                                                                |                                                                                                                                                   
gctCTTCCGGCCACTGTCCTTTGGTTGCGGGTTCATAAGCCAGACCTTTACCTTTCGACATGCCGTAGGTGGTTTTCATACCAGGGATTCCGCACTGCACGCGGATGGCTTTGAATCCAAGCTCTTGATGACGGGCATAATCATCCana                                                                                                                                                 <  3:169514/149‑1 (MQ=255)
gctCTTCCGGCCACTGTCCTTTGGTTGCGGGTTCATAAGCCAGACCTTTACCTTTCGACATGCCGTAGGTGGTTTTCATACCAGGGATTCCGCACTGCACGCGGATGGCTTTGAATCCAAGCTCTTGATGACGGGCATAATCATCCaga                                                                                                                                                 <  3:205801/149‑1 (MQ=255)
  tctTCCGGCCACTGTCCTTTGGTTGCGGGTTCATAAGCCAGACCTTTACCTTTCGACATGCCGTAGGTGGTTTTCATACCAGGGATTCCGCACTGCACGCGGATGGCTTTGAATCCAAGCTCTTGATGACGGGCATAATCATCCAGAGc                                                                                                                                               >  2:160100/1‑149 (MQ=255)
  tctTCCGGCCACTGTCCTTTGGTTGCGGGTTCATAAGCCAGACCTTTACCTTTCGACATGCCGTAGGTGGTTTTCATACCAGGGATTCCGCACTGCACGCGGATGGCTTTGAATCCAAGCTCTTGATGACGGGCATAATCATCCAGAGc                                                                                                                                               >  3:230895/1‑149 (MQ=255)
   ctTCCGGCCACTGTCCTTTGGTTGCGGGTTCATAAGCCAGACCTTTACCTTTCGACATGCCGTAGGTGGTTTTCATACCAGGGATTCCGCACTGCACGCGGATGGCTTTGAATCCAAGCTCTTGATGACGGGCATAATCATCCAGAGCt                                                                                                                                              >  4:413372/1‑149 (MQ=255)
    ttCCGGCCACTGTCCTTTGGTTGCGGGTTCATAAGCCAGACCTTTACCTTTCGACATGCCGTAGGTGGTTTTCATACCAGGGATTCCGCACTGCACGCGGATGGCTTTGAATCCAAGCTCTTGATGACGGGCATAATCATCCAGAGCtt                                                                                                                                             >  1:181514/1‑149 (MQ=255)
    ttCCGGCCACTGTCCTTTGGTTGCGGGTTCATAAGCCAGACCTTTACCTTTCGACATGCCGTAGGTGGTTTTCATACCAGGGATTCCGCACTGCACGCGGATGGCTTTGAATCCAAGCTCTTGATGACGGGCATAATCATCCAGAGCtt                                                                                                                                             >  1:431033/1‑149 (MQ=255)
     tCCGGCCACTGTCCTTTGGTTGCGGGTTCATAAGCCAGACCTTTACCTTTCGACATGCCGTAGGTGGTTTTCATACCAGGGATTCCGCACTGCACGCGGATGGCTTTGAATCCAAGCTCTTGATGACGGGCATAATCATCCAGAGCTtc                                                                                                                                            <  2:308573/149‑1 (MQ=255)
       cGGCCACTGTCCTTTGGTTGCGGGTTCATAAGCCAGACCTTTACCTTTCGACATGCCGTAGGTGGTTTTCATACCAGGGATTCCGCACTGCACGCGGATGGCTTTGAATCCAAGCTCTTGATGACGGGCATAATCATCCAGAGCTtcat                                                                                                                                          <  1:32283/149‑1 (MQ=255)
       cGGCCACTGTCCTTTGGTTGCGGGTTCATAAGCCAGACCTTTACCTTTCGACATGCCGTAGGTGGTTTTCATACCAGGGATTCCGCACTGCACGCGGATGGCTTTGAATCCAAGCTCTTGATGACGGGCATAATCATCCAGAGCTtcat                                                                                                                                          >  3:162538/1‑149 (MQ=255)
        ggCCACTGTCCTTTGGTTGCGGGTTCATAAGCCAGACCTTTACCTTTCGACATGCCGTAGGTGGTTTTCATACCAGGGATTCCGCACTGCACGCGGATGGCTTTGAATCCAAGCTCTTGATGACGGGCATAATCATCCAGAGCTtcatc                                                                                                                                         <  2:167440/149‑1 (MQ=255)
            aCTGTCCTTTGGTTGCGGGTTCATAAGCCAGACCTTTACCTTTCGACATGCCGTAGGTGGTTTTCATACCAGGGATTCCGCACTGCACGCGGATGGCTTTGAATCCAAGCTCTTGATGACGGGCATAATCATCCAGAGCTTCATCAATa                                                                                                                                     <  4:162538/149‑1 (MQ=255)
            aCTGTCCTTTGGTTGCGGGTTCATAAGCCAGACCTTTACCTTTCGACATGCCGTAGGTGGTTTTCATACCAGGGATTCCGCACTGCACGCGGATGGCTTTGAATCCAAGCTCTTGATGACGGGCATAATCATCCAGAGCTTCATCAATa                                                                                                                                     >  1:562637/1‑149 (MQ=255)
                 ccTTTGGTTGCGGGTTCATAAGCCAGACCTTTACCTTTCGACATGCCGTAGGTGGTTTTCATACCAGGGATTCCGCACTGCACGCGGATGGCTTTGAATCCAAGCTCTTGATGACGGGCATAATCATCCAGAGCTTCATCAATACtgtg                                                                                                                                <  2:181514/149‑1 (MQ=255)
                 ccTTTGGTTGCGGGTTCATAAGCCAGACCTTTACCTTTCGACATGCCGTAGGTGGTTTTCATACCAGGGATTCCGCACTGCACGCGGATGGCTTTGAATCCAAGCTCTTGATGACGGGCATAATCATCCAGAGCTTCATCAATACtgtg                                                                                                                                <  2:395091/149‑1 (MQ=255)
                  cTTTGGTTGCGGGTTCATAAGCCAGACCTTTACCTTTCGACATGCCGTAGGTGGTTTTCATACCAGGGATTCCGCACTGCACGCGGATGGCTTTGAATCCAAGCTCTTGATGACGGGCATAATCATCCAGAGCTTCATCAATACTGTGa                                                                                                                               >  2:69270/1‑149 (MQ=255)
                   tttGGTTGCGGGTTCATAAGCCAGACCTTTACCTTTCGACATGCCGTAGGTGGTTTTCATACCAGGGATTCCGCACTGCACGCGGATGGCTTTGAATCCAAGCTCTTGATGACGGGCATAATCATCCAGAGCTTCATCAATACTGTGAc                                                                                                                              >  1:379911/1‑149 (MQ=255)
                   tttGGTTGCGGGTTCATAAGCCAGACCTTTACCTTTCGACATGCCGTAGGTGGTTTTCATACCAGGGATTCCGCACTGCACGCGGATGGCTTTGAATCCAAGCTCTTGATGACGGGCATAATCATCCAGAGCTTCATCAATACTGTGAc                                                                                                                              <  4:21306/149‑1 (MQ=255)
                   tttGGTTGCGGGTTCATAAGCCAGACCTTTACCTTTCGACATGCCGTAGGTGGTTTTCATACCAGGGATTCCGCACTGCACGCGGATGGCTTTGAATCCAAGCTCTTGATGACGGGCATAATCATCCAGAGCTTCATCAATACTGTGAc                                                                                                                              >  1:257185/1‑149 (MQ=255)
                    ttGGTTGCGGGTTCATAAGCCAGACCTTTACCTTTCGACATGCCGTAGGTGGTTTTCATACCAGGGATTCCGCACTGCACGCGGATGGCTTTGAATCCAAGCTCTTGATGACGGGCATAATCATCCAGAGCTTCATCAATACTGTGAcc                                                                                                                             >  1:492305/1‑149 (MQ=255)
                     tGGTTGCGGGTTCATAAGCCAGACCTTTACCTTTCGACATGCCGTAGGTGGTTTTCATACCAGGGATTCCGCACTGCACGCGGATGGCTTTGAATCCAAGCTCTTGATGACGGGCATAATCATCCAGAGCTTCATCAATACTGTGACCg                                                                                                                            <  4:503830/149‑1 (MQ=255)
                        ttGCGGGTTCATAAGCCAGACCTTTACCTTTCGACATGCCGTAGGTGGTTTTCATACCAGGGATTCCGCACTGCACGCGGATGGCTTTGAATCCAAGCTCTTGATGACGGGCATAATCATCCAGAGCTTCATCAATACTGTGACCggtg                                                                                                                         >  1:237582/1‑149 (MQ=255)
                               ttCATAAGCCAGACCTTTACCTTTCGACATGCCGTAGGTGGTTTTCATACCAGGGATTCCGCACTGCACGCGGATGGCTTTGAATCCAAGCTCTTGATGACGGGCATAATCATCCAGAGCTTCATCAATACTGTGACCGGTGGTATGGc                                                                                                                  <  4:457866/149‑1 (MQ=255)
                                  aTAAGCCAGACCTTTACCTTTCGACATGCCGTAGGTGGTTTTCATACCAGGGATTCCGCACTGCACGCGGATGGCTTTGAATCCAAGCTCTTGATGACGGGCATAATCATCCAGAGCTTCATCAATACTGTGACCGGTGGTATGGCAAt                                                                                                               >  3:532142/1‑149 (MQ=255)
                                  aTAAGCCAGACCTTTACCTTTCGACATGCCGTAGGTGGTTTTCATACCAGGGATTCCGCACTGCACGCGGATGGCTTTGAATCCAAGCTCTTGATGACGGGCATAATCATCCAGAGCTTCATCAATACTGTGACCGGTGGTATGGCAAt                                                                                                               >  3:481319/1‑149 (MQ=255)
                                  aTAAGCCAGACCTTTACCTTTCGACATGCCGTAGGTGGTTTTCATACCAGGGATTCCGCACTGCACGCGGATGGCTTTGAATCCAAGCTCTTGATGACGGGCATAATCATCCAGAGCTTCATCAATACTGTGACCGGTGGTATGGCAAt                                                                                                               >  4:90645/1‑149 (MQ=255)
                                   tAAGCCAGACCTTTACCTTTCGACATGCCGTAGGTGGTTTTCATACCAGGGATTCCGCACTGCACGCGGATGGCTTTGAATCCAAGCTCTTGATGACAGGCATAATCATCCAGAGCTTCATCAATACTGTGACCGGTGGTATGGCAATa                                                                                                              >  2:480472/1‑149 (MQ=255)
                                     aGCCAGACCTTTACCTTTCGACATGCCGTAGGTGGTTTTCATACCAGGGATTCCGCACTGCACGCGGATGGCTTTGAATCCAAGCTCTTGATGACGGGCATAATCATCCAGAGCTTCATCAATACTGTGACCGGTGGTATGGCAATaaa                                                                                                            <  1:95104/149‑1 (MQ=255)
                                     aGCCAGACCTTTACCTTTCGACATGCCGTAGGTGGTTTTCATACCAGGGATTCCGCACTGCACGCGGATGGCTTTGAATCCAAGCTCTTGATGACGGGCATAATCATCCAGAGCTTCATCAATACTGTGACCGGTGGTATGGCAATaaa                                                                                                            <  1:482457/149‑1 (MQ=255)
                                      gCCAGACCTTTACCTTTCGACATGCCGTAGGTGGTTTTCATACCAGGGATTCCGCACTGCACGCGGATGGCTTTGAATCCAAGCTCTTGATGACGGGCATAATCATCCAGAGCTTCATCAATACTGTGACCGGTGGTATGGCAATAAAc                                                                                                           >  3:106983/1‑149 (MQ=255)
                                        cAGACCTTTACCTTTCGACATGCCGTAGGTGGTTTTCATACCAGGGATTCCGCACTGCACGCGGATGGCTTTGAATCCAAGCTCTTGATGACGGGCATAATCATCCAGAGCTTCATCAATACTGTGACCGGTGGTATGGCAATAAACCa                                                                                                         >  4:409100/1‑149 (MQ=255)
                                             ctttacctttCGACATGCCGTAGGTGGTTTTCATACCAGGGATTCCGCACTGCACGCGGATGGCTTTGAATCCAAGCTCTTGATGACGGGCATAATCATCCAGAGCTTCATCAATACTGTGACCGGTGGTATGGCAATAAACCATCAcc                                                                                                    >  2:188366/1‑149 (MQ=255)
                                             ctttacctttCGACATGCCGTAGGTGGTTTTCATACCAGGGATTCCGCACTGCACGCGGATGGCTTTGAATCCAAGCTCTTGATGACGGGCATAATCATCCAGAGCTTCATCAATACTGTGACCGGTGGTATGGCAATAAACCATCAcc                                                                                                    >  1:499213/1‑149 (MQ=255)
                                                    tttCGACATGCCGTAGGTGGTTTTCATACCAGGGATTCCGCACTGCACGCGGATGGCTTTGAATCCAAGCTCTTGATGACGGGCATAATCATCCAGAGCTTCATCAATACTGTGACCGGTGGTATGGCAATAAACCATCACCCCTTCag                                                                                             >  1:417944/1‑148 (MQ=255)
                                                      tCGACATGCCGTAGGTGGTTTTCATACCAGGGATTCCGCACTGCACGCGGATGGCTTTGAATCCAAGCTCTTGATGACGGGCATAATCATCCAGAGCTTCATCAATACTGTGACCGGTGGTATGGCAATAAACCATCACCCCTTCACga                                                                                           <  1:570397/149‑1 (MQ=255)
                                                        gACATGCCGTAGGTGGTTTTCATACCAGGGATTCCGCACTGCACGCGGATGGCTTTGAATCCAAGCTCTTGATGACGGGCATAATCATCCAGAGCTTCATCAATACTGTGACCGGTGGTATGGCAATAAACCATCACCCCTTCACgaga                                                                                         <  2:477894/149‑1 (MQ=255)
                                                           aTGCCGTAGGTGGTTTTCATACCAGGGATTCCGCACTGCACGCGGATGGCTTTGAATCCAAGCTCTTGATGACGGGCATAATCATCCAGAGCTTCATCAATACTGTGACCGGTGGTATGGCAATAAACCATCACCCCTTCACGAGAggc                                                                                      >  2:189506/1‑149 (MQ=255)
                                                                 tAGGTGGTTTTCATACCAGGGATTCCGCACTGCACGCGGATGGCTTTGAATCCAAGCTCTTGATGACGGGCATAATCATCCAGAGCTTCATCAATACTGTGACCGGTGGTATGGCAATAAACCATCACCCCTTCACGAGACGCgccgcc                                                                                >  4:406124/1‑149 (MQ=255)
                                                                 tAGGTGGTTTTCATACCAGGGATTCCGCACTGCACGCGGATGGCTTTGAATCCAAGCTCTTGATGACGGGCATAATCATCCAGAGCTTCATCAATACTGTGACCGGTGGTATGGCAATAAACCATCACCCCTTCACGAGACGCgccgcc                                                                                >  3:415452/1‑149 (MQ=255)
                                                                  aGGTGGTTTTCATACCAGGGATTCCGCACTGCACGCGGATGGCTTTGAATCCAAGCTCTTGATGACGGGCATAATCATCCAGAGCTTCATCAATACTGTGACCGGTGGTATGGCAATAAACCATCACCCCTTCACGAGACGCGccgccg                                                                               <  4:291427/149‑1 (MQ=255)
                                                                   ggtggtTTTCATACCAGGGATTCCGCACTGCACGCGGATGGCTTTGAATCCAAGCTCTTGATGACGGGCATAATCATCCAGAGCTTCATCAATACTGTGACCGGTGGTATGGCAATAAACCATCACCCCTTCACGAGACGCGCCGCCGa                                                                              <  4:418236/149‑1 (MQ=255)
                                                                   ggtggtTTTCATACCAGGGATTCCGCACTGCACGCGGATGGCTTTGAATCCAAGCTCTTGATGACGGGCATAATCATCCAGAGCTTCATCAATACTGTGACCGGTGGTATGGCAATAAACCATCACCCCTTCACGAGACGCGCCGCCGa                                                                              <  3:399876/149‑1 (MQ=255)
                                                                   ggtggtTTTCATACCAGGGATTCCGCACTGCACGCGGATGGCTTTGAATCCAAGCTCTTGATGACGGGCATAATCATCCAGAGCTTCATCAATACTGTGACCGGTGGTATGGCAATAAACCATCACCCCTTCACGAGACGCGCCGCCGa                                                                              >  1:48128/1‑149 (MQ=255)
                                                                    gtggtTTTCATACCAGGGATTCCGCACTGCACGCGGATGGCTTTGAATCCAAGCTCTTGATGACGGGCATAATCATCCAGAGCTTCATCAATACTGTGACCGGTGGTATGGCAATAAACCATCACCCCTTCACGAGACGCGCCGCCGAg                                                                             >  4:173704/1‑149 (MQ=255)
                                                                      ggtTTTCATACCAGGGATTCCGCACTGCACGCGGATGGCTTTGAATCCAAGCTCTTGATGACGGGCATAATCATCCAGAGCTTCATCAATACTGTGACCGGTGGTATGGCAATAAACCATCACCCCTTCACGAGACGCGCCGCCGAGTa                                                                           >  2:579877/1‑149 (MQ=255)
                                                                        ttttCATACCAGGGATTCCGCACTGCACGCGGATGGCTTTGAATCCAAGCTCTTGATGACGGGCATAATCATCCAGAGCTTCATCAATACTGTGACCGGTGGTATGGCAATAAACCATCACCCCTTCACGAGACGCGCCGCCGAGTAAc                                                                         <  3:90645/149‑1 (MQ=255)
                                                                        ttttCATACCAGGGATTCCGCACTGCACGCGGATGGCTTTGAATCCAAGCTCTTGATGACGGGCATAATCATCCAGAGCTTCATCAATACTGTGACCGGTGGTATGGCAATAAACCATCACCCCTTCACGAGACGCGCCGCCGAGTAAc                                                                         <  1:261530/149‑1 (MQ=255)
                                                                                 cAGGGATTCCGCACTGCACGCGGATGGCTTTGAATCCAAGCTCTTGATGACGGGCATAATCATTCAGAGCTtc                                                                                                                                            <  4:314517/73‑1 (MQ=255)
                                                                                 cAGGGATTCCGCACTGCACGCGGATGGCTTTGAATCCAAGCTCTTGATGACGGGCATAATCATTCAGAGCTtc                                                                                                                                            >  1:234604/1‑73 (MQ=255)
                                                                                 cAGGGATTCCGCACTGCACGCGGATGGCTTTGAATCCAAGCTCTTGATGACGGGCATAATCATTCAGAGCTtc                                                                                                                                            >  3:314517/1‑73 (MQ=255)
                                                                                 cAGGGATTCCGCACTGCACGCGGATGGCTTTGAATCCAAGCTCTTGATGACGGGCATAATCATTCAGAGCTtc                                                                                                                                            <  2:234604/73‑1 (MQ=255)
                                                                                 cAGGGATTCCGCACTGCACGCGGATGGCTTTGAATCCAAGCTCTTGATGACGGGCATAATCATCCAGAGCTTCATCAATACTGTGACCGGTGGTATGGCAATAAACCATCACCCCTTCACGAGACGCGCCGCCGAGTAACTGGTAAAGc                                                                <  2:519206/149‑1 (MQ=255)
                                                                                 cAGGGATTCCGCACTGCACGCGGATGGCTTTGAATCCAAGCTCTTGATGACGGGCATAATCATCCAGAGCTTCATCAATACTGTGACCGGTGGTATGGCAATAAACCATCACCCCTTCACGAGACGCGCCGCCGAGTAACTGGTAAAGc                                                                <  2:377901/149‑1 (MQ=255)
                                                                                      aTTCCGCACTGCACGCGGATGGCTTTGAATCCAAGCTCTTGATGACGGGCATAATCATCCAGAGCTTCATCAATACTGTGACCGGTGGTATGGCAATAAACCATCACCCCTTCACGAGACGCGCCGCCGAGTAACTGGTAAAGCGGCAt                                                           <  1:562480/149‑1 (MQ=255)
                                                                                           gCACTGCACGCGGATGGCTTTGAATCCAAGCTCTTGATGACGGGCATAATCATCCAGAGCTTCATCAATACTGTGACCGGTGGTATGGCAATAAACCATCACCCCTTCACGAGACGCGCCGCCGAGTAACTTGTAAAGCGGCATGTTgg                                                      <  4:214083/149‑1 (MQ=255)
                                                                                                    gcgGATGGCTTTGAATCCAAGCTCTTGATGACGGGCATAATCATCCAGAGCTTCATCAATACTGTGACCGGTGGTATGGCAATAAACCATCACCCCTTCACGAGACGCGCCGCCGAGTAACTGGTAAAGCGGCATGTTGGCAGCTTTgg                                             <  4:396357/149‑1 (MQ=255)
                                                                                                       gATGGCTTTGAATCCAAGCTCTTGATGACGGGCATAATCATCCAGAGCTTCATCAATACTGTGACCGGTGGTATGGCAATAAACCATCACCCCTTCACGAGACGCGCCGCCGAGTAACTGGTAAAGCGGCATGTTGGCAGCTTTGGCtt                                          <  1:154538/149‑1 (MQ=255)
                                                                                                                  aTCCAAGCTCTTGATGACGGGCATAATCATCCAGAGCTTCATCAATACTGTGACCGGTGGTATGGCAATAAACCATCACCCCTTCACGAGACGCGCCGCCGAGTAACTGGTAAAGCGGCATGTTGGCAGCTTTGGCTTTAATATCCcac                               >  1:179333/1‑149 (MQ=255)
                                                                                                                                tgaCGGGCATAATCATCCAGAGCTTCATCAATACTGTGACCGGTGGTATGGCAATAAACCATCACCCCTTCACGAGACGCGCCGCCGAGTAACTGGTAAAGCGGCATGTTGGCAGCTTTGGCTTTAATATCCCACAGCGCCATATCAAc                 <  2:38672/149‑1 (MQ=255)
                                                                                                                                  aCGGGCATAATCATCCAGAGCTTCATCAATACTGTGACCGGTGGTATGGCAATAAACCATCACCCCTTCACGAGACGCGCCGCCGAGTAACTGGTAAAGCGGCATGTTGGCAGCTTTGGCTTTAATATCCCACAGCGCCATATCAACCg               <  4:375235/149‑1 (MQ=255)
                                                                                                                                    gggCATAATCATCCAGAGCTTCATCAATACTGTGACCGGTGGTATGGCAATAAACCATCACCCCTTCACGAGACGCGCCGCCGAGTAACTGGTAAAGCGGCATGTTGGCAGCTTTGGCTTTAATATCCCACAGCGCCATATCAACCGCt             >  4:108284/1‑149 (MQ=255)
                                                                                                                                           atcatcCAGAGCTTCATCAATACTGTGACCGGTGGTATGGCAATAAACCATCACCCCTTCACGAGACGCGCCGCCGAGTAACTGGTAAAGCGGCATGTTGGCAGCTTTGGCTTTAATATCCCACAGCGCCATATCAACCGCTGAAATgg      <  3:173704/149‑1 (MQ=255)
                                                                                                                                            tcatcCAGAGCTTCATCAATACTGTGACCGGTGGTATGGCAATAAACCATCACCCCTTCACGAGACGCGCCGCCGAGTAACTGGTAAAGCGGCATGTTGGCAGCTTTGGCTTTAATATCCCACAGCGCCATATCAACCGCTGAAATGGc     <  3:8783/149‑1 (MQ=255)
                                                                                                                                            tcatcCAGAGCTTCATCAATACTGTGACCGGTGGTATGGCAATAAACCATCACCCCTTCACGAGACGCGCCGCCGAGTAACTGGTAAAGCGGCATGTTGGCAGCTTTGGCTTTAATATCCCACAGCGCCATATCAACCGCTGAAATGGc     >  1:456738/1‑149 (MQ=255)
                                                                                                                                            tcatcCAGAGCTTCATCAATACTGTGACCGGTGGTATGGCAATAAACCATCACCCCTTCACGAGACGCGCCGCCGAGTAACTGGTAAAGCGGCATGTTGGCAGCTTTGGCTTTAATATCCCACAGCGCCATATCAACCGCTGAAATGGc     <  1:227284/149‑1 (MQ=255)
                                                                                                                                             catcCAGAGCTTCATCAATACTGTGACCGGTGGTATGGCAATAAACCATCACCCCTTCACGAGACGCGCCGCCGAGTAACTGGTAAAGCGGCATGTTGGCAGCTTTGGCTTTAATATCCCACAGCGCCATATCAACCGCTGAAATGGcc    <  4:346096/149‑1 (MQ=255)
                                                                                                                                               tcCAGAGCTTCATCAATACTGTGACCGGTGGTATGGCAATAAACCATCACCCCTTCACGAGACGCGCCGCCGAGTAACTGGTAAAGCGGCATGTTGGCAGCTTTGGCTTTAATATCCCACAGCGCCATATCAACCGCTGAAATGGCCGa  <  3:206360/149‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                   
GCTCTTCCGGCCACTGTCCTTTGGTTGCGGGTTCATAAGCCAGACCTTTACCTTTCGACATGCCGTAGGTGGTTTTCATACCAGGGATTCCGCACTGCACGCGGATGGCTTTGAATCCAAGCTCTTGATGACGGGCATAATCATCCAGAGCTTCATCAATACTGTGACCGGTGGTATGGCAATAAACCATCACCCCTTCACGAGACGCGCCGCCGAGTAACTGGTAAAGCGGCATGTTGGCAGCTTTGGCTTTAATATCCCACAGCGCCATATCAACCGCTGAAATGGCCGA  >  NC_000913/1654601‑1654892

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