Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,363,448 A→G 5.2% D62D (GAT→GAC) ‡ yfaY ← inactive PncC family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,363,4480AG 5.2% 327.0 / 5.9 96D62D (GAT→GAC) ‡yfaYinactive PncC family protein
Reads supporting (aligned to +/- strand):  ref base A (41/50);  new base G (2/3);  total (43/53)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.99e-01

CATTACTCGTCCACGTTCGTGGAAATAGCGTTCCATCTCTTTGAGCCAGGCTTCATGCAGCACCAGGCCTTCACCTTTTGCTGTCGCAGCGGCGAGTGCGCTTAAATCATCGCTGGTCGGTCCCAGCCCGCCGTTAACGATCAGCACATCGGCGTGCTGACTACGTTCGCGAAGAATGGTGACTAAGTCATCAAGGTTATCCCCCACCGTATTGCGGCGAGATAATGGCAACCCCTGATGAAAGAAAAAATCGGCCAGCCAGGCAGCGTTAGTGTCAACGATTTGCCCGTGTAACAC  >  NC_000913/2363301‑2363597
                                                                                                                                                   |                                                                                                                                                     
cATTACTCGTCCACGTTCGTGGAAATAGCGTTCCATCTCTTTGAGCCAGGCTTCATGCAGCACCAGGCCTTCACCTTTTGCTGTCGCAGCGGCGAGTGCGCTTAAATCATCGCTGGTCGGTCCCAGCCCGCCGTTAACGATCAGCACATc                                                                                                                                                     >  2:166266/1‑150 (MQ=255)
  ttACTCGTCCACGTTCGTGGAAATAGCGTTCCATCTCTTTGAGCCAGGCTTCATGCAGCACCAGGCCTTCACCTTTTGCTGTCGCAGCGGCGAGTGCGCTTAAATCATCGCTGGTCGGTCCCAGCCCGCCGTTAACGATCAGCACATCgg                                                                                                                                                   >  2:1114639/1‑150 (MQ=255)
  ttACTCGTCCACGTTCGTGGAAATAGCGTTCCATCTCTTTGAGCCAGGCTTCATGCAGCACCAGGCCTTCACCTTTTGCTGTCGCAGCGGCGAGTGCGCTTAAATCATCGCTGGTCGGTCCCAGCCCGCCGTTAACGATCAGCACATCgg                                                                                                                                                   >  2:1114607/1‑150 (MQ=255)
   tACTCGTCCACGTTCGTGGAAATAGCGTTCCATCTCTTTGAGCCAGGCTTCATGCAGCACCAGGCCTTCACCTTTTGCTGTCGCAGCGGCGAGTGCGCTTAAATCATCGCTGGTCGGTCCCAGCCCGCCGTTAACGATCAGCACATCGGc                                                                                                                                                  <  1:1026155/150‑1 (MQ=255)
      tCGTCCACGTTCGTGGAAATAGCGTTCCATCTCTTTGAGCCAGGCTTCATGCAGCACCAGGCCTTCACCTTTTGCTGTCGCAGCGGCGAGTGCGCTTAAATCATCGCTGGTCGGTCCCAGCCCGCCGTTAACGATCAGCACATCGGCGTg                                                                                                                                               >  2:1399400/1‑150 (MQ=255)
      tCGTCCACGTTCGTGGAAATAGCGTTCCATCTCTTTGAGCCAGGCTTCATGCAGCACCAGGCCTTCACCTTTTGCTGTCGCAGCGGCGAGTGCGCTTAAATCATCGCTGGTCGGTCCCAGCCCGCCGTTAACGATCAGCACATCGGCGTg                                                                                                                                               >  1:994413/1‑150 (MQ=255)
      tCGTCCACGTTCGTGGAAATAGCGTTCCATCTCTTTGAGCCAGGCTTCATGCAGCACCAGGCCTTCACCTTTTGCTGTCGCAGCGGCGAGTGCGCTTAAATCATCGCTGGTCGGTCCCAGCCCGCCGTTAACGATCAGCACATCGGCGTg                                                                                                                                               <  1:854454/150‑1 (MQ=255)
      tCGTCCACGTTCGTGGAAATAGCGTTCCATCTCTTTGAGCCAGGCTTCATGCAGCACCAGGCCTTCACCTTTTGCTGTCGCAGCGGCGAGTGCGCTTAAATCATCGCTGGTCGGTCCCAGCCCGCCGTTAACGATCAGCACATCGGCGTg                                                                                                                                               >  1:1241798/1‑150 (MQ=255)
        gTCCACGTTCGTGGAAATAGCGTTCCATCTCTTTGAGCCAGGCTTCATGCAGCACCAGGCCTTCACCTTTTGCTGTCGCAGCGGCGAGTGCGCTTAAATCATCGCTGGTCGGTCCCAGCCCGCCGTTAACGATCAGCACATCGGCGTGCt                                                                                                                                             <  1:143557/150‑1 (MQ=255)
          ccACGTTCGTGGAAATAGCGTTCCATCTCTTTGAGCCAGGCTTCATGCAGCACCAGGCCTTCACCTTTTGCTGTCGCAGCGGCGAGTGCGCTTAAATCATCGCTGGTCGGTCCCAGCCCGCCGTTAACGATCAGCACATCGGCGTGCTGa                                                                                                                                           <  2:823247/150‑1 (MQ=255)
           cACGTTCGTGGAAATAGCGTTCCATCTCTTTGAGCCAGGCTTCATGCAGCACCAGGCCTTCACCTTTTGCTGTCGCAGCGGCGAGTGCGCTTAAATCATCGCTGGTCGGTCCCAGCCCGCCGTTAACGATCAGCACATCGGCGTGCTGAc                                                                                                                                          >  1:665800/1‑150 (MQ=255)
            aCGTTCGTGGAAATAGCGTTCCATCTCTTTGAGCCAGGCTTCATGCAGCACCAGGCCTTCACCTTTTGCTGTCGCAGCGGCGAGTGCGCTTAAATCATCGCTGGTCGGTCCCAGCCCGCCGTTAACGATCAGCACGCCGATGTGCTGatc                                                                                                                                         <  1:1168460/150‑3 (MQ=255)
            aCGTTCGTGGAAATAGCGTTCCATCTCTTTGAGCCAGGCTTCATGCAGCACCAGGCCTTCACCTTTTGCTGTCGCAGCGGCGAGTGCGCTTAAATCATCGCTGGTCGGTCCCAGCCCGCCGTTAACGATCAGCACGCCGATGTGCTGatc                                                                                                                                         <  1:1168457/150‑3 (MQ=255)
            aCGTTCGTGGAAATAGCGTTCCATCTCTTTGAGCCAGGCTTCATGCAGCACCAGGCCTTCACCTTTTGCTGTCGCAGCGGCGAGTGCGCTTAAATCATCGCTGGTCGGTCCCAGCCCGCCGTTAACGATCAGCACATCGGCGTGCTGACt                                                                                                                                         <  2:860210/150‑1 (MQ=255)
                tCGTGGAAATAGCGTTCCATCTCTTTGAGCCAGGCTTCATGCAGCACCAGGCCTTCACCTTTTGCTGTCGCAGCGGCGAGTGCGCTTAAATCATCGCTGGTCGGTCCCAGCCCGCCGTTAACGATCAGCACGCCGATGTGCTGatcgtta                                                                                                                                     <  2:348882/150‑7 (MQ=255)
                  gTGGAAATAGCGTTCCATCTCTTTGAGCCAGGCTTCATGCAGCACCAGGCCTTCACCTTTTGCTGTCGCAGCGGCGAGTGCGCTTAAATCATCGCTGGTCGGTCCCAGCCCGCCGTTAACGATCAGCACATCGGCGTGCTGACTACGTTc                                                                                                                                   >  2:350061/1‑150 (MQ=255)
                  gTGGAAATAGCGTTCCATCTCTTTGAGCCAGGCTTCATGCAGCACCAGGCCTTCACCTTTTGCTGTCGCAGCGGCGAGTGCGCTTAAATCATCGCTGGTCGGTCCCAGCCCGCCGTTAACGATCAGCACATCGGCGTGCTGACTACGTTc                                                                                                                                   >  2:582752/1‑150 (MQ=255)
                     gAAATAGCGTTCCATCTCTTTGAGCCAGGCTTCATGCAGCACCAGGCCTTCACCTTTTGCTGTCGCAGCGGCGAGTGCGCTTAAATCATCGCTGGTCGGTCCCAGCCCGCCGTTAACGATCAGCACATCGGCGTGCTGACTACGTTcgcg                                                                                                                                <  1:565586/150‑1 (MQ=255)
                     gAAATAGCGTTCCATCTCTTTGAGCCAGGCTTCATGCAGCACCAGGCCTTCACCTTTTGCTGTCGCAGCGGCGAGTGCGCTTAAATCATCGCTGGTCGGTCCCAGCCCGCCGTTAACGATCAGCACATCGGCGTGCTGACTACGTTcgcg                                                                                                                                <  1:741764/150‑1 (MQ=255)
                          agggTTCCATCTCTTTGAGGGAGGCTTCTTGCAGCACCCGGCGGTCACCTTTTGCTGTCGCAGCGGCGAGTGCGCTTAAATCCTCGCTGGTCGGTCCCAGCCCGCCGTTAACGCTCTGCACATCGGCGTGGTGACTACGTTCGCgaagaa                                                                                                                           <  2:1235182/147‑1 (MQ=255)
                          aGCGTTCCATCTCTTTGAGCCAGGCTTCATGCAGCACCAGGCCTTCACCTTTTGCTGGCGCAGCGGCGAGTGCGCTTAAATCATCGCTGGTCGGTCCCAGCCCGCCGTTAACGATCAGCACATCGGCGTGCTGACTACGTTCGCgaagaa                                                                                                                           <  2:1235181/150‑1 (MQ=255)
                            cGTTCCATCTCTTTGAGCCAGGCTTCATGCAGCACCAGGCCTTCACCTTTTGCTGTCGCAGCGGCGAGTGCGCTTAAATCATCGCTGGTCGGTCCCAGCCCGCCGTTAACGATCAGCACATCGGCGTGCTGACTACGTTCGCGAAGAATg                                                                                                                         >  1:1168967/1‑150 (MQ=255)
                              ttCCATCTCTTTGTGACTGGCTTCTTGCTGCACCAGGCCTTCACCTTTTGCTGTCGCAGCGGCGTGTGCGCTTAAATAATCGCTGGTCGGTCCCAGCCCGCCGTTAACGATCAGCACATCGGCGTTCTGACTACGTTCGCGAAGAATGGt                                                                                                                       <  2:1392887/150‑1 (MQ=255)
                              ttCCATCTCTTTGAGCCAGGCTTCATGCAGCACCAGGCCTTCACCTTTTGCTGTCGCAGCGGCGAGTGCGCTTAAATCATCGCTGGTCGGTCCCAGCCCGCCGTTAACGATCAGCACATCGGCGTGCTGACTACGTTCGCGAAGAATGGt                                                                                                                       <  2:520665/150‑1 (MQ=255)
                              ttCCATCTCTTTGAGCCAGGCTTCATGCAGCACCAGGCCTTCACCTTTTGCTGTCGCAGCGGCGAGTGCGCTTAAATCATCGCTGGTCGGTCCCAGCCCGCCGTTAACGATCAGCACATCGGCGTGCTGACTACGTTCGCGAAGAATGGt                                                                                                                       <  1:825332/150‑1 (MQ=255)
                              ttCCATCTCTTTGAGCCAGGCTTCATGCAGCACCAGGCCTTCACCTTTTGCTGTCGCAGCGGCGAGTGCGCTTAAATCATCGCTGGTCGGTCCCAGCCCGCCGTTAACGATCAGCACATCGGCGTGCTGACTACGTTCGCGAAGAATGGt                                                                                                                       <  2:11674/150‑1 (MQ=255)
                               tCCATCTCTTTGAGCCAGGCTTCATGCAGCACCAGGCCTTCACCTTTTGCTGTCGCAGCGGCGAGTGCGCTTAAATCATCGCTGGTCGGTCCCAGCCCGCCGTTAACGATCAGCACATCGGCGTGCTGACTACGTTCGCGAAGAATGGTg                                                                                                                      >  2:731164/1‑150 (MQ=255)
                               tCCATCTCTTTGAGCCAGGCTTCATGCAGCACCAGGCCTTCACCTTTTGCTGTCGCAGCGGCGAGTGCGCTTAAATCATCGCTGGTCGGTCCCAGCCCGCCGTTAACGATCAGCACATCGGCGTGCTGACTACGTTCGCGAAGAATGGTg                                                                                                                      >  2:379870/1‑150 (MQ=255)
                                ccATCTCTTTGAGCCAGGCTTCATGCAGCACCCGGCCTTCACCTTTTGCTGTCGCAGCGGCGAGTGCGCTTAAATCATCGCTGGTCGGTCCCAGCCCGCCGTTAACGATCAGCACATCGGCGTGCTGACTACGTTCGCGAAGAATGGTGa                                                                                                                     <  2:915593/150‑1 (MQ=255)
                                ccATCTCTTTGAGCCAGGCTTCATGCAGCACCAGGCCTTCACCTTTTGCTGTCGCAGCGGCGAGTGCGCTTAAATCATCGCTGGTCGGTCCCAGCCCGCCGTTAACGATCAGCACATCGGCGTGCTGACTACGTTCGCGAAGAATGGTGa                                                                                                                     <  2:590028/150‑1 (MQ=255)
                                ccATCTCTTTGAGCCAGGCGTCATGCAGCACCAGGCCTTCACCTTTTGCTGTCGCAGCGGCGAGTGCGCTTCAATCATCGCTGGTCGTTCCCATCCCGCCGTTAACGATCAGCACATCGGCGTGCTGACTACGTTCGCGTAGAATGGTGa                                                                                                                     <  2:775599/150‑1 (MQ=255)
                                   tctcTTTGAGCCAGGCTTCATGCAGCACCAGGCCTTCACCTTTTGCTGTCGCAGCGGCGAGTGCGCTTAAATCATCGCTGGTCGGTCCCAGCCCGCCGTTAACGATCAGCACATCGGCGTGCTGACTACGTTCGCGAAGAATGGTGACTa                                                                                                                  >  1:527910/1‑150 (MQ=255)
                                       tttGAGCCAGGCTTCATGCAGCACCAGGCCTTCACCTTTTGCTGTCGCAGCGGCGAGTGCGCTTAAATCATCGCTGGTCGGTCCCAGCCCGCCGTTAACGATCAGCACATCGGCGTGCTGACTACGTTCGCGAAGAATGGTGACTAAGtc                                                                                                              >  2:808440/1‑150 (MQ=255)
                                       tttGAGCCAGGCTTCATGCAGCACCAGGCCTTCACCTTTTGCTGTCGCAGCGGCGAGTGCGCTTAAATCATCGCTGGTCGGTCCCAGCCCGCCGTTAACGATCAGCACATCGGCGTGCTGACTACGTTCGCGAAGAATGGTGACTAAGtc                                                                                                              <  2:567138/150‑1 (MQ=255)
                                       tttGAGCCAGGCTTCATGCAGCACCAGGCCTTCACCTTTTGCTGTCGCAGCGGCGAGTGCGCTTAAATCATCGCTGGTCGGTCCCAGCCCGCCGTTAACGATCAGCACATCGGCGTGCTGACTACGTTCGCGAAGAATGGTGACTAAGtc                                                                                                              >  1:720437/1‑150 (MQ=255)
                                        ttGAGCCAGGCTTCATGCAGCACCAGGCCTTCACCTTTTGCTGTCGCAGCGGCGAGTGCGCTTAAATCATCGCTGGTCGGTCCCAGCCCGCCGTTAACGATCAGCACATCGGCGTGCTGACTACGTTCGCGAAGAATGGTGACTAAGtca                                                                                                             <  2:34870/150‑1 (MQ=255)
                                          gcgCCAGGCTTCCTGCCGCCCCCGGCCTGCACCTCGCGCTGGCGCAGCGGCGAGTGCGCATAACCCATCGCTGGACGGTCCCAGCCCGCCGAGACCGACCAGCACAGCGGCGTGCTGACAACGTTCGCGAAGAATGGTGACTAAGtcatc                                                                                                           <  2:1052265/148‑1 (MQ=255)
                                           aGCCAGGCTTCATGCAGCACCAGGCCTTCACCTTTTGCTGTCGCAGCGGCGAGTGCGCTTAAATCATCGCTGGTCGGTCCCAGCCCGCCGTTAACGATCAGCACATCGGCGTGCTGACTACGTTCGCGAAGAATGGTGACTAAGtcatca                                                                                                          >  1:1367445/1‑150 (MQ=255)
                                             ccAGGCTTCATGCAGCACCAGGCCTTCACCTTTTGCTGTCGCAGCGGCGAGTGCGCTTAAATCATCGCTGGTCGGTCCCAGCCCGCCGTTAACGATCAGCACATCGGCGTGCTGACTACGTTCGCGAAGAATGGTGACTAAGTCATCAAg                                                                                                        <  1:855756/150‑1 (MQ=255)
                                              cAGGCTTCATGCAGCACCAGGCCTTCACCTTTTGCTGTCGCAGCGGCGAGTGCGCTTAAATCATCGCTGGTCGGTCCCAGCCCGCCGTTAACGATCAGCACATCGGCGTGCTGACTACGTTCGCGAAGAATGGTGACTAAGTCATCAAgg                                                                                                       >  2:104605/1‑150 (MQ=255)
                                              cAGGCTTCATGCAGCACCAGGCCTTCACCTTTTGCTGTCGCAGCGGCGAGTGCGCTTAAATCATCGCTGGTCGGTCCCAGCCCGCCGTTAACGATCAGCACATCGGCGTGCTGACTACGTTCGCGAAGAATGGTGACTAAGTCATCAAgg                                                                                                       >  2:1133360/1‑150 (MQ=255)
                                              cAGGCTTCATGCAGCACCAGGCCTTCACCTTTTGCTGTCGCAGCGGCGAGTGCGCTTAAATCATCGCTGGTCGGTCCCAGCCCGCCGTTAACGATCAGCACATCGGCGTGCTGACTACGTTCGCGAAGAATGGTGACTAAGTCATCAAgg                                                                                                       >  2:448726/1‑150 (MQ=255)
                                              cAGGCTTCATGCAGCACCACGCCTTCACCTTTTGCTGTCGCAGCGGCGAGTGCGCTTAAATCATCGCTGGTCGGTCCCAGCCCGCCGTTAACGATCAGCACATCGGCGTGCTGACTACGTTCGCGAAGAATGGTGACTAAGTCATCAAgg                                                                                                       <  1:396888/150‑1 (MQ=255)
                                               aGGCTTCATGCAGCACCAGGCCTTCACCTTTTGCTGTCGCAGCGGCGAGTGCGCTTAAATCATCGCTGGTCGGTCCCAGCCCGCCGTTAACGATCAGCACATCGGCGTGCTGACTACGTTCGCGAAGAATGGTGACTAAGTCATCAAGGt                                                                                                      >  1:442416/1‑150 (MQ=255)
                                               aGGCTTCATGCAGCACCAGGCCTTCACCTTTTGCTGTCGCAGCGGCGAGTGCGCTTAAATCATCGCTGGTCGGTCCCAGCCCGCCGTTAACGATCAGCACATCGGCGTGCTGACTACGTTCGCGAAGAATGGTGACTAAGTCATCAAGGt                                                                                                      >  1:442380/1‑150 (MQ=255)
                                                   ttCATGCAGCACCAGGCCTTCACCTTTTGCTGTCGCAGCGGCGAGTGCGCTTAAATCATCGCTGGTCGGTCCCAGCCCGCCGTTAACGATCAGCACATCGGCGTGCTGACTACGTTCGCGAAGAATGGTGACTAAGTCATCAAGGTTATc                                                                                                  <  1:299907/150‑1 (MQ=255)
                                                     cATGCAGCACCAGGCCTTCACCTTTTGCTGTCGCAGCGGCGAGTGCGCTTAAATCATCGCTGGTCGGTCCCAGCCCGCCGTTAACGATCAGCACATCGGCGTGCTGACTACGTTCGCGAAGAATGGTGACTAAGTCATCAAGGTTATccc                                                                                                <  1:135384/150‑1 (MQ=255)
                                                     cATGCAGCACCAGGCCTTCACCTTTTGCTGTCGCAGCGGCGAGTGCGCTTAAATCATCGCTGGTCGGTCCCAGCCCGCCGTTAACGATCAGCACATCGGCGTGCTGACTACGTTCGCGAAGAATGGTGACTAAGTCATCAAGGTTATccc                                                                                                <  1:135321/150‑1 (MQ=255)
                                                        gcagcaCCAGGCCTTCACCTTTTGCTGTCGCAGCGGCGAGTGCGCTTAAATCATCGCTGGTCGGTCCCAGCCCGCCGTTAACGATCAGCACATCGGCGTGCTGACTACGTTCGCGAAGAATGGTGACTAAGTCATCAAGGTTATCCCCCa                                                                                             <  2:51981/150‑1 (MQ=255)
                                                                aGGCCTTCACCTTTTGCTGTCGCAGCGGCGAGTGCGCTTAAATCATCGCTGGTCGGTCCCAGCCCGCCGTTAACGATCAGCACATCGGCGTGCTGACTACGTTCGCGAAGAATGGTGACTAAGTCATCAAGGTTATCCCCCACCGTATTg                                                                                     >  1:938893/1‑150 (MQ=255)
                                                                  gCCTTCACCTTTTGCTGTCGCAGCGGCGAGTGCGCTTAAATCATCGCTGGTCGGTCCCAGCCCGCCGTTAACGATCAGCACATCGGCGTGCTGACTACGTTCGCGAAGAATGGTGACTAAGTCATCAAGGTTATCCCCCACCGTATTgcg                                                                                   >  2:1189958/1‑150 (MQ=255)
                                                                   ccTTCACCTTTTGCTGTCGCAGCGGCGAGTGCGCTTAAATCATCGCTGGTCGGTCCCAGCCCGCCGTTAACGATCAGCACATCGGCGTGCTGACTACGTTCGCGAAGAATGGTGACTAAGTCATCAAGGTTATCCCCCACCGTATTgcgg                                                                                  <  1:1160985/150‑1 (MQ=255)
                                                                   ccTTCACCTTTTGCTGTCGCAGCGGCGAGTGCGCTTAAATCATCGCTGGTCGGTCCCAGCCCGCCGTTAACGATCAGCACATCGGCGTGCTGACTACGTTCGCGAAGAATGGTGACTAAGTCATCAAGGTTATCCCCCACCGTATTgcgg                                                                                  <  1:1160931/150‑1 (MQ=255)
                                                                    cTTCACCTTTTGCTGTCGCAGCGGCGAGTGCGCTTAAATCATCGCTGGTCGGTCCCAGCCCGCCGTTAACGATCAGCACATCGGCGTGCTGACTACGTTCGCGAAGAATGGTGACTAAGTCATCAAGGTTATCCCCCACCGTATTgcggc                                                                                 >  2:848208/1‑150 (MQ=255)
                                                                    cTTCACCTTTTGCTGTCGCAGCGGCGAGTGCGCTTAAATCATCGCTGGTCGGTCCCAGCCCGCCGTTAACGATCAGCACATCGGCGTGCTGACTACGTTCGCGAAGAATGGTGACTAAGTCATCAAGGTTATCCCCCACCGTATTgcggc                                                                                 >  1:1042251/1‑150 (MQ=255)
                                                                    cTTCACCTTTTGCTGTCGCAGCGGCGAGTGCGCTTAAATCATCGCTGGTCGGTCCCAGCCCGCCGTTAACGATCAGCACATCGGCGTGCTGACTACGTTCGCGAAGAATGGTGACTAAGTCATCAAGGTTATCCCCCACCGTATTgcggc                                                                                 <  1:339699/150‑1 (MQ=255)
                                                                    cTTCACCTTTTGCTGTCGCAGCGGCGAGTGCGCTTAAATCATCGCTGGTCGGTCCCAGCCCGCCGTTAACGATCAGCACATCGGCGTGCTGACTACGTTCGCGAAGAATGGTGACTAAGTCATCAAGGTTATCCCCCACCGTATTgcggc                                                                                 <  1:504653/150‑1 (MQ=255)
                                                                    cTTCACCTTTTGCTGTCGCAGCGGCGAGTGCGCTTAAATCATCGCTGGTCGGTCCCAGCCCGCCGTTAACGATCAGCACATCGGCGTGCTGACTACGTTCGCGAAGAATGGTGACTAAGTCATCAAGGTTATCCCCCACCGTATTgcggc                                                                                 >  1:584008/1‑150 (MQ=255)
                                                                     ttaaCCTTTTGCTGTCGCAGCGGCGAGTGCGCTTAAATCATCGCTGGTCGGTCCCAGCCCGACGTTAACGATCAGCACATCGGCGTGCTGACTCCGTTCGCGAAGAATGGTGACTAAGTCATCAAGGTTATCCCCCACCGTATTgaggcg                                                                                >  1:1266772/4‑150 (MQ=255)
                                                                     ttCACCTTTTGCTGTCGCAGCGGCGAGTGCGCTTAAATCATCGCTGGTCGGTCCCAGCCCGCCGTTAACGATCAGCACATCGGCGTGCTGACTACGTTCGCGAAGAATGGTGACTAAGTCATCAAGGTTATCCCCCACCGTATTgcggcg                                                                                >  1:1266775/1‑150 (MQ=255)
                                                                     ttCACCTTTTGCTGTCGCAGCGGCGAGTGCGCTTAAATCATCGCTGGTCGGTCCCAGCCCGCCGTTAACGATCAGCACATCGGCGTGCTGACTACGTTCGCGAAGAATGGTGACTAAGTCATCAAGGTTATCCCCCACCGTATTgcggcg                                                                                >  1:528212/1‑150 (MQ=255)
                                                                     ttCACCTTTTGCTGTCGCAGCGGCGAGTGCGCTTAAATCATCGCTGGTCGGTCCCAGCCCGCCGTTAACGATCAGCACATCGGCGTGCTGACTACGTTCGCGAAGAATGGTGACTAAGTCATCAAGGTTATCCCCCACCGTATTgcggcg                                                                                >  2:225994/1‑150 (MQ=255)
                                                                            tttGCTGTCGCAGCGGCGAGTGCGCTTAAATCATCGCTGGTCGGTCCCAGCCCGCCGTTAACGATCAGCACATCGGCGTGCTGACTACGTTCGCGAAGAATGGTGACTAAGTCATCAAGGTTATCCCCCACCGTATTGCGGCGAGATAAt                                                                         <  1:797783/150‑1 (MQ=255)
                                                                                  gTCGCAGCGGCGAGTGCGCTTAAATCATCGCTGGTCGGTCCCAGCCCGCCGTTAACGATCAGCACATCGGCGTGCTGACTACGTTCGCGAAGAATGGTGACTAAGTCATCAAGGTTATCCCCCACCGTATTGCGGCGAGATAATGGCAAc                                                                   <  1:641663/150‑1 (MQ=255)
                                                                                  gTCGCAGCGGCGAGTGCGCTTAAATCATCGCTGGTCGGTCCCAGCCCGCCGTTAACGATCAGCACATCGGCGTGCTGACTACGTTCGCGAAGAATGGTGACTAAGTCATCAAGGTTATCCCCCACCGTATTGCGGCGAGATAATGGCAAc                                                                   <  1:641594/150‑1 (MQ=255)
                                                                                  gTCGCAGCGGCGAGTGCGCTTAAATCATCGCTGGTCGGTCCCAGCCCGCCGTTAACGATCAGCACATCGGCGTGCTGACTACGTTCGCGAAGAATGGTGACTAAGTCATCAAGGTTATCCCCCACCGTATTGCGGCGAGATAATGGCAAc                                                                   <  1:641575/150‑1 (MQ=255)
                                                                                      cAGCGGCGAGTGCGCTTAAATCATCGCTGGTCGGTCCCAGCCCGCCGTTAACGATCAGCACATCGGCGTGCTGACTACGTTCGCGAAGAATGGTGACTAAGTCATCAAGGTTATCCCCCACCGTATTGCGGCGAGATAATGGCAACCCCt                                                               <  1:509872/150‑1 (MQ=255)
                                                                                      cAGCGGCGAGTGCGCTTAAATCATCGCTGGTCGGTCCCAGCCCGCCGTTAACGATCAGCACATCGGCGTGCTGACTACGTTCGCGAAGAATGGTGACTAAGTCATCAAGGTTATCCCCCACCGTATTGCGGCGAGATAATGGCAACCCCt                                                               <  1:509913/150‑1 (MQ=255)
                                                                                        gcggcgAGTGCGCTTAAATCATCGCTGGTCGGTCCCAGCCCGCCGTTAACGATCAGCACATCGGCGTGCTGACTACGTTCGCGAAGAATGGTGACTAAGTCATCAAGGTTATCCCCCACCGTATTGCGGCGAGATAATGGCAACCCCtga                                                             >  1:974584/1‑150 (MQ=255)
                                                                                        gcggcgAGTGCGCTTAAATCATCGCTGGTCGGTCCCAGCCCGCCGTTAACGATCAGCACATCGGCGTGCTGACTACGTTCGCGAAGAATGGTGACTAAGTCATCAAGGTTATCCCCCACCGTATTGCGGCGAGATAATGGCAACCCCtga                                                             <  2:448035/150‑1 (MQ=255)
                                                                                             gAGTGCGCTTATATCATCGCTGGTCGGTCCCAGCGCGCCGTTAACGATCAGCTCATCGGCGTGCTGACTACGTTCGCGAATAATGGTGACTAAGTCATCTAGGTTATCCCCCACCTTATTGCGGCGATATAATGGCAACCCCTGATgaaa                                                        <  2:793551/150‑1 (MQ=255)
                                                                                             gAGTGCGCTTAAATCATCGCTGGTCGGTCCCAGCCCGCCGTTAACGATCAGCACATCGGCGTGCTGACTACGTTCGCGAAGAATGGTGACTAAGTCATCAAGGTTATCCCCCACCGTATTGCGGCGAGATAATGGCAACCCCTGATgaaa                                                        <  1:574012/150‑1 (MQ=255)
                                                                                             gAGTGCGCTTAAATCATCGCTGGTCGGTCCCAGCCCGCCGTTAACGATCAGCACATCGGCGTGCTGACTACGTTCGCGAAGAATGGTGACTAAGTCATCAAGGTTATCCCCCACCGTATTGCGGCGAGATAATGGCAACCCCTGATgaaa                                                        <  1:236516/150‑1 (MQ=255)
                                                                                             gAGTGCGCTTAAATCATCGCTGGTCGGTCCCAGCCCGCCGTTAACGATCAGCACATCGGCGTGCTGACTACGTTCGCGAAGAATGGTGACTAAGTCATCAAGGTTATCCCCCACCGTATTGCGGCGAGATAATGGCAACCCCTGATgaaa                                                        <  1:1226747/150‑1 (MQ=255)
                                                                                             gAGTGCGCTTAAATCATCGCTGGTCGGTCCCAGCCCGCCGTTAACGATCAGCACATCGGCGTGCTGACTACGTTCGCGAAGAATGGTGACTAAGTCATCAAGGTTATCCCCCACCGTATTGCGGCGAGATAATGGCAACCCCTGATgaaa                                                        <  2:1129755/150‑1 (MQ=255)
                                                                                             gAGTGCGCTTAAATCATCGCTGGTCGGTCCCAGCCCGCCGTTAACGATCAGCACATCGGCGTGCTGACTACGTTCGCGAAGAATGGTGACTAAGTCATCAAGGTTATCCCCCACCGTATTGCGGCGAGATAATGGCAACCCCTGATgaaa                                                        <  1:874063/150‑1 (MQ=255)
                                                                                                      tAAATCATCGCTGGTCGGTCCCAGCCCGCCGTTAACGATCAGCACATCGGCGTGCTGACTACGTTCGCGAAGAATGGTGACTAAGTCATCAAGGTTATCCCCCACCGTATTGCGGCGAGATAATGGCAACCCCTGATGAAAGAAAAAATc                                               <  1:1206813/150‑1 (MQ=255)
                                                                                                       aaaTCATCGCTGGTCGGTCCCAGCCCGCCGTTAACGATCAGCACATCGGCGTGCTGACTACGTTCGCGAAGAATGGTGACTAAGTCATCAAGGTTATCCCCCACCGTATTGCGGCGAGATAATGGCAACCCCTGATGAAAGAAAAAATCg                                              <  2:759619/150‑1 (MQ=255)
                                                                                                        aaTCATCGCTGGTCGGTCCCAGCCCGCCGTTAACGATCAGCACATCGGCGTGCTGACTACGTTCGCGAAGAATGGTGACTAAGTCATCAAGGTTATCCCCCACCGTATTGCGGCGAGATAATGGCAACCCCTGATGAAAGAAAAAATCgg                                             <  2:1313329/150‑1 (MQ=255)
                                                                                                            atcGCTGGTCGGTCCCAGCCCGCCGTTAACGATCAGCACATCGGCGTGCTGACTACGTTCGCGAAGAATGGTGACTAAGTCATCAAGGTTATCCCCCACCGTATTGCGGCGAGATAATGGCAACCCCTGATGAAAGAAAAAATCGgccag                                         >  1:1258251/1‑150 (MQ=255)
                                                                                                               gCTGGTCGGTCCCAGCCCGCCGTTAACGATCAGCACATCGGCGTGCTGACTACGTTCGCGAAGAATGGTGACTAAGTCATCAAGGTTATCCCCCACCGTATTGCGGCGAGATAATGGCAACCCCTGATGAAAGAAAAAATCGgccagcca                                      <  2:1128063/150‑1 (MQ=255)
                                                                                                                    tcggtcCCAGCCCGCCGTTAACGATCAGCACATCGGCGTGCTGACTACGTTCGCGAAGAATGGTGACTAAGTCATCAAGGTTATCCCCCACCGTATTGCGGCGAGATAATGGCAACCCCTGATGAAAGAAAAAATCGGCCAGCCAGGCAg                                 >  2:1124956/1‑150 (MQ=255)
                                                                                                                    tcggtcCCAGCCCGCCGTTAACGATCAGCACATCGGCGTGCTGACTACGTTCGCGAAGAATGGTGACTAAGTCATCAAGGTTATCCCCCACCGTATTGCGGCGAGATAATGGCAACCCCTGATGAAAGAAAAAATCGGCCAGCCAGGCAg                                 >  2:227570/1‑150 (MQ=255)
                                                                                                                      ggtcCCTGCCCGCCGTTAACGATCAGCACATCGGCGTGCTGACTAGGTTCGCTATTATTGGTGACTAAGTCATCAAGGTTATCCCCCACCTTTTTGCGGCGATTTAATGGCAACCCCTTATGAAATAAAAAATCGGCCAGCCAGGCAGCg                               <  1:897027/150‑1 (MQ=255)
                                                                                                                         cccAGCCCGCCGTTAACGATCAGCACATCGGCGTGCTGACTACGTTCGCGAAGAATGGTGACTAAGTCATCAAGGTTATCCCCCACCGTATTGCGGCGAGATAATGGCAACCCCTGATGAAAGAAAAAATCGGCCAGCCAGGCAGCGTTa                            >  2:731377/1‑150 (MQ=255)
                                                                                                                          ccAGCCCGCCGTTAACGATCAGCACATCGGCGTGCTGACTACGTTCGCGAAGAATGGTGACTAAGTCATCAAGGTTATCCCCCACCGTATTGCGGCGAGATAATGGCAACCCCTGATGAAAGAAAAAATCGGCCAGCCAGGCAGCGTTAg                           >  1:1236040/1‑150 (MQ=255)
                                                                                                                          ccAGCCCGCCGTTAACGATCAGCACATCGGCGTGCTGACTACGTTCGCGAAGAATGGTGACTAAGTCATCAAGGTTATCCCCCACCGTATAGCGGCGAGATAATGGCAACCCCTGATGAAAGAAAAAATCGGCCAGCCAGGCAGCGTTAg                           >  1:1236050/1‑150 (MQ=255)
                                                                                                                           cAGCCCGCCGTTAACGATCAGCACATCGGCGTGCTGACTACGTTCGCGAAGAATGGTGACTAAGTCATCAAGGTTATCCCCCACCGTATTGCGGCGAGATAATGGCAACCCCTGATGAAAGAAAAAATCGGCCAGCCAGGCAGCGTTAgt                          <  2:146587/150‑1 (MQ=255)
                                                                                                                              cccGCCGTTAACGATCAGCACATCGGCGTGCTGACTACGTTCGCGAAGAATGGTGACTAAGTCATCAAGGTTATCCCCCACCGTATTGCGGCGAGATAATGGCAACCCCTGATGAAAGACAAAATCGGCCAGCCAGGCAGCGTTAGTGTc                       >  2:729875/1‑150 (MQ=255)
                                                                                                                                     ttAACGATCAGCACATCGGCGTGCTGACTACGTTCGCGAAGAATGGTGACTAAGTCATCAAGGTTATCCCCCACCGTATTGCGGCGAGATAATGGCAACCCCTGATGAAAGAAAAAATCGGCCAGCCAGGCAGCGTTAGTGTCAACGagt                >  2:1060865/1‑148 (MQ=255)
                                                                                                                                       aaCGATCAGCACATCGGCGTGCTGACTACGTTCGCGAAGAATGGTGACTAAGTCATCAAGGTTATCCCCCACCGTATTGCGGCGAGATAATGGCAACCCCTGATGAAAGAAAAAATCGGCCAGCCAGGCAGCGTTAGTGTCAACGATTTg              >  1:748539/1‑150 (MQ=255)
                                                                                                                                         tagtCAGCACGCCGATGTGCTGACTACGTTCGCGAAGAATGGTGACTAAGTCATCAAGGTTATCCCCCACCGTATTGCGGCGAGATAATGGCAACCCCTGATGAAAGAAAAAATCGGCCAGCCAGGCAGCGTTAGTGTCAACGATTTGcc            >  2:597425/4‑150 (MQ=255)
                                                                                                                                           gtCAGCACGCCGATGTGCTGACTACGTTCGCGAAGAATGGTGACTAAGTCATCAAGGTTATCCCCCACCGTATTGCGGCGAGATAATGGCAACCCCTGATGAAAGAAAAAATCGGCCAGCCAGGCAGCGTTAGTGTCAACGATTTGCCCg          >  1:498339/2‑150 (MQ=255)
                                                                                                                                           aTCAGCACATCGGCGTGCTGACTACGTTCGCGAAGAATGGTGACTAAGTCATCAAGGTTATCCCCCACCTTATTGCGGCGAGATATTGGAAACCCATGATGAAAGAAAAAATCGGCCAGCCAGGCAGCGTTAGTGTCAACGATTTGCCCg          <  1:416598/150‑1 (MQ=255)
                                                                                                                                            tccgCACATGGGCGTGCTGAGTCGGTTCGCGACGAAGGGTGACTAAGTCATCTCGGTTATCCCCCACGGTATTGGGGCGTGTTAATGGCAACCCCTTATTAAAGAAAAAATCGGCCAGCCAGGCAGCGTTAGTGTCAACGATTTGCCCgt         <  1:551194/147‑1 (MQ=255)
                                                                                                                                               gCACATCGGCGTGCTGACTACGTTCGCGAAGAATGGTGACTAAGTCATCAAGGTTATCCCCCACCGTATTGCGGCGAGATAATGGCAACCCCTGATGAAAGAAAAAATCGGCCAGCCAGGCAGCGTTAGTGTCAACGATTTGCCCGTGTa      >  1:1375087/1‑150 (MQ=255)
                                                                                                                                                cacaTCGGCGTGCTGACTACGTTCGCGAAGAATGGTGACTAAGTCATCAAGGTTATCCCCCACCGTATTGCGGCGAGATAATGGCAACCCCTGATGAAAGAAAAAATCGGCCAGCCAGGCAGCGTTAGTGTCAACGATTTGCCCGTGTaa     >  1:1144033/1‑150 (MQ=255)
                                                                                                                                                  caTCGGCGTGCTGACTACGTTCGCGAAGAATGGTGACTAAGTCATCAAGGTTATCCCCCACCGTATTGCGGCGAGATAATGGCAACCCCTGATGAAAGAAAAAATCGGCCAGCCAGGCAGCGTTAGTGTCAACGATTTGCCCGTGTAaca   <  2:183249/150‑1 (MQ=255)
                                                                                                                                                   aTCGGCGTGCTGACTACGTTCGCGAAGAATGGTGACTAAGTCATCAAGGTTATCCCCCACCGTATTGCGGCGAGATAATGGCAACCCCTGATGAAAGAAAAAATCGGCCAGCCAGGCAGCGTTAGTGTCAACGATTTGCCCGTGTAacac  <  2:1291705/150‑1 (MQ=255)
                                                                                                                                                   aTCGGCGTGCTGACGACGTTCGCGAAGACTGGTGACTAAGTCATCAAGGTTATCCCCCACCGTATTGCGGCGAGATAATGGCTACCCCTGATGAAAGAAAAAATCGGCCAGCCAGGCAGCGTTAGTGTCAACGTTTTGCACGTGTAacac  <  2:1291684/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
CATTACTCGTCCACGTTCGTGGAAATAGCGTTCCATCTCTTTGAGCCAGGCTTCATGCAGCACCAGGCCTTCACCTTTTGCTGTCGCAGCGGCGAGTGCGCTTAAATCATCGCTGGTCGGTCCCAGCCCGCCGTTAACGATCAGCACATCGGCGTGCTGACTACGTTCGCGAAGAATGGTGACTAAGTCATCAAGGTTATCCCCCACCGTATTGCGGCGAGATAATGGCAACCCCTGATGAAAGAAAAAATCGGCCAGCCAGGCAGCGTTAGTGTCAACGATTTGCCCGTGTAACAC  >  NC_000913/2363301‑2363597

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: