Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 354,617 T→C 9.4% intergenic (+25/‑305) prpE → / → codB propionate‑‑CoA ligase/cytosine transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913354,6170TC 9.4% 445.8 / 31.0 149intergenic (+25/‑305)prpE/codBpropionate‑‑CoA ligase/cytosine transporter
Reads supporting (aligned to +/- strand):  ref base T (63/72);  new base C (3/11);  total (66/83)
Fisher's exact test for biased strand distribution p-value = 9.17e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.31e-01

Only 100 of 156 total aligned reads displayed.

CGCGATCCGGAAAAATGCTGCGCCGCACGATCCAGGCGATTTGCGAAGGACGCGATCCTGGAGATCTGACGACCATTGATGATCCTGCGTCGTTGGATCAGATCCGCCAGGCGATGGAAGAGTAGGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGATGCGACGCTGACGCGTTTTATCATGCCTACGGACCTGAACCGTAGGTCGGATAAGGCGCTCGCGTCGCATCCGACACCATGCTCAGATGCCTGATGCGACGCTGACGCGTCTTATCAGGCC  >  NC_000913/354468‑354766
                                                                                                                                                     |                                                                                                                                                     
cgcgATCCGGAAAAATGCTGCGCCGCACGATCCAGGCGATTTGCGAAGGACGCGATCCTGGAGATCTGACGACCATTGATGATCCTGCGTCGTTGGATCAGATCCGCCAGGCGATGGAAGAGTAGGTTATTGTCGGATGCGTCGCGCGGt                                                                                                                                                       <  2:90080/150‑1 (MQ=255)
cgcgATCCGGAAAAATGCTGCGCCGCACGATCCAGGCGATTTGCGAAGGACGCGATCCTGGAGATCTGACGACCATTGATGATCCTGCGTCGTTGGATCAGATCCGCCAGGCGATGGAAGAGTAGGTTATTGTCGGATGCGTCGCGCGGt                                                                                                                                                       <  2:410/150‑1 (MQ=255)
cgcgATCCGGAAAAATGCTGCGCCGCACGATCCAGGCGATTTGCGAAGGACGCGATCCTGGAGATCTGACGACCATTGATGATCCTGCGTCGTTGGATCAGATCCGCCAGGCGATGGAAGAGTAGGTTATTGTCGGATGCGTCGCGCGGt                                                                                                                                                       <  1:1439242/150‑1 (MQ=255)
    aTCCGGAAAAATGCTGCGCCGCACGATCCAGGCGATTTGCGAAGGACGCGATCCTGGAGATCTGACGACCATTGATGATCCTGCGTCGTTGGATCAGATCCGCCAGGCGATGGAAGAGTAGGTTATTGTCGGATGCGTCGCGCGGTGCAt                                                                                                                                                   >  2:366766/1‑150 (MQ=255)
     tCCGGAAAAATGCTGCGCCGCACGATCCAGGCGATTTGCGAAGGACGCGATCCTGGAGATCTGACGACCATTGATGATCCTGCGTCGTTGGATCAGATCCGCCAGGCGATGGAAGAGTAGGTTATTGTCGGATGCGTCGCGCGGTGCATc                                                                                                                                                  >  1:2171355/1‑150 (MQ=255)
     tCCGGAAAAATGCTGCGCCGCACGATCCAGGCGATTTGCGAAGGACGCGATCCTGGAGATCTGACGACCATTGATGATCCTGCGTCGTTGGATCAGATCCGCCAGGCGATGGAAGAGTAGGTTATTGTCGGATGCGTCGCGCGGTGCATc                                                                                                                                                  <  1:1296026/150‑1 (MQ=255)
     tCCGGAAAAATGCTGCGCCGCACGATCCAGGCGATTTGCGAAGGACGCGATCCTGGAGATCTGACGACCATTGATGATCCTGCGTCGTTGGATCAGATCCGCCAGGCGATGGAAGAGTAGGTTATTGTCGGATGCACCGCGCGACGCATc                                                                                                                                                  <  2:1728489/150‑1 (MQ=255)
     tCCGGAAAAATGCTGCGCCGCACGATCCAGGCGATTTGCGAAGGACGCGATCCTGGAGATCTGACGACCATTGATGATCCTGCGTCGTTGGATCAGATCCGCCAGGCGATGGAAGAGTAGGTTATTGTCGGATGCACCGCGCGACGCATc                                                                                                                                                  <  2:1683812/150‑1 (MQ=255)
     tCCGGAAAAATGCTGCGCCGCACGATCCAGGCGATTTGCGAAGGACGCGATCCTGGAGATCTGACGACCATTGATGATCCTGCGTCGTTGGATCAGATCCGCCAGGCGATGGAAGAGTAGGTTATTGTCGGATGCACCGCGCGACGCATc                                                                                                                                                  <  2:2649129/150‑1 (MQ=255)
     tCCGGAAAAATGCTGCGCCGCACGATCCAGGCGATTTGCGAAGGACGCGATCCTGGAGATCTGACGACCATTGATGATCCTGCGTCGTTGGATCAGATCCGCCAGGCGATGGAAGAGTAGGTTATTGTCGGATGCACCGCGCGACGCATc                                                                                                                                                  <  2:2409271/150‑1 (MQ=255)
     tCCGGAAAAATGCTGCGCCGCACGATCCAGGCGATTTGCGAAGGACGCGATCCTGGAGATCTGACGACCATTGATGATCCTGCGTCGTTGGATCAGATCCGCCAGGCGATGGAAGAGTAGGTTATTGTCGGATGCACCGCGCGACGCATc                                                                                                                                                  <  2:436560/150‑1 (MQ=255)
     tCCGGAAAAATGCTGCGCCGCACGATCCAGGCGATTTGCGAAGGACGCGATCCTGGAGATCTGACGACCATTGATGATCCTGCGTCGTTGGATCAGATCCGCCAGGCGATGGAAGAGTAGGTTATTGTCGGACGCGTCGCGCGGTGCATc                                                                                                                                                  <  1:38706/150‑1 (MQ=255)
     tCCGGAAAAATGCTGCGCCGCACGATCCAGGCGATTTGCGAAGGACGCGATCCTGGAGATCTGACGACCATTGATGATCCTGCGTCGTTGGATCAGATCCGCCAGGCGATGGAAGAGTAGGTTATTGCCGGATGCACCGCGCGACGCATc                                                                                                                                                  <  2:37845/150‑1 (MQ=255)
       cGGAAAAATGCTGCGCCGCACGATCCAGGCGATTTGCGAAGGACGCGATCCTGGAGATCTGACGACCATTGATGATCCTGCGTCGTTGGATCAGATCCGCCAGGCGATGGAAGAGTAGGTTATTGTCGGATGCGTCGCGCGGTGCATCCg                                                                                                                                                >  1:1173391/1‑150 (MQ=255)
        ggAAAAATGCTGCGCCGCACGATCCAGGCGATTTGCGAAGGACGCGATCCTGGAGATCTGACGACCATTGATGATCCTGCGTCGTTGGATCAGATCCGCCAGGCGATGGAAGAGTAGGTTATTGTCGGATGCGTCGCGCGGTGCATCCgg                                                                                                                                               <  1:1642078/150‑1 (MQ=255)
        ggAAAAATGCTGCGCCGCACGATCCAGGCGATTTGCGAAGGACGCGATCCTGGAGATCTGACGACCATTGATGATCCTGCGTCGTTGGATCAGATCCGCCAGGCGATGGAAGAGTAGGTTATTGTCGGATGCGTCGCGCGGTGCATCCgg                                                                                                                                               <  1:2273490/150‑1 (MQ=255)
        ggAAAAATGCTGCGCCGCACAATCCAGGCGATTTGCGAAGGACGCGATCCTGGAGATCTGACGACCATTGATGATCCTGCGTCGTTGGATCAGATCCGCCAGGCGATGGAAGAGTAGGTTATTGTCGGATGCGTCGCGCGGTGCATCCgg                                                                                                                                               <  1:1464402/150‑1 (MQ=255)
              aTGCTGCGCCGCACGATCCAGGCGATTTGCGAAGGACGCGATCCTGGAGATCTGACGACCATTGATGATCCTGCGTCGTTGGATCAGATCCGCCAGGCGATGGAAGAGTAGGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACtgt                                                                                                                                         <  2:2160752/150‑1 (MQ=255)
               tgctgcGCCGCACGATCCTGGCGATTGGCGAAGGACGCTATCCGGGAGATCTGACGACCATTGTTGATCCTGCGTGTTTGGATCAGATCCGCCTGGCGATGGAAGAGTAGGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACtgtg                                                                                                                                        <  1:1570736/150‑1 (MQ=255)
                    cgccgcACGATCCAGGCGATTTGCGAAGGACGCGATCCTGGAGATCTGACGACCATTGATGATCCTGCGTCGTTGGATCAGATCCGCCAGGCGATGGAAGAGTAGGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCg                                                                                                                                   >  1:32434/1‑150 (MQ=255)
                          aCGATCCAGGCGATTTGCGAAGGACGCGATCCTGGAGATCTGACGACCATTGATGATCCTGCGTCGTTGGATCAGATCCGCCAGGCGATGGAAGAGTAGGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCt                                                                                                                             <  2:1177337/150‑1 (MQ=255)
                          aCGATCCAGGCGATTTGCGAAGGACGCGATCCTGGAGATCTGACGACCATTGATGATCCTGCGTCGTTGGATCAGATCCGCCAGGCGATGGAAGAGTAGGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCt                                                                                                                             <  2:811025/150‑1 (MQ=255)
                              tCCCGGCTATTTGCGAAGGACGCGATCCTGGAGATCTGACGACCATTGATGATCCTGCGTCGTTGGATCAGATCCGCCAGGCGATGGAAGAGTTGGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGATg                                                                                                                         <  1:980231/150‑1 (MQ=255)
                                   gCGATTTGCGAAGGACGCGATCCTGGAGATCTGACGACCATTGATGATCCTGCGTCGTTGGATCAGATCCGCCAGGCGATGGAAGAGTAGGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGATGCGACg                                                                                                                    >  2:1285393/1‑150 (MQ=255)
                                   gCGATTTGCGAAGGACGCGATCCTGGAGATCTGACGACCATTGATGATCCTGCGTCGTTGGATCAGATCCGCCAGGCGATGGAAGAGTAGGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGATGCGACg                                                                                                                    >  2:1475289/1‑150 (MQ=255)
                                   gCGATTTGCGAAGGACGCGATCCTGGAGATCTGACGACCATTGATGATCCTGCGTCGTTGGATCAGATCCGCCAGGCGATGGAAGAGTAGGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGATGCGACg                                                                                                                    >  2:503174/1‑150 (MQ=255)
                                        ttGCGTAGGACGCGATCCTGGAGATCTGACGACCATTGATGATCCTGCGTCGTTGGATCAGATCCGCCAGGCGATGGAAGAGTAGGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGATGCGACGCTGAc                                                                                                               <  2:148425/150‑1 (MQ=255)
                                        ttGCGAAGGACGCGATCCTGGAGATCTGACGACCATTGATGATCCTGCGTCGTTGGATCAGATCCGCCAGGCGATGGAAGAGTAGGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGATGCGACGCTGAc                                                                                                               <  2:1321899/150‑1 (MQ=255)
                                              tggACGCGTTCCTGGAGATCTGACGACCTTTGTTGATCCTGCGTGTTTGTATCAGTTCCGCCAGGCGATGGAAGAGTAGTTTTTTGTCGGTTGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGATGCGACGCTGACGCGttt                                                                                                         <  1:1619920/149‑1 (MQ=255)
                                              aGGACGCGATCCTGGAGATCTGACGACCATTGATGATCCTGCGTCGTTGGATCAGATCCGCCAGGCGATGGAAGAGTAGGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGATGCGACGCTGACGCGttt                                                                                                         >  2:1902418/1‑150 (MQ=255)
                                              aGGACGCGATCCTGGAGATCTGACGACCATTGATGATCCTGCGTCGTTGGATCAGATCCGCCAGGCGATGGAAGAGTAGGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGATGCGACGCTGACGCGttt                                                                                                         >  2:1008943/1‑150 (MQ=255)
                                              aGGACGCGATCCTGGAGATCTGACGACCATTGATGATCCTGCGTCGTTGGATCAGATCCGCCAGGCGATGGAAGAGTAGGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGATGCGACGCTGACGCGttt                                                                                                         >  2:1636375/1‑150 (MQ=255)
                                              aGGACGCGATCCTGGAGATCTGACGACCATTGATGATCCTGCGTCGTTGGATCAGATCCGCCAGGCGATGGAAGAGTAGGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGATGCGACGCTGACGCGttt                                                                                                         >  2:2338971/1‑150 (MQ=255)
                                                   gcgATCCTGGAGATCTGACGACCATTGATGATCCTGCGTCGTTGGATCAGATCCGCCAGGCGATGGAAGAGTAGGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGATGCGACGCTGACGCGTTTTATCa                                                                                                    >  1:1826185/1‑150 (MQ=255)
                                                   gcgATCCTGGAGATCTGACGACCATTGATGATCCTGCGTCGTTGGATCAGATCCGCCAGGCGATGGAAGAGTAGGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGATGCGACGCTGACGCGTTTTATCa                                                                                                    >  1:1826245/1‑150 (MQ=255)
                                                     gATCCTGGAGATCTGACGACCATTGATGATCCTGCGTCGTTGGATCAGATCCGCCAGGCGATGGAAGAGTAGGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGATGCGACGCTGACGCGTTTTATCATg                                                                                                  <  2:1895798/150‑1 (MQ=255)
                                                        ccTGGATATCTCACGACCATTGATGATCCTGCGTCGTTGGATCAGAACCGCCAGGCGATGGAAGAGTAGGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTTATGCGACGCTGACGCGATTTATCATGCCt                                                                                               <  1:680430/150‑1 (MQ=255)
                                                        ccTGGAGATCTGACGACCATTGATGATCCTGCGTCGTTGGATCAGATCCGCCAGGCGATGGAAGAGTAGGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGATGCGACGCTGACGCGTTTTATCATGCCt                                                                                               <  1:1331639/150‑1 (MQ=255)
                                                          tGGAGATCTGACGACCATTGATGATCCTGCGTCGTTGGATCAGATCCGCCAGGCGATGGAAGAGTAGGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGATGCGACGCTGACGCGTTTTATCATGCCTAc                                                                                             >  2:622432/1‑150 (MQ=255)
                                                          tGGAGATCTGACGACCATTGATGATCCTGCGTCGTTGGATCAGATCCGCCAGGCGATGGAAGAGTAGGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGATGCGACGCTGACGCGTTTTATCATGCCTAc                                                                                             >  2:1239355/1‑150 (MQ=255)
                                                                tCTGACGACCATTGATGATCCTGCGTCGTTGGATCAGATCCGCCAGGCGATGGAAGAGTTGGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGATGCGACGCTGACGCGTTTTATCATGCCTACGGACCt                                                                                       <  1:1520040/150‑1 (MQ=255)
                                                                tCTGACGACCATTGATGATCCTGCGTCGTTGGATCAGATCCGCCAGGCGATGGAAGAGTAGGTTATTGTCGGTTGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGATGCGACGCTGACGCGTTTTATCATGCCTACGGACCt                                                                                       <  1:899157/150‑1 (MQ=255)
                                                                tCTGACGACCATTGATGATCCTGCGTCGTTGGATCAGATCCGCCAGGCGATGGAAGAGTAGGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGATGCGACGCTGACGCGTTTTATCATGCCTACGGACCt                                                                                       <  2:58745/150‑1 (MQ=255)
                                                                tCTGACGACCATTGATGATCCTGCGTCGTTGGATCAGATCCGCCAGGCGATGGAAGAGTAGGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGATGCGACGCTGACGCGTTTTATCATGCCTACGGACCt                                                                                       <  2:114336/150‑1 (MQ=255)
                                                                     cgacCATTGATGATCCTGCGTCGTTGGATCAGATCCGCCAGGCGATGGAAGAGTAGGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGATGCGACGCTGACGCGTTTTATCATGCCTACGGACCTGAAcc                                                                                  <  2:1890047/150‑1 (MQ=255)
                                                                     cgacCATTGATGATCCTGCGTCGTTGGATCAGATCCGCCAGGCGATGGAAGAGTAGGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGATGCGACGCTGACGCGTTTTATCATGCCTACGGACCTGAAcc                                                                                  <  2:735894/150‑1 (MQ=255)
                                                                     cgacCATTGATGATCCTGCGTCGTTGGATCAGATCCGCCAGGCGATGGAAGAGTAGGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGATGCGACGCTGACGCGTTTTATCATGCCTACGGACCTGAAcc                                                                                  <  1:1797538/150‑1 (MQ=255)
                                                                     cgacCATTGATGATCCTGCGTCGTTGGATCAGATCCGCCAGGCGATGGAAGAGTAGGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGATGCGACGCTGACGCGTTTTATCATGCCTACGGACCTGAAcc                                                                                  <  2:621580/150‑1 (MQ=255)
                                                                      gacCATTGATGATCCTGCGTCGTTGGATCAGATCCGCCAGGCGATGGAAGAGTAGGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGATGCGACGCTGACGCGTTTTATCATGCCTACGGACCTGAACCg                                                                                 >  1:1978359/1‑150 (MQ=255)
                                                                      gacCATTGATGATCCTGCGTCGTTGGATCAGATCCGCCAGGCGATGGAAGAGTAGGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGATGCGACGCTGACGCGTTTTATCATGCCTACGGACCTGAACCg                                                                                 >  1:1276049/1‑150 (MQ=255)
                                                                      gacCATTGATGATCCTGCGTCGTTGGATCAGATCCGCCAGGCGATGGAAGAGTAGGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGATGCGACGCTGACGCGTTTTATCATGCCTACGGACCTGAACCg                                                                                 >  1:2582789/1‑150 (MQ=255)
                                                                       acCATTGATGATCCTGCGTCGTTGGATCAGATCCGCCAGGCGATGGAAGAGTAGGTTTTTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGATGCGACGCTGACGCGTTTTATCATGCCTACGGACCTGAACCGt                                                                                <  1:2211582/150‑1 (MQ=255)
                                                                       acCATTGATGATCCTGCGTCGTTGGATCAGATCCGCCAGGCGATGGAAGAGTAGGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGATGCGACGCTGACGCGTTTTATCATGCCTACGGACCTGAACCGt                                                                                <  1:1481569/150‑1 (MQ=255)
                                                                          aTTGATGATCCTGCGTCGTTGGATCAGATCCGCCAGGCGAGGGAAGAGTAGGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGATGCGACGCTGACGCGTTTTATCATGCCTACGGACCTGAACCGTAgg                                                                             <  2:479203/150‑1 (MQ=255)
                                                                           ttGATGATCCTGCGTCGTTGGATCAGATCCGCCAGGCGATGGAAGAGTAGGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGATGCGACGCTGACGCGTTTTATCATGCCTACGGACCTGAACCGTAGGt                                                                            >  2:1635752/1‑150 (MQ=255)
                                                                           ttGATGATCCTGCGTCGTTGGATCAGATCCGCCAGGCGATGGAAGAGTAGGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGATGCGACGCTGACGCGTTTTATCATGCCTACGGACCTGAAACGTAGGt                                                                            >  2:2271345/1‑150 (MQ=255)
                                                                                   ccTGCGTCGTTGGATCAGATCCGCCTGGCGATGGAAGAGTAGGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGATGCGACGCTGACGCGTTTTATCATGCCTACGGACCTGAACCGTAGGTCGGATAAg                                                                    <  2:2464054/150‑1 (MQ=255)
                                                                                   ccTGCGTCGTTGGATCAGATCCGCCAGGCGATGGAAGAGTAGGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGATGCGACGCTGACGCGTTTTATCATGCCTACGGACCTGAACCGTAGGTCGGATAAg                                                                    <  2:743105/150‑1 (MQ=255)
                                                                                        gtcgtTGGATCAGATCCGCCAGGCGATGGAAGAGTAGGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGATGCGACGCTGACGCGTTTTATCATGCCTACGGACCTGAACCGTAGGTCGGATAAGGCGCt                                                               >  2:1089213/1‑150 (MQ=255)
                                                                                        gtcgtTGGATCAGATCCGCCAGGCGATGGAAGAGTAGGTTATTGTCGGATGCGTAGCGCGGTGCATACGGCACTGTGTGCCGATGCCTGATGCGACGCTGACGCGTTTTATCATGCCTACGGACCTGAACCGTAGGTCGGATAAGGCGCt                                                               >  1:859863/1‑150 (MQ=255)
                                                                                         tcgtTGGATCAGATCCGCCAGGCGATGGAAGAGTAGGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGATGCGACGCTGACGCGTTTTATCATGCCTACGGACCTGAACCGTAGGTCGGATAAGGCGCTc                                                              >  2:852570/1‑150 (MQ=255)
                                                                                         tcgtTGGATCAGATCCGCCAGGCGATGGAAGAGTAGGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGATGCGACGCTGACGCGTTTTATCATGCCTACGGACCTGAACCGTAGGTCGGATAAGGCGCTc                                                              >  2:2695755/1‑150 (MQ=255)
                                                                                         tcgtTGGATCAGATCCGCCAGGCGATGGAAGAGTAGGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGATGCGACGCTGACGCGTTTTATCATGCCTACGGACCTGAACCGTAGGTCGGAGAAGGCGCTc                                                              >  2:2695716/1‑150 (MQ=255)
                                                                                             tGGATCAGATCCGCCAGGCGATGGAAGAGTAGGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGATGCGACGCTGACGCGTCTTATCATGCCTACGGACCTGAACCGTAGGTCGGATAAGGCGCTCGCGt                                                          <  2:1903623/150‑1 (MQ=255)
                                                                                               gATCAGATCCGCCAGGCGATGGAAGAGTAGGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGATGCGACGCTGACGCGTTTTATCATGCCTACGGACCTGAACCGTAGGTCGGATAAGGCGCTCGCGTCg                                                        >  2:2180062/1‑150 (MQ=255)
                                                                                               gATCAGATCCGCCAGGCGATGGAAGAGTAGGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGATGCGACGCTGACGCGTTTTATCATGCCTACGGACCTGAACCGTAGGTCGGATAAGGCGCTCGCGTCg                                                        >  2:2179780/1‑150 (MQ=255)
                                                                                               gATCAGATCCGCCAGGCGATGGAAGAGTAGGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGATGCGACGCTGACGCGTTTTATCATGCCTACGGACCTGAACCGTAGGTCGGATAAGGCGCTCGCGTCg                                                        >  2:1571366/1‑150 (MQ=255)
                                                                                               gATCAGATCCGCCAGGCGATGGAAGAGTAGGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGATGCGACGCTGACGCGTTTTATCATGCCTACGGACCTGAACCGTAGGTCGGATAAGGCGCTCGCGTCg                                                        >  1:1572450/1‑150 (MQ=255)
                                                                                               gATCAGATCCGCCAGGCGATGGAAGAGTAGGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGATGCGACGCTGACGCGTTTTATCATGCCTACGGACCTGAACCGTAGGTCGGATAAGGCGCTCGCGTCg                                                        >  1:1826394/1‑150 (MQ=255)
                                                                                               gATCAGATCCGCCAGGCGATGGAAGAGTAGGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGATGCGACGCTGACGCGTTTTATCATGCCTACGGACCTGAACCGTAGGTCGGATAAGGCGCTCGCGTCg                                                        >  1:805989/1‑150 (MQ=255)
                                                                                               gATCAGATCCGCCAGGCGATGGAAGAGTAGGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGATGCGACGCTGACGCGTTTTATCATGCCTACGGACCTGAACCGTAGGTCGGATAAGGCGCTCGCGTCg                                                        >  1:1847635/1‑150 (MQ=255)
                                                                                               gATCAGATCCGCCAGGCGATGGAAGAGTAGGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGATGCGACGCTGACGCGTTTTATCATGCCTACGGACCTGAACCGTAGGTCGGATAAGGCGCTCGCGTCg                                                        >  2:149455/1‑150 (MQ=255)
                                                                                                 tCAGATCCGCCAGGCGATGGAAGAGTAGGTTATTGTCGGATTCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGATGCGACGCTGACGCGTTTTATCATGCCTACGGACCTGAACCGTAGGTCGGATAATGCGCTCGCGTCGCa                                                      >  1:668853/1‑150 (MQ=255)
                                                                                                 tCAGATCCGCCAGGCGATGGAAGAGTAGGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGATGCGACGCTGACGCGTTTTATCATGCCTACGGACCTGAACCGTAGGTCGGATAAGGCGCTCGCGTCGCa                                                      >  1:2363284/1‑150 (MQ=255)
                                                                                                 tCAGATCCGCCAGGCGATGGAAGAGTAGGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGATGCGACGCTGACGCGTTTTATCATGCCTACGGACCTGAACCGTAGGTCGGATAAGGCGCTCGCGTCGCa                                                      >  1:603713/1‑150 (MQ=255)
                                                                                                 tCAGATCCGCCAGGCGATGGAAGAGTAGGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGATGCGACGCTGACGCGTTTTATCATGCCTACGGACCTGAACCGTAGGTCGGATAAGGCGCTCGCGTCGCa                                                      >  2:2075611/1‑150 (MQ=255)
                                                                                                 tCAGATCCGCCAGGCGATGGAAGAGTAGGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGATGCGACGCTGACGCGTTTTATCATGCCTACGGACCTGAACCGTAGGTCGGATAAGGCGCTCGCGTCGCa                                                      >  1:2029181/1‑150 (MQ=255)
                                                                                                 tCAGATCCGCCAGGCGATGGAAGAGGAGGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGAAGCCTGATGCGACGCTGACGCGTTTTATCATGCCTACGGACCTGAACCGTAGGTCGGATAAGGCGCTCGCGTCGCa                                                      >  2:1574992/1‑150 (MQ=255)
                                                                                                  cAGATCCGCCAGGCGATGGAAGAGTAGGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGATGCGACGCTGACGCGTTTTATCATGCCTACGGACCTGAACCGTAGGTCGGATAAGGCGCTCGCGTCGCAt                                                     >  2:1898991/1‑150 (MQ=255)
                                                                                                  cAGATCCGCCAGGCGATGGAAGAGTAGGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGATGCGACGCTGACGCGTTTTATCATGCCTACGGACCTGAACCGTAGGTCGGATAAGGCGCTCGCGTCGCAt                                                     >  2:31227/1‑150 (MQ=255)
                                                                                                   aGATCCGCCAGGCGATGGAAGAGTAGGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGATGCGACGCTGACGCGTTTTATCATGCCTACGGACCTGAACCGTAGGTCGGATAAGGCGCTCGCGTCGCATc                                                    >  2:1694156/1‑150 (MQ=255)
                                                                                                   aGATCCGCCAGGCGATGGAAGAGTAGGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGATGCGACGCTGACGCGTTTTATCATGCCTACGGACCTGAACCGTAGGTCGGATAAGGCGCTCGCGTCGCATc                                                    <  1:2661621/150‑1 (MQ=255)
                                                                                                   aGATCCGCCAGGCGATGGAAGAGTAGGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGATGCGACGCTGACGCGTTTTATCATGCCTACGGACCTGAACCGTAGGTCGGATAAGGCGCTCGCGTCGCATc                                                    <  2:358665/150‑1 (MQ=255)
                                                                                                   aGATCCGCCAGGCGATGGAAGAGTAGGTTATTGTCGGATGCGGCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGATGCGACGCTGACGCGTTTTATCATGCCTACGGACCTGAACCGTAGGTCGGATAAGGCGCTCGCGTCGCATc                                                    >  2:2504654/1‑150 (MQ=255)
                                                                                                         gCCAGGCGATGGAAGAGTAGGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGATGCGACGCTGACGCGTTTTATCATGCCTACGGACCTGAACCGTAGGTCGGATAAGGCGCTCGCGTCGCATCCGacac                                              <  2:2582789/150‑1 (MQ=255)
                                                                                                          ccAGGCGATGGAAGAGTAGGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGATGCGACGCTGACGCGTTTTATCATGCCTACGGACCTGAACCGTAGGTCGGATAAGGCGCTCGCGTCGCATCCGACAcc                                             >  2:261181/1‑150 (MQ=255)
                                                                                                             ggCGATGGAAGAGTAGGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGATGCGACGCTGACGCGTTTTATCATGCCTACGGACCTGAACCGTAGGTCGGATAAGGCGCTCGCGTCGCATCCGACACCATg                                          <  2:1173391/150‑1 (MQ=255)
                                                                                                             ggCGATGGAAGAGTAGGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGATGCGACGCTGACGCGTTTTATCATGCCTACGGACCTGAACCGTAGGTCGGATAAGGCGCTCGCGTCGCATCCGACACCATg                                          <  2:1037791/150‑1 (MQ=255)
                                                                                                             ggCGATGGAAGAGTAGGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGATGCGACGCTGACGCGTTTTATCATGCCTACGGACCTGAACCGTAGGTCGGATAAGGCGCTCGCGTCGCATCCGACACCATg                                          <  2:1887570/150‑1 (MQ=255)
                                                                                                                 aTGGAAGAGTAGGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGATGCGACGCTGACGCGTTTTATCATGCCTACGGACCTGAACCGTAGGTCGGATAAGGCGCTCGCGTCGCATCCGACACCATGCTCa                                      >  1:870973/1‑150 (MQ=255)
                                                                                                                               tATTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGATGCGACGCTGACGCGTTTTATCATGCCTACGGACCTGAACCGTAGGTCGGATAAGGCGCTCGCGTCGCATCCGACACCATGCTCAGATGCCTGATGCGa                        <  1:79650/150‑1 (MQ=255)
                                                                                                                                aTTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGATGCGACGCTGACGCGTTTTATCATGCCTACGGACCTGAACCGTAGGTCGGATAAGGCGCTCGCGTCGCATCCGACACCATGCTCAGATGCCTGATGCGAc                       <  2:287274/150‑1 (MQ=255)
                                                                                                                                   gTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGATGCGACGCTGACGCGTTTTATCATGCCTACGGACCTGAACCGTAGGTCGGATAAGGCGCTCGCGTCGCATCCGACACCATGCTCAGATGCCTGATGCGACGCt                    >  2:1224453/1‑150 (MQ=255)
                                                                                                                                    tCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGATGCGACGCTGACGCGTTTTATCATGCCTACGGACCTGAACCGTAGGTCGGATAAGGCGCTCGCGTCGCATCCGACACCATGCTCAGATGCCTGATGCGACGCTg                   >  2:2548056/1‑150 (MQ=255)
                                                                                                                                    ccGGATGCACCGCGCGACGCATCCGGCACTGTGTGCCGATGCCTGATGCGACGCTGACGCGTTTTATCATGCCTACGGACCTGAACCGTAGGTCGGATAAGGCGCTCGCGTCGCATCCGACACCATGCTCAGATGCCTGATGCGACGCTg                   >  1:1204843/2‑150 (MQ=255)
                                                                                                                                                 gcgGTGCATCCGGCACTGTGTGCCGATGCCTGATGCGACGCTGACGCGTTTTATCATGCCTACGGACCTGAACCGTAGGTCGGATAAGGCGCTCGCGTCGCATCCGACACCATGCTCAGATGCCTGATGCGACGCTGACGCGTCTTATCa      <  2:1734566/150‑1 (MQ=255)
                                                                                                                                                 gcgGTGCATCCGGCACTGTGTGCCGATGCCTGATGCGACGCTGACGCGTTTTATCATGCCTACGGACCTGAACCGTAGGTCGGATAAGGCGCTCGCGTCGCATCCGACACCATGCTCAGATGCCTGATGCGACGCTGACGCGTCTTATCa      <  2:2162581/150‑1 (MQ=255)
                                                                                                                                                 gcgGTGCATCCGGCACTGTGTGCCGATGCCTGATGCGACGCTGACGCGTTTTATCATGCCTACGGACCTGAACCGTAGGTCGGATAAGGCGCTCGCGTCGCATCCGACACCATGCTCAGATGCCTGATGCGACGCTGACGCGTCTTATCa      <  2:2162573/150‑1 (MQ=255)
                                                                                                                                                   ggTGCATCCGGCACTGTGTGCCGATGCCTGATGCGACGCTGACGCGTTTTATCATGCCTACGGACCTGAACCGTAGGTCGGATAAGGCGCTCGCGTCGCATCCGACACCATGCTCAGATGCCTGATGCGACGCTGACGCGTCTTATCAgg    <  2:2057221/150‑1 (MQ=255)
                                                                                                                                                     tGCATCCGGCACTGTGTGCCGATGCCTGATGCGACGCTGACGCGTTTTATCATGCCTACGGACCTGAACCGTAGGTCGGATCAGGCGCTCGCGTCGCATCCGACACCATGCTCAGATGCCTGATGCGACGCTGACGCGTCTTATCAGGcc  >  2:1862537/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CGCGATCCGGAAAAATGCTGCGCCGCACGATCCAGGCGATTTGCGAAGGACGCGATCCTGGAGATCTGACGACCATTGATGATCCTGCGTCGTTGGATCAGATCCGCCAGGCGATGGAAGAGTAGGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGATGCGACGCTGACGCGTTTTATCATGCCTACGGACCTGAACCGTAGGTCGGATAAGGCGCTCGCGTCGCATCCGACACCATGCTCAGATGCCTGATGCGACGCTGACGCGTCTTATCAGGCC  >  NC_000913/354468‑354766

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