Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,982,021 A→T 6.6% noncoding (28/207 nt) PP_mr43 ← Cobalamin regulatory region

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,982,0210AT 6.6% 211.6 / 4.8 72noncoding (28/207 nt)PP_mr43Cobalamin regulatory region
Reads supporting (aligned to +/- strand):  ref base A (33/34);  new base T (3/2);  total (36/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.58e-01

TGGCTTGCACCGTTCCATGTCCCCGTTCGGGCTGGCGGGATGACGCCTTCCCGCGCGTTGGCGCAGTGGCTGTGTCGATCCGTTGTCGGTGCTTACAGTTGCGGGGGCAGCCGCGGCTTGTACCGCGTTCCCGTCTTAGCTTCGGCCAGGCCGAAGAACCTCGAAGTGGCAAGGCTACGCAGTGGCTGGCGGGTGGTCAACCATTGCTGCGGGGCCGGTCGACAACGCTGGGCGCAAGTGCGTTACAGTGGCGCCTTATCCACGGTCAAGGACGCACCATGCCCCTCGCC  >  NC_002947/3981874‑3982163
                                                                                                                                                   |                                                                                                                                              
tGGCTTGCACCGTTCCATGTCCCCGTTCGGGCTGGCGGGATGACGCCTTCCCGCGCGATTGCGCAGTGGCTGTGTCGATCCGTTGTCGGTGCTTAAAAGCCTTGGGGCAGCCGCGGCTTGTACCGCGTTCCCGTCTTAGCTTCGGCCAgg                                                                                                                                              >  1:1310518/1‑150 (MQ=255)
    ttGCACCGTTCCATGTCCCCGTTCGGGCTGGCGGGATGACGCCTTCCCGCGCGTTGGCGCAGTGGCTGTGTCGATCCGTTGTCGGTGCTTACAGTTGCGGGGGCAGCCGCGGCTTGTACCGCGTTCCCGTCTTAGCTTCGGCCTGGCCga                                                                                                                                          >  2:1579722/1‑150 (MQ=255)
        aCCGTTCCATGTCCCCGTTCGGGCTGGCGGGATGACGCCTTCCCGCGCGTTGGCGCAGTGGCTGTGTCGATCCGTTGTCGGTGCTTACAGTTGCGGGGGCAGCCGCGGCTTGTACCGCGTTCCCGTCTTAGCTTCGGCCAGGCCgaagaa                                                                                                                                      >  1:433069/1‑150 (MQ=255)
         ccGTTCCATGTCCCCGTTCGGGCTGGCGGGATGACGCCTTCCCGCGCGTTGGCGCAGTGGCTGTGTCGATCCGTTGTCGGTGCTTACAGTTGCGGGGGCAGCCGCGGCTTGTACCGCGTTCCCGTCTTAGCTTCGGCCTGGCCgaagcta                                                                                                                                     <  1:1579722/150‑4 (MQ=255)
          cGTTCCATGTCCCCGTTCGGGCTGGCGGGATGACGCCTTCCCGCGCGGTGGCGCAGTGGCTGTGTCGATCCGTTGTCGGTGCTTACAGTTGCGGGGGCAGCCGCGGCTTGTACCGCGTTCCCGTCTTAGCTTCGGCCAGGCCGAAGAAcc                                                                                                                                    >  1:398307/1‑150 (MQ=255)
               cATGTCCCCGTTCGGGCTGGCGGGATGACGCCTTCCCGCGCGTTGGCGCAGTGGCTGTGTCGATCCGTTGTCGGTGCTTACAGTTGCGGGGGCAGCCGCGGCTTGTACCGCGTTCCCGTCTTAGCTTCGGCCAGGCCGAAGAACCTCGaa                                                                                                                               >  1:57566/1‑150 (MQ=255)
                  ggccccGTTCGGGCTGGCGGGAGGACGCCTTCCCGCGCGTTGGCGCAGTGGCTGTGTCGATCCGTTGTCGGTGCTTACAGTTGCGGGGGCAGCCGCGGCTTGTACCGCGTTCCCGTCTTAGCTTCGGCCAGGCCGAAGAACCTCGAAGTg                                                                                                                            <  1:2231642/148‑1 (MQ=255)
                     cccGTTCGGGCTGGCGGGATGACGCCTTCCCGCGCGTTGGCGCAGTGGCTGTGTCGATCCGTTGTCGGTGCTTACAGTTGCGGGGGCAGCCGCGGCTTGTACCGCGTTCCCGTCTTAGCTTCGGCCAGGCCGAAGAACCTCGAAGTGGCa                                                                                                                         >  1:1468274/1‑150 (MQ=255)
                     cccGTTCGGGCTGGCGGGATGACGCCTTCCCGCGCGTTGGCGCAGCGGCTGTGTCGATCCGTTGTCGGTGCTTACAGTTGCGGGGGCAGCCGCGGCTTGTACCGCGTTCCCGTCTTAGCTTCGGCCAGGCCGAAGAACCTCGAAGTGGCa                                                                                                                         >  1:1272829/1‑150 (MQ=255)
                      ccGTTCGGGCTGGCGGGATGACGGCTTCTCGCGCGTTGGCGCAGTGGCTGTGTCGATCCGTTGTCGGTGCTTACAGTTGCGGGGGCAGCCGCGGCTTGTACCGCGTTCCCGTCTTAGCTTCGGCCAGGCCGAAGAACCTCGATTTGGCaa                                                                                                                        <  2:1310518/150‑1 (MQ=255)
                      ccGTTCGGGCTGGCGGGATGACGCCTTCCCGCGCGTTGGCGCAGTGGCTGTGTCGATCCGTTGTCGGTGCTTACAGTTGCGGGGGCAGCCGCGGCTTGTACCGCGTTCCCGTCTTAGCTTCGGCCAGGCCGAAGAACCTCGAAGTGGCaa                                                                                                                        <  2:1952164/150‑1 (MQ=255)
                      ccGTTCGGGCTGGCGGGATGACGCCTTCCCGCGCGTTGGCGCAGTGGCTGTGTCGATCCGTTGTCGGTGCTTACAGTTGCGGGGGCAGCCGCGGCTTGTACCGCGTTCCCGTCTTAGCTTCGGCCAGGCCGAAGAACCTCGAAGTGGCaa                                                                                                                        <  1:1823283/150‑1 (MQ=255)
                       cGTTCGGGCTGGCGGGATGACGCCTTCCCGCGCGTTGGCGCAGTGGCTGTGTCGATCCGTTGTCGGTGCTTACAGTTGCGGGGGCAGCCGCGGCTTGTACCGCGTTCCCGTCTTAGCTTCGGCCAGGCCGAAGAACCTCGAAGTGGCAAg                                                                                                                       <  2:947892/150‑1 (MQ=255)
                             ggCTGGCGGGATGACGCCTTCCCGCGCGTTGGCGCAGTGGCTGTGTCGATCCGTTGTCGGTGCTTACAGTTGCGGGGGCAGCCGCGGCTTGTACCGCGTTCCCGTCTTAGCTTCGGCCAGGCCGAAGAACCTCGCAGTGGc                                                                                                                          >  1:1991306/1‑141 (MQ=255)
                             ggCTGGCGGGATGACGCCTTCCCGCGCGTTGGCGCAGTGGCTGTGTCGATCCGTTGTCGGTGCTTACAGTTGCGGGGGCAGCCGCGGCTTGTACCGCGTTCCCGTCTTAGCTTCGGCCAGGCCGAAGAACCTCGAAGTGGc                                                                                                                          <  2:1991306/141‑1 (MQ=255)
                                 ggCGGGATGACGCCTTCCCGCGCGTTGGCGCAGTGGCTGTGTCGATCCGTTGTCGGTGCTTACAGTTGCGGGGGCAGCCGCGGCTTGTACCGCGTTCCCGTCTTAGCTTCGGCCAGGCCGAAGAACCTCGAAGTGGCAAGGCTACGCAGt                                                                                                             <  1:2311766/150‑1 (MQ=255)
                                  gCGGGATGACGCCTTCCCGCGCGTTGGCGCAGTGGCTGTGTCGATCCGTTGTCGGTGCTTACAGTTGCGGGGGCAGCCGCGGCTTGTACCGCGTTCCCGTCTTAGCTTCGGCCAGGCCGAAGAACCTCGAAGTGGCAAGGCTACGCAGtg                                                                                                            <  2:433069/150‑1 (MQ=255)
                                    gggATGACGCCTTCCCGCGCGTTGGCGCAGTGGCTGTGTCGATCCGTTGTCGGTGCTTACAGTTGCGGGGGCAGCCGCGGCTTGTACCGCGTTCCCGTCTTAGCTTCGGCCAGGCCGAAGAACCTCGAAGTGGCaa                                                                                                                        >  1:1209630/1‑136 (MQ=255)
                                    gggATGACGCCTTCCCGCGCGTTGGCGCAGTGGCTGTGTCGATCCGTTGTCGGTGCTTACAGTTGCGGGGGCAGCCGCGGCTTGTACCGCGTTCCCGTCTTAGCTTCGGCCAGGCCGAAGAACCTCGAAGTGGCaa                                                                                                                        <  2:1209630/136‑1 (MQ=255)
                                    gggATGACGCCTTCCCGCGCGTTGGCGCAGTGGCTGTGTCGATCCGTTGTCGGTGCTTACAGTTGCGGGGGCAGCCGCGGCTTGTACCGCGTTCCCGTCTTAGCTTCGGCCAGGCCGAAGAACCTCGAAGTGGCAAGGCTACGCAGtggc                                                                                                          >  1:1292452/1‑150 (MQ=255)
                                         gACGCCTTCCCGCGCGTTGGCGCAGTGGCTGTGTCGATCCGTTGTCGGTGCTTACAGTTGCGGGGGCAGCCGCGGCTTGTACCGCGTTCCCGTCTTAGCTTCGGCCAGGCCGAAGAACCTCGAAGTGGCAAGGCTACGCAg                                                                                                              >  2:1838218/1‑141 (MQ=255)
                                         gACGCCTTCCCGCGCGTTGGCGCAGTGGCTGTGTCGATCCGTTGTCGGTGCTTACAGTTGCGGGGGCAGCCGCGGCTTGTACCGCGTTCCCGTCTTAGCTTCGGCCAGGCCGAAGAACCTCGAAGTGGCAAGGCTACGCAg                                                                                                              <  1:1838218/141‑1 (MQ=255)
                                         gACGCCTTCCCGCGCGTTGGCGCAGTGGCTGTGTCGATCCGTTGTCGGTGCTTACAGTTGCGGGGGCAGCCGCGGCTTGTACCGCGTTCCCGTCTTAGCTTCGGCCAGGCCGAAGAACCTCGAAGTGGCAAGGCTACGCAGTGGCTGGCg                                                                                                     >  2:1986774/1‑150 (MQ=255)
                                           cGCCTTCCCGCGCGTTGGCGCAGTGGCTGTGTCGATCCGTTGTCGGTGCTTACAGTTGCGGGGGCAGCCGCGGCTTGTACCGCGTTCCCGTCTTAGCTTCGGCCAGGCCGAAGAACCTCGAAGTGGCAAGGCTACGCAGTGGCTGGCggg                                                                                                   <  1:1718532/150‑1 (MQ=255)
                                             ccTTCCCGCGCGTTGGCGCAGTGGCTGTGTCGATCCGTTGTCGGTGCTTACAGTTGCGGGGGCAGCCGCGGCTTGTACCGCGTTCCCGTCTTAGCTTCGGCCAGGCCGAAGAACCTCGAAGTGGCAAGGCTACGCAGTGGCTGGCGggtg                                                                                                 <  1:1754634/150‑1 (MQ=255)
                                              cTTCCCGCGCGTTGGCGCAGTGGCTGTGTCGATCCGTTGTCGGTGCTTACAGTTGCGGGGGCAGCCGCGGCTTGTACCGCGTTCCCGTCTTAGCTTCGGCCAGGCCGAAGAAcc                                                                                                                                    <  2:508109/114‑1 (MQ=255)
                                              cTTCCCGCGCGTTGGCGCAGTGGCTGTGTCGATCCGTTGTCGGTGCTTACAGTTGCGGGGGCAGCCGCGGCTTGTACCGCGTTCCCGTCTTAGCTTCGGCCAGGCCGAAGAAcc                                                                                                                                    >  1:508109/1‑114 (MQ=255)
                                              cTTCCCGCGCGTTGGCGCAGTGGCTGTGTCGATCCGTTGTCGGTGCTTACAGTTGCGGGGGCAGCCGCGGCTTGTACCGCGTTCCCGTCTTAGCTTCGGCCAGGCCGAAGAACCTCGAAGTGGCAAGGCTACGCAGTGGCTGGCGggtgg                                                                                                <  2:1513241/150‑1 (MQ=255)
                                              cTTCCCGCGCGTTGGCGCAGTGGCTGTGTCGATCCGTTGTCGGTGCTTACAGTTGCGGGGGCAGCCGCGGCTTGTACCGCGTTCCCGTCTTAGCTTCGGCCAGGCCGAAGAACCTCGAAGTGGCAAGGCTACGCAGTGGCTGGCGggtgg                                                                                                >  2:2047837/1‑150 (MQ=255)
                                               ttCCCGCGCGTTGGCGCAGTGGCTGTGTCGATCCGTTGTCGGCGCTTACAGTTGCGGGGGCAGCCGCGGCTTGTACCGCGTTCCCGTCTTAGCTTCGGCCAGGCCGAAGAACCTCGAAGTGGCAAGGCTACGCAGTGGCTGGCGggtggt                                                                                               >  1:314417/1‑150 (MQ=255)
                                                      gtgTTGGCGCAGTGGCTGTGTCGATCCGTTGTCGGTGCTTACAGTTGCGGGGGCAGCCGCGGCTTGTACCGCGTTCCCGTCTTAGCTTCGGCCAGGCCGAAGAACCTCGAAGTGGCAAGGCTACGCAGTGGCTGGCGGGTGGTCAACCAt                                                                                        <  1:1492543/148‑1 (MQ=255)
                                                       cgTTGGCGCAGTGGCTGTGTCGATCCGTTGTCGGTGCTTACAGTTGCGGGGGCAGCCGCGGCTTGTACCGCGTTCCCGTCTTAGCTTCGGCCAGGCCGAAGAACCTCGAAGTGGCAAGGCTACGCAGTGGCTGGCGGGTGGTCAACCAtt                                                                                       <  2:57566/150‑1 (MQ=255)
                                                        gTTGGCGCAGTGGCTGTGTCGATCCGTTGTCGGTGCTTACAGTTGCGGGGGCAGCCGCGGCTTGTACCGCGTTCCCGTCTTAGCTTCGGCCAGGCCGAAGAACCTCGAAGTGGCAAGGCTACGCAGTGGCTGGCGGGTGGTCAACCATtg                                                                                      <  1:1986774/150‑1 (MQ=255)
                                                            gcgcAGTGGCTGTGTCGATCCGTTGTCGGTGCTTACAGTTGCGGGGGCAGCCGCGGCTTGTACCGCGTTCCCGTCTTAGCTTCGGCCAGGCCGAAGAACCTCGAAGTGGCAAGGCTACGCAGTGGCTGGCGGGTGGTCAACCATtgctgc                                                                                  >  2:1376699/1‑150 (MQ=255)
                                                              gcAGTGGCTGTGTCGATCCGTTGTCGGCGCTTACAGTTGCGGGGGCAGCCGCGGCTTGTACCGCGTTCCCGTCTTAGCTTCGGCCAGGCCGAAGAACCTCGAAGTGGCAAGGCTACGCAGTGGCTGGCGGGTGGTCAACCATTGCTGCgg                                                                                <  2:314417/150‑1 (MQ=255)
                                                                    gCTGTGTCGATCCGTTGTCGGTGCTTACAGTTGCGGGGGCAGCCGCGGCTTGTACCGCGTTCCCGTCTTAGCTTCGGCCAGGCCGAAGAACCTCGAAGTGGCAAGGCTACGCAGTGGCTGGCGGGTGGTCAACCATTGCTGCGGGGCCgg                                                                          >  2:1485616/1‑150 (MQ=255)
                                                                         gtCGATCCGTTGTCGGTGCTTACAGTTGCGGGGGCAGCCGCGGCTTGTACCGCGTTCCCGTCTTAGCTTCGGCCAGGCCGAAGAACCTCGAAGTGGCAAGGCTACGCAGTGGCTGGCGTGTGGTCAACCATTGCTGCGGGGCCGGTCGAc                                                                     <  1:1033852/150‑1 (MQ=255)
                                                                            gATCCGTTGTCGGTGCTTACAGTTGCGGGGGCAGCCGCGGCTTGTACCGCGTTCCCGTCTTAGCTTCGGCCAGGCCGAAGAACCTCGAAGTGGCAAGGCTACGCAGTGGCTGGCGGGTGGTCAACCATTGCTGCGGGGCCGGTCGACAAc                                                                  >  2:1158753/1‑150 (MQ=255)
                                                                                 gTTGTCGGTGCTTACAGTTGCGGGGGCAGCCGCGGCTTGTACCGCGTTCCCGTCTTAGCTTCGGCCAGGCCGAAGAACCTCGAAGTGGCAAGGCTACGCAGTGGCTGGCGGGTGGTCAACCATTGCTGCGGGGCCGGTCGACAACGCTgg                                                             <  1:1502642/150‑1 (MQ=255)
                                                                                       ggTGCTTACAGTTGCGGTGTCAGCCGCGGTTTGTACCGCGCTCCCGTCTTAGCTTCCGCCAGGCCGAAGAACCTCGAAGTGGCAAGGCTACGCCGTGGCTGGCAGGTGGTCAACCATTGCTGCGGGGCCGGTCGACAACGCTGGGCACaa                                                       <  2:66495/150‑1 (MQ=255)
                                                                                        gTGCTTACAGTTGCGGGGGCAGCCGCGGCTTGTACCGCGTTCCCGTCTTAGCTTCGGCCAGGCCGAAGAACCTCGAAGTGGCAAGGCTACGCAGTGGCTGGCGGGTGGTCAACCATTGCTGCGGGGCCGGTCGACAACGCTGGGCGCAAg                                                      >  2:414262/1‑150 (MQ=255)
                                                                                            ttACAGTTGCGGGGGCAGCCGCGGCTTGTACCGCGTTCCCGTCTTAGCTTCGGCCAGGCCGAAGAACCTCGCAGTGTGAAGGCTACGCAGTGGCTGGCGAGTGGTCAAACATTGCTGGGGGGCCGGTCGACAACGCTGGGAGCAAGTGCg                                                  >  2:1141241/1‑150 (MQ=255)
                                                                                             tACAGTTGCGGGGGCAGCCGCGGCTTGTACCGCGTTCCCGTCTTAGCTTCGGCCAGGCCGAAGAACCTCGAAGTGGCAAGGCTACGCAGTGGCTGGCGGGTGGTCAACCATTGCTGCGGGGCCGGTCGACAACGCTGGGCGCAAGTGCGt                                                 <  1:2138346/150‑1 (MQ=255)
                                                                                                         ggCAGCCGCGGCTTGTACCGCGTTCCCGTCTTAGCTTCGGCCAGGCCGAAGAACCTCGAAGTGGCAAGGCTACGCAGTGGCTGGCGGGTGGTCAACCATTGCTGCGGGGCCGGTCGACAACGCTGGGCGCAAGTGCGTTACAGTGGCGcc                                     <  1:569302/150‑1 (MQ=255)
                                                                                                            aGCCGCGGCTTGTACCGCGTTCCCGTCTTAGCTTCGGCCAGGCCGAAGAACCTCGAAGTGGCAAGGCTACGCAGTGGCTGGCGGGTGGTCAACCATTGCTGCGGGGCCGGTCGACAACGCTGGGCGCAAGTGCGTTACAGTGGCGCCTTa                                  <  2:1272829/150‑1 (MQ=255)
                                                                                                                   gCTTGTACCGCGTTCCCGTCTTAGCTTCGGCCAGGCCGAAGAACCTCGAAGTGGCAAGGCTACGCAGTGGCTGGCGGGTGGTCAACCATTGCTGCGGGGCCGGTCGACAACGCTGGGCGCAAGTGCGTTACAGTGGCGCCTTATCCACgg                           <  1:1485616/150‑1 (MQ=255)
                                                                                                                           cgcgTTCCCGTCTTAGCTTCGGCCAGGCCGAAGAACCTCGAAGTGGCAAGGCTACGCAGTGGCTGGCGGGTGGTCAAc                                                                                           >  2:1735193/1‑78 (MQ=255)
                                                                                                                           cgcgTTCCCGTCTTAGCTTCGGCCAGGCCGAAGAACCTCGAAGTGGCAAGGCTACGCAGTGGCTGGCGGGTGGTCAAc                                                                                           <  1:1735193/78‑1 (MQ=255)
                                                                                                                           cgcgTTCCCGTCTTAGCTTCGGCCAGGCCGAAGAACCTCGAAGTGGCAAGGCTACGCAGTGGCTGGCGGGTGGTCAACCATTGCTGCGGGGCCGGTCGACAACGCTGGGCGCAAGTGCGTTACAGt                                           <  1:426165/126‑1 (MQ=255)
                                                                                                                           cgcgTTCCCGTCTTAGCTTCGGCCAGGCCGAAGAACCTCGAAGTGGCAAGGCTACGCAGTGGCTGGCGGGTGGTCAACCATTGCTGCGGGGCCGGTCGACAACGCTGGGCGCAAGTGCGTTACAGt                                           >  2:426165/1‑126 (MQ=255)
                                                                                                                           cgcgTTCCCGTCTTAGCTTCGGCCAGGCCGAAGAACCTCGAAGTGGCAAGGCTACGCAGTGGCTGGCGGGTGGTCAACCATTGCTGCGGGGCCGGTCGACAACGCTGGGCGCAAGTGCGTTACAGTGGCGCCTTATCCACGGTCAAGGAc                   >  1:1741340/1‑150 (MQ=255)
                                                                                                                            gcgTTCCCGTCTTAGCTTCGGCCAGGCCGAAGAACCTCGAAGTGGCAAGGCTACGCAGTGGCTGGCGGGTGGTCAACCATTGCTGCGGGGCCGGTCGACAACGCTGGGCGCAAGTGCGTTACAGTGGCGCCTTATCCACGGTCAAGGACg                  >  2:1713537/1‑150 (MQ=255)
                                                                                                                                 cccGTCTTAGCTTCGGCCAGGCCGAAGAACCTCGAAGTGGCAAGGCTACGCAGTGGCTGGCGGGTGGTCAACCATTGCTGCGGGGCCGGTCGACAACGCTGGGCGCAAGTGCGTTACAGTGGCGcc                                     <  1:26463/126‑1 (MQ=255)
                                                                                                                                 cccGTCTTAGCTTCGGCCAGGCCGAAGAACCTCGAAGTGGCAAGGCTACGCAGTGGCTGGCGGGTGGTCAACCATTGCTGCGGGGCCGGTCGACAACGCTGGGCGCAAGTGCGTTACAGTGGCGcc                                     >  2:26463/1‑126 (MQ=255)
                                                                                                                                  ccGTCTTAGCTTCGGCCAGGCCGAAGAACCTCGAAGTGGCAAGGCTACGCAGTGGCTGGCGGGTGGTCAACCATTGCTGCGGGGCCGGTCGACAACGCTgg                                                             >  2:444111/1‑101 (MQ=255)
                                                                                                                                  ccGTCTTAGCTTCGGCCAGGCCGAAGAACCTCGAAGTGGCAAGGCTACGCAGTGGCTGGCGGGTGGTCAACCATTGCTGCGGGGCCGGTCGACAACGCTgg                                                             <  1:444111/101‑1 (MQ=255)
                                                                                                                                    gTCTTAGCTTCGGCCAGGCCGAAGAACCTCGAAGTGGCAAGGCTACGCAGTGGCTGGCGGGTGGTCAACCATTGCTGCGGGGCCGGTCGACAACGCTGGGCGCAAGTGCGTTACa                                             >  1:2315189/1‑115 (MQ=255)
                                                                                                                                    gTCTTAGCTTCGGCCAGGCCGAAGAACCTCGAAGTGGCAAGGCTACGCAGTGGCTGGCGGGTGGTCAACCATTGCTGCGGGGCCGGTCGACAACGCTGGGCGCAAGTGCGTTACa                                             <  2:2315189/115‑1 (MQ=255)
                                                                                                                                       ttAGCTTCGGCCAGGCCGAAGAACCTCGAAGTGGCAAGGCTACGCAGTGGCTGGCGGGTGGTCAACCATTGCTGCGGGGCCGGTCGACAACGCTGGGCGCAAGTGCGTTACAGTGGCGCCTTATCCACGGTCAAGGACGCACCATGcccc       >  1:1254356/1‑150 (MQ=255)
                                                                                                                                        tAGCTTCGGCCAGGCCGAAGAACCTCGAAGTGGCAAGGCTACGCAGTGGCTGGCGGGTGGTCAACCATTGCTGCGGGGCCGGTCGACAACGCTGGGCGCAAGTGCGTTACAGTGGCGCCTTATCCACGGTCAAGGACGCACCATGCCCCt      >  2:1016491/1‑150 (MQ=255)
                                                                                                                                        gttcTTCGGCCTGGCCGAAGAACCTCGAAGTGGCAAGGCTACGCAGTGGCTGGCGGGTGGTCAACCATTGCTGCGGGGCCGGTCGACAACGCTGGGCGCAAGTGCGTTACAGTGGCGCCTTATCCACGGTCAAGGACGCACCATGCCCCt      >  2:948310/4‑150 (MQ=255)
                                                                                                                                         aGCTTCGGCCAGGCCGAAGAACCTCGAAGTGGCAAGGCTACGCAGTGGCTGGCGGGTGGTCAACCATTGCTGCGGGGCCGGt                                                                         >  2:40265/1‑82 (MQ=255)
                                                                                                                                         aGCTTCGGCCAGGCCGAAGAACCTCGAAGTGGCAAGGCTACGCAGTGGCTGGCGGGTGGTCAACCATTGCTGCGGGGCCGGt                                                                         <  1:40265/82‑1 (MQ=255)
                                                                                                                                         aGCTTCGGCCAGGCCGAAGAACCTCGAAGTGGCAAGGCTACGCAGTGGCTGGCGGGTGGTCAACCATTGCTGCGGGGCCGGTCGACAACGCTGGGCGCAAGTGCGTTACAGTGGCGCCTTATCCACGGTCAAGGACGCACCATGCCCCTc     <  2:1468274/150‑1 (MQ=255)
                                                                                                                                           cTTCGGCCAGGCCGAAGAACCTCGAAGTGGCAAGGCTACGCAGTGGCTGGCGGGTGGTCAACCATTGCTGCGGGGCCGGTCGACAACGCTGGGCGCAAGTGCGTTACAGTGGCGCCTTATCCACGGTCaa                       >  2:1955246/1‑130 (MQ=255)
                                                                                                                                           cTTCGGCCAGGCCGAAGAACCTCGAAGTGGCAAGGCTACGCAGTGGCTGGCGGGTGGTCAACCATTGCTGCGGGGCCGGTCGACAACGCTGGGCGCAAGTGCGTTACAGTGGCGCCTTATCCACGGTCaa                       <  1:1955246/130‑1 (MQ=255)
                                                                                                                                            ttCGGCCTGGCCGAAGAACCTCGAAGTGGCAAGGCTACGCAGTGGCTGGCGGGTGGTCAACCATTGCTGCGGGGCCGGTCGACAACGCTGGgcgc                                                         <  1:702434/95‑1 (MQ=255)
                                                                                                                                            ttCGGCCTGGCCGAAGAACCTCGAAGTGGCAAGGCTACGCAGTGGCTGGCGGGTGGTCAACCATTGCTGCGGGGCCGGTCGACAACGCTGGgcgc                                                         >  2:702434/1‑95 (MQ=255)
                                                                                                                                            ttCGGCCAGGCCGAAGAACCTCGAAGTGGCAAGGCTACGCAGTGGCTGGCGGGTGGTCAACCATTGCTGCGGGGCCGGTCGACAACGCTGGGCGCAAGTGCGTTACAGTGGCGCCTTATCCACGGTCAAGGACGCACCATGCCCCTCGcc  >  2:2156594/1‑150 (MQ=255)
                                                                                                                                            ttCGGCCAGGCCGAAGAAACGCGAAGTGGCAAGGCTACGCGGTGGCTGGCGAGGGGGCAACCATTGCTTCGGGGCCGGTCGACAACGCTGGGCGCAAGTGCGTTACAGTGGCGCCTTATCCACGGTCAAGGACGCACCATGCCCCTCGcc  >  2:326889/1‑150 (MQ=255)
                                                                                                                                                gCCAGGCCGAAGAACCTCGAAGTGGCAAGGCTACGCAGTGGCTGGCGGGTGGTCAACCATTGCTGCGGGGCCGGTCGACAACGCTGGGCGCAAGTGCGTTACAGTGGCGCCTTATCCACGGTCaa                       <  2:55527/125‑1 (MQ=255)
                                                                                                                                                gCCAGGCCGAAGAACCTCGAAGTGGCAAGGCTACGCAGTGGCTGGCGGGTGGTCAACCATTGCTGCGGGGCCGGTCGACAACGCTGGGCGCAAGTGCGTTACAGTGGCGCCTTATCCACGGTCaa                       >  1:55527/1‑125 (MQ=255)
                                                                                                                                                   |                                                                                                                                              
TGGCTTGCACCGTTCCATGTCCCCGTTCGGGCTGGCGGGATGACGCCTTCCCGCGCGTTGGCGCAGTGGCTGTGTCGATCCGTTGTCGGTGCTTACAGTTGCGGGGGCAGCCGCGGCTTGTACCGCGTTCCCGTCTTAGCTTCGGCCAGGCCGAAGAACCTCGAAGTGGCAAGGCTACGCAGTGGCTGGCGGGTGGTCAACCATTGCTGCGGGGCCGGTCGACAACGCTGGGCGCAAGTGCGTTACAGTGGCGCCTTATCCACGGTCAAGGACGCACCATGCCCCTCGCC  >  NC_002947/3981874‑3982163

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