Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,560,807 C→G 8.2% A41A (GCG→GCC ddpB ← putative D,D‑dipeptide ABC transporter membrane subunit DdpB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,560,8070CG 8.2% 309.8 / 14.9 97A41A (GCG→GCCddpBputative D,D‑dipeptide ABC transporter membrane subunit DdpB
Reads supporting (aligned to +/- strand):  ref base C (48/41);  new base G (3/5);  total (51/46)
Fisher's exact test for biased strand distribution p-value = 4.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.00e-01

Only 100 of 104 total aligned reads displayed.

CCAGCACCGGACGCCCGGTACGAATGGATGTTCCCAGGTCACCATGAAACAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGTCCAGTCCCAGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTAATCACGCAGACGCCCGCCACCACGAGCACCAGCCCCCAGCAGCGTTGGCGTAAAATACTCCAGAAGGTCACGTTATTTACTCATGGTATTGATATTGA  >  NC_000913/1560659‑1560956
                                                                                                                                                    |                                                                                                                                                     
ccAGCACCGGACGCCCGGTACGAATGGATGTTCCCAGGTCACCATGAAACAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGTCCAGTCCCAGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGTCACCCg                                                                                                                                                      >  1:921845/1‑150 (MQ=255)
ccAGCACCGGACGCCCGGTACGAATGGATGTTCCCAGGTCACCATGAAACAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGTCCAGTCCCAGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGTCACCCg                                                                                                                                                      >  2:1346290/1‑150 (MQ=255)
  aGCACCGGACGCCCGGTACGAATGGATGTTCCCAGGTCACCATGAAACAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGTCCAGTCCCAGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGTCACCCGcc                                                                                                                                                    <  1:1346290/150‑1 (MQ=255)
    cACCGGACGCCCGGTACGAATGGATGTTCCCAGGTCACCATGAAACAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGTCCAGTCCCAGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGTCACCCGCCAg                                                                                                                                                  >  1:507075/1‑150 (MQ=255)
    cACCGGACGCCCGGTACGAATGGATGTTCCCAGGTAACCATGAAACAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGTCCAGTCCCAGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGTCACCCGCCAg                                                                                                                                                  <  2:1442942/150‑1 (MQ=255)
       cGGACGCCCGGTACGAATGGATGTTCCCAGGTCACCATGAAACAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGTCCAGTCCCAGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGTCACCCGCCAGTaa                                                                                                                                               >  2:546509/1‑150 (MQ=255)
        ggACGCCCGGTACGAATGGATGTTCCCAGGTCACCATGAAACAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGTCCAGTCCCAGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGTCACCCGCCAGTAAc                                                                                                                                              >  2:67165/1‑150 (MQ=255)
         gACGCCCGGTACGAATGGATGTTCCCAGGTCACCATGAAACAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGTCCAGTCCCAGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGTCACCCGCCAGTAAcc                                                                                                                                             <  1:1311885/150‑1 (MQ=255)
               cGGTACGAATGGATGTTCCCAGGTCACCATGAAACAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGTCCAGTCCCAGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGTTACTGGCGGGTGACCGCGcca                                                                                                                                       <  2:453871/150‑2 (MQ=255)
               cGGTACGAATGGATGTTCCCAGGTCACCATGAAACAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGTCCAGTCCCAGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGTTACTGGCGGGTGACCGCGcca                                                                                                                                       <  1:639901/150‑2 (MQ=255)
                ggTACGAATGGATGTTCCCAGGTCACCATGAAACAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGTCCAGTCCCAGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGTTACTGGCGGGTGACCGCGccag                                                                                                                                      <  1:327831/150‑3 (MQ=255)
                ggTACGAATGGATGTTCCCAGGTCACCATGAAACAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGTCCAGTCCCAGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCgg                                                                                                                                      <  1:1379685/150‑1 (MQ=255)
                 gTACGAATGGATGTTCCCCTGTCACCATGAAACAGGTCGCTGACGTAGCGGTAAAACTTTACATACAGTGGCTGGGCCAGTCCCAGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGTTACTGGCGGGTGACCGCGccagc                                                                                                                                     <  2:467647/150‑4 (MQ=255)
                  tACGAATGGATGTTCCCAGGTCACCATGAAACAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGTCCAGTCCCAGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGTTACTGGCGGGTGACCGCGccagng                                                                                                                                    <  1:320752/150‑5 (MQ=255)
                     gAATGGATGTTCCCAGGTCACCATGAAACAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGTCCAGTCCCAGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCAc                                                                                                                                 <  2:1414295/150‑1 (MQ=255)
                             gTTCCCAGGTCACCATGAAACAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGTCCAGTCCCAGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATc                                                                                                                         <  2:924584/150‑1 (MQ=255)
                              ttCCCAGGTCACCATGAAACAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGTCCAGTCCCAGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCa                                                                                                                        <  1:526691/150‑1 (MQ=255)
                               tCCCAGGTCACCATGAAACAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGTCCAGTCCCAGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAg                                                                                                                       <  1:981656/150‑1 (MQ=255)
                               tCCCAGGTCACCATGAAACAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGTCCAGTCCCAGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAg                                                                                                                       >  2:640887/1‑150 (MQ=255)
                               tCCCAGGTCACCATGAAACAGGTCGCTGACGTAGCGGTAAAAATGTACATACAGTGGCTGGTCCAGTCCCAGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAg                                                                                                                       >  2:646278/1‑150 (MQ=255)
                                cccAGGTCACCATGAAACAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGTCCAGTCCCAGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGa                                                                                                                      >  1:1393677/1‑150 (MQ=255)
                                cccAGGTCACCATGAAACAGGTCGCTGACGTAGCGGTAAAACTCTACATACAGTGGCTGGTCCAGTCCCAGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGa                                                                                                                      >  1:918345/1‑150 (MQ=255)
                                  cAGGTCACCATGAAACAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGTCCAGTCCCAGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATg                                                                                                                    <  2:819321/150‑1 (MQ=255)
                                  cAGGTCACCATGAAACAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGTCCAGTCCCAGATGCTGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATg                                                                                                                    >  1:1441582/1‑150 (MQ=255)
                                   aGGTCACCATGAAACAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGTCCAGTCCCAGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGc                                                                                                                   <  2:698807/150‑1 (MQ=255)
                                       cACCATGAAACAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGTCCAGTCCCAGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCgaga                                                                                                               >  2:1273463/1‑150 (MQ=255)
                                        aCCATGAAACAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGTCCAGTCCCAGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGat                                                                                                              <  1:559723/150‑1 (MQ=255)
                                        aCCATGAAACAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGTCCAGTCCCAGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGat                                                                                                              <  1:1597432/150‑1 (MQ=255)
                                        aCCATGAAACAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGTCCAGTCCCAGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGat                                                                                                              <  2:1056007/150‑1 (MQ=255)
                                          cATGAAACAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGTCCAGTCCCAGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGataa                                                                                                            <  1:1034965/150‑1 (MQ=255)
                                          cATGAAACAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGTCCAGTCCCAGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGataa                                                                                                            >  1:1374831/1‑150 (MQ=255)
                                                 cAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGTCCAGTCCCAGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGt                                                                                                     >  1:940485/1‑150 (MQ=255)
                                                 cAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGTCCAGTCCCAGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGt                                                                                                     <  1:67165/150‑1 (MQ=255)
                                                   ggTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGTCCAGTCCCAGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCGGATc                                                                                                                                   <  1:315593/118‑1 (MQ=255)
                                                   ggTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGTCCAGTCCCAGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTaa                                                                                                   >  1:371655/1‑150 (MQ=255)
                                                    gTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGTCCAGTCCCAGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTAAt                                                                                                  <  1:913742/150‑1 (MQ=255)
                                                       gCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGTCCAGTCCCAGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTAATCAc                                                                                               >  2:1528124/1‑150 (MQ=255)
                                                          gACGTAGCGGTAAAACTGTACATACAGTGGCTGGTCCAGTCCCAGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTAATCACGCa                                                                                            >  1:517886/1‑150 (MQ=255)
                                                          gACGTAGCGGTAAAACTGTACATACAGTGGCTGGTCCAGTCCCAGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTAATCACGCa                                                                                            <  1:607655/150‑1 (MQ=255)
                                                                gCGGTAAAACTGTACATACAGTGGCTGGTCCAGTCCCAGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTAATCACGCAGAcgcc                                                                                      >  1:1476146/1‑150 (MQ=255)
                                                                  ggTAAAACTGTACATACAGTGGCTGGTCCAGTCCCAGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTAATCACGCAGAcgcccg                                                                                    >  1:352110/1‑150 (MQ=255)
                                                                    tAAAACTGTACATACAGTGGCTGGTCCAGTCCCAGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTAATCACGCAGAcgcccgca                                                                                  <  2:940485/150‑2 (MQ=255)
                                                                      aaaCTGTACATACAGTGGCTGGTCCAGTCCCAGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTAATCACGCAGACGCCCGccac                                                                                >  2:61214/1‑150 (MQ=255)
                                                                       aaCTGTACATACAGTGGCTGGTCCAGTCCCAGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTAATCACGCAGACGCCCGccacc                                                                               >  1:703155/1‑150 (MQ=255)
                                                                       aaCTGTACATACAGTGGCTGGTCCAGTCCCAGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTAATCACGCAGACGCCCGccacc                                                                               >  1:798246/1‑150 (MQ=255)
                                                                            tacatacaGTGGCTGGTCCAGTCCCAGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTAATCACGCAGACGCCCGCCACCACGAg                                                                          <  2:507075/150‑1 (MQ=255)
                                                                              catacaGTGGCTGGTCCAGTCCCAGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTAATCACGCAGACGCCCGCCACCACGAGcc                                                                        >  1:483627/1‑149 (MQ=255)
                                                                              catacaGTGGCTGGTCCAGTCCCAGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTAATCACGCAGACGCCCGCCACCACGAGCa                                                                        <  2:250879/150‑1 (MQ=255)
                                                                                tacaGTGGCTGGTCCAGTCCCAGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTAATCACGCAGACGCCCGCCACCACGAGCAcc                                                                      >  1:1526783/1‑150 (MQ=255)
                                                                                    gTGGCTGGTCCAGTCCCAGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTAATCACGCAGACGCCCGCCACCACGAGCACCAGcc                                                                  <  2:918345/150‑1 (MQ=255)
                                                                                     tGGCTGGTCCAGTCCCAGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTAATCACGCAGACGCCCGCCACCACGAGCACCAGccc                                                                 >  2:411071/1‑150 (MQ=255)
                                                                                     tGGCTGGTCCAGTCCCAGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTAATCACGCAGACGCCCGCCACCACGAGCACCAGccc                                                                 <  1:1133641/150‑1 (MQ=255)
                                                                                      ggCTGGTCCAGTCCCAGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTAATCACGCAGACGCCCGCCACCACGAGCACCAGcccc                                                                >  1:1386828/1‑150 (MQ=255)
                                                                                          gggccAGTCCCAGTTGCTGGCGTATATTTTCCACGATTGCATCGCTGGCGCGGTCACCCGCCAGTAACCGCGGCGGATCACTATGGATCAGATGCGAGATTATAAACGTAATCTAGCAGAGGCTCGCCACCATGAGCAGCAGACGCcagc                                                            <  2:1441582/147‑1 (MQ=255)
                                                                                             ccAGTCCCAGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTAATCACGCAGACGCCCGCCACCACGAGCACCAGCCCCCagcagc                                                         <  1:156998/150‑1 (MQ=255)
                                                                                             ccAGTCCCAGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTAATCACGCAGACGCCCGCCACCACGAGCACCAGCCCCCagcagc                                                         >  1:645995/1‑150 (MQ=255)
                                                                                                   ccAGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTAATCACGCAGACGCCCGCCACCACGAGCACCCGCCCCCAGCAGCGTTGGc                                                   >  1:1042167/1‑150 (MQ=255)
                                                                                                   ccAGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTAATCACGCAGACGCCCGCCACCACGAGCACCAGCCCCCAGCAGCGTTGGc                                                   <  2:376606/150‑1 (MQ=255)
                                                                                                    cAGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTAATCACGCAGACGCCCGCCACCACGAGCACCAGCCCCCAGCAGCGTTGGCg                                                  <  1:448197/150‑1 (MQ=255)
                                                                                                        tGCTGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTAATCACGCAGACGCCCGCCACCACGAGCACCAGCCCCCAGCAGCGTTGGCGTaaa                                              >  1:589377/1‑150 (MQ=255)
                                                                                                          cTGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTAATCACGCAGACGCCCGCCACCACGAGCACCAGCCCCCAGCAGCGTTGGCGTAAAAt                                            >  2:1619344/1‑150 (MQ=255)
                                                                                                            ggCGAATATTTTCCACGATAGCATCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTAATCACGCAGACGCCCGCCACCACGAGCACCAGCCCCCAGCAGCGTTGGCGTAAAATAc                                          >  2:506040/1‑150 (MQ=255)
                                                                                                             gCGAATATTTTCCACGATAGCATCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTAATCACGCAGACGCCCGCCACCAGGAGCACCAGCCCCCAGCAGCGTTGGCGTAAAATACt                                         >  1:17975/1‑150 (MQ=255)
                                                                                                             gCGAATATTTTCCACGATAGCATCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTAATCACGCAGACGCCCGCCACCACGAGCACCAGCCCCCAGCAGCGTTGGCGTAAAATACt                                         <  2:1073018/150‑1 (MQ=255)
                                                                                                              cGAATATTTTCCACGATAGCATCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTAATCACGCAGACGCCCGCCACCACGAGCACCAGCCCCCAGCAGCGTTGGCGTAAAATACTc                                        >  1:1484444/1‑150 (MQ=255)
                                                                                                              cGAATATTTTCCACGATAGCATCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTAATCACGCAGACGCCCGCCACCACGAGCACCAGCCCCCAGCAGCGTTGGCGTAAAATACTc                                        >  2:1023351/1‑150 (MQ=255)
                                                                                                                    ttttCCACGATAGCATCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTAATCACGCAGACGCCCGCCACCACGAGCACCAGCCCCCAGCAGCGTTGGc                                                   <  1:1194119/133‑1 (MQ=255)
                                                                                                                    ttttCCACGATAGCATCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTAATCACGCAGACGCCCGCCACCACGAGCACCAGCCCCCAGCAGCGTTGGc                                                   >  2:1194119/1‑133 (MQ=255)
                                                                                                                        ccACGATAGCATCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTAATCACGCAGACGCCCGCCACCACGAGCACCAGCCCCCAGCAGCGTTGGCGTAAAATACTCCAGAAGGTCa                              <  2:798246/150‑1 (MQ=255)
                                                                                                                        ccACGATAGCATCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTAATCACGCAGACGCCCGCCACCACGAGCACCAGCCCCCAGCAGCGTTGGCGTAAAATACTCCAGAAGGTCa                              <  2:703155/150‑1 (MQ=255)
                                                                                                                         cACGATAGCATCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTAATCACGCAGACGCCCGCCACCACGAGCACCAGCCCCCAGCAGCGTTGGCGTAAAATACTCCAGAAGGTCAc                             >  2:1438503/1‑150 (MQ=255)
                                                                                                                                gCATCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTAATCACGCAGACGCCCGCCACCACGAGCACCAGCCCCCAGCAGCGTTGGCGTAAAATACTCCAGAAGGTCACGTTAttt                      >  2:675430/1‑150 (MQ=255)
                                                                                                                                gCATCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTAATCACGCAGACGCCCGCCACCACGAGCACCAGCCCCCAGCAGCGTTGGCGTAAAATACTCCAGAAGGTCACGTTAtt                       <  2:1538864/149‑1 (MQ=255)
                                                                                                                                gCATCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTAATCACGCAGACGCCCGCCACCACGAGCACCAGCCCCCAGCAGCGTTGGCGTAAAATACTCCAGAAGGTCACGTTAtt                       >  1:1538864/1‑149 (MQ=255)
                                                                                                                                 cATCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTAATCACGCAGACGCCCGCCACCACGAGCACCAGCCCCCAGCAGCGTTGGCGTAAAATACTCCAGAAGGTCACGTTATTTa                     <  2:921845/150‑1 (MQ=255)
                                                                                                                                   tCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTAATCACGCAGACGCCCGCCACCACGAGCACCAGCCCCCAGCAGCGTTGGCGTAAAATACTCCa                                      <  2:314024/131‑1 (MQ=255)
                                                                                                                                   tCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTAATCACGCAGACGCCCGCCACCACGAGCACCAGCCCCCAGCAGCGTTGGCGTAAAATACTCCa                                      >  1:314024/1‑131 (MQ=255)
                                                                                                                                   tCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTAATCACGCAGACGCCCGCCACCACGAGCACCAGCCCCCAGCAGCGTTGGCGTAAAATACTCCAGAAGGTCACGTTATTTACt                   >  1:339571/1‑150 (MQ=255)
                                                                                                                                     tccggCGCGGTTACTGGCGGGTGACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTAATCACGCAGACGCCCGCCACCACGAGCACCAGCCCCCAGCAGCGTTGGCGTAAAATACTCCAGAAGGTCACGTTATTTACTCa                 >  2:1376562/4‑150 (MQ=255)
                                                                                                                                     gCTGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTAATCACGCAGACGCCCGCCACCACGAGCACCAGCCCCCAGCAGCGTTGGCGTAAAATACTCCAGAAGGTCACGTTATTTACTCa                 <  2:498348/150‑1 (MQ=255)
                                                                                                                                     gCTGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTAATCACGCAGACGCCCGCCACCACGAGCACCAGCCCCCAGCAGCGTTGGCGTAAAATACTCCAGAAGGTCACGTTATTTACTCa                 >  2:1528773/1‑150 (MQ=255)
                                                                                                                                      cTGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTAATCACGCAGACGCCCGCCACCACGAGCACCAGCCCCCAGCAGCGTTGGCGTAAAATACTCCAGAAGGTCACGTTATTTACTCAt                >  2:1097616/1‑150 (MQ=255)
                                                                                                                                        ggCGCGGTTACTGGCGGGTGACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTAATCACGCAGACGCCCGCCACCACGAGCACCAGCCCCCAGCAGCGTTGGCGTAAAATACTCCAGAAGGTCACGTTATTTACTCATgg              >  1:596275/1‑150 (MQ=255)
                                                                                                                                        ggCGCGGTTACTGGCGGGTGACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTAATCACGCAGACGCCCGCCACCACGAGCACCAGCCCCCAGCAGCGTTGGCGTAAAATACTCCAGAAGGTCACGTTATTTACTCATgg              >  2:801408/1‑150 (MQ=255)
                                                                                                                                            cgGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTAATCACGCAGACGCCCGCCACCACGAGCACCAGCCCCCAGCAGCGTTGGCGTAAAATACTCCAGAAGGTCACGTTATTTACTCATGGtatt          <  1:675430/150‑1 (MQ=255)
                                                                                                                                            cgGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTAATCACGCAGACGCCCGCCACCACGAGCACCAGCCCCCAGCAGCGTTGGCGTAAAATACTCCAGAAGGTCACGTTATTTACTCATGGtatt          >  2:1303217/1‑150 (MQ=255)
                                                                                                                                             ggTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTAATCACGCAGACGCCCGCCACCACGAGCACCAGCCCCCAGCAGCGTTGGCGTAAAATACTCCGGAAGGTCACGTTATTTAGTCATGGaattt         >  2:1228705/1‑149 (MQ=255)
                                                                                                                                             ggTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTAATCACGCAGACGCCCGCCACCACGAGCACCAGCCCCCAGCAGCGTTGGCGTAAAATACTCCAGAAGGTCACGTTATTTACTCATGGtattg         <  1:1303217/150‑1 (MQ=255)
                                                                                                                                             ggTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTAATCACGCAGACGCCCGCCACCACGAGCACCAGCCCCCAGCAGCGTTGGCGTAAAATACTCCAGAAGGTCACGTTATTTACTCATGGtattg         <  1:646278/150‑1 (MQ=255)
                                                                                                                                             ggTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTAATCACGCAGACGCCCGCCACCACGAGCACCAGCCCCCAGCAGCGTTGGCGTAAAATACTCCAGAAGGTCACGTTATTTACTCATGGtattg         <  1:640887/150‑1 (MQ=255)
                                                                                                                                             ggTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTAATCACGCAGACGCCCGCCACCACGAGCACCAGCCCCCAGCAGCGTTGGCGTAAAATACTCCAGAAGGTCACGTTATTTACTCATGGtattg         >  2:1404694/1‑150 (MQ=255)
                                                                                                                                               tCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTAATCACGCAGACGCCCGCCACCACGAGCACCAGCCCCCAGCAGCGTTGGCGTAAAATACTCCAGAAGGTCACGTTATTTACTCATGGtattgat       <  2:324660/150‑1 (MQ=255)
                                                                                                                                                cACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTAATCACGCAGACGCCCGCCACCACGAGCACCAGCCCCCAGCAGCGTTGGCGTAAAATACTCCAGAAGGTCACGTTATTTACTCATGGtattgata      >  2:287494/1‑150 (MQ=255)
                                                                                                                                                  cccGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTAATCACGCAGACGCCTGCCACCACGAGCACCAGCCCCCAGCAGCGTTGGCGTAAAATACTCCAGAAGGTCACGTTATTTACTCATGGtattgatatt    >  1:1115881/1‑150 (MQ=255)
                                                                                                                                                  cccGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTAATCACGCAGACGCCCGCCACCCCGAGCACCAGCCCCCCGCAGCGTTGGCGTAAAATACTCCAGAAGGTCACGTTATTTACTCATGGtattgatatt    >  2:1460190/1‑150 (MQ=255)
                                                                                                                                                  cccGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTAATCACGCAGACGCCCGCCACCACGAGCACCAGCCCCCAGCAGCGTTGGCGTAAAATACTCCAGAAGGTCACGTTATTTACTCATGGtattgatatt    >  2:120826/1‑150 (MQ=255)
                                                                                                                                                  cccGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTAATCACGCAGACGCCCGCCACCACGAGCACCAGCCCCCAGCAGCGTTGGCGTAAAATACTCCAGAAGGTCACGTTATTTACTCATGGtattgatatt    >  1:203225/1‑150 (MQ=255)
                                                                                                                                                   ccGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTAATCACGCAGACGCCCGCCACCACGAGCACCAGCCCCCAGCAGCGTTGGCGTAAAATACTCCAGAAGGTCACGTTATTTACTCATGGtattgatattg   <  2:645995/150‑1 (MQ=255)
                                                                                                                                                    cGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTAATCACGCAGACGCCCGCCACCACGAGCACCAGCCCCCAGCAGCGTTGGCGTAAAATACTCCAGAAGGTCACGTTATTTACTCATGGtattgatattga  <  2:608041/150‑1 (MQ=255)
                                                                                                                                                    cGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTAATCACGCAGACGCCCGCCACCACGAGCACCAGCCCCCAGCAGCGTTGGCGTAAAATACTCCAGAAGGTCACGTTATTTACTCATGGtattgatattga  >  1:1604211/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CCAGCACCGGACGCCCGGTACGAATGGATGTTCCCAGGTCACCATGAAACAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGTCCAGTCCCAGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTAATCACGCAGACGCCCGCCACCACGAGCACCAGCCCCCAGCAGCGTTGGCGTAAAATACTCCAGAAGGTCACGTTATTTACTCATGGTATTGATATTGA  >  NC_000913/1560659‑1560956

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