Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,560,811 A→G 5.2% L40P (CTG→CCG) ‡ ddpB ← D,D‑dipeptide ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,560,8110AG 5.2% 361.1 / 6.4 116L40P (CTG→CCG) ‡ddpBD,D‑dipeptide ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base A (56/54);  new base G (3/3);  total (59/57)
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.17e-01

Only 100 of 118 total aligned reads displayed.

GCACCGGACGCCCGGTACGAATGGATGTTCCCAGGTCACCATGAAACAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGTCCAGTCCCAGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTAATCACGCAGACGCCCGCCACCACGAGCACCAGCCCCCAGCAGCGTTGGCGTAAAATACTCCAGAAGGTCACGTTATTTACTCATGGTATTGATATTGAAG  >  NC_000913/1560662‑1560958
                                                                                                                                                     |                                                                                                                                                   
gCACCGGACGCCCGGTACGAATGGATGTTCCCAGGTCACCATGAAACAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGTCCAGTCCCAGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGTCACCCGCCa                                                                                                                                                     >  1:606257/1‑150 (MQ=255)
gCACCGGACGCCCGGTACGAATGGATGTTCCCAGGTCACCATGAAACAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGTCCAGTCCCAGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGTCACCCGCCa                                                                                                                                                     <  2:1090617/150‑1 (MQ=255)
 cACCGGACGCCCGGTACGAATGGATGTTCCCAGGTCACCATGAAACAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGTCCAGTCCCAGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGTCACCCGCCAg                                                                                                                                                    <  2:1547563/150‑1 (MQ=255)
     ggACGCCCGGTACGAATGGATGTTCCCAGGTCACCATGAAACAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGTCCAGTCCCAGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGTCACCCGCCAGTAAc                                                                                                                                                >  1:1506062/1‑150 (MQ=255)
          cccGGTACGAATGGATGTTCCCAGGTCACCATGAAACAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGTCCAGTCCCAGTTGCTGGCGAATATTTTCCACGATTGCATCGCTGGCGCGGTCACCCGCCAGTAACCgcgc                                                                                                                                           >  1:967158/1‑150 (MQ=255)
          cccGGTACGAATGGATGTTCCCAGGTCACCATGAAACAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGTCCAGTCCCAGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGTTACTGGCGGGTGACCgcgc                                                                                                                                           <  2:485346/150‑1 (MQ=255)
          cccGGTACGAATGGATGTTCCCAGGTCACCATGAAACAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGTCCAGTCCCAGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGTCACCCGCCAGTAACCgcgc                                                                                                                                           <  2:847928/150‑1 (MQ=255)
          cccGGTACGAATGGATGTTCCCAGGTCACCATGAAACAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGTCCAGTCCCAGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGTCACCCGCCAGTAACCgcgc                                                                                                                                           >  1:1583210/1‑150 (MQ=255)
          cccGGAACGAATGGATGTTCCCAGGTCACCATGAAACAGGTCGCTGACGGAGCGGTAAAACTGTACATACAGTGGCTGGTCCAGTCCCAGTAGCTGGCGAATATTATACACGAAAGCATCGCTGGCGCGGACACCCGCCACTAACCgcgc                                                                                                                                           >  1:1884966/1‑150 (MQ=255)
              gTACGAATGGATGTTCCCAGGTCACCATGAAACAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGTCCAGTCCCAGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCGGa                                                                                                                                       >  1:644262/1‑150 (MQ=255)
                aCGAATGGATGTTCCCAGGTCACCATGAAACAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGTCCAGTCCCAGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCGGATc                                                                                                                                     <  2:1357953/150‑1 (MQ=255)
                  gAATGGATGTTCCCAGGTCACCATGAAACAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGTCCAGTCCCAGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGTTACTGGCGGGTGACCGCGccagcgatg                                                                                                                                   <  2:394680/150‑8 (MQ=255)
                  gAATGGATGTTCCCAGGTCACCATGAAACAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGTCCAGTCCCAGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCAc                                                                                                                                   >  2:144821/1‑150 (MQ=255)
                  gAATGGATGTTCCCAGGTCACCATGAAACAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGTCCAGTCCCAGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCAc                                                                                                                                   >  2:1295105/1‑150 (MQ=255)
                     tGGATGTTCCCAGGTCACCATGAAACAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGTCCAGTCCCAGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCACCAg                                                                                                                                <  2:379282/150‑1 (MQ=255)
                       gATGGTCCCAGGTCACCATGAAACAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGTCCAGTCCCAGGTGCTGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCACCAggg                                                                                                                              <  1:1017025/150‑1 (MQ=255)
                         tGTTCCCAGGTCACCATGAAACAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGTCCAGTCCCAGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGAt                                                                                                                            <  2:1447484/150‑1 (MQ=255)
                            tCCCAGGTCACCATGAAACAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGTCCAGTCCCAGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAg                                                                                                                         >  2:748596/1‑150 (MQ=255)
                             cccAGGTCACCATGAAACAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGTCCAGGCCCAGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGa                                                                                                                        <  1:1643475/150‑1 (MQ=255)
                             cccAGGTCACCATGAAACAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCGGGTCCAGTCCCAGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGGCGCGCTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGa                                                                                                                        >  2:276935/1‑150 (MQ=255)
                              ccAGGTCACCATGAAACAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGTCCAGTCCCAGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGAt                                                                                                                       >  2:917182/1‑150 (MQ=255)
                               cAGGTCACCATGAAACAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGTCCAGTCCCAGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATg                                                                                                                      >  1:699042/1‑150 (MQ=255)
                               cAGGTCACCATGAAACAGGTCGCTGACGTAGCGGGAAAACTGTACATACAGTGGCTGGTCCAGTCACAGTTGCTGGCGAATATTTTCCACGATCGCATCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATg                                                                                                                      <  1:1161875/150‑1 (MQ=255)
                                 ggTCACCATGAAACAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGTCCAGTCCCAGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCACCCGGGATCAGATGCg                                                                                                                    >  2:244309/1‑150 (MQ=255)
                                 ggTCACCATGAAACAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGTCCAGTCCCAGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCg                                                                                                                    <  2:639321/150‑1 (MQ=255)
                                  gTCACCATGAAACAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGTCCAGTCCCAGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCga                                                                                                                   <  2:1795077/150‑1 (MQ=255)
                                  gTCACCATGAAACAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGTCCAGTCCCAGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCga                                                                                                                   >  2:395478/1‑150 (MQ=255)
                                   tCACCATGAAACAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGTCCAGTCCCAGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCgag                                                                                                                  >  1:793221/1‑150 (MQ=255)
                                   tCACCATGAAACAGGGCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGTCCAGTCCCAGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCgag                                                                                                                  <  2:1484045/150‑1 (MQ=255)
                                     cccATGAACCATGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGTCCAGTCCCAGTTTCTGGCGAATATTTTTCACGATAGCATCGCTGGCGCGGTCAGCCGACAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGat                                                                                                                <  1:795328/149‑1 (MQ=255)
                                      ccATGAAACAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGTCCAGTCCCAGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGAGAGata                                                                                                               >  2:1015782/1‑150 (MQ=255)
                                      ccATGAAACAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGTCCAGTCCCAGTTGCTGGCGAATATTTTACACGATAGCATCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCCCCAGGGATCAGATGCgagaga                                                                                                               >  2:1773504/1‑148 (MQ=255)
                                       cATGAAACAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGTCCAGTCCCAGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGataa                                                                                                              >  1:924617/1‑150 (MQ=255)
                                       cATGAAACAGGTCACTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGTCCAGTCCCAGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGataa                                                                                                              <  1:1707097/150‑1 (MQ=255)
                                        aTGAAACAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGTCCAGTCCCAGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGataat                                                                                                             >  1:1125830/1‑150 (MQ=255)
                                         tGAAACAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGTCCAGTCCCAGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGataata                                                                                                            <  2:91600/150‑1 (MQ=255)
                                          gAAACAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGTCCAGTCCCAGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGAtaataa                                                                                                           >  2:1471363/1‑150 (MQ=255)
                                            aaCAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGTCCAGTCCCAGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGCTAATAAAc                                                                                                         >  1:1312665/1‑150 (MQ=255)
                                             aCAGGTCGCTGACGTAGGGGTAAAACTGTACATACAGTGGCTGGTCCAGTCCCAGTGGCTGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCTGATGCGAGATAATAAACg                                                                                                        <  2:231338/150‑1 (MQ=255)
                                                  tCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGTCCAGTCCCAGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTAATc                                                                                                   <  1:18102/150‑1 (MQ=255)
                                                  tCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGTCCAGTCCCAGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTAATc                                                                                                   <  1:10852/150‑1 (MQ=255)
                                                   cGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGTCCAGTCCCTGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGTCACCCGCCAGGAACCGCGCCGGATCACCAGGGATCAGATGCGTGATAATAAACGTAATCa                                                                                                  >  1:794886/1‑150 (MQ=255)
                                                   cGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGTCCAGTCCCAGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTAATCa                                                                                                  <  2:1245330/150‑1 (MQ=255)
                                                    gCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGTCCAGTCCCAGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTAATCAc                                                                                                 >  1:536147/1‑150 (MQ=255)
                                                               ggTAAAACTGTACATACAGTGGCTGGTCCAGTCCCAGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGGCGTGGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTAATCACGCAGAcgcccg                                                                                      <  1:1162841/150‑1 (MQ=255)
                                                                    aaCTGTACATACAGTGGCTGGTCCAGTCCCAGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATATTAAACGTATTCACGCAGACGCCCGccacc                                                                                 >  1:923908/1‑150 (MQ=255)
                                                                           catagaGTGGCTGGTCGAGTCCCAGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATTATAAACGTAATCACGCAGTCGCCCGCCACCTCGAGCa                                                                          <  2:1309880/150‑1 (MQ=255)
                                                                           catacaGTGGCTGGTCCAGTCCCAGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTAATCACGCAGACGCCCGCCACCACGAGCa                                                                          >  1:528395/1‑150 (MQ=255)
                                                                                aGTGGCTGGTCCAGTCCTAGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTAATCACGCAGACGCCCGCCACCACGAGCACCAGc                                                                     >  1:620356/1‑150 (MQ=255)
                                                                                aGTGGCTGGTCCAGTCCCAGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGGCGCAGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTAATCACGCAGACGCCCGCCACCACGAGCACCAGc                                                                     <  1:573040/150‑1 (MQ=255)
                                                                                     cTGGTCCAGTCCCAGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTAATCACGCAGACGCCCGCCACCACGAGCCCCAGcccccc                                                                >  2:1823292/1‑149 (MQ=255)
                                                                                     cTGGTCCAGTCCCAGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTAATCACGCAGACGCCCGCCACCACGAGCACCAGCCCCca                                                                <  1:1391144/150‑1 (MQ=255)
                                                                                     cTGGTCCAGTCCCAGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTAATCACGCAGACGCCCGCCACCACGAGCACCAGCCCCca                                                                >  2:1220170/1‑150 (MQ=255)
                                                                                           cAGTCCCAGTTGCTGGCGAATATTTTCCACGATAGCAGCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTAATCACGCAGACGCCCGCCACCACGAGCACCAGCCCCCAGCAGCg                                                          <  2:726724/150‑1 (MQ=255)
                                                                                            aGTCCCAGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTATTCACGCAGACGCCCGCCACCACGAGCACCAGCCCCCAGCAGCGt                                                         <  2:89459/150‑1 (MQ=255)
                                                                                            aGTCCCAGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTAATCACGCAGACGCCCGCCACCACGAGCACCAGCCCCCAGCAGCGt                                                         >  1:795545/1‑150 (MQ=255)
                                                                                              tCCCAGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTAATCACGCAGACGCCCGCCACCACGAGCACCAGCCCCCAGCAGCGTTg                                                       <  2:1281600/150‑1 (MQ=255)
                                                                                               cccAGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTAATCACGCAGACGCCCGCCACCACGAGCACCAGCCCCCAGCAGCGTTgg                                                      >  1:1375307/1‑150 (MQ=255)
                                                                                               cccAGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTAATCACGCAGACGCCCGCCACCACGAGCACCAGCCCCCAGCAGCGTTgg                                                      <  2:1331979/150‑1 (MQ=255)
                                                                                                  aGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTAATCACGCAGACGCCCGCCACCACGAGCACCAGCCCCCAGCAGCGTTGGCGt                                                   >  1:586057/1‑150 (MQ=255)
                                                                                                   gTTGCTGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTAATCACGCAGACGCCCGCCACCACGAGCACCAGCCCCCAGCAGCGTTGGCGTa                                                  >  1:832181/1‑150 (MQ=255)
                                                                                                       cTGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTAATCACGCAGACGCCCGCCACCACGAGCACCAGCCCCCAGCAGCGTTGGCGTAAAAt                                              >  1:558274/1‑150 (MQ=255)
                                                                                                         ggCGAATATTTTCCACGATAGCATCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTAATCACGCAGACGCCCGCCACCACGCGCACCAGCCCCCAGCAGCGTTGGCGTAAAATAc                                            >  1:1728649/1‑150 (MQ=255)
                                                                                                         ggCGAATATTTTCCACGATAGCATCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTAATCACGCAGACGCCCGCCACCACGAGCACCAGCCCCCAGCAGCGTTGGCGTAAAATAc                                            <  2:40114/150‑1 (MQ=255)
                                                                                                           cGAATATTTTCCACGATAGCATCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTAATCACGCAGACGCCCGCCACCACGAGCACCAGCCCCCAGCAGCGTTGGCGTAAAATACTc                                          >  2:632216/1‑150 (MQ=255)
                                                                                                           cGAATATTTTCCACGATAGCATCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTAATCACGCAGACGCCCGCCACCACGAGCACCAGCCCCCAGCAGCGTTGGCGTAAAATACTc                                          <  2:438317/150‑1 (MQ=255)
                                                                                                             aaTATTTTCCAAGATAGCATCGCTGGCGCGGTCAACCGCCACTAACCGCGCCGGATCACCAGGGAACAGATGCGAGATAATAAACGTAATCACGCAGACGCCCGCCACCACGAGACCCAGCCCCCAGCAGCGTTGGCGTAAAATACTCCa                                        >  2:263196/1‑150 (MQ=255)
                                                                                                                atTTTCCACGATAGCATCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTAATCACGCAGACGCCCGCCACCACGAGCACCAGCCCCCAGCAGCGTTGGCGTAAAATACTCCAGaa                                     <  1:395478/150‑1 (MQ=255)
                                                                                                                 ttttCCACGATAGCATCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTAATCACGCAGACGCCCGCCACCACGAGCACCAGCCCCCAGCAGCGTTGGCGTAAAATACTCCAGatg                                    >  1:1143534/1‑148 (MQ=255)
                                                                                                                   ttCCACGATAGCATCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTAATCACGCAGACGCCCGCCACCACGAGCCCCAGCCCCCAGCAGCGTTGGCGTAAAATACTCCAGAAGGt                                  >  2:529715/1‑150 (MQ=255)
                                                                                                                    tCCACGATAGCATCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTAATCACGCAGACGCCCGCCACCACGAGCACCAGCCCCCAGCAGCGTTGGCGTAAAATACTCCAGAAGGTc                                 >  1:892169/1‑150 (MQ=255)
                                                                                                                       aCGATAGCATCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTAATCACGCAGACGCCCGCCACCACGAGCACCAGCCCCCAGCAGCGTTGGCGTAAAATACTCCAGAAGGTCACg                              <  1:638541/150‑1 (MQ=255)
                                                                                                                        cGATAGCATCGCTGGCGCGTTCACCCGCCAGTAACCGAGCCGGATCAGCAGGGATCAGATGCGAGATAATAAACGTAATCACGCAGACGCCCGCCACCACGAGCACCAGCCCCCAGCAGCGTTGGCGTAAAATACTCCAGAAGGTCACGt                             <  2:192527/150‑1 (MQ=255)
                                                                                                                            aGCATCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTAATCACGCAGACGCCCGCCACCACGAGCACCAGCCCCCAGCAGCGTTGGCGTAAAATACTCCAGAAGGTCACGTTAtt                         >  2:1516297/1‑150 (MQ=255)
                                                                                                                            aGCATCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTAATCACGCAGACGCCCGCCACCACGAGCACCAGCCCCCAGCAGCGTTGGCGTAAAATACTCCAGAAGGTCACGTTAtt                         <  2:727630/150‑1 (MQ=255)
                                                                                                                             gCATCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTAATCACGCAGACGCCCGCCACCACGAGCACCAGCCCCCAGCAGCGTTGGCGTAAAATACTCCAGAAGGTCACGTTAttt                        <  2:1863767/150‑1 (MQ=255)
                                                                                                                             gCATCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTAATCACGCAGACGCCCGCCACCACGAGCACCAGCCCCCAGCAGCGTTGGCGTAAAATACTCCAGAAGGTCACGTTAttt                        <  1:467015/150‑1 (MQ=255)
                                                                                                                               aTCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTATTCACGCAGACGCCCGCCACCCCGCGCCCCAGCCCCCAGCAGCGTTGGCGTAAAATACTCCAGAAGGTCACGGTATTTAc                      >  1:1245839/1‑150 (MQ=255)
                                                                                                                                   cTGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTAATCACGCAGACGCCCGCCACCACGAGCACCAGCCCCCAGCAGCGTTGGCGTAAAATACTCCAGAAGGTCACGTTATTTACTCAt                  >  2:1698740/1‑150 (MQ=255)
                                                                                                                                   cTGGCGCGGTCACCCGCCAGTAACCGCGCAGGATCACCAGCGATCAGATGCGAGATAATAAACGTAATCACGCAGACGCCCGCCACCACGATCACCAGCCCCCAGCAGCGTTGGCGTAAAATACTCCAGAAGGTCACGTTATTTACTCAt                  <  1:1637573/150‑1 (MQ=255)
                                                                                                                                    cggCGCGGTTACTGGCGGGTGACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTAATCACGCAGACGCCCGCCACCACGAGCACCAGCCCCCAGCAGCGTTGGCGTAAAATACTCCAGAAGGTCACGTTATTTACTCATg                 >  1:416671/2‑150 (MQ=255)
                                                                                                                                     ggCGCGGTTACTGGCGGGTGACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTAATCACGCAGACGCCCGCCACCACGAGCACCAGCCCCCAGCAGCGTTGGCGGAAAATACTCCAGAAGGTCACGTTATTTACTCATgg                >  2:502155/1‑150 (MQ=255)
                                                                                                                                     ggCGCGGTCACCCGACAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTAATCACGCAGACGCCCGCCCCCACGAGCACCAGCCCCCAGCAGCGTTGGCGTAAAATACTCCAGAAGGTCACGTTATTTACTCATgg                >  2:1732904/1‑150 (MQ=255)
                                                                                                                                      gcgcGGTTACTGGCGGGTGACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTAATCACGCAGACGCCCGCCACCACGAGCACCAGCCCCCAGCAGCGTTGGCGTAAAATACTCCAGAAGGTCACGTTATTTACTCATGGt               >  1:1556186/1‑150 (MQ=255)
                                                                                                                                      gcgcGGTCACCCGCCAGTAACCGCGCCGGATCACCCGGGATCAGATGCGAGATAATAAACGTAATCACGCAGACGCCCGCCACCACGAGCACCAGCCCCCAGCAGCGTTGGCGTAAAATACTCCAGAAGGTCACGTTATTTACTCATGGt               <  1:1516297/150‑1 (MQ=255)
                                                                                                                                      gcgcGGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTAATCACGCAGACGCCCGCCACCACGAGCACCAGCCCCCAGCAGCGTTGGCGTAAAATACTCCAGAAGGTCACGTTATTTACTCATGGt               <  1:1471363/150‑1 (MQ=255)
                                                                                                                                        gcgGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTAATCACGCAGACGCCCGCCACCACGAGCACCAGCCCCCCGCAGCGTTGGCGTAAAATACTCCAGAAGGTCACGTTATTTACTCATGGtat             >  2:738088/1‑150 (MQ=255)
                                                                                                                                         cgGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTAATCACGCAGACGCCCGCCACCACGAGCACCAGCCCCCAGCAGCGTTGGCGTAAAATACTCCAGAAGGTCACGTTATTTACTCATGGtatt            >  2:856077/1‑150 (MQ=255)
                                                                                                                                         cgGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTAATCACGCAGACGCCCGCCACCACGAGCACCAGCCCCCAGCAGCGTTGGCGTAAAATACTCCACAAGGTCACGTTATTTACTCATGGtatt            <  2:1583210/150‑1 (MQ=255)
                                                                                                                                           ggcACCCGCCAGTAACGGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTAATCACGCAGACGCCCGCCACCACGAGCACCAGCCCCCATCAGCGTTGGCGTAAAATACTCCAGAAGGTCACGTTATTTACTCATGGtattga          <  2:1362342/148‑1 (MQ=255)
                                                                                                                                           ggcACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTAATCACGCAGACGCCCGCCACCACGAGCACCAGCCCCCAGCAGCGTTGGCGTAAAATACTCCAGAAGGTCACGTTATTTACTCATGGtattga          <  1:1431110/148‑1 (MQ=255)
                                                                                                                                             cACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTAATCACGCAGACGCCCGCCCCCACGAGCCCCAGCCCCCAGCAGCGTTGGCGTAAAATCCTCCAGAAGGTCACGTTATTTACTCATGGtattgata        >  1:1273036/1‑150 (MQ=255)
                                                                                                                                               cccGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTAATCACGCAGACGCCCGCCACCACGAGCACCAGCCCCCAGCAGCGTTGGCGTAAAATACTCCAGAAGGTCACGTTATTTACTCATGGtattgatatt      >  1:435007/1‑150 (MQ=255)
                                                                                                                                               cccGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTAATCACGCAGACGCCCGCCACCACGAGCACCAGCCCCCAGCAGCGTTGGCGGAAAATACTCCAGAAGGTCACGTTATTTACTCATGGtattgatatt      >  1:276822/1‑150 (MQ=255)
                                                                                                                                                ccGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTAATCACGCAGACGCCCGCCACCACGAGCACCAGCCCCCAGCAGCGTTGGCGTAAAATACTCCAGAAGGTCACGTTATTTACTCATGGtattgatattg     <  2:644262/150‑1 (MQ=255)
                                                                                                                                                ccGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTAATCACGCAGACGCCCGCCACCACGAGCACCAGCCCCCAGCAGCGTTGGCGTAAAATACTCCAGAAGGTCACGTTATTTACTCATGGtattgatattg     <  2:606257/150‑1 (MQ=255)
                                                                                                                                                ccGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTAATCACGCAGACGCCCGCCACCACGAGCACCAGCCCCCAGCAGCGTTGGCGTAAAATACTCCAGAAGGTCACGTTATTTACTCATGGtattgatattg     <  2:1506062/150‑1 (MQ=255)
                                                                                                                                                 cGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTAATCACGCAGACGCCCGCCACCACGAGCACCAGCCCCCAGCAGCGTTGGCGTAAAATACTCCAGAAGGTCACGTTATTTACTCATGGtattgatattga    >  1:977303/1‑150 (MQ=255)
                                                                                                                                                  gCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTAATCACGCAGACGCCCGCCACCACGAGCACCAGCCCCCAGCAGCGTTGGCGTAAAATACTCCAGAAGGTCACGTTATTTACTCATGGTATTGATATTGaa   <  2:536147/150‑1 (MQ=255)
                                                                                                                                                   ccAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTAATCACGCAGACGCCCGCCACCACGAGCACCAGCCCCCAGCAGCGTTGGCGTAAAATACTCCAGAAGGTCACGTTATTTACTCATGGTATTGATATTGAAg  <  1:422416/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
GCACCGGACGCCCGGTACGAATGGATGTTCCCAGGTCACCATGAAACAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGTCCAGTCCCAGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTAATCACGCAGACGCCCGCCACCACGAGCACCAGCCCCCAGCAGCGTTGGCGTAAAATACTCCAGAAGGTCACGTTATTTACTCATGGTATTGATATTGAAG  >  NC_000913/1560662‑1560958

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