Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,987,175 C→G 5.3% M236I (ATG→ATC hemY ← protein HemY

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,987,1750CG 5.3% 505.2 / 12.0 148M236I (ATG→ATChemYprotein HemY
Reads supporting (aligned to +/- strand):  ref base C (67/73);  new base G (4/4);  total (71/77)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.86e-01

Only 100 of 152 total aligned reads displayed.

ATAATTTGCTGGGCAGTATCATGATCGTCACATTCAATAAGATGTTCCGCCATTGCCACCTGCAACGCTACCTGATGACGCGTTTTCCGGCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTCGCCATTGATGGGATAATATCCAGCAGCGAACTCCATGCACCTGTGCGGATATACGCCTGTTCCGCCAGACGTAATACTTCCG  >  NC_000913/3987026‑3987323
                                                                                                                                                     |                                                                                                                                                    
atAATTTGCTGGGCAGTATCATGATCGTCACATTCAATAAGATGTTCCGCCATTGCCACCTGCAACGCTACCTGATGACGCGTTTTCCGGCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATcc                                                                                                                                                      <  2:1189207/150‑1 (MQ=255)
atAATTTGCTGGGCAGTATCATGATCGTCACATTCAATAAGATGTTCCGCCATTGCCACCTGCAACGCTACCTGATGACGCGTTTTCCGGCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATcc                                                                                                                                                      >  2:60449/1‑150 (MQ=255)
atAATTTGCTGGGCAGTATCATGATCGTCACATTCAATAAGATGTTCCGCCATTGCCACATGCAACGCTACCTGATGACGCGTTTTCCTGGTTTGGTTTTTCGACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGGGCCTGATcc                                                                                                                                                      <  2:1106914/150‑1 (MQ=255)
 tAATTTGCTGGGCAGTATCATGATCGTCACATTCAATAAGATGTTCCGCCATTGCCACCTGCAACGCTACCTGATGACGCGTTTTCCGGCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCa                                                                                                                                                     >  1:954615/1‑150 (MQ=255)
   aTTTGCTGGGCAGTATCATGATCGTCACATTCAATAAGATGTTCCGCCATTGCCACCTGCAACGCTACCTGATGACGCGTTTTCCGGCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATc                                                                                                                                                   >  1:154070/1‑150 (MQ=255)
     ttGCTGGGCAGTATCATGATCGTCACATTCAATAAGATGTTCCGCCATTGCCACCTGCAACGCTACCTGATGACGCGTTTTCCGGCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAg                                                                                                                                                 <  1:839974/150‑1 (MQ=255)
     ttGCTGGGCAGTATCATGATCGTCACATTCAATAAGATGTTCCGCCATTGCCACCTGCAACGCTACCTGATGACGCGTTTTCCGGCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAg                                                                                                                                                 <  2:843368/150‑1 (MQ=255)
        cTGGGCAGTATCATGATCGTCACATTCAATAAGATGTTCCGCCATTGCCACCTGCAACGCTACCTGATGACGCGTTTTCCGGCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGcc                                                                                                                                              >  2:1825536/1‑150 (MQ=255)
        cTGGGCAGTATCATGATCGTCACATTCAATAAGATGTTCCGCCATTGCCACCTGCAACGCTACCTGATGACGCGTTTTCCGGCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATACATCAGGcc                                                                                                                                              >  2:1502924/1‑150 (MQ=255)
         tGGGCAGTATCATGATCGTCACATTCAATAAGATGTTCCGCCATTGCCACCTGCAACGCTACCTGATGACGCGTTTTCCGGCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCa                                                                                                                                             <  2:627451/150‑1 (MQ=255)
         tGGGCAGTATCATGATCGTCACATTCAATAAGATGTTCCGCCATTGCCACCTGCAACGCTACCTGATGACGCGTTTTCCGGCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCa                                                                                                                                             >  2:1801642/1‑150 (MQ=255)
          gggCAGTATCATGATCGTCACATTCAATAAGATGTTCCGCCATTGCCACCTGCAACGCTACCTGATGACGCGTTTTCCGGCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATGGATCAGGcgcg                                                                                                                                            <  2:65160/150‑4 (MQ=255)
          gggCAGTATCATGATCGTCACATTCAATAAGATGTTCCGCCATTGCCACCTGCAACGCTACCTGATGACGCGTTTTCCGGCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATGGATCAGGcgcg                                                                                                                                            <  1:1722807/150‑4 (MQ=255)
           ggCAGTATCATGATCGTCACATTCAATAAGATGTTCCGCCATTGCCACCTGCAACGCTACCTGATGACGCGTTTTCCGGCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAAt                                                                                                                                           <  2:448700/150‑1 (MQ=255)
           ggCAGTATCATGATCGTCACATTCAATAAGATGTTCCGCCATTGCCACCTGCAACGCTACCTGATGACGCGTTTTCCGGCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAAt                                                                                                                                           >  1:832479/1‑150 (MQ=255)
             cAGTATCATGATCGTCACATTCAATAAGATGTTCCGCCATTGCCCCCTGCAACGCTACCTGATGACGCGTTTTCCGGCTGTGCGTTGTCCCCCGGATACGCAAACCTTCGCTACCGTGATCGGCACGCGCCAGATCCATCAGGCCAATcc                                                                                                                                         >  2:1363759/1‑150 (MQ=255)
             cAGTATCATGATCGTCACATTCAATAAGATGTTCCGCCATTGCCACCTGCAACGCTACCTGATGACGCGTTTTCCGGCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATGGATCAGGcgcgtgc                                                                                                                                         <  1:711970/150‑7 (MQ=255)
                tATCATGATCGTCACATTCAATAAGATGTTCCGCCATTGCCACCTGCAACGCTACCTGATGACGCGTTTTCCGGCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATg                                                                                                                                      <  2:2296788/150‑1 (MQ=255)
                 aTCATGATCGTCACATTCAATAAGATGTTCCGCCATTGCCACCTGCAACGCTACCTGATGACGCGTTTTCCGGCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGc                                                                                                                                     >  2:2250835/1‑150 (MQ=255)
                  tCATGATCGTCACATTCAATAAGATGTTCCGCCATTGCCACCTGCAACGCTACCTGATGACGCGTTTTCCGGCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGcc                                                                                                                                    >  1:531390/1‑150 (MQ=255)
                  tCATGATCGTCACATTCAATAAGATGTTCCGCCATTGCCACCTGCAACGCTACCTGATGACGCGTTTTCCGGCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGcc                                                                                                                                    <  2:1457544/150‑1 (MQ=255)
                        tCGTCACATTCAATAAGATGTTCCGCCATTGCCACCTGCAACGCTACCTGATGACGCGTTTTCCGGCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCtgttgt                                                                                                                              >  1:676476/1‑150 (MQ=255)
                             acaTTCAATAAGATGTTCCGCCATTGCCACCTGCAACGCTACCTGATGACGCGTTTTCCGGCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAg                                                                                                                         <  2:1089651/150‑1 (MQ=255)
                             acaTTCAATAAGATGTTCCGCCATTGCCACCTGCAACGCTACCTGATGACGCGTTTTCCGGCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAg                                                                                                                         <  1:1859317/150‑1 (MQ=255)
                              caTTCAATAAGATGTTCCGCCATTGCCACCTGCAACGCTACCTGATGACGCGTTTTCCGGCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGc                                                                                                                        >  1:253539/1‑150 (MQ=255)
                                 tCAATAAGATGTTCCGCCATTGCCACCTGCAACGCTACCTGATGACGCGTTTTCCGGCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCAtt                                                                                                                     <  1:1314317/150‑1 (MQ=255)
                                  cAATAAGATGTTCCGCCATTGCCACCTGCAACGCTACCTGATGACGCGTTTTCCGGCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTg                                                                                                                    >  1:1471178/1‑150 (MQ=255)
                                    aTAAGATGTTCCGCCATTGCCACCTGCAACGCTACCTGATGACGCGTTTTCCGGCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCa                                                                                                                  >  1:984728/1‑150 (MQ=255)
                                       aGATGTTCCGCCATTGCCACCTGCAACGCTACCTGATGACGCGTTTTCCGGCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGa                                                                                                               >  2:18076/1‑150 (MQ=255)
                                           gTTCCGCCATTGCCACCTGCAACGCTACCTGATGACGCGTTTTCCGGCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGtt                                                                                                           >  1:1822230/1‑150 (MQ=255)
                                             tCCGCCATTGCCACCTGCAACGCTACCTGATGACGCGTTTTCCGGCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGttct                                                                                                         <  1:1210763/150‑1 (MQ=255)
                                                gCCATTGCCACCTGCAACGCTACCTGATGACGCGTTTTCCGGCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTtca                                                                                                      <  2:1422824/150‑1 (MQ=255)
                                                  cATTGCCACCTGCAACGCTACCTGATGACGCGTTTTCCGGCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTtcatc                                                                                                    <  1:934803/150‑1 (MQ=255)
                                                      gCCACCTGCAACGCTACCTGATGACGCGTTTTCCGGCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCa                                                                                                <  2:1620557/150‑1 (MQ=255)
                                                          ccTGCAACGCTACCTGATGACGCGTTTTCCGGCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACAt                                                                                            >  1:1640103/1‑150 (MQ=255)
                                                            tGCAACGCTACCTGATGACGCGTTTTCCGGCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATgg                                                                                          >  1:2236638/1‑150 (MQ=255)
                                                             gCAACGCTACCTGATGACGCGTTTTCCGGCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATggg                                                                                         <  2:1034339/150‑1 (MQ=255)
                                                             gCAACGCTACCTGATGACGCGTTTTCCGGCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATggg                                                                                         >  2:778631/1‑150 (MQ=255)
                                                              cAACGCTACCTGATGACGCGTTTTCCGGCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGGCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGc                                                                                        <  2:134236/150‑1 (MQ=255)
                                                               aaCGCTACCTGATGACGCGTTTTCCGGCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCt                                                                                       <  2:1756127/150‑1 (MQ=255)
                                                                    tACCTGATGACGCGTTTTCCGGCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTCGc                                                                                  >  1:368035/1‑150 (MQ=255)
                                                                        tgatgaCGCGTTTTCCGGCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTCGCCAtt                                                                              <  2:253539/150‑1 (MQ=255)
                                                                         gatgaCGCGTTTTCCGGCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTCGCCATTg                                                                             <  1:406961/150‑1 (MQ=255)
                                                                         gatgaCGCGTTTTCCGGCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTCGCCATTg                                                                             >  2:461385/1‑150 (MQ=255)
                                                                           tgaCGCGTTTTCCGGCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTCGCCATTGAt                                                                           >  2:577410/1‑150 (MQ=255)
                                                                            gaCGCGTTTTCCGGCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTCGCCATTGATg                                                                          >  1:2150442/1‑150 (MQ=255)
                                                                                cgTTTTCCGGCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTCGCCATTGATGGGat                                                                      <  2:1272181/150‑1 (MQ=255)
                                                                                cgTTTTCCGGCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTCGCCATTGATGGGat                                                                      >  2:936376/1‑150 (MQ=255)
                                                                                cgTTTTCCGGCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTCGCCATTGATGGGat                                                                      <  1:60449/150‑1 (MQ=255)
                                                                                  ttttCCGGCTTTTGTTTTTCCTCCAGTTAGGCAAACCTTCGATACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATGCATGCCTGTTGTTCCAGCATTGCACGATGTTATTCATCACCAACATGGGCTTTCGCCATTGATGGTataa                                                                    <  2:1185964/150‑1 (MQ=255)
                                                                                  ttttCCGGCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTCGCCATTGATGGGataa                                                                    <  1:103099/150‑1 (MQ=255)
                                                                                  ttttCCGGCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTCGCCATTGATGGGataa                                                                    <  1:121640/150‑1 (MQ=255)
                                                                                   tttCCGGCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTCGCCATTGATGGGataat                                                                   <  2:954615/150‑1 (MQ=255)
                                                                                    ttCCGGCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTCGCCATTGATGGGataata                                                                  <  1:75390/150‑1 (MQ=255)
                                                                                      ccGGCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTCGCCATTGATGGGATAATATc                                                                >  2:1711496/1‑150 (MQ=255)
                                                                                        ggCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTCGCCATTGATGGGATAATATcaa                                                              >  2:2295500/1‑148 (MQ=255)
                                                                                        ggCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTCGCCATTGATGGGATAATATCca                                                              >  2:1281621/1‑150 (MQ=255)
                                                                                        ggCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTCGCCATTGATGGGATAATATCca                                                              >  1:1716320/1‑150 (MQ=255)
                                                                                        ggCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTCGCCATTGATGGGATAATATCca                                                              >  1:209048/1‑150 (MQ=255)
                                                                                          cTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTCGCCATTGATGGGATAATATCcagc                                                            >  2:2043309/1‑150 (MQ=255)
                                                                                          cTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTCGCCATTGATGGGATAATATCcagc                                                            >  1:550607/1‑150 (MQ=255)
                                                                                          cTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGACACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTCGCCATTGATGGGATAATATCcagc                                                            >  2:1586939/1‑150 (MQ=255)
                                                                                            ttGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTCGCCATTGATGGGATAATATCcagcag                                                          >  2:1839877/1‑150 (MQ=255)
                                                                                              ggTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTCGCCATTGATGGGATAATATCCAGCAGCg                                                        <  1:1822524/150‑1 (MQ=255)
                                                                                               gTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTCGCCATTGATGGGATAATATCCAGCAGCGa                                                       <  1:2423015/150‑1 (MQ=255)
                                                                                               gTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTCGCCATTGATGGGATAATATCCAGCAGCGa                                                       >  1:1372330/1‑150 (MQ=255)
                                                                                                    tCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTCGCCATTGATGGGATAATATCCAGCAGCGAACTcc                                                  <  2:902243/150‑1 (MQ=255)
                                                                                                     ccaccaGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTCGCCATTGATGGGATAATATCCAGCAGCGAACTCCa                                                 >  2:1469049/1‑150 (MQ=255)
                                                                                                     ccaccaGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTCGCCATTGATGGGATAATATCCAGCAGCGAACTCCa                                                 <  2:26396/150‑1 (MQ=255)
                                                                                                      caccaGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTCGCCATTGATGGGATAATATCCAGCAGCGAACTCCAt                                                >  2:1871870/1‑150 (MQ=255)
                                                                                                        ccaGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTCGCCATTGATGGGATAATATCCAGCAGCGAACTCCATGc                                              >  2:1511182/1‑150 (MQ=255)
                                                                                                         caGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTCGCCATTGATGGGATAATATCCAGCAGCGAACTCCATGCa                                             <  1:2434047/150‑1 (MQ=255)
                                                                                                          aGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTCGCCATTGATGGGATAATATCCAGCAGCGAACTCCATGCAc                                            >  1:885140/1‑150 (MQ=255)
                                                                                                           tttACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTCGCCATTGATGGGATAATATCCAGCAGCGAACTCCATGCAcc                                           <  2:376242/149‑1 (MQ=255)
                                                                                                           gTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTCGCCATTGATGAGATAATATCCAGCAGCGAACTCCATGCAcc                                           <  2:980226/150‑1 (MQ=255)
                                                                                                              aCGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTCGCCATTGATGGGATAATATCCAGCAGCGAACTCCATGCACCtgt                                        >  1:1424176/1‑150 (MQ=255)
                                                                                                                 cAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTCGCCATTGATGGGATAATATCCAGCAGCGAACTCCATGCACCTGTGCg                                     <  1:1032291/150‑1 (MQ=255)
                                                                                                                 cAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTCGCCATTGATGGGATAATATCCAGCAGCGAACTCCATGCACCTGTGCg                                     >  1:791776/1‑150 (MQ=255)
                                                                                                                    aCCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTCGCCATTGATGGGATAATATCCAGCAGCGAACTCCATGCACCTGTGCGGat                                  >  1:888576/1‑150 (MQ=255)
                                                                                                                     ccTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTCGCCATTGATGGGATAATATCCAGCAGCGAACTCCATGCACCTGTGCGGata                                 <  1:282144/150‑1 (MQ=255)
                                                                                                                         cGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTCGCCATTGATGGGATAATATCCAGCAGCGAACTCCATGCACCTGTGCGGATATACg                             <  1:1132087/150‑1 (MQ=255)
                                                                                                                           cTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTCGCCATTGATGGGATAATATCCAGCAGCGAACTCCATGCACCTGTGCGGATATACGcc                           <  2:2236638/150‑1 (MQ=255)
                                                                                                                            tACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTCGCCATTGATGGGATAATATCCAGCAGCGAACTCCATGCACCTGTGCGGATATACGCCt                          <  1:150846/150‑1 (MQ=255)
                                                                                                                                  tATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTCGCCATTGATGGGATAATATCCAGCAGCGAACTCCATGCACCTGTGCGGATATACGCCTGTTCCg                    <  2:1936970/150‑1 (MQ=255)
                                                                                                                                     cGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTCGCCATTGATGGGATAATATCCAGCAGCGAACTCCATGCACCTGTGCGGATATACGCCTGTTCCGCCa                 >  1:663046/1‑150 (MQ=255)
                                                                                                                                      ggCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTCGCCATTGATGGGATAATATCCAGCAGCGAACTCCATGCACCTGTGCGGATATACGCCTGTTCCGCCAg                <  1:2332607/150‑1 (MQ=255)
                                                                                                                                      ggCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTCGCCATTGATGGGATAATATCCAGCAGCGAACTCCATGCACCTGTGCGGATATACGCCTGTTCCGCCAg                >  1:2044136/1‑150 (MQ=255)
                                                                                                                                       ggattggcCTGATGGATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTCGCCATTGATGGGATAATATCCAGCAGCGAACTCCATGCACCTGTGCGGATATACGCCTGTTCCGCCAGa               >  2:82260/7‑150 (MQ=255)
                                                                                                                                       gCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTCGCCATTGATGGGATAATATCCAGCAGCGAACTCCATGCACCTGTGCGGATATACGCCTGTTCCGCCAGa               <  1:1316261/150‑1 (MQ=255)
                                                                                                                                        cACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTCGCCATTGATGGGATAATATCCAGCAGCGAACTCCATGCACCTGTGCGGATATACGCCTGTTCCGCCAGAc              >  1:1195092/1‑150 (MQ=255)
                                                                                                                                        cACGCGCCTGATCCATCAGGCAAATACATGCCTGTTGTTCCAGCATTGCACGATGTTCTTAGTCACAAACATGGGCTTTCGCCATTGATGGGATAATATACAGCAGCGAACTCCATGCACCTGTGCGGATATACGCCTGTTCCGCAAGAc              >  2:274800/1‑150 (MQ=255)
                                                                                                                                         attggcCTGATGGATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTATTCATCACCAACATGGGCTTTCGCCATTGATGGGATAATATACAGCAGCGAACTCCATGCACCTGTGCGGATATACGCCTGTTCCGCCAGaag             >  2:833741/5‑148 (MQ=255)
                                                                                                                                           tggcCTGATGGATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTCGCCATTGATGGGATAATATCCAGCAGCGAACTCCATGCACCTGTGCGGATATACGCCTGTTCCGCCAGACGTa           >  2:1853604/3‑150 (MQ=255)
                                                                                                                                           gcgcCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTCGCCATTGATGGGATAATATCCAGCAGCGAACTCCATGCACCTGTGCGGATATACGCCTGTTCCGCCAGACGTa           <  2:154070/150‑1 (MQ=255)
                                                                                                                                            cgcCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTCGCCATTGATGGGATAATATCCAGCAGCGAACTCCATGCACCTGTGCGGATATACGCCTGTTCCGCCAGACGTaa          >  1:1849630/1‑150 (MQ=255)
                                                                                                                                             gcCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTCGCCATTGATGGGATAATATCCAGCAGCGAACTCCATGCACCTGTGCGGATATACGCCTGTTCCGCCAGACGTAAt         <  1:988740/150‑1 (MQ=255)
                                                                                                                                               cTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTCGCCATTGATGGGATAATATCCAGCAGCGAACTCCATGCACCTGTGCGGATATACGCCTGTTCCGCCAGACGTAATAc       >  1:1809587/1‑150 (MQ=255)
                                                                                                                                                tGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTCGCCATTGATGGGATAATATCCAGCAGCGAACTCCATGCACCTGTGCGGATATACGCCTGTTCCGCCAGACGTAATACt      <  2:676476/150‑1 (MQ=255)
                                                                                                                                                   tCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTCGCCATTGATGGGATAATATCCAGCAGCGAACTCCATGCACCTGTGCGGATATACGCCTGTTCCGCCAGACGTAATACTTcc   <  1:18076/150‑1 (MQ=255)
                                                                                                                                                    ccATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTCGCCATTGATGGGATAATATCCAGCAGCGAACTCCATGCACCTGTGCGGATATACGCCTGTTCCGCCAGACGTAATACTTCCg  <  2:2041284/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
ATAATTTGCTGGGCAGTATCATGATCGTCACATTCAATAAGATGTTCCGCCATTGCCACCTGCAACGCTACCTGATGACGCGTTTTCCGGCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTCGCCATTGATGGGATAATATCCAGCAGCGAACTCCATGCACCTGTGCGGATATACGCCTGTTCCGCCAGACGTAATACTTCCG  >  NC_000913/3987026‑3987323

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