Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 851,020 G→C 6.0% intergenic (‑600/‑205) pcs ← / → hemA phosphatidylcholine synthase/glutamyl‑tRNA reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947851,0200GC 6.0% 943.6 / 36.2 283intergenic (‑600/‑205)pcs/hemAphosphatidylcholine synthase/glutamyl‑tRNA reductase
Reads supporting (aligned to +/- strand):  ref base G (180/86);  new base C (14/3);  total (194/89)
Fisher's exact test for biased strand distribution p-value = 2.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.82e-01

Only 100 of 287 total aligned reads displayed.

GCTTCTAAGCGCGCGGTTGATCAGGCGCTGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCGGCAAGACCTGAAAACACCGCTAGGGCAGGCTACGAAAAGGCGATAGCCGCTGCGCTATCAGCACAATAGCGAACGGTGGTTGTCC  >  NC_002947/850872‑851169
                                                                                                                                                    |                                                                                                                                                     
gCTTCTAAGCGCGCGGTTGATCAGGCGCTGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGc                                                                                                                                                      >  2:2462654/1‑150 (MQ=255)
gCTTCTAAGCGCGCGGTTGATCAGGCGCTGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCGTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGc                                                                                                                                                      <  1:4693149/150‑1 (MQ=255)
cttttgatctggcttttgatcctgatcttgatcttcgcgacttcaggaggccgagcgtagggCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGc                                                                                                                                                      >  2:1920593‑M1/63‑150 (MQ=255)
  tttgatctggcttttgatcctgatcttgatcttcgcgacttcaggaggccgaacgtagggCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCt                                                                                                                                                    >  2:5190597‑M1/61‑150 (MQ=255)
  ttCTAAGCGCGCGGTTGATCAGGCGCTGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCt                                                                                                                                                    >  1:847066/1‑150 (MQ=255)
  ttCTAAGCGCGCGGTTGATCAGGCGCTGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCt                                                                                                                                                    >  2:2962588/1‑150 (MQ=255)
   tCTAAGCGCGCGGTTGATCAGGCGCTGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTg                                                                                                                                                   >  2:5496216/1‑150 (MQ=255)
    cTAAGCGCGCGGTTGATCAGGCGCTGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGc                                                                                                                                                  >  2:2061037/1‑150 (MQ=255)
     tAAGCGCGCGGTTGATCAGGCGCTGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCg                                                                                                                                                 >  1:4027271/1‑150 (MQ=255)
       aGCGCGCGGTTGATCAGGCGCTGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGtt                                                                                                                                               <  2:6094452/150‑1 (MQ=255)
        gcgcgcGGTTGATCAGGCGCTGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTc                                                                                                                                              <  1:1528938/150‑1 (MQ=255)
          gcgcgGTTGATCAGGCGCTGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCgg                                                                                                                                            >  2:6577499/1‑150 (MQ=255)
           cgcgGTTGATCAGGCGCTGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGc                                                                                                                                           <  1:4102020/150‑1 (MQ=255)
              ggTTGATCAGGCGCTGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGcctc                                                                                                                                        <  2:6007619/150‑1 (MQ=255)
                gtGATCAGGCGCTGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGcctcct                                                                                                                                      <  1:1715632/149‑1 (MQ=255)
                  gATCAGGCGCTGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGa                                                                                                                                    >  2:967318/1‑150 (MQ=255)
                    tCAGGCGCTGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAg                                                                                                                                  >  2:3552087/1‑150 (MQ=255)
                       ggCGCTGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTcg                                                                                                                               >  2:1497157/1‑150 (MQ=255)
                        gcgcTGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTcgc                                                                                                                              >  2:385980/1‑150 (MQ=255)
                          gcTGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGa                                                                                                                            <  2:252746/150‑1 (MQ=255)
                              cAGAGTGCGACTTCTGGAGGCCGCGCGTAGGGCTTGCTCCGGGGAGGACCGCGCTGCCGGCCCCATCCTCGGGCACAGCGCTGGACGGGCAGCCATGCCCGTCACCTCCCTGCGCACGGCCTGCGTGCGGCCTCCTGAAGTCGCGaagat                                                                                                                        <  2:3707992/150‑1 (MQ=21)
                                       aCTTCAGGAGGCCGAGCGTAGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAg                                                                                                               >  2:1038398/1‑150 (MQ=255)
                                        cTTCAGGAGGCCGAGCGTAGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTCGGGCACAGCGCTTGACCGGCATCCATGCCCGGCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAgg                                                                                                              >  1:6690210/1‑150 (MQ=255)
                                         ttCAGGAGGCCGAGCGTAGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGCAGTCGCGAAGATCAAGATCAGGa                                                                                                             >  2:6556284/1‑150 (MQ=255)
                                            aggaggccgaacgcaggccttgcgcagggaggtgacgggcatggatgcccgtcaagcgctgtgcccgaggatggggcctgcagcgcggtcctcccgggagcaagCCCTACGCTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCa                                                                                                          >  1:174159‑M2/105‑150 (MQ=255)
                                            aggaggCCGAGCGTAGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCa                                                                                                          >  2:6677909/1‑150 (MQ=255)
                                             ggaggCCGAGCGTAGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCaa                                                                                                         >  1:2838372/1‑150 (MQ=255)
                                               cggCCGCGCGTCGGGCGTGCTCCCGGGCGGACCGCGCTGCAGGCCCCATCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGGCACCTCCCTGCGCAAGGCCTGCGTCGGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCaaaa                                                                                                       <  2:4935582/149‑1 (MQ=18)
                                                   cGAGCGTAGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCAACCTCGGGCACCGCGCTTGCCGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCa                                                                                                   <  1:2106891/150‑1 (MQ=255)
                                                    gAGCGTAGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAg                                                                                                  >  1:5379786/1‑150 (MQ=255)
                                                     aacgcaggccttgcgcagggaggtgacgggcatggatgcccgtcaagcgctgtgcccgaggatggggcctgcagcgcggtcctcccgggagcaagCCCTACGCTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCa                                                                                                   <  1:4122892‑M2/53‑1 (MQ=255)
                                                      gCGTAGGGCTTGCCCCCGGGAGGACCGCGCTGCAGGCCCCTTCGTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGAt                                                                                                <  1:1987024/150‑1 (MQ=255)
                                                        gTAGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCTGGATCAAAAGCCAGATCa                                                                                              >  1:6635843/1‑150 (MQ=255)
                                                         tAGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCaa                                                                                             >  2:5276198/1‑150 (MQ=255)
                                                          aggccttgcgcagggaggtgacgggcatggatgcccgtcaagcgctgtgcccgaggatggggcctgcagcgcggtcctcccgggagcaagCCCTACGCTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCaaa                                                                                            >  2:6819417‑M2/91‑150 (MQ=255)
                                                           gggCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCaaaa                                                                                           >  1:1046535/1‑150 (MQ=255)
                                                           gggCTTGCTCCCGGGAGGACCGCGCGGCAGGCCCCATCCTCGGGCACAGCGCTTGACGGGCATCCAGGCCCGTCACCGCCGTGCGCATGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCaaaa                                                                                           <  2:5451046/150‑1 (MQ=18)
                                                             gCTTGCTCCCGGTAGGACCGCGCTGCAGGCCCCATCCTCGGGCTCAGCGCTTGACGGTCATCCATGTCCTTCACATCCCTGCGCAAGGCCTGCTTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGc                                                                                         <  2:1567526/150‑1 (MQ=17)
                                                               TTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCg                                                                                       >  2:4562686‑M1/1‑150 (MQ=255)
                                                                   tCCCGGGAGGACCGCGCTGCAGGCCCCATCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCGGCaa                                                                                   >  1:531414/1‑150 (MQ=255)
                                                                   gcagggaggtgacgggcatggatgcccgtcaagcgctgtgcccgaggatggggcctgcagcgcggtcctcccgggagcaagCCCTACGCTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCGGCaa                                                                                   >  2:5122178‑M2/82‑150 (MQ=255)
                                                                   gcagggaggtgacgggcatggatgcccgtcaagcgctgtgcccgaggatggggcctgcagcgcggtcctcccgggagcaagCCCTACGCTCGGCCTCCTGAAGTCGCAAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCGGCaa                                                                                   >  1:3467961‑M2/82‑150 (MQ=255)
                                                                    cccGGGAGGACCGCGCTGCAGGCCCCATCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCGGCAAg                                                                                  <  1:1907730/150‑1 (MQ=255)
                                                                     gggggAGGACCGCGCTGCAGGCCCCATCCTCGGGCACAGCGGTTGCCGGGCACCCATGCCCGTAACGTCCGTGCGCAAGGCCGGCGTTCGGCCTCCTTATGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCGGCAAGa                                                                                 <  1:1232045/148‑1 (MQ=21)
                                                                     cggggAGGACCGCGCTGCCGGCCCCATCCTGGGGCACTGCGGTTGACGGGCATCCATGCCGGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCGGCAAGa                                                                                 <  1:6624021/148‑1 (MQ=18)
                                                                       gggAGGACCGCGCTGCAGGCCCCATCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCGGCAAGAcc                                                                               <  1:5519813/150‑1 (MQ=255)
                                                                        gaggtgacgggcatggatgcccgtcaagcgctgtgcccgaggatggggcctgcagcgcggtcctcccgggagcaagCCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCGGCAAGACCt                                                                              >  1:3598479‑M2/77‑150 (MQ=255)
                                                                          ggtgacgggcatggatgcccgtcaagcgctgtgcccgaggatggggcctgcagcgcggtcctcccgggagcaagGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCGGCAAGACCTGa                                                                            >  1:2281750‑M2/75‑150 (MQ=255)
                                                                                gcgcTGCCGGCCCCATCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCATGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCGGCAAGACCTGAAAacac                                                                      <  1:1327494/150‑1 (MQ=255)
                                                                                    tGCAGGCCCCATCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCGGCAAGACCTGAAAACACCGCt                                                                  >  2:1008658/1‑150 (MQ=255)
                                                                                    tGCAGGCCCCATCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCGGCAAGACCTGAAAACACCGCt                                                                  >  2:6343939/1‑150 (MQ=255)
                                                                                      ggatgcccgtcaagcgctgtgcccgaggatggggcctgcagcgcggtcctcccgggagcaagCCCTACGCTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCGGCAAGACCTGAAAACACCGCTAg                                                                >  1:6308355‑M2/63‑150 (MQ=255)
                                                                                      ggatgcccgtcaagcgctgtgcccgaggatggggcctgcagcgcggtcctcccgggagcaagCCCTACGCTCGGCCTCCTGAAGTCGCAAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCGGCAAGACCTGAAAACACCGCTAg                                                                >  1:2042657‑M2/63‑150 (MQ=255)
                                                                                        ggCCCCATCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCGGCAAGACCTGAAAACACCGCTAggg                                                              >  2:6205387/1‑150 (MQ=255)
                                                                                         gCCCCATCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCGGCAAGACCTGAAAACACCGCTAGGGc                                                             >  1:5810214/1‑150 (MQ=255)
                                                                                          ccccATCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCGGCAAGACCTGAAAACACCGCTAGGGCa                                                            <  2:4781131/150‑1 (MQ=255)
                                                                                            ccATCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGGTCAAGATCAGGATCAAAAGCCAGATCAAAAGCCGGCAAGACCTGAAAACACCGCTAGGGCAgg                                                          >  1:6471168/1‑150 (MQ=255)
                                                                                            ccATCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCGGCAAGACCTGAAAACACCGCTAGGGCAgg                                                          >  1:5791009/1‑150 (MQ=255)
                                                                                            ccATCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCGGCAAGACCTGAAAACACCGCTAGGGCAgg                                                          >  2:4444225/1‑150 (MQ=255)
                                                                                            ccATCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCGGCAAGACCTGAAAACACCGCTAGGGCAgg                                                          >  2:407109/1‑150 (MQ=255)
                                                                                                 cTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCGGCAAGACCTGAAAACACCGCTAGGGCAGGCTACg                                                     >  1:3151991/1‑150 (MQ=255)
                                                                                                    gggCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCGGCAAGACCTGAAAACACCGCTAGGGCAGGCTACGaaa                                                  >  1:1472797/1‑150 (MQ=255)
                                                                                                      gCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGa                                                                                                             >  2:712121/1‑89 (MQ=255)
                                                                                                          aGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCGGCAAGACCTGAAAACACCGCTAGGGCAGGCTACGAAAAGGCGa                                            >  2:1138065/1‑150 (MQ=255)
                                                                                                          aGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCGGCAAGACCTGAAAACACCGCTAGGGCAGGCTACGAAAAGGCGa                                            >  1:4350075/1‑150 (MQ=255)
                                                                                                          aGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCGGCAAGACCTGAAAACACCGCTAGGGCAGGCTACGAAAAGGCGa                                            >  1:2671594/1‑150 (MQ=255)
                                                                                                          aGCGCTGGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCGGCAAGACCTGAAAACACCGCTAGGGCAGGCTACGAAAAGGCGa                                            <  2:4325912/150‑1 (MQ=255)
                                                                                                               ttGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCGGCAAGACCTGAAAACACCGCTAGGGCAGGCTACGAAAAGGCGATAGcc                                       >  1:2186159/1‑150 (MQ=255)
                                                                                                               ttGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCGGCAAGACCTGAAAACACCGCTAGGGCAGGCTACGAAAAGGCGATAGcc                                       >  2:1783664/1‑150 (MQ=255)
                                                                                                               gtGACGGGCATGGATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCGGCAAGACCTGAAAACACCGCTAGGGCAGGCTACGAAAAGGCGATAGcc                                       >  2:2621482/2‑150 (MQ=255)
                                                                                                                tGACGGGCATGGATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCGGCAAGACCTGAAAACACCGCTAGGGCAGGCTACGAAAAGGCGATAGCCg                                      >  2:5210241/1‑150 (MQ=255)
                                                                                                                    gggCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCGGCAAGACCTGAAAACACCGCTAGGGCAGGCTACGAAAAGGCGATAGCCGCTgc                                  >  2:5906805/1‑150 (MQ=255)
                                                                                                                    gggCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCGGCAAGACCTGAAAACACCGCTAGGGCAGGCTACGAAAAGGCGATAGCCGCTgc                                  >  2:5766627/1‑150 (MQ=255)
                                                                                                                    gggCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCGGCAAGACCTGAAAACACCGCTAGGGCAGGCTACGAAAAGGCGATAGCCGCTgc                                  >  1:3860979/1‑150 (MQ=255)
                                                                                                                         cccATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCGGCAAGACCTGAAAACACCGCTAGGGCAGGCTACGAAAAGGCGATAGCCGCTGCGCTAt                             >  2:4653087/2‑150 (MQ=255)
                                                                                                                          ccATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCGGCAAGACCTGAAAACACCGCTAGGGCAGGCTACGAAAAGGCGATAGCCGCTGCGCTATc                            >  1:6009070/1‑150 (MQ=255)
                                                                                                                           cATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCGGCAAGACCTGAAAACACCGCTAGGGCAGGCTACGAAAAGGCGATAGCCGCTGCGCTATCa                           >  2:169490/1‑150 (MQ=255)
                                                                                                                           cATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCGGCAAGACCTGAAAACACCGCTAGGGCAGGCTACGAAAAGGCGATAGCCGCTGCGCTATCa                           <  1:482786/150‑1 (MQ=255)
                                                                                                                             tGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCGGCAAGACCTGAAAACACCGCTAGGGCAGGCTACGAAAAGGCGATAGCCGCTGCGCTATCAGc                         <  2:3133649/150‑1 (MQ=255)
                                                                                                                              gCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCGGCAAGACCTGAAAACACCGCTAGGGCAGGCTACGAAAAGGCGATAGCCGCTGCGCTATCAGca                        <  1:4245977/150‑1 (MQ=255)
                                                                                                                               cccGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCGGCAAGACCTGAAAACACCGCTAGGGCAGGCTACGAAAAGGCGATAGCCGCTGCGCTATCAGcac                       >  2:4601241/1‑150 (MQ=255)
                                                                                                                                 cGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCGGCAAGACCTGAAAACACCGCTAGGGCAGGCTACGAAAAGGCGATAGCCGCTGCGCTATCAGCACaa                     <  2:5605231/150‑1 (MQ=255)
                                                                                                                                   tCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCGGCAAGACCTGAAAACACCGCTAGGGCAGGCTACGAAAAGGCGATAGCCGCTGCGCTATCAGCACAATa                   >  2:332778/1‑150 (MQ=255)
                                                                                                                                    tcctcccgggagcaagGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCGGCAAGACCTGAAAACACCGCTAGGGCAGGCTACGAAAAGGCGATAGCCGCTGCGCTATCAGCACAATAg                  >  1:710316‑M2/17‑150 (MQ=255)
                                                                                                                                     aCCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCGGCAAGACCTGAAAACACCGCTAGGGCAGGCTACGAAAAGGCGATAGCCGCTGCGCTATCAGCACAATAGc                 >  1:5045462/1‑150 (MQ=255)
                                                                                                                                     aCCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCGGCAAGACCTGAAAACACCGCTAGGGCAGGCTACGAAAAGGCGATAGCCGCTGCGCTATCAGCACAATAGc                 >  1:2389087/1‑150 (MQ=255)
                                                                                                                                       cTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCGGCAAGACCTGAAAACACCGCTAGGGCAGGCTACGAAAAGGCGATAGCCGCTGCGCTATCAGCACAATAGCGa               >  1:5183586/1‑150 (MQ=255)
                                                                                                                                       cTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCGGCAAGACCTGAAAACACCGCTAGGGCAGGCTACGAAAAGGCGATAGCCGCTGCGCTATCAGCACAATAGCGa               >  2:1973093/1‑150 (MQ=255)
                                                                                                                                        tCCCTGCGCACGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCCAAAGCCAGATCAAAAGCCGGCCAGACCTGAAAACACCGCTAGGGCTGGCTACGACATGGCGATAGCCGCTGCGCTATCAGCACACTTGCGac              >  1:5940033/1‑149 (MQ=37)
                                                                                                                                         cccTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCGGCAAGACCTGAAAACACCGCTAGGGCAGGCTACGAAAAGGCGATAGCCGCTGCGCTATCAGCACAATAGCGAAc             >  2:4180591/1‑150 (MQ=255)
                                                                                                                                         cccTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCGGCAAGACCTGAAAACACCGCTAGGGCAGGCTACGAAAAGGCGATAGCCGCTGCGCTATCAGCACAATAGCGAAc             >  1:6107497/1‑150 (MQ=255)
                                                                                                                                          ccTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCGGCAAGACCTGAAAACACCGCTAGGGCAGGCTACGAAAAGGCGATAGCCGCTGCGCTATCAGCACAATAGCGAACg            >  1:446023/1‑150 (MQ=255)
                                                                                                                                          ccTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCGGCAAGACCTGAAAACACCGCTAGGGCAGGCTACGAAAAGGCGATAGCCGCTGCGCTATCAGCACAATAGCGAACg            >  2:4259846/1‑150 (MQ=255)
                                                                                                                                          ccTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCGGCAAGACCTGAAAACACCGCTAGGGCAGGCTACGAAAAGGCGATAGCCGCTGCGCTATCAGCACAATAGCGAACg            >  2:794995/1‑150 (MQ=255)
                                                                                                                                            ggagcaagCCCTACGCTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCGGCAAGACCTGAAAACACCGCTAGGGCAGGCTACGAAAAGGCGATAGCCGCTGCGCTATCAGCACAATAGCGAACggt          >  1:5102521‑M2/9‑150 (MQ=255)
                                                                                                                                             gcgcAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCGGCAAGACCTGAAAACACCGCTAGGGCAGGCTACGAAAAGGCGATAGCCGCTGCGCTATCAGCACAATAGCGAACggtg         >  2:3478692/1‑150 (MQ=255)
                                                                                                                                             gcgcAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCGGCAAGACCTGAAAACACCGCTAGGGCAGGCTACGAAAAGGCGATAGCCGCTGCGCTATCAGCACAATAGCGAACggtg         >  1:605000/1‑150 (MQ=255)
                                                                                                                                             gcgaaagGCCTGCGTTCGGCCTCCTTAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCGGCAAGACCTGAAAACACCGCTAGGGCAGGCTACGAAAAGGCGATAGCCGCTGCGCTATCAGCACAATAGCGAACggtg         <  1:2462654‑M2/143‑1 (MQ=255)
                                                                                                                                              cgcAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCGGCAAGACCTGAAAACACCGCTAGGGCAGGCTACGAAAAGGCGATAGCCGCTGCGCTATCAGCACAATAGCGAACggtgg        <  2:2343865/150‑1 (MQ=255)
                                                                                                                                               gcaagGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCGGCAAGACCTGAAAACACCGCTAGGGCAGGCTACGAAAAGGCGATAGCCGCTGCGCTATCAGCACAATAGCGAACggtggt       >  2:2219621‑M2/6‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GCTTCTAAGCGCGCGGTTGATCAGGCGCTGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCGGCAAGACCTGAAAACACCGCTAGGGCAGGCTACGAAAAGGCGATAGCCGCTGCGCTATCAGCACAATAGCGAACGGTGGTTGTCC  >  NC_002947/850872‑851169

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