Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 88,649 T→A 5.6% A223A (GCA→GCT)  betC ← choline‑sulfatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00294788,6490TA 5.6% 1211.4 / 39.4 356A223A (GCA→GCT) betCcholine‑sulfatase
Reads supporting (aligned to +/- strand):  ref base T (134/201);  new base A (9/11);  total (143/213)
Fisher's exact test for biased strand distribution p-value = 6.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.74e-01

Only 100 of 363 total aligned reads displayed.

AGGGTCTGCAGCAATTGGCCGATGTTGTCATCGATATAGCTGCACGCGCCGAAGTAGGCGCGGCGGGCATCGCGAACCTTATCCACAGGCAGTGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTGGCTGAACTCGGCACGGGGCATGGGGATATCCACACCTTCGTAGCGGTCCCAGTAGCGTTTGGCAATGGTGTAGGGGTCATGCGGGTGGGTCATCGAAACGGTCAGGCAAAATGGCCGGCCATCGTTTTCGCGCACATGGTCGTA  >  NC_002947/88500‑88798
                                                                                                                                                     |                                                                                                                                                     
aGGGTCTGCAGCAATTGGCCGATGTTGTCATCGATATAGCTGCACGCGCCGAAGTAGGCGCGGCGGGCATCGCGAACCTTATCCACAGGCAGTGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCt                                                                                                                                                       <  1:3307948/150‑1 (MQ=255)
 gggTCTGCAGCAATTGGCCGATGTTGTCATCGATATAGCTGCACGCGCCGAAGTAGGCGCGGCGGGCATCGCGAACCTTATCCACAGGCAGTGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTg                                                                                                                                                      <  1:1126439/150‑1 (MQ=255)
 gggTCTGCAGCAATTGGCCGATGTTGTCATCGATATAGCTGCACGCGCCGAAGTAGGCGCGGCGGGCATCGCGAACCTTATCCACAGGCAGTGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTg                                                                                                                                                      <  1:2320713/150‑1 (MQ=255)
    tCTGCAGCAATTGGCCGATGTTGTCATCGATATAGCTGCACGCGCCGAAGTAGGCGCGGCGGGCATCGCGAACCTTATCCACAGGCAGTGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCt                                                                                                                                                   <  2:7973158/150‑1 (MQ=255)
    tCTGCAGCAATTGGCCGATGTTGTCATCGATATAGCTGCACGCGCCGAAGTAGGCGCGGCGGGCATCGCGAACCTTATCCACAGGCAGTGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCt                                                                                                                                                   <  1:8525544/150‑1 (MQ=255)
        cagcaATTGGCCGATGTTGTCATCGATATAGCTGCACGCGCCGAAGTAGGCGCGGCGGGCATCGCGAACCTTATCCACAGGCAGTGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTGGCt                                                                                                                                               <  2:7630902/150‑1 (MQ=255)
        cagcaATTGGCCGATGTTGTCATCGATATAGCTGCACGCGCCGAAGTAGGCGCGGCGGGCATCGCGAACCTTATCCACAGGCAGTGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTGGCt                                                                                                                                               >  1:6744218/1‑150 (MQ=255)
          gcaATTGGCCGATGTTGTCATCGATATAGCTGCACGCGCCGAAGTAGGCGCGGCGGGCATCGCGAACCTTATCCACAGGCAGTGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTGGCTGa                                                                                                                                             >  1:5171110/1‑150 (MQ=255)
            aaTTGGCCGATGTTGTCATCGATATAGCTGCACGCGCCGAAGTAGGCGCGGCGGGCATCGCGAACCTTATCCACAGGCAGTGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTGGCTGAAc                                                                                                                                           >  1:5880666/1‑150 (MQ=255)
              ttGGCCGATGTTGTCATCGATATAGCTGCACGCGCCGAAGTAGGCGCGGCGGGCATCGCGAACCTTATCCACAGGCAGTGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTGGCTGAACTc                                                                                                                                         <  2:4005081/150‑1 (MQ=255)
               tGGCCGATGTTGTCATCGATATAGCTGCACGCGCCGAAGTAGGCGCGGCGGGCATCGCGAACCTTATCCACAGGCAGTGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTGGCTGAACTCg                                                                                                                                        <  1:1910722/150‑1 (MQ=255)
                ggCCGATGTTGTCATCGATATAGCTGCACGCGCCGAAGTAGGCGCGGCGGGCATCGCGAACCTTATCCACAGGCAGTGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTGGCTGAACTCgg                                                                                                                                       <  1:7233416/150‑1 (MQ=255)
                 gCCGATGTTGTCATCGATATAGCTGCACGCGCCGAAGTAGGCGCGGCGGGCATCGCGAACCTTATCCACAGGCAGTGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTgg                                                                                                                                                 <  1:4189115/139‑1 (MQ=255)
                  ccGATGTTGTCATCGATATAGCTGCACGCGCCAAAGTAGGCGCGGCGGGCATCGCGAACCTTATCCACAGGCAGTGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGTGCCGAGTTCAGCCAGGCAGAACTCgacc                                                                                                                                     <  2:4558741/150‑4 (MQ=255)
                   cGATGTTGTCATCGATATAGCTGCACGCGCCGAAGTAGGCGCGGCGGGCATCGCGAACCTGATCCACAGGCAGTGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGTGCCGAGTTCAGCCAGGCAGAACTCgaccc                                                                                                                                    <  2:6809720/150‑5 (MQ=255)
                    gATGTTGTCATCGATATAGCTGCACGCGCCGAAGTAGGCGCGGCGGGCATCGCGAACCTTATCCACAGGCAGTGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTGGCTGAACTCGGCACg                                                                                                                                   >  2:4918260/1‑150 (MQ=255)
                        ttgtCATCGATATAGCTGCACGCGCCGAAGTAGGCGCGGCGGGCATCGCGAACCTTATCCACAGGCAGTGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTGGCTGAACTCGGCACGGGGc                                                                                                                               <  2:8286185/150‑1 (MQ=255)
                           tCATCGATATAGCTGCACGCGCCGAAGTAGGCGCGGCGGGCATCGCGAACCTTATCCACAGGCAGTGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTGGCTGAACTCGGCACGGGGCATg                                                                                                                            >  1:7097286/1‑150 (MQ=255)
                             aTCGATCTAGCTGCACGCGCCGAAGTAGGCGCGGCGGGCATCGCGAACCTTATCCACAGGCAGTGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTGGCTGAACTCGGCACGGGGCATggg                                                                                                                          >  2:5107187/1‑150 (MQ=255)
                              tCGATATAGCTGCACGCGCCGAAGTAGGCGCGGCGGGCATCGCGAACCTTATCCACAGGCAGTGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTGGCTGAACTCGGCACGGGGCATgggg                                                                                                                         <  1:36474/150‑1 (MQ=255)
                              tCGATATAGCTGCACGCGCCGAAGTAGGCGCGGCGGGCATCGCGAACCTTATCCACAGGCAGTGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTGGCTGAACTCGGCACGGGGCATgggg                                                                                                                         >  1:4778782/1‑150 (MQ=255)
                                gATATAGCTGCACGCGCCGAAGTAGGCGCGGCGGGCATCGCGAACCTTATCCACAGGCAGTGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGTGTCGAGTTCTGCCTGGCTGAACTCGGCACGGGGCATGGGGat                                                                                                                       <  2:6118998/150‑1 (MQ=255)
                                gATATAGCTGCACGCGCCGAAGTAGGCGCGGCGGGCATCGCGAACCTTATCCACAGGCAGTGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTGGCTGAACTCGGCACGGGGCATGGGGat                                                                                                                       <  2:2464890/150‑1 (MQ=255)
                                 atatAGCTGCACGCGCCGAAGTAGGCGCGGCGGGCATCGCGAACCTTATCCACAGGCAGTGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTGGCTGAACTCGGCACGGGGCATGGGGata                                                                                                                      >  1:1357969/1‑150 (MQ=255)
                                   ataGCTGCACGCGCCGAAGTAGGCGCGGCGGGCATCGCGAACCTTATCCACAGGCAGTGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTGGCTGAACTCGGCACGGGGCATggggcaccg                                                                                                                    >  1:413815/1‑145 (MQ=255)
                                   ataGCTGCACGCGCCGAAGTAGGCGCGGCGGGCATCGCGAACCTTATCCACAGGCAGTGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTGGCTGAACTCGGCACGGGGCATggggcaacg                                                                                                                    >  2:2626784/1‑145 (MQ=255)
                                   ataGCTGCACGCGCCGAAGTAGGCGCGGCGGGCATCGCGAACCTTATCCACAGGCAGTGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTGGCTGAACTCGGCACGGGGCATGGGGATATc                                                                                                                    <  2:1779226/150‑1 (MQ=255)
                                   ataGCTGCACGCGCCGAAGTAGGCGCGGCGGGCATCGCGAACCTTATCCACAGGCAGTGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTGGCTGAACTCGGCACGGGGCATGGGGATATc                                                                                                                    >  1:4338399/1‑150 (MQ=255)
                                         gCACGCGCCGAAGTAGGCGCGGCGGGCATCGCGAACCTTATCCACAGGCAGTGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTGGCTGAACTCGGCACGGGGCATGGGGATATCCACAcc                                                                                                              >  1:3178257/1‑150 (MQ=255)
                                            cgcgCCGAAGTAGGCGCGGCGGGCATCGCGAACCTTATCCACAGGCAGTGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTGGCTGAACTCGGCACGGGGCATGGGGATATCCACACCTTc                                                                                                           <  1:509660/150‑1 (MQ=255)
                                             gcgcCGAAGTAGGCGCGGCGGGCATCGCGAACCTTATCCACAGGCAGTGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTGGCTGAACTCGGCACGGGGCATGGGGATATCCACACCTTCg                                                                                                          <  2:7088778/150‑1 (MQ=255)
                                             gcgcCGAAGTAGGCGCGGCGGGCATCGCGAACCTTATCCACAGGCAGTGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTGGCTGAACTCGGCACGGGGCATGGGGATATCCACACCTTCg                                                                                                          >  1:6687360/1‑150 (MQ=255)
                                               gcCGAAGTAGGCGCGGCGGGCATCGCGAACCTTATCCACAGGCAGTGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTGGCTGAACTCGGCACGGGGCATGGGGATATCCACACCTTCGTa                                                                                                        >  2:6248975/1‑150 (MQ=255)
                                               gcCGAAGTAGGCGCGGCGGGCATCGCGAACCTTATCCACAGGCAGTGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTGGCTGAACTCGGCACGGGGCATGGGGATATCCACACCTTCGTa                                                                                                        <  1:6982099/150‑1 (MQ=255)
                                               gcCGAAGTAGGCGCGGCGGGCATCGCGAACCTTATCCACAGGCAGTGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTGGCTGAACTCGGCACGGGGCATGGGGATATCCACACCTTCGTa                                                                                                        <  1:6047941/150‑1 (MQ=255)
                                                  gAAGTAGGCGCGGCGGGCATCGCGAACCTTATCCACAGGCAGTGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTGGCTGAACTCGGCACGGGGCATGGGGATATCCACACCTTCGTAGCg                                                                                                     <  2:6902713/150‑1 (MQ=255)
                                                    aGTAGGCGCGGCGGGCATCGCGAACCTTATCCACAGGCAGTGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTGGCTGAACTCGGCACGGGGCATGGGGATATCCACACCTTCGTAGCg                                                                                                     <  2:2112833/148‑1 (MQ=255)
                                                    aGTAGGCGCGGCGGGCATCGCGAACCTTATCCACAGGCAGTGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTGGCTGAACTCGGCACGGGGCATGGGGATATCCACACCTTCGTAGCg                                                                                                     >  1:2112833/1‑148 (MQ=255)
                                                      tAGGCGCGGCGGGCATCGCGAACCTTATCCACAGGCAGTGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTGGCTGAACTCGGCACGGGGCATGGGGATATCCACACCTTCGTAGCGGTcc                                                                                                 >  1:6226431/1‑150 (MQ=255)
                                                      tAGGCGCGGCGGGCATCGCGAACCTTATCCACAGGCAGTGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTGGCTGAACTCGGCACGGGGCATGGGGATATCCACACCTTCGTAGCGGTcc                                                                                                 <  1:4921691/150‑1 (MQ=255)
                                                      tAGGCGCGGCGGGCATCGCGAACCTTATCCACAGGCAGTGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCATGGCTGAACTCGGCACGGGGCATGGGGATATCCACACCTTCGTAGCGGTcc                                                                                                 <  2:5650728/150‑1 (MQ=255)
                                                        ggCGCGGCGGGCATCGCGAACCTTATCCACAGGCAGTGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTGGCTGAACTCGGCACGGGGCATGGGGATATCCACACCTTCGTAGCGGTCCCa                                                                                               <  1:7999648/150‑1 (MQ=255)
                                                         gcgcGGCGGGCATCGCGAACCTTATCCACAGGCAGTGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTGGCTGAACTCGGCACGGGGCATGGGGATATCCACACCTTCGTAGCGGTCCCAg                                                                                              >  2:1444150/1‑150 (MQ=255)
                                                           gcggcgGGCATCGCGAACCTTATCCACAGGCAGTGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTGGCTGAACTCGGCACGGGGCATGGGGATATCCACACCTTCGTAGCGGTCCCAGTa                                                                                            >  1:1771757/1‑150 (MQ=255)
                                                                       gcgAACCTTATCCACAGGCAGTGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTGGCTGAACTCGGCACGGGGCATGGGGATATCCACACCTTCGTAGCGGTCCCAGTAGCGTTTGGCAAt                                                                                >  1:95291/1‑150 (MQ=255)
                                                                          aaCCTTATCCACAGGCAGTGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTGGCTGAACTCGGCACGGGGCATGGGGATATCCACACCTTCGTAGCGGTCCCAGTAGCGTTTGGCAATGgt                                                                             <  2:4193100/150‑1 (MQ=255)
                                                                          aaCCTTATCCACAGGCAGTGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTGGCTGAACTCGGCACGGGGCATGGGGATATCCACACCTTCGTAGCGGTCCCAGTAGCGTTTGGCAATGgt                                                                             <  1:7905003/150‑1 (MQ=255)
                                                                          aaCCTTATCCACAGGCAGTGGCTTGTTCCAAACGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTGGCTGAACTCGGCACGGGGCATGGGGATATCCACACCTTCGTAGCGGTCCCAGTAGCGTTTGGCAATGgt                                                                             <  1:8258121/150‑1 (MQ=255)
                                                                            ccTTATCCACAGGCAGTGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTGGCTGAACTCGGCACGGGGCATGGGGATATCCACACCTTCGTAGCGGTCCCAGTAGCGTTTGGCAATGgtgt                                                                           >  2:5365452/1‑150 (MQ=255)
                                                                              ttATCCACAGGCAGTGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTGGCTGAACTCGGCACGGGGCATGGGGATATCCACACCTTCGTAGCGGTCCCAGTAGCGTTTGGCAATGGTGTAg                                                                         >  1:2640912/1‑150 (MQ=255)
                                                                              ttATCCACAGGCAGTGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTGGCTGAACTCGGCACGGGGCATGGGGATATCCACACCTTCGTAGCGGTCCCAGTAGCGTTTGGCAATGGTGTAg                                                                         <  2:52768/150‑1 (MQ=255)
                                                                              ttATCCACAGGCAGTGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCCGTGAATGCGGGTCGAGTTCTGCCTGGCTGAACTCGGCACGGGGCATGGGGATATCCACACCTTCGTAGCGGTCCCAGTAGCGTTTGGCAATGGTGTAg                                                                         >  2:7609488/1‑150 (MQ=255)
                                                                               tATCCACAGGCAGTGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTGGCTGAACTCGGCACGGGGCATGGGGATATCCACACCTTCGTAGCGGTCCCAGTAGCGTTTGGCAATGGTGTAgg                                                                        <  2:6575808/150‑1 (MQ=255)
                                                                                 tCCACAGGCAGTGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTGGCTGAACTCGGCACGGGGCATGGGGATATCCACACCTTCGTAGCGGTCCCAGTAGCGTTTGGCAATGGTGTAgggg                                                                      >  1:8762796/1‑150 (MQ=255)
                                                                                             ggCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTGGCTGAACTCGGCACGGGGCATGGGGATATCCACACCTTCGTAGCGGTCCCAGTAGCGTTTGGCAATGGTGTAGGGGTCATGCgggtgg                                                          <  2:1126070/150‑1 (MQ=255)
                                                                                             ggCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTGGCTGAACTCGGCACGGGGCATGGGGATATCCACACCTTCGTAGCGGTCCCAGTAGCGTTTGGCAATGGTGTAGGGGTCATGCgggtgg                                                          <  1:2799897/150‑1 (MQ=255)
                                                                                             ggCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTGGCTGAACTCGGCACGGGGCATGGGGATATCCACACCTTCGTAGCGGTCCCAGTAGCGTTTGGCAATGGTGTAGGGGTCATGCgggtgg                                                          <  1:6433185/150‑1 (MQ=255)
                                                                                             ggCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTGGCTGAACTCGGCACGGGGCATGGGGATATCCACACCTTCGTAGCGGTCCCAGTAGCGTTTGGCAATGGTGTAGGGGTCATGCgggtgg                                                          <  1:7391521/150‑1 (MQ=255)
                                                                                             ggCTTGTTCCAAAGGTCATAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTGGCTGAACTCGGCACGGGGCATGGGGATATCCACACCTTCGTAGCGGTCCCAGTAGCGTTTGGCAATGGTGTAGGGGTCATGCgggtgg                                                          <  1:2405207/150‑1 (MQ=255)
                                                                                              gCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTGGCTGAACTCGGCACGGGGCATGGGGATATCCACACCTTCGTAGCGGTCCCAGTAGCGTTTGGCAATGGTGTAGGGGTCATGCgggtggg                                                         <  1:4013629/150‑1 (MQ=255)
                                                                                                ttGTTCCAATGGTGGTAGTCCTTCAGCAGGCGGTGTGAATGCGGGTCGTGTTCTGGGTGGCTGAACTCGGCAGGGGGGATGGGGATATCGACAGGTTCGTAGCGGTCGCTGTAGCGTTTGGCTATGGTGTAGGGGTCATGCGGGTGGGTc                                                       >  2:7415188/1‑150 (MQ=255)
                                                                                                 tGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTGGCTGAACTCGGCACGGGGCATGGGGATATCCACACCTTCGTAGCGGTCCCAGTAGCGTTTGGCAATGGTGTAGGGGTCATGCGGGTGGGTCa                                                      >  1:7662435/1‑150 (MQ=255)
                                                                                                     ccAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTGGCTGAACTCGGCACGGGGCATGGGGATATCCACACCTTCGTAGCGGTCCCAGTAGCGTTTGGCAATGGTGTAGGGGTCATGCGGGTGGGTCATCGa                                                  >  2:3193343/1‑150 (MQ=255)
                                                                                                      cAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTGGCTGAACTCGGCACGGGGCATGGGGATATCCACACCTTCGTAGCGGTCCCAGTAGCGTTTGGCAATGGGGTAGGGGTCATGCGGGTGGGTCATCGaa                                                 >  1:2169283/1‑150 (MQ=255)
                                                                                                       aaaGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTGGCTGAACTCGGCACGGGGCATGGGGATATCCACACCTTCGTAGCGGTCCCAGTAGCGTTTGGCAATGGTGTAGGGGTCATGCGGGTGGGTCATCGaaa                                                >  2:2716727/1‑150 (MQ=255)
                                                                                                       aaaGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTGGCTGAACTCGGCACGGGGCATGGGGATATCCACACCTTCGTAGCGGTCCCAGTAGCGTTTGGCAATGGTGTAGGGGTCATGCGGGTGGGTCATCGaaa                                                <  2:382624/150‑1 (MQ=255)
                                                                                                          ggTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTGGCTGAACTCGGCACGGGGCATGGGGATATCCACACCTTCGTAGCGGTCCCAGTAGCGTTTGGCAATGGTGTAGGGGTCATGCGGGTGGGTCATCGAAACgg                                             >  1:7570012/1‑150 (MQ=255)
                                                                                                            tCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTGGCTGAACTCGGCACGGGGCATGGGGATATCCACACCTTCGTAGCGGTCCCAGTAGCGTTTGGCAATGGTGTAGGGGTCATGCGGGTGGGTCATCGAAACGGTc                                           >  2:2344486/1‑150 (MQ=255)
                                                                                                             cGTAGGCCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTGGCTGAACTCGGCACGGGGCATGGGGATATCCACACCTTCGTAGCGGTCCCAGTAGCGTTTGGCAATGGTGTAGGGGTCATGCGGGTGGGTCATCGAAACGGTCa                                          <  2:5711641/150‑1 (MQ=255)
                                                                                                             cGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTGGCTGAACTCGGCACGGGGCATGGGGATATCCACACCTTCGTAGCGGTCCCAGTAGCGTTTGGCAATGGTGTAGGGGTCATGCGGGTGGGTCATCGAAACGGTCa                                          <  2:6326116/150‑1 (MQ=255)
                                                                                                             cGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTGGCTGAACTCGGCACGGGGCATGGGGATATCCACACCTTCGTAGCGGTCCCAGTAGCGTTTGGCAATGGTGTAGGGGTCATGCGGGTGGGTCATCGAAACGGTCa                                          <  2:3077259/150‑1 (MQ=255)
                                                                                                             cGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTGGCTGAACTCGGCACGGGGCATGGGGATATCCACACCTTCGTAGCGGTCCCAGTAGCGTTTGGCAATGGTGTAGGGGTCATGCGGGTGGGTCATCGAAACGGTCa                                          <  1:164451/150‑1 (MQ=255)
                                                                                                              gTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTGGCTGAACTCGGCACGGGGCATGGGGATATCCACACCTTCGTAGCGGTCCCAGTAGCGTTTGGCAATGGTGTAGGGGTCATGCGGGTGGGTCATCGAAACGGTCAg                                         >  2:3945705/1‑150 (MQ=255)
                                                                                                              gTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTGGCTGAACTCGGCACGGGGCATGGGGATATCCACACCTTCGTAGCGGTCCCAGTAGCGTTTGGCAATGGTGTAGGGGTCATGCGGGTGGGTCATCGAAACGGTCAg                                         <  2:8444546/150‑1 (MQ=255)
                                                                                                              gTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTGGCTGAACTCGGCACGGGGCATGGGGATATCCACACCTTCGTAGCGGTCCCAGTAGCGTTTGGCAATGGTGTAGGGGTCATGCGGGTGGGTCATCGAAACGGTCAg                                         <  1:2951318/150‑1 (MQ=255)
                                                                                                               tAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTGGCTGAACTCGGCACGGGGCATGGGGATATCCACACCTTCGTAGCGGTCCCAGTAGCGTTTGGCAATGGTGTAGGGGTCATGCGGGTGGGTCATCGAAACGGTCAgg                                        <  2:2502898/150‑1 (MQ=255)
                                                                                                                aGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTGGCTGAACTCGGCACGGGGCATGGGGATATCCACACCTTCGTAGCGGTCCCAGTAGCGTTTGGCAATGGTGTAGGGGTCATGCgg                                                              >  2:4396137/1‑127 (MQ=255)
                                                                                                                aGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTGGCTGAACTCGGCACGGGGCATGGGGATATCCACACCTTCGTAGCGGTCCCAGTAGCGTTTGGCAATGGTGTAGGGGTCATGCgg                                                              <  1:4396137/127‑1 (MQ=255)
                                                                                                                 gACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTGGCTGAACTCGGCACGGGGCATGGGGATATCCACACCTTCGTAGCGGTCCCAGTAGCGTTTGGCAATGGTGTAGGGGTCATGCGGGTGGGTCATCGAAACGGTCAGGCa                                      >  2:180217/1‑150 (MQ=255)
                                                                                                                  aCCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTGGCTGAACTCGGCACGGGGCATGGGGATATCCACACCTTCGTAGCGGTCCCAGTAGCGTTTGGCAATGGTGTAGGGGTCATGCGGGTGGGTCATCGAAACGGTCAGGCaa                                     >  2:8753426/1‑150 (MQ=255)
                                                                                                                     ttCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTGGCTGAACTCGGCACGGGGCATGGGGATATCCACACCTTCGTAGCGGTCCCAGTAGCGTTTGGCAATGGTGTAGGGGTCATGCgggtgg                                                          >  2:6035667/1‑126 (MQ=255)
                                                                                                                         gcagGCGCTGTGAATGCGGGTCGAGTTCTGCCTGGCTGAACTCGGCACGGGGCATGGGGATATCCACACCTTCGTAGCGGTCCCAGTAGCGTTTGGCAATGGTGTAGGGGTCATGCGGGTGGGTCATCGAAACGGTCAGGCAAAATggcc                              <  2:780950/150‑1 (MQ=255)
                                                                                                                         gcagGCGCTGTGAATGCGGGTCGAGTTCTGCCTGGCTGAACTCGGCACGGGGCATGGGGATATCCACACCTTCGTAGCGGTCCCAGTAGCGTTTGGCAATGGTGTAGGGGTCATGCGGGTGGGTCATCGAAACGGTCAGGCAAAATggcc                              <  2:6445872/150‑1 (MQ=255)
                                                                                                                           tggCGCTGTGACTGCGGGTCGAGTTCTGCCTGGCTGTTCTCGGCAGGGGGGATGGGGATATCCAGTCCTTGGTAGCGGTCCCTGTAGCGTTTGGGATTGGTGTTGGGGTCATGCGGGTGGGTCATCGAATGGGTCTGGGAAAATggccgg                            <  2:1888777/149‑1 (MQ=255)
                                                                                                                           agGCGCTGTGAATGCGGGTCGAGTTCTGCCTGGCTGAACTCGGCACGGGGCATGGGGATATCCACACCTTCGTAGCGGTCCCAGTAGCGTTTGGCAATGGTGTAGGGGTCATGCGGGTGGGTCATCGAAACGGTCAGGCAAAATggccgg                            >  1:4589677/1‑150 (MQ=255)
                                                                                                                           agGCGCTGTGAATGCGGGTCGAGTTCTGCCTGGCTGAACTCGGCACGGGGCATGGGGATATCCACACCTTCGTAGCGGTCCCAGTAGCGTTTGGCAATGGTGTAGGGGTCATGCGGGTGGGTCATCGAAACGGTCAGGCAAAATggccgg                            <  2:7363632/150‑1 (MQ=255)
                                                                                                                                    gAATGCGGGTCGAGTTCTGCCTGGCTGAACTCGGCACGGGGCATGGGGATATCCACACCTTCGTAGCGGTCCCAGTAGCGTTTGGCAATGGTGTAGGGGTCATGCGGGTGGGTCATCGAAACGGTCAGGCAAAATGGCCGGCCATCGttt                   >  2:8041823/1‑150 (MQ=255)
                                                                                                                                        ccgtgccGAGTTCAGCCAGGCAGAACTCGACCCGGGGCATGGGGATATCCACACCTTCGTAGCGGTCCCAGTAGCGTTTGGCAATGGTGTAGGGGTCATGCGGGTGGGTCATCGAAACGGTCAGGCAAAATGGCCGGCCATCGTTTTcgc               >  1:3784570/7‑150 (MQ=255)
                                                                                                                                          gtgTCGAGTTCTGCCTGGCTGAACTCGGCACGGGGCATGGGGATATCCACACCTTCGTAGCGGTCCCATTAGCGTTTGGCAATGGTGTAGGGGTCATGCGGGTGGGTCATCGAAACGGTCAGGCAAAATGGCCGGCCATCGTTTTCgcgc             >  2:2796636/3‑150 (MQ=255)
                                                                                                                                          gggTCGAGTTCTGCCTGGCTGAACTCGGCACGGGGCATGGGGATATCCACACCTTCGTAGCGGTCCCAGTAGCGTTTGGCAATGGTGTAGGGGTCATGCGGGTGGGTCATCGAATCGGTCAGGCAAAATGGCCGGCCATCGTTTTCgcgc             >  2:5190205/1‑150 (MQ=255)
                                                                                                                                          gggTCGAGTTCTGCCTGGCTGAACTCGGCACGGGGCATGGGGATATCCACACCTTCGTAGCGGTCCCAGTAGCGTTTGGCAATGGTGTAGGGGTCATGCGGGTGGGTCATCGAAACGGTCAGGCAAAATGGCCGGCCATCGTTTTCgcgc             <  1:4673916/150‑1 (MQ=255)
                                                                                                                                           ggTCGAGTTCTGCCTGGCTGAACTCGGCACGGGGCATGGGGATATCCACACCTTCGTAGCGGTCCCAGTAGCGTTTGGCAATGGTGTAGGGGTCATGCGGGTGGGTCATCGAAACGGTCAGGCAAAATGGCCGGCCATCGTTTTc                 <  1:1590/145‑1 (MQ=255)
                                                                                                                                           ggTCGAGTTCTGCCTGGCTGAACTCGGCACGGGGCATGGGGATATCCACACCTTCGTAGCGGTCCCAGTAGCGTTTGGCAATGGTGTAGGGGTCATGCGGGTGGGTCATCGAAACGGTCAGGCAAAATGGCCGGCCATCGTTTTCGCGca            <  1:8219931/150‑1 (MQ=255)
                                                                                                                                            gTCGAGTTCTGCCTGGCTGAACTCGGCACGGGGCATGGGGATATCCACACCTTCGTAGCGGTCCCAGTAGCGTTTGGCAATGGTGTAGGGGTCATGCGGGTGGGTCATCGAAACGGTCAGGCAAAATGGCCGGCCATCGTTTTCGCGcac           >  1:5439193/1‑150 (MQ=255)
                                                                                                                                              cGAGTTCAGCCAGGCAGAACTCGGCACGGGGCATGGGGATATCCACACCTTCGTAGCGGTCCCAGTAGCGTTTGGCAATGGTGTAGGGGTCATGCGGGTGGGTCATCGAAACGGTCAGGCAAAATGGCCGGCCATCGTTTTCGCGCACAt         >  1:5651803/1‑150 (MQ=255)
                                                                                                                                              cGAGTTCAGCCAGGCAGAACTCGGCACGGGGCATGGGGATATCCACACCTTCGTAGCGGTCCCAGTAGCGTTTGGCAATGGTGTAGGGGTCATGCGGGTGGGTCATCGAAACGGTCAGGCAAAATGGCCGGCCATCGTTTTCGCGCACAt         >  2:8146877/1‑150 (MQ=255)
                                                                                                                                              cGAGTTCAGCCAGGCAGAACTCGGCACGGGGCATGGGGATATCCACACCTTCGTAGCGGTCCCAGTAGCGTTTGGCAAGGGTGTAGGGGTCATGCGGGTGGGTCATCGAAACGGTCAGGCAAAATGGCCGGCCATCGTTTTCGCGCACAt         >  2:4992124/1‑150 (MQ=255)
                                                                                                                                               gAGGTCGGCCTGGCTGCACTCGGCACGGGGCATGGGGAGCTCCACACCTTCGTAGCGGTCCCAGTCGCGTTTGGCCATGGTGTAGGGGTCAGGCGGGTGGGTCATCGCAACGGTCAGGCACAATGGCCGGCCATCGTTTTCgcgccccgg        <  2:1661280/150‑6 (MQ=255)
                                                                                                                                                  ttCTGCCTGGCTGAACTCGGCACGGGGCATGGGGATATCCACACCTTCGTAGCGGTCCCAGTAGCGTTTGGCAATGGTGTAGGGGTCATGCGGGTGGGTCATCGAAACGGTCAGGCAAAATGGCCGGCCATCGTTTTCGCGCACATGGTc     >  2:1563685/1‑150 (MQ=255)
                                                                                                                                                     tGCCTGGCTGAACTCGGCACGGGGCATGGGGATATCCACACCTTCGTAGCGGTCCCAGTAGCGTTTGGCAATGGTGTAGGGGTCATGCGGGTGGGTCATCGAAACGGTCAGGCAAAATGGCCGGCCATCGTTTTCGCGCACATGGTCGTa  >  1:1532384/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
AGGGTCTGCAGCAATTGGCCGATGTTGTCATCGATATAGCTGCACGCGCCGAAGTAGGCGCGGCGGGCATCGCGAACCTTATCCACAGGCAGTGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTGGCTGAACTCGGCACGGGGCATGGGGATATCCACACCTTCGTAGCGGTCCCAGTAGCGTTTGGCAATGGTGTAGGGGTCATGCGGGTGGGTCATCGAAACGGTCAGGCAAAATGGCCGGCCATCGTTTTCGCGCACATGGTCGTA  >  NC_002947/88500‑88798

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: ‑