Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,982,021 A→T 7.2% 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 7.2% 250.1 / 8.7 83noncoding (28/207 nt)PP_mr43Cobalamin regulatory region
Reads supporting (aligned to +/- strand):  ref base A (44/33);  new base T (2/4);  total (46/37)
Fisher's exact test for biased strand distribution p-value = 4.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.86e-01

CTGGCTTGCACCGTTCCATGTCCCCGTTCGGGCTGGCGGGATGACGCCTTCCCGCGCGTTGGCGCAGTGGCTGTGTCGATCCGTTGTCGGTGCTTACAGTTGCGGGGGCAGCCGCGGCTTGTACCGCGTTCCCGTCTTAGCTTCGGCCAGGCCGAAGAACCTCGAAGTGGCAAGGCTACGCAGTGGCTGGCGGGTGGTCAACCATTGCTGCGGGGCCGGTCGACAACGCTGGGCGCAAGTGCGTTACAGTGGCGCCTTATCCACGGTCAAGGACGCACCATGCCCCTCGCCCCCG  >  NC_002947/3981873‑3982167
                                                                                                                                                    |                                                                                                                                                  
cTGGCTTGCACCGTTCCATGTCCCCGTTCGGGCTGGCGGGATGACGCCTTCCCGCGCGTTGGCGCAGTGGCTGTGTCGATCCGTTGTCAGTGCTTACAGTTGCGGGGGCAGCCGCGGCTTGTACCGCGTTCCCGTCTTAGCTTCGGCCAg                                                                                                                                                   >  2:2141656/1‑150 (MQ=255)
 tGGCTTGCACCGTTCCATGTCCCCGTTCGGGCTGGCGGGATGACGCCTTCCCGCGCGTTGGCGCAGTGGCTGTGTCGATCCGTTGTCGGTGCTTACAGTTGCGGGGGCAGCCGCGGCTTGTACCGCGTTCCCGTCTTAGCTTCGGCCAgg                                                                                                                                                  >  1:1437773/1‑150 (MQ=255)
   gCTTGCACCGTTCCATGTCCCCGTTCGGGCTGGCGGGATGACGCCTTCCCGCGCGTTGGCGCAGTGGCTGTGTCGATCCGTTGTCGGTGCTTACAGTTGCGGGGGCAGCCGCGGCTTGTACCGCGTTCCCGTCTTAGCTTCGGCCAGGcc                                                                                                                                                >  2:423601/1‑150 (MQ=255)
     ttGCACCGTTCCATGTCCCCGTTCGGGCTGGCGGGATGACGCCTTCCCGCGCGTTGGCGCAGTGGCTGTGTCGATCCGTTGTCGGTGCTTACAGTTGCGGGGGCAGCCGCGGCTTGTACCGCGTTCCCGTCTTAGCTTCGGCCAGGCCga                                                                                                                                              >  2:1812219/1‑150 (MQ=255)
       gCACCGTTCCATGTCCCCGTTCGGGCTGGCGGGATGACGCCTTCCCGCGCGTTGGCGCAGTGGCTGTGTCGATCCGTTGTCGGTGCTTACAGTTGCGGGGGCAGCCGCGGCTTGTACCGCGTTCCCGTCTTAGCTTCGGCCTGGCCgaag                                                                                                                                            <  2:2083981/150‑1 (MQ=255)
       gCACCGTTCCATGTCCCCGTTCGGGCTGGCGGGATGACGCCTTCCCGCGCGTTGGCGCAGTGGCTGTGTCGATCCGTTGTCGGTGCTTACAGTTGCGGGGGCAGCCGCGGCTTGTACCGCGTTCCCGTCTTAGCTTCGGCCTGGCCgaag                                                                                                                                            <  1:1015511/150‑1 (MQ=255)
       gCACCGTTCCATGTCCCCGTTCGGGCTGGCGGGATGACGCCTTCCCGCGCGTTGGCGCAGTGGCTGTGTCGATCCGTTGTCGGTGCTTACAGTTGCGGGGGCAGCCGCGGCTTGTACCGCGTTCCCGTCTTAGCTTCGGCCAGGCCgaag                                                                                                                                            >  2:1057907/1‑150 (MQ=255)
         aCCGTTCCATGTCCCCGTTCGGGCTGGCGGGATGACGCCTTCCCGCGCGTTGGCGCAGTGGCTGTGTCGATCCGTTGTCGGTGCTTACAGTTGCGGGGGCAGCCGCGGCTTGTACCGCGTTCCCGTCTTAGCTTCGGCCAGGCCgaagaa                                                                                                                                          >  1:555333/1‑150 (MQ=255)
           cGTTCCATGTCCCCGTTCGGGCTGGCGGGATGACGCCTTCCCGCGCGTTGGCGCAGTGGCTGTGTCGATCCGTTGTCGGTGCTTACAGTTGCGGGGGCAGCCGCGGCTTGTACCGCGTTCCCGTCTTAGCTTCGGCCTGGCCgaagctaa                                                                                                                                        <  2:370195/150‑5 (MQ=255)
             ttCCATGTCCCCGTTCGGGCTGGCGGGATGACGCCTTCCCGCGCGTTGGCGCAGTGGCTGTGTCGATCCGTTGTCGGTGCTTACAGTTGCGGGGGCAGCCGCGGCTTGTACCGCGTTCCCGTCTTAGCTTCGGCCAGGCCGAAGAACCTc                                                                                                                                      >  2:1064919/1‑150 (MQ=255)
              tCCATGTCCCCGTTCGGGCTGGCGGGATGACGCCTTCCCGCGCGTTGGCGCAGTGGCTGTGTCGATCCGTTGTCGGTGCTTACAGTTGCGGGGGCAGCCGCGGCTTGTACCGCGTTCCCGTCTTAGCTTCGGCCAGGCCGAAGAACCTCg                                                                                                                                     <  2:2348025/150‑1 (MQ=255)
                cATGTCCTCGTTCGGGCTGGCGGGATGACGCCTTCCCGCGCGTTGGCGCAGTGGCTGTGTCGATCCGTTGTCGGTGCTTACAGTTGCGGGGGCAGCCGCGGCTTGCACCGCTTCCCCGTCTTAGCTTCGGCCAGGCCGAAGAACCTCGaa                                                                                                                                   <  2:1668043/150‑1 (MQ=255)
                cATGTCCCCGTTCGGGCTGGCGGGATGACGCCTTCCCGCGCGTTGGCGCAGTGGCTGTGTCGATCCGTTGTCGGTGCTTACAGTTGCGGGGGCAGCCGCGGCTTGTACCGCGTTCCCGTCTTAGCTTCGGCCAGGCCGAAGAACCTCGaa                                                                                                                                   <  1:1056565/150‑1 (MQ=255)
                 aTGTCCCCGTTCGGGCTGGCGGGATGACGCCTTCCCGCGCGTTGGCGCAGTGGCTGTTTCGATCCGTTGTCGGTGCTTACAGTTGCGGGGGCAGCCGCGGCTTGTACCGCGTTCCCGTCTTAGCTTCGGCCAGGCCGAAGAACCTCGAAg                                                                                                                                  >  2:823092/1‑150 (MQ=255)
                    tCCCCGTTCGGGCTGGCGGGATGACGCCTTCCCGCGCGTTGGCGCAGTGGCTGTGTCGATCCGTTGTCGGTGCTTACAGTTGCGGGGGCAGCCGCGGCTTGTACCGCGTTCCCGTCTTAGCTTCGGCCAGGCCGAAGAACCTCGAAGTgg                                                                                                                               >  1:870716/1‑150 (MQ=255)
                    tCCCCGTTCGGGCTGGCGGGATGACGCCTTCCCGCGCGTTGGCGCAGTGGCTGTGTCGATCCGTTGTCGGTGCTTACAGTTGCGGGGGCAGCCGCGGCTTGTACCGCGTTCCCGTCTTAGCTTCGGCCAGGCCGAAGAACCTCGAAGTgg                                                                                                                               <  1:131111/150‑1 (MQ=255)
                     ccccGTTCGGGCTGGCGGGATGACGCCTTCCCGCGCGTTGGCGCAGTGGCTGTGTCGATCCGTTGTCGGTGCTTACAGTTGCGGGGGCAGCCGCGGCTTGTACCGCGTTCCCGTCTTAGCTTCGGCCAGGCCGAAGAACCTCGAAGTGGc                                                                                                                              >  1:1728738/1‑150 (MQ=255)
                      cccGTTCGGGCTGGCGGGATGACGCCTTCCCGCGCGTTGGCGCAGTGGCTGTGTCGATCCGTTGTCGGTGCTTACAGTTGCGGGGGCAGCCGCGGCTTGTACCGCGTTCCCGTCTTAGCTTCGGCCAGGCCGAAGAACCTCGCAGTGGCa                                                                                                                             >  1:1119063/1‑150 (MQ=255)
                                   gCGGGATGACGCCTTCCCGCGCGTTGGCGCAGTGGCTGTGTCGATCCGTTGTCGGTGCTTACAGTTGCGGGGGCAGCCGCGGCTTGTACCGCGTTCCCGTCTTAGCTTCGGCCAGGCCGAAGAACCTCGAAGTGGCAAGGCTACGCAGtg                                                                                                                <  2:768076/150‑1 (MQ=255)
                                            cGCCTTCCCGCGCGTTGGCGCAGTGGCTGTGTCGATCCGTTGTCGGTGCTTACAGTTGCGGGGGCAGCCGCGGCTTGTACCGCGTTCCCGTCTTAGCTTCGGCCAGGCCGAAGAACCTCGAAGTGGCAAGGCTACGCAGTGGCTGGCggg                                                                                                       >  2:1400216/1‑150 (MQ=255)
                                             gCCTTCCCGCGCGTTGGCGCAGTGGCTGTGTCGATCCGTTGTCGGTGCTTACAGTTGCGGGGGCAGCCGCGGCTTGTACCGCGTTCCCGTCTTAGCTTCGGCCAGGCCGAAGAACCTCGAAGTGGCaa                                                                                                                            >  2:416597/1‑128 (MQ=255)
                                             gCCTTCCCGCGCGTTGGCGCAGTGGCTGTGTCGATCCGTTGTCGGTGCTTACAGTTGCGGGGGCAGCCGCGGCTTGTACCGCGTTCCCGTCTTAGCTTCGGCCAGGCCGAAGAACCTCGAAGTGGCaa                                                                                                                            <  1:416597/128‑1 (MQ=255)
                                             gCCTTCCCGCGCGTTGGCGCAGTGGCTGTGTCGATCCGTTGTCGGTGCTTACAGTTGCGGGGGCAGCCGCGGCTTGTACCGCGTTCCCGTCTTAGCTTCGGCCAGGCCGAAGAACCTCGAAGTGGCAAGGCTACGCAGTGGCTGGCGggt                                                                                                      >  1:381547/1‑150 (MQ=255)
                                              ccTTCCCGCGCGTTGGCGCAGTGGCTGTGTCGATCCGTTGTCGGTGCTTACAGTTGCGGGGGCAGCCGCGGCTTGTACCGCGTTCCCGTCTTAGCTTCGGCCAGGCCGAAGAACCTCGAAGTGGCAAGGCTACGCAGTGGCTGGCGggtg                                                                                                     <  2:919661/150‑1 (MQ=255)
                                              ccTTCCCGCGCGTTGGCGCAGTGGCTGTGTCGATCCGTTGTCGGTGCTTACAGTTGCGGGGGCAGCCGCGGCTTGTACCGCGTTCCCGTCTTAGCTTCGGCCAGGCCGAAGAACCTCGAAGTGGCAAGGCTACGCAGTGGCTGGCGggtg                                                                                                     <  2:381547/150‑1 (MQ=255)
                                               cTTCCCGCGCGTTGGCGCAGTGGCTGTGTCGATCCGTTGTCGGTGCTTACAGTTGCGGGGGCAGCCGCGGCTTGTACCACGTTCCCGTCTTAGCTTCGGCCAGGCCGAAGAACCTCGAAGTGGCAAGGCTACGCAGTGGCTGGCGggtgg                                                                                                    >  2:805506/1‑150 (MQ=255)
                                                 tCCCGCGCGTTGGCGCAGTGGCTGTGTCGATCCGTTGTCGGTGCTTACAGTTGCGGGGGCAGCCGCGGCTTGTACCGCGTTCCCGTCTTAGCTTCGGCCAGGCCGAAGAACCTCGAAGTGGCAAGGCTACGCAGTGGCTGGCGGGTGGTc                                                                                                  <  1:1064919/150‑1 (MQ=255)
                                                   ccGCGCGTTGGCGCAGTGGCTGTGTCGATCCGTTGTCGGTGCTTACAGTTGCGGGGGCAGCCGCGGCTTGTACCGCGTTCCCGTCTTAGCTTCGGCCAGGCCGAAGAACCTCGAAGTGGCAAGGCTACGCAGTGGCTGGCGGGTGGTCaa                                                                                                >  1:1573018/1‑150 (MQ=255)
                                                   ccGCGCGTTGGCGCAGTGGCTGTGTCGATCCGTTGTCGGTGCTTACAGTTGCGGGGGCAGCCGCGGCTTGTACCGCGTTCCCGTCTTAGCTTCGGCCAGGCCGAAGAACCTCGAAGTGGCAAGGCTACGCAGTGGCTGGCGGGTGGTCaa                                                                                                >  1:1092930/1‑150 (MQ=255)
                                                   ccGCGCGTTGGCGCAGTGGCTGTGTCGATCCGTTGTCGGTGCTTACAGTTGCGGGGGCAGCCGCGGCTTGTACCGCGTTCCCGTCTTAGCTTCGGCCAGGCCGAAGAACCTCGAAGTGGCAAGGCTACGCAGTGGCTGGCGGGTGGTCaa                                                                                                >  1:2250785/1‑150 (MQ=255)
                                                    cgcgcgTTGGCGCAGTGGCTGTGTCGATCCGTTGTCGGTGCTTACAGTTGCGGGGGCAGCCGCGGCTTGTACCGCGTTCCCGTCTTAGCTTCGGCCAGGCCGAAGAACCTCGAAGTGGCAAGGCTACGCAGTGGCTGGCGGGTGGTCAAc                                                                                               <  1:311074/150‑1 (MQ=255)
                                                        cgTTGGCGCAGTGGCTGTGTCGATCCGTTGTCGGTGCTTACAGTTGCGGGGGCAGCCGCGGCTTGTACCGCGTTCCCGTCTTAGCTTCGGCCAGGCCGAAGAACCTCGAAGTGGCAAGGCTACGCAGTGGCTGGCGGGTGGTCAACCAtt                                                                                           <  1:1997016/150‑1 (MQ=255)
                                                             gcgcAGTGGCTGTGTCGATCCGTTGTCGGTGCTTACAGTTGCGGGGGCAGCCGCGGCTTGTACCGCGTTCCCGTCTTAGCTTCGGCCAGGCCGAAGAACCTCGAAGTGGCAAGGCTACGCAGtggctgg                                                                                                           >  1:2284038/1‑129 (MQ=255)
                                                             gcgcAGTGGCTGTGTCGATCCGTTGTCGGTGCTTACAGTTGCGGGGGCAGCCGCGGCTTGTACCGCGTTCCCGTCTTAGCTTCGGCCAGGCCGAAGAACCTCGAAGTGGCAAGGCTACGCAGtggctgg                                                                                                           <  2:2284038/129‑1 (MQ=255)
                                                               gcAGTGGCTGTGTCGATCCGTTGTCGGTGCTTACAGTTGCGGGGGCAGCCGCGGCTTGTACCGCGTTCCCGTCTTAGCTTCGGCCAGGCCGAAGAACCTCGAAGTGGCAAGGCTACGCAGTGGCTGGCGGGTGGTCAACCATTGCTGCgg                                                                                    <  1:423601/150‑1 (MQ=255)
                                                               gcAGTGGCTGTGTCGATCCGTTGTCGGTGCTTACAGTTGCGGGGGCAGCCGCGGCTTGTACCGCGTTCCCGTCTTAGCTTCGGCCAGGCCGAAGAACCTCGAAGTGGCAAGGCTACGCAGTGGCTGGCGGGTGGTCAACCATTGCTGCgg                                                                                    >  2:1326653/1‑150 (MQ=255)
                                                                cAGTGGCTGTGTCGATCCGTTGTCGGTGCTTACAGTTGCGGGGGCAGCCGCGGCTTGTACCGCGTTCCCGTCTTAGCTTCGGCCTGGCCgaagc                                                                                                                                           <  1:1701536/94‑2 (MQ=255)
                                                                cAGTGGCTGTGTCGATCCGTTGTCGGTGCTTACAGTTGCGGGGGCAGCCGCGGCTTGTACCGCGTTCCCGTCTTAGCTTCGGCCTGGCCgaagc                                                                                                                                           >  2:1701536/1‑93 (MQ=255)
                                                                 aGTGGCTGTGTCGATCCGTTGTCGGTGCTTACAGTTGCGGGGGCAGCCGCGGCTTGTACCGCGTTCCCGTCTTAGCTTCGGCCAGGCCGAAGAACCTCGAAGTGGCAAGGCTACGCAGTGGCTGGCGGGTGGTCAACCATTGCTGCgggg                                                                                  >  1:454201/1‑150 (MQ=255)
                                                                     gCTGTGTCGATCCGTTGTCGGTGCTTACAGTTGCGGGGGCAGCCGCGGCTTGTACCGCGTTCCCGTCTTAGCTTCGGCCAGGCCGAAGAACCTCGAAGTGGCAAGGCTACGCAGTGGCTGGCGGGTGGTCAACCATTGCTGCGGGGCCgg                                                                              >  1:1793520/1‑150 (MQ=255)
                                                                      cTGTGTCGATCCGTTGTCGGTGCTTACAGTTGCGGGGGCAGCCGCGGCTTGTACCGCGTTCCCGTCTTAGCTTCGGCCAGGCCGAAGAACCTCGAAGTGGCAAGGCTACGCAGt                                                                                                                 <  1:974162/114‑1 (MQ=255)
                                                                      cTGTGTCGATCCGTTGTCGGTGCTTACAGTTGCGGGGGCAGCCGCGGCTTGTACCGCGTTCCCGTCTTAGCTTCGGCCAGGCCGAAGAACCTCGAAGTGGCAAGGCTACGCAGt                                                                                                                 >  2:974162/1‑114 (MQ=255)
                                                                         tgtCGATCCGTTGTCGGTGCTTACAGTTGCGGGGGCAGCCGCGGCTTGTACCGCGTTCCCGTCTTAGCTTCGGCCAGGCCGAAGAACCTCGAAGTGGCAAGGCTACGCAGTGGCTGGCGGGTGGTCAACCATTGCTGCGGGGCCGGTCGa                                                                          <  2:1437773/150‑1 (MQ=255)
                                                                         tgtCGATCCGTTGTCGGTGCTTACAGTTGCGGGGGCAGCCGCGGCTTGTACCGCGTTCCCGTCTTAGCTTCGGCCAGGCCGAAGAACCTCGAAGTGGCAAGGCTACGCAGTGGCTGGCGGGTGGTCAACCATTGCTGCGGGGCCGGTCGa                                                                          <  2:1573018/150‑1 (MQ=255)
                                                                             gATCCGTTGTCGGTGCTTACAGTTGCGGGGGCAGCCGCGGCTTGTACCGCGTTCCCGTCTTAGCTTCGGCCAGGCCGAAGAAc                                                                                                                                         >  2:2257199/1‑83 (MQ=255)
                                                                             gATCCGTTGTCGGTGCTTACAGTTGCGGGGGCAGCCGCGGCTTGTACCGCGTTCCCGTCTTAGCTTCGGCCAGGCCGAAGAAc                                                                                                                                         <  1:2257199/83‑1 (MQ=255)
                                                                             gATCCGTTGTCGGTGCTTACAGTTGCGGGGGCAGCCGCGGCTTGTACCGCGTTCCCGTCTTAGCTTCGGCCAGGCCGAAGAACCTCGAAGTGGCAAGGCTACGCAGTGGCTGGCGGGTGGTCAACCATTGCTGCGGGGCCGGTCGACAAc                                                                      >  1:7006/1‑150 (MQ=255)
                                                                               tacGTTGTCGGTGCTTACAGTTGCGGGGTCAGCCGCGGCTTGTACCGCGTTCCCGTCTTAGCTTCGGCCAGGCCGAAGAACCTCGCAGCTGCGAAACAACGCAGAACCGTCCGAGGGGCCAACCATTGCGTCGGGGCCGGTCTACAACGc                                                                    >  2:2330200/3‑150 (MQ=255)
                                                                                 cGTTGTCGGTGCTTACAGTTGCGGGGGCAGCCGCGGCTTGTACCGCGTTCCCGTCTTAGCTTCGGCCAGGCCGAAGAACCTCGAAGTGGCAAGGCTACGCAGTGGCTGGCGGGTGGTCAACCATTGCTGCGGGGCCGGTCGACAACGCTg                                                                  <  2:2250785/150‑1 (MQ=255)
                                                                                  gTTGTCGGTGCTTCCACATGCAGTTGCGGCCGCGGCTTGTAGCGCGTTCCCGTCTTAGCTTCGGCCAGGCCGAAGAACCTCGAAGTGGCAAGGCTACGCAGTGGCTGGCGGGTGGTCAACCATTGCTGCGGGGCCGGTCGACAACGCTgg                                                                 >  2:594039/1‑150 (MQ=255)
                                                                                   ttGTCGGTGCTTACAGTTGCGGGGGCAGCCGCGGCTTGTACCGCGTTCCCGTCTTAGCTTCGGCCAGGCCGAAGAACCTCGAAGTGGCAAGGCTACGCAGTGGCTGGCGggtg                                                                                                     >  2:1857252/1‑113 (MQ=255)
                                                                                   ttGTCGGTGCTTACAGTTGCGGGGGCAGCCGCGGCTTGTACCGCGTTCCCGTCTTAGCTTCGGCCAGGCCGAAGAACCTCGAAGTGGCAAGGCTACGCAGTGGCTGGCGggtg                                                                                                     <  1:1857252/113‑1 (MQ=255)
                                                                                        ggTGCTTACAGTTGCGGGGGCAGCCGCGGCTTGTACCGCGTTCCCGTCTTAGCTTCGGCCAGGCCGAAGAACCTCGAAGTGGCAAGGCTACGCAGTGGCTGGCGGGTGGTCAACCATTGCTGCGGGGCCGGTCGACAACGCTGGGCGCaa                                                           >  1:1765626/1‑150 (MQ=255)
                                                                                         gTGCTTACAGTTGCGGGGGCAGCCGCGGCTTGTACCGCGTTCCCGTCTTAGCTTCGGCCAGGCCGAAGAACCTCGAAGTGGCAAGGCTACGCAGTGGCTGGCGGGTGGTCAACCATTGCTGCGGGGCCGGTCGACAACGCTGGGCGCAAg                                                          >  1:680508/1‑150 (MQ=255)
                                                                                              tACAGTTGCGGGGGCAGCCGCGGCTTGTACCGCGTTCCCGTCTTAGCTTCGGCCAGGCCGAAGAACCTCGAAGTGGCAAGGCTACGCAGTGGCTGGCGGGTGGTCAACCATTGCTGCGGGGCCGGTCGACAACGCTGGGCGCAAGTGCGt                                                     <  2:680508/150‑1 (MQ=255)
                                                                                                  gTTGCGGGGGCAGCCGCGGCTTGTACCGCGTTCCCGTCTTAGCTTCGGCCAGGCCGAAGAACCTCGAAGTGGCAAGGCTACGCAGTGGCTGGCGGGTGGTCAACCATTGCTGCGGGGCCGGTCGACAACGCTGGGCGCAAGTGCGTTACa                                                 <  1:1400216/150‑1 (MQ=255)
                                                                                                  gTTGCGGGGGCAGCCGCGGCTTGTACCGCGTTCCCGTCTTAGCTTCGGCCAGGCCGAAGAACCTCGAAGTGGCAAGGCTACGCAGTGGCTGGCGGGTGGTCAACCATTGCTGCGGGGCCGGTCGACAACGCTGGGCGCAAGTGCGTTACa                                                 <  2:1350487/150‑1 (MQ=255)
                                                                                                       gggggCAGCCGCGGCTTGTACCGCGTTCCCGTCTTAGCTTCGGCCAGGCCGAAGAACCTCGAAGTGGCAAGGCTACGCAGTGGCTGGCGGGTGGTCAACCATTGCTGCGGGGCCGGTCGACAACGCTGGGCGCAAGTGCGTTACAGt                                               <  2:409408/147‑1 (MQ=255)
                                                                                                       gggggCAGCCGCGGCTTGTACCGCGTTCCCGTCTTAGCTTCGGCCAGGCCGAAGAACCTCGAAGTGGCAAGGCTACGCAGTGGCTGGCGGGTGGTCAACCATTGCTGCGGGGCCGGTCGACAACGCTGGGCGCAAGTGCGTTACAGt                                               >  1:409408/1‑147 (MQ=255)
                                                                                                                  cgGCTTGTACCGCGTTCCCGTCTTAGCTTCGGCCAGGCCGAAGAACCTCGAAGTGGCAAGGCTACGCAGTGGCTGGCGGGTGGTCAACCATTGCTGCGGGGCCGGTCGACAACGCTGGGCGCAAGTGCGTTACAGTGGCGCCTTATCCAc                                 <  1:1057907/150‑1 (MQ=255)
                                                                                                                  cgGCTTGTACCGCGTTCCCGTCTTAGCTTCGGCCAGGCCGAAGAACCTCGAAGTGGCAAGGCTACGCAGTGGCTGGCGGGTGGTCAACCATTGCTGCGGGGCCGGTCGACAACGCTGGGCGCAAGTGCGTTACAGTGGCGCCTTATCCAc                                 <  2:1092930/150‑1 (MQ=255)
                                                                                                                     cTTGTACCGCGTTCCCGTCTTAGCTTCGGCCAGGCCGAAGAACCTCGAAGTGGCAAGGCTACGCAGTGGCTGGCGGGTGGTCAACCATTGCTGCGGGGCCGGTCGACAACGCTGGGCGCAAGTGCGTTACAGTGGCGCCTTGTCCACGGt                              >  1:833367/1‑150 (MQ=255)
                                                                                                                         tACCGCGTTCCCGTCTTAGCTTCGGCCAGGCCGAAGAACCTCGAAGTGGCAAGGCTACGCAGTGGCTGGCGGGTGGTCAACCATTGCTGCGGGGCCGGTCGACAACGCTGGGCGCAAGTGCGTTACAGTGGCGCCTTATCCACGGTCAAg                          >  2:1248266/1‑150 (MQ=255)
                                                                                                                          aCCGCGTTCCCGTCTTAGCTTCGGCCAGGCCGAAGAACCTCGAAGTGGCAAGGCTACGCAGTGGCTGGCGGGTGGTCAACCATTGCTGCGGGGCCGGTCGACAACGCTGGGCGCAAGTGCGTTACAGTGGCGCCTTATCCACGGTCAAgg                         >  2:458305/1‑150 (MQ=255)
                                                                                                                           ccGCGTTCCCGTCTTAGCTTCGGCCAGGCCGAAGAACCTCGAAGTGGCAAGGCTACGCAGTGGCTGGCGGGTGGTCAACCATTGCTGCgg                                                                                    >  2:1641040/1‑90 (MQ=255)
                                                                                                                           ccGCGTTCCCGTCTTAGCTTCGGCCAGGCCGAAGAACCTCGAAGTGGCAAGGCTACGCAGTGGCTGGCGGGTGGTCAACCATTGCTGCGGGGCCGGTCGACAACGCTGGGCGCAAGTGCGTTACAGTGGCGCCTTATCCACGGTCaa                           >  1:1348998/1‑147 (MQ=255)
                                                                                                                           ccGCGTTCCCGTCTTAGCTTCGGCCAGGCCGAAGAACCTCGAAGTGGCAAGGCTACGCAGTGGCTGGCGGGTGGTCAACCATTGCTGCGGGGCCGGTCGACAACGCTGGGCGCAAGTGCGTTACAGTGGCGCCTTATCCACGGTCaa                           <  2:1348998/147‑1 (MQ=255)
                                                                                                                           ccGCGTTCCCGTCTTAGCTTCGGCCAGGCCGAAGAAACTCGAAGTGGCAAGGCTACGCAGTGGCTGGCGGGTGGTCAACCATTGCTGCgg                                                                                    <  1:1641040/90‑1 (MQ=255)
                                                                                                                            cgcgTTCCCGTCTTAGCTTCGGCCAGGCCGAAGAACCTCGAAGTGGCAAGGCTACGCAGTGGCTGGCGGGTGGTCAACCATTGCTGCGGGGCCGGTCGACAACGCTGGGCGCAAGTGCGTTACAGTGGCGCCTTATCCACGGTCAAGGAc                       >  1:843387/1‑150 (MQ=255)
                                                                                                                            cgcgTTCCCGTCTTAGCTTCGGCCAGGCCGAAGAACCTCGAAGTGGCAAGGCTACGCAGTGGCTGGCGGGTGGTCAACCATTGCTGCGGGGCCGGTCGACAACGCTGGGCGCAAGTGCGTTACAGTGGCGCCTTATCCACGGTCAAGGAc                       <  1:823092/150‑1 (MQ=255)
                                                                                                                              cgTTCCCGTCTTAGCTTCGGCCAGGCCGAAGAACCTCGAAGTGGCAAGGCTACGCAGTGGCTGGCGGGTGGTCAACCATTGCTGCGGGGCCGGTCGACAACGCTGGGCGCAAGTGCGTTACAGTGGCGCCTTATCCACGGTCAAGGACGc                     <  1:2330200/150‑1 (MQ=255)
                                                                                                                               gTTCCCGTCTTAGCTTCGGCCAGGCCGAAGAACCAAGAAGTGGGCGGGCTACGCGGTCTCTGGCGGGTGGTCAACCATTGCTGCGGGGCCGGTCGACAACGCTGGGAGCACGAGCGTTACAGTGGCGCTTTATCTACTGTCAAggtctca                    <  2:1728738/150‑6 (MQ=255)
                                                                                                                                   ccGTCTTAGCTTCGGCCAGGCCGAAGAACCTCGAAGTGGCAAGGCTACGCAGTGGCTGGCGGGTGGTCAACCATTGCTGCGGGGCCGGTCGACAACGCTGGGCGCAAGTGCGTTACAGTGGCGCCTTATCCACGGTCAAGGACGCACCAt                >  1:2179395/1‑150 (MQ=255)
                                                                                                                                   ccGTCTTAGCTTCGGCCAGGCCGAAGAACCTCGAAGTGGCAAGGCTACGCAGTGGCTGGCGGGTGGTCAACCATTGCTGCGGGGCCGGTCGACAACGCTGGGCGCAAGTGCGTTACAGTGGCGCCTTATCCACGGTCAAGGACGCACCAt                >  1:527743/1‑150 (MQ=255)
                                                                                                                                    cGTCTTAGCTTCGGCCAGGCCGAAGAACCTCGAAGTGGCAAGGCTACGCAGTGGCTGGCGGGTGGTCAACCATTGCTGCGGGGCCGGTCGACAACGCTGGGCGCAAGTGCGTTACACTGGCGCCTTATCCACGGTCAAGGACGCACCATg               >  1:2109586/1‑150 (MQ=255)
                                                                                                                                        ttAGCTTCGGCCAGGCCGAAGAACCTCGAAGTGGCAAGGCTACGCAGTGGCTGGCGGGTGGTCAACCATTGCTGCGGGGCCGGTCGACAACGCTGGGCGCAAGTGCGTTACAGTGGCGCCTTATCCACGGTCAAGGACGCACCATGcccc           <  2:1765626/150‑1 (MQ=255)
                                                                                                                                          ttcTTCGGCCTGGCCGAAGAACCTCGAAGTGGCAAGGCTACGCAGTGGCTGGCGGGTGGTCAACCATTGCTGCGGGGCCGGTCGACAACGCTGGGCGCAAGTGCGTTACAGTGGCGCCTTATCCACGGTCAAGGACGCACCATGCCCCTc         >  2:2029746/3‑150 (MQ=255)
                                                                                                                                          aGCTTCGGCCAGGCCGAAGAACCTCGAAGTGGCAAGGCTACGCAGTGGCTGGCGGGTGGTCAACCATTGCTGCGGGGCCGGTCGACAACGCTGGGCGCAAGTGCGTTACAGTGGCGCCTTATCCACGGTCAAGGACGCACCATGCCCCTc         >  1:429354/1‑150 (MQ=255)
                                                                                                                                          aGCTTCGGCCAGGCCGAAGAACCTCGAAGTGGCAAGGCTACGCAGTGGCTGGCGGGTGGTCAACCATTGCTGCGGGGCCGGTCGACAACGCTGGGCGCAAGTGCGTTACAGTGGCGCCTTATCCACGGTCAAGGACGCACCATGCCCCTc         <  2:843387/150‑1 (MQ=255)
                                                                                                                                            cTTCGGCCAGGCCGAAGAACCTCGAAGTGGCAAGGCTACGCAGTGGCTGGCGGGTGGTCAACCATTGCTGCGGGGCCGGTCGACAACGCTGGGCGCAAGTGCGTTACAGTGGCGCCTTATCCACGGTCAAGGACGCACCATGCCCCTCGc       >  1:1276847/1‑150 (MQ=255)
                                                                                                                                                ggCCAGGCCGAAGAACCTCGAAGTGGCAATGCTACGCAGTGGCTGGCGGGTGGTCAACCATTGCTGCGGGGCCGGTCGACAACGCTGGGCGCAAGTGCGTTACAGTGGCGCCTTATCCACGGTCAAGGACGCACCATGCCCCTCGccccc   <  1:1812219/150‑1 (MQ=255)
                                                                                                                                                 gCCAGGCCGAAGAACCTCGAAGTGGCAAGGCTACGCAGTGGCTGGCGGGTGGTCAACCATTGCTGCGGGGCCGGTCGACAACGCTGGGCGCAAGTGCGTTACAGTGGCGCCTTATCCACGGTCAAGGACGCACCATGCCCCTCGCCCCCg  >  1:395950/1‑150 (MQ=255)
                                                                                                                                                 gCCAGGCCGAAGAACCTCGAAGTGGCAAGGCTACGCAGTGGCTGGCGGGTGGTCAACCATTGCTGCGGGGCCGGTCGACAACGCTGGGCGCAAGTGCGTTACAGTGGCGCCTTATCCACGGTCAAGGACGCACCATGCCCCTCGCCCCCg  >  1:937792/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
CTGGCTTGCACCGTTCCATGTCCCCGTTCGGGCTGGCGGGATGACGCCTTCCCGCGCGTTGGCGCAGTGGCTGTGTCGATCCGTTGTCGGTGCTTACAGTTGCGGGGGCAGCCGCGGCTTGTACCGCGTTCCCGTCTTAGCTTCGGCCAGGCCGAAGAACCTCGAAGTGGCAAGGCTACGCAGTGGCTGGCGGGTGGTCAACCATTGCTGCGGGGCCGGTCGACAACGCTGGGCGCAAGTGCGTTACAGTGGCGCCTTATCCACGGTCAAGGACGCACCATGCCCCTCGCCCCCG  >  NC_002947/3981873‑3982167

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