Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,865,807 G→C 5.7% L79L (CTC→CTG)  ligA ← DNA ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,865,8070GC 5.7% 1415.8 / 51.2 417L79L (CTC→CTG) ligADNA ligase
Reads supporting (aligned to +/- strand):  ref base G (201/192);  new base C (16/8);  total (217/200)
Fisher's exact test for biased strand distribution p-value = 1.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 419 total aligned reads displayed.

GCTGGCCGTCGCGGTACAGCAGGCTCACGGCCAGGCCATCGAGCTTGGGTTCGCAGCTGTAGTCGACCGCGCCCGGCTGATCGAGGCCTTCTACCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTCGGGGGTTACCAGGTGCGGGTGTTCAGCTTCCAGCGCCTTGAGCTCGTTGAACAGGCGGTCGTATTCGGCATCAG  >  NC_002947/4865659‑4865955
                                                                                                                                                    |                                                                                                                                                    
  tGGCCGTCGCGGTACAGCGGGCTCACGGCCAGGCCATCGAGCTTGGGTTCGCAGCTGTAGTCGACCGCGCCCGGCTGATCGAGGCCTTCTACCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGc                                                                                                                                                   <  1:1656663/150‑1 (MQ=255)
     ccGTCGCGGTACAGCAGGCTCACGGCCAGGCCATCGAGCTTGGGTTCGCAGCTGTAGTCGACCGCGCCCGGCTGATCGAGGCCTTCTACCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATg                                                                                                                                                <  2:5802090/150‑1 (MQ=255)
     ccGTCGCGGTACAGCAGGCTCACGGCCAGGCCATCGAGCTTGGGTTCGCAGCTGTAGTCGACCGCGCCCGGCTGATCGAGGCCTTCTACCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATg                                                                                                                                                <  1:5628741/150‑1 (MQ=255)
      cGTCGCGGTACAGCAGGCTCACGGCCAGGCCATCGAGCTTGGGTTCGCAGCTGTAGTCTACCGCGCCCGGCTGATCGAGGCCTTCTACCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATgg                                                                                                                                               <  2:7254133/150‑1 (MQ=255)
          gcgGTACAGCAGGCTCACGGCCAGGCCATCGAGCTTGGGTTCGCAGCTGTAGTCGACCGCGCCCGGCTGATCGAGGCCTTCTACCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCATGCTCAGCCTGGGCaa                                                                                                                                           <  2:87916/150‑2 (MQ=255)
                        tCACGGCCAGGCCATCGAGCTTGGGTTCGCAGCTGTAGTCGACCGCGCCCGGCTGATCGAGGCCTTCTACCACGCGACGGCCGAACTCGCGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCt                                                                                                                             >  2:7029855/1‑150 (MQ=255)
                        tCACGGCCAGGCCATCGAGCTTGGGTTCGCAGCTGTAGTCGACCGCGCCCGGCTGATCGAGGCCTTCTACCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCt                                                                                                                             >  1:3985858/1‑150 (MQ=255)
                          aCGGCCAGGCCATCGAGCTTGGGTTCGCAGCTGTAGTCGACCGTGCCCGGCTGATCGAGGCCTTCTACCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTgg                                                                                                                           >  2:2921688/1‑150 (MQ=255)
                            ggccaggccaTCGAGCTTGGGTTCGCAGCTGTAGTCGACCGCGCCCGGCTGATCGAGGCCTTCTACCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCt                                                                                                                         <  2:7837571/150‑1 (MQ=255)
                               caggccaTCGAGCTTGGGTTCGCAGCTGTAGTCGACCGCGCCCGGCTGATCGAGGCCTTCTACCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGaa                                                                                                                      >  2:4883397/1‑150 (MQ=255)
                               caggccaTCGAGCTTGGGTTCGCAGCTGTAGTCGACCGCGCCCGGCTGATCGAGGCCTTCTACCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGaa                                                                                                                      >  1:6546103/1‑150 (MQ=255)
                                 ggccaTCGAGCTTGGGTTCGCAGCTGTAGTCGACCGCGCCCGGCTGATCGAGGCCTTCTACCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGACcg                                                                                                                    <  2:751991/150‑1 (MQ=255)
                                  gccaTCGAGCTTGGGTTCGCAGCTGTAGTCGACCGCGCCCGGCTGATCGAGGCCTTCTACCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAAcgc                                                                                                                   <  2:4203569/150‑1 (MQ=255)
                                    caTCGAGCTTGGGTTCGCAGCTGTAGTCGACCGCGCCCGGCTGATCGAGGCCTTCTACCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAAcgccg                                                                                                                 >  2:5263874/1‑150 (MQ=255)
                                      tCGAGCTTGGGTTCGCAGCTGTAGTCGACCGCGCCCGGCTGATCGAGGCCTTCTACCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACgccgcc                                                                                                               >  2:4447213/1‑150 (MQ=255)
                                        gAGCTTGGGTTCGCAGCTGTAGTCGACCGCGCCCGGCTGATCGAGGCCTTCTACCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCaa                                                                                                             <  2:6760244/150‑1 (MQ=255)
                                        gAGCTTGGGTTCGCAGCTGTAGTCGACCGCGCCCGGCTGATCGAGGCCTTCTACCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCaa                                                                                                             <  1:7187058/150‑1 (MQ=255)
                                         aGCTTGGGTTCGCAGCTGTAGTCGACCGCGCCCGGCTGATCGAGGCCTTCTACCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAg                                                                                                            <  2:4048000/150‑1 (MQ=255)
                                          gCTTGGGTTCGCAGCTGTAGTCGACCGCGCCCGGCTGATCGAGGCCTTCTACCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAgg                                                                                                           >  1:5417343/1‑150 (MQ=255)
                                            ttGGGTTCGCAGCTGTAGTCGACCGCGCCCGGCTGATCGAGGCCTTCTACCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCa                                                                                                         >  2:7584670/1‑150 (MQ=255)
                                               ggTTCGCAGCTGTAGTCGACCGCGCCCGGCTGATCGAGGCCTTCTACCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCt                                                                                                      <  1:6024774/150‑1 (MQ=255)
                                                gTTCGCAGCTGTAGTCGACCGCGCCCGGCTGATCGAGGCCTTCTACCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTc                                                                                                     >  1:442710/1‑150 (MQ=255)
                                                    gCAGCTGTAGTCGACCGCGCCCGGCTGATCGAGGCCTTCTACCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTccacc                                                                                                 >  1:362525/1‑150 (MQ=255)
                                                    gCAGCTGTAGTCGACCGCGCCCGGCTGATCGAGGCCTTCTACCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTccacc                                                                                                 >  1:4606289/1‑150 (MQ=255)
                                                      aGCTGTAGTCGACCGCGCCCGGCTGATCGAGGCCTTCTACCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTccaccca                                                                                               >  1:4105746/1‑150 (MQ=255)
                                                         tGTAGTCGACCGCGCCCGGCTGATCGAGGCCTTCTACCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCacac                                                                                            >  1:2090483/1‑150 (MQ=255)
                                                           tAGTCGACCGCGCCCGGCTGATCGAGGCCTTCTACCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGc                                                                                          >  1:6132957/1‑150 (MQ=255)
                                                             gTCGACCGCGCCCGGCTGATCGAGGCCTTCTACCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCtg                                                                                        >  2:3414541/1‑150 (MQ=255)
                                                                  ccGCGCCCGGCTGATCGAGGCCTTCTACCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCtgggtgg                                                                                   <  1:1841211/150‑1 (MQ=255)
                                                                   cgcgCCCGGCTGATCGAGGCCTTCTACCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCtgggtggg                                                                                  <  1:6743708/150‑1 (MQ=255)
                                                                      gcCCGGCTGCTCGAGGCCTTCGCCCCCGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTGCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACCCGCTGGGTGGGCGa                                                                               <  2:7409319/150‑1 (MQ=255)
                                                                      gcCCGGCTGATCGAGGCCTTCTACCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGa                                                                               >  2:6017064/1‑150 (MQ=255)
                                                                       cccGGCTGATCGAGGCCTTCTACCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAg                                                                              >  1:930118/1‑150 (MQ=255)
                                                                       cccGGCTGATCGAGGCCTTCTACCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAg                                                                              >  1:6954195/1‑150 (MQ=255)
                                                                        ccGGCTGATCGAGGCCTTCTACCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGt                                                                             >  2:3104419/1‑150 (MQ=255)
                                                                         cGGCTGATCGAGGCCTTCTACCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTc                                                                            >  2:4644637/1‑150 (MQ=255)
                                                                           gCTGATCGAGGCCTTCTACCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTCgg                                                                          >  2:7350678/1‑150 (MQ=255)
                                                                            cTGATCGAGGCCTTCTACCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCAGCCAAGGCAGATCCACCCACACGCTGGGTGGGCGAGTCggg                                                                         <  2:1743333/150‑1 (MQ=255)
                                                                            cTGATCGAGGCCTTCGACCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTCggg                                                                         <  2:773242/150‑1 (MQ=255)
                                                                             tGATCGAGGCCTTCTACCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTCgggg                                                                        <  1:6366587/150‑1 (MQ=255)
                                                                              gATCGAGGCCTTCTACCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTCggggg                                                                       <  1:5216076/150‑1 (MQ=255)
                                                                              gATCGAGGCCTTCTACCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCATTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTCggggg                                                                       >  2:8385230/1‑150 (MQ=255)
                                                                               aTCGAGGCCTTCTACCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTCGGGGGt                                                                      >  1:3889466/1‑150 (MQ=255)
                                                                                tCGAGGCCTTCTACCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTCGGGGGtt                                                                     <  2:6898305/150‑1 (MQ=255)
                                                                                  gAGGCCTTCTACCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTCGGGGGTTAc                                                                   >  2:1350817/1‑150 (MQ=255)
                                                                                     gCCTTCTACCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTATCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTCGGGGGTTACCAg                                                                <  2:6137896/150‑1 (MQ=255)
                                                                                       cTTCTACCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTCGGGGGTTACCAGGt                                                              >  2:2047426/1‑150 (MQ=255)
                                                                                        ttCTACCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTCGGGGGTTACCAGGTg                                                             >  1:7797600/1‑150 (MQ=255)
                                                                                        ttCTACCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTCGGGGGTTACCAGGTg                                                             >  2:5866569/1‑150 (MQ=255)
                                                                                           tACCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTCGGGGGTTACCAGGTGCgg                                                          <  2:5141543/150‑1 (MQ=255)
                                                                                            aCCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTCGGGGGTTACCAGGTGCggg                                                         >  1:8383543/1‑150 (MQ=255)
                                                                                               aCGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTCGGGGGTTACCAGGTGCGGgtgt                                                      <  1:2724635/150‑1 (MQ=255)
                                                                                                    aCGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGTGCGAGTCGGGGGTTACCAGGTGCGGGTGTTCAGc                                                 <  2:4699563/150‑1 (MQ=255)
                                                                                                       gCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTCGGGGGTTACCAGGTGCGGGTGTTCAGCTTc                                              >  2:8639122/1‑150 (MQ=255)
                                                                                                       gCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTCGGGGGTTACCAGGTGCGGGTGTTCAGCTTc                                              >  1:6511373/1‑150 (MQ=255)
                                                                                                        ccGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTCGGGGGTTACCAGGTGCGGGTGTTCAGCTTcc                                             >  2:2820699/1‑150 (MQ=255)
                                                                                                        ccGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTCGGGGGTTACCAGGTGCGGGTGTTCAGCTTcc                                             >  2:3835090/1‑150 (MQ=255)
                                                                                                         cGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTCGGGGGTTACCAGGTGCGGGTGTTCAGCTTCCa                                            >  1:29573/1‑150 (MQ=255)
                                                                                                          gAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTCGGGGGTTACCAGGTGCGGGTGTTCAGCTTCCAg                                           >  2:499013/1‑150 (MQ=255)
                                                                                                          gAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTCGGGGGTTACCAGGTGCGGGTGTTCAGCTTCCAg                                           >  1:4841590/1‑150 (MQ=255)
                                                                                                              tCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTCGGGGGTTACCAGGTGCGGGTGTTCAGCTTCCa                                            >  2:1180389/1‑145 (MQ=255)
                                                                                                              tCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTCGGGGGTTACCAGGTGCGGGTGTTCAGCTTCCa                                            <  1:1180389/145‑1 (MQ=255)
                                                                                                              tCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTCGGGGGTTACCAGGTGCGGGTGTTCAGCTTCCAGCGcc                                       >  2:2634573/1‑150 (MQ=255)
                                                                                                              tCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTCGGGGGTTACCAGGTGCGGGTGTTCAGCTTCCAGCGcc                                       <  1:2300130/150‑1 (MQ=255)
                                                                                                                aCGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTCGGGGGTTACCAGGTGCGGGTGTTCAGCTTCCAGCGCCtt                                     >  2:6807753/1‑150 (MQ=255)
                                                                                                                aCGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTCGGGGGTTACCAGGTGCGGGTGTTCAGCTTCCAGCGCCtt                                     >  2:939176/1‑150 (MQ=255)
                                                                                                                aCGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTCGGGGGTTACCAGGTGCGGGTGTTCAGCTTCCAGCGCCtt                                     >  2:3247700/1‑150 (MQ=255)
                                                                                                                aCGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTCGGGGGTTACCAGGTGCGGGTGTTCAGCTTCCAGCGCCtt                                     >  1:7097021/1‑150 (MQ=255)
                                                                                                                 cGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTCGGGGGTTACCAGGTGCGGGTGTTCAGCTTCCAGCGCCTTg                                    >  1:6227953/1‑150 (MQ=255)
                                                                                                                  gCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTCGGGGGTTACCAGGTGCGGGTGTTCAGCTTCCAGCGCCTTGa                                   <  2:8281740/150‑1 (MQ=255)
                                                                                                                  gCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTCGGGGGTTACCAGGTGCGGGTGTTCAGCTTCCAGCGCCTTGa                                   <  2:802531/150‑1 (MQ=255)
                                                                                                                   cAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTCGGGGGTTACCAGGTGCGGGTGTTCAGCTTCCAGCGCCTTGAg                                  <  2:2090483/150‑1 (MQ=255)
                                                                                                                    aGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTCGGGGGTTACCAGGTGCGGGTGTTCAGCTTCCAGCGCCTTGAGc                                 >  2:2339905/1‑150 (MQ=255)
                                                                                                                     ggTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTCGGGGGTTACCAGGTGCGGGTGTTCAGCTTCCAGCGCCTTGAGCt                                >  2:4251786/1‑150 (MQ=255)
                                                                                                                       tCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTCGGGGGTTACCAGGTGCGGGTGTTCAGCTTCCAGCGCCTTGAGCTCg                              >  2:2196426/1‑150 (MQ=255)
                                                                                                                         ggCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTCGGGGGTTACCAGGTGCGGGTGTTCAGCTTCCAGCGCCTTGAGCTCGtt                            <  1:3510969/150‑1 (MQ=255)
                                                                                                                         ggCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTCGGGGGTTACCAGGTGCGGGTGTTCAGCTTCCAGCGCCTTGAGCTCGtt                            <  2:5677093/150‑1 (MQ=255)
                                                                                                                         ggCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTCGGGGGTTACCAGGTGCGGGTGTTCAGCTTCCAGCGCCTTGAGCTCGtt                            >  1:8605897/1‑150 (MQ=255)
                                                                                                                            ttcttcGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTCGGGGGTTACCAGGTGCGGGTGTTCAGCTTCCAGCGCCTTGAGCTCGTTGaa                         >  1:103686/1‑150 (MQ=255)
                                                                                                                            ttcttcGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTCGGGGGTTACCAGGTGCGGGTGTTCAGCTTCCAGCGCCTTGAGCTCGTTGaa                         >  1:3387592/1‑150 (MQ=255)
                                                                                                                               ttcGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTCGGGGGTTACCAGGTGCGGGTGTTCAGCTTCCAGCGCCTTGAGCTCGTTGAACAg                      >  1:2040476/1‑150 (MQ=255)
                                                                                                                               ttcGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTCGGGGGTTACCAGGTGCGGGTGTTCAGCTTCCAGCGCCTTGAGCTCGTTGAACAg                      >  1:4183991/1‑150 (MQ=255)
                                                                                                                                  gAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTCGGGGGTTACCAGGTGCGGGTGTTCAGCTTCCAGCGCCTTGAGCTCGTTGAACAGGCg                   >  2:1125772/1‑150 (MQ=255)
                                                                                                                                  gAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTCGGGGGTTACCAGGTGCGGGTGTTCAGCTTCCAGCGCCTTGAGCTCGTTGAACAGGCg                   >  2:8110032/1‑150 (MQ=255)
                                                                                                                                   cacgaagtGCCCATGCTCAGCCTGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTCGGGGGTTACCAGGTGCGGGTGTTCAGCGTCCAGCGCCTTGAGCTCGTTGAACAGGCgg                  >  2:3694934/8‑150 (MQ=255)
                                                                                                                                   aaGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTCGGGGGTTACCAGGTGCGGGTGTTCAGCTTCCAGCGCCTTGAGCTCGTTGAACAGGCgg                  >  1:7469000/1‑150 (MQ=255)
                                                                                                                                    acgaagtGCCCATGCTCAGCCTGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTCGGGGGTTACCAGGTGCGGGTGTTCAGCTTCCAGCGCCTTGAGCTCGTTGAACAGGCGGt                 >  1:3075958/7‑150 (MQ=255)
                                                                                                                                     ggCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTCGGGGGTTACCAGGTGCGGGTGTTCAGCTTCCAGCGCCTTGAGCTCGTTGAACAGGCGGTc                <  1:4950468/150‑1 (MQ=255)
                                                                                                                                     ggCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTCGGGGGTTACCAGGTGCGGGTGTTCAGCTTCCAGCGCCTTGAGCTCGTTGAACAGGCGGTc                <  2:5387061/150‑1 (MQ=255)
                                                                                                                                      gCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTCGGGGGTTACCAGGTGCGGGTGTTCAGCTTCCAGCGCCTTGAGCTCGTTGAACAGGCGGTCg               <  2:7162176/150‑1 (MQ=255)
                                                                                                                                         ttGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTCGGGGGTTACCAGGTGCGGGTGTTCAGCTTCCAGCGCCTTGAGCTCGTTGAACAGGCGGTCGTAt            <  1:7952203/150‑1 (MQ=255)
                                                                                                                                           gCCCAGGCTGCGCATGGGCACTGCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTCGGGGGTTACCAGGTGCGGGTGTTCAGCTTCCAGCGCCTTGAGCTCGTTGAACAGGCGGTCGTATTc          <  2:5655584/150‑1 (MQ=255)
                                                                                                                                           gCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTCGGGGGTTACCAGGTGCGGGTGTTCAGCTTCCAGCGCCTTGAGCTCGTTGAACAGGCGGTCGTATTc          <  2:6479597/150‑1 (MQ=255)
                                                                                                                                             ccATGCTCAGCCTGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTCGGGGGTTACCAGGTGCGGGTGTTCAGCTTCCAGCGCCTTGAGCTCGTTGAACAGGCGGTCGTATTCgg        >  1:7288133/1‑150 (MQ=255)
                                                                                                                                             ccAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTCGGGGGTTACCAGGTGCGGGTGTTCAGCTTCCAGCGCCTTGAGCTCGTTGAACAGGCGGTCGTATTCgg        <  1:3252255/150‑1 (MQ=255)
                                                                                                                                                  cTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTCGGGGGTTACCAGGTGCGGGTGTTCAGCTTCCAGCGCCTTGAGCTCGTTGAACAGGCGGTCGTATTCGGCATCa   <  2:3641714/150‑1 (MQ=255)
                                                                                                                                                  cTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTCGGGGGTTACCAGGTGCGGGTGTTCAGCTTCCAGCGCCTTGAGCTCGTTGAACAGGCGGTCGTATTCGGCATCa   <  1:1544901/150‑1 (MQ=255)
                                                                                                                                                  cTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTCGGGGGTTACCAGGTGCGGGTGTTCAGCTTCCAGCGCCTTGAGCTCGTTGAACAGGCGGTCGTATTCGGCATCa   >  2:3943446/1‑150 (MQ=255)
                                                                                                                                                   tGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTCGGGGGTTACCAGGTGCGGGTGTTCAGCTTCCAGCGCCTTGAGCTCGTTGAACAGGCGGTCGTATTCGGCATCAg  <  1:4533941/150‑1 (MQ=255)
                                                                                                                                                   tGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCACCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTCGGGGGTTACCAGGTGCGGGTGTTCAGCTTCCAGCGCCTTGAGCTCGTTGAACAGGCGGTCGTATTCGGCATCAg  <  1:5893946/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GCTGGCCGTCGCGGTACAGCAGGCTCACGGCCAGGCCATCGAGCTTGGGTTCGCAGCTGTAGTCGACCGCGCCCGGCTGATCGAGGCCTTCTACCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTCGGGGGTTACCAGGTGCGGGTGTTCAGCTTCCAGCGCCTTGAGCTCGTTGAACAGGCGGTCGTATTCGGCATCAG  >  NC_002947/4865659‑4865955

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