Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,125,923 A→C 5.7% Q221P (CAG→CCG)  PP_5374 → choline/carnitine/betaine transporter family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029476,125,9230AC 5.7% 398.8 / 14.8 140Q221P (CAG→CCG) PP_5374choline/carnitine/betaine transporter family protein
Reads supporting (aligned to +/- strand):  ref base A (68/64);  new base C (4/4);  total (72/68)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.87e-01

Only 100 of 142 total aligned reads displayed.

CATGCGCTCGATCCTTTACCCGCTGATTGGTGAGCGCATCTATGGGCCCATCGGACATGTGGTGGACATCCTTGCCGTGGTGATTACCGCGTTCGGTGTATCCCAGTCGCTGGGCCTGGGCGTGGTGCAAATGAACACCGGGCTGAGCCAGGTGTTTGACCTGCCCATCAGCCTGGGGGTGCAGATCACCTTGATCGTGCTCATTACCCTCGTTACTACCGTATCGGTGATGGCCGGCGTGTCGCGTGGCATGAAGCGCCTTTCGGAATGGAATATGCTGTTGTCGGTGGTGCTGGTGG  >  NC_002947/6125774‑6126072
                                                                                                                                                     |                                                                                                                                                     
cATGCGCTCGATCCTTTACCCGCTGATTGGTGAGCGCATCTATGGGCCCATCGGACATGTGGTGGACATCCTTGCCGTGGTGATTACCGCGTTCGGTGTATCCCAGTCGCTGGGCCTGGGCGTGGTGCAAATGAACACCGGGCTGAGCCa                                                                                                                                                       >  1:923292/1‑150 (MQ=255)
 ctGCGCTCGATACTTTACCCGCTGATTGGTGAGCGCATCTATGGGCCCATCGGACATGCGGTGGACATTCTTGCCGTGGTGGTTCGCGCCTTCGGTGTATCCCAGTCGCTGGGCCTGGGCGTGGTGCAAATGAACACAGGGCTGAGCCAg                                                                                                                                                      <  2:198035/149‑1 (MQ=255)
   gcgcTCGATCCTTTACCCGCTGATTGGTGAGCGCATCTATGGGCCCATCGGACATGTGGTGGACATCCTTGCCGTGGTGATTACCGCGTTCGGTGTATCCCAGTCGCTGTGCCTGGGCGTGGTGCAAATGAACACCGGGCTGAGCCAg                                                                                                                                                      >  1:854543/1‑148 (MQ=255)
   gcgcTCGATCCTTTACCCGCTGATTGGTGAGCGCATCTATGGGCCCATCGGACATGTGGTGGACATCCTTGCCGTGGTGATTACCGCGTTCGGTGTATCCCAGTCGCTGGGCCTGGGCGTGGTGCAAATGAACACCGGGCTGAGCCAg                                                                                                                                                      <  2:854543/148‑1 (MQ=255)
    cgcTCGATCCTTTACCCGCTGATTGGTGAGCGCATCTATGGGCCCATCGGACATGTGGTGGACATCCTTGCCGTGGTGATTACCGCGTTCGGTGTATCCCAGTCGCTGGGCCTGGGCGTGGTGCAAATGAACACCTGGCTCAGCCCGgtg                                                                                                                                                   <  1:1770287/150‑1 (MQ=255)
     gcTCGATCCTTTCCCAGCTGATTGGTGAGCGCATCTATGGGCCCATCGGACATGTGGTGGACATCCTTGCCGTGGTGATTACCGCGTTCGGTGTATCCCAGTCGCTGGGCCTGGGCGTGGTGCAACTGAACACCGGGCTGAGCCAGgtgt                                                                                                                                                  <  2:277131/150‑1 (MQ=255)
     gcTCGATCCTTTACCCGCTGATTGGTGAGCGCATCTATGGGCCCATCGGACATGTGGTGGACATCCTTGCCGTGGTGATTACCGCGTTCGGTGTATCCCAGTCGCTGGGCCTGGGCGTGGTGCAAATGAACACCGGGCTGAGCCAGgtgt                                                                                                                                                  >  1:368755/1‑150 (MQ=255)
       tCGATCCTTTACCCGCTGATTGGTGAGCGCATCTATGGGCCCATCGGACATGTGGTGGACATCCTTGCCGTGGTGATTACCGCGTTCGGTGTATCCCAGTCGCTGGGCCTGGGCGTGGTGCAAATGAACACCGGGCTGAGCCAGGTGttt                                                                                                                                                >  1:420036/1‑150 (MQ=255)
         gATCCTTTACCCGCTGATTGGTGAGCGCATCTATGGGCCCATCGGACATGTGGTGGACATCCTTGCCGTGGTGATTACCGCGTTCGGTGTATCCCAGTCGCTGGGCCTGGGCGTGGTGCAAATGAACCCCGGGCTGAGCCAGGTGTTTGa                                                                                                                                              <  2:2348786/150‑1 (MQ=255)
              tttACCCGCTGATTGGTGAGCGCATCTATGGGCCCATCGGACATGTGGTGGACATCCTTGCCGTGGTGATTACCGCGTTCGGTGTATCCCAGTCGCTGGGCCTGGGCGTGGTGCAAATGAACACCGGGCTGAGCCAGGTGTTTGACCTGc                                                                                                                                         >  1:1551672/1‑150 (MQ=255)
              tttACCCGCTGATTGGTGAGCGCATCTATGGGCCCATCGGACATGTGGTGGACATCCTTGCCGTGGTGATTACCGCGTTCGGTGTATCCCAGTCGCTGGGCCTGGGCGTGGTGCAAATGAACACCGGGATGAGCCAGGAGTTTGACCTGc                                                                                                                                         >  1:1204131/1‑150 (MQ=255)
                   ccGCTGATTGGTGAGCGCATCTATGGGCCCATCGGACATGTGGTGGACATCCTTGCCGTGGTGATTACCGCGTTCGGTGTATCCCAGTCGCTGGGCCTGGGCGTGGTGCAAATGAACACCGGGCTGAGCCAg                                                                                                                                                      >  2:481436/1‑132 (MQ=255)
                   ccGCTGATTGGTGAGCGCATCTATGGGCCCATCGGACATGTGGTGGACATCCTTGCCGTGGTGATTACCGCGTTCGGTGTATCCCAGTCGCTGGGCCTGGGCGTGGTGCAAATGAACACCGGGCTGAGCCAg                                                                                                                                                      <  1:481436/132‑1 (MQ=255)
                     gCTGATTGGTGAGCGCATCTATGGGCCCATCGGACATGTGGTGGACATCCTTGCCGTGGTGATTACCGCGTTCGGTGTATCCCAGTCGCTGGGCCTGGGCGTGGTGCAAATGAACACCGGGCTGAGCCAGGTGTTTGACCTGCCCATCAg                                                                                                                                  >  2:969385/1‑150 (MQ=255)
                     gCTGATTGGTGAGCGCATCTATGGGCCCATCGGACATGTGGTGGACATCCTTGCCGTGGTGATTACCGCGTTCGGTGTATCCCAGTCGCTGGGCCTGGGCGTGGTGCAAATGAACACCGGGCTGAGCCAGGTGTTTGACCTGCCCATCAg                                                                                                                                  >  1:2167928/1‑150 (MQ=255)
                       tGATTGGTGAGCGCATCTATGGGCCCATCGGACATGTGGTGGACATCCTTGCCGTGGTGATTACCGCGTTCGGTGTATCCCAGTCGCTGGGCCTGGGCGTGGTGCAAATGAACACCGGGCTGAGCCAGGTGTTTGACCTGCCCATCAGcc                                                                                                                                <  2:2995440/150‑1 (MQ=255)
                         aTTGGTGAGCGCATCTATGGGCCCATCGGACATGTGGTGGACATCCTTGCCGTGGTGATTACCGCGTTCGGTGTATCCCAGTCGCTGGGCCTGGGCGTGGTGCAAATGAACACCGGGCTGAGCCAGGTGTTTGACCTGCCCATCAGCCTg                                                                                                                              >  1:327699/1‑150 (MQ=255)
                          ttGGTGAGCGCATCTATGGGCCCATCGGACATGTGGTGGACATCCTTGCCGTGGTGATTACCGCGTTCGGTGTATCCCAGTCGCTGGGCCTGGGCGTGGTGCAAATGAACACCGGGCTGAGCCAGGTGTTTGACCTGCCCATCAGCCTgg                                                                                                                             <  2:1599380/150‑1 (MQ=255)
                           tGGTGAGCGCATCTATGGGCCCATCGGACATGTGGTCGGCATCCTTGCCGTGGTGATTACCGCGTTCGGTGTATCCCAGTCGCTGGGCCTGGGCGTGGTGCAAATGAACACCGGGCTGAGCCAGGTGTTTGACCTGCCCATCAGCCTggg                                                                                                                            <  1:1411828/150‑1 (MQ=255)
                            ggTGAGCGCATCTATGGGCCCATCGGACATGTGGTGGACATCCTTGCCGTGGTGATTACCGCGTTCGGTGTATCCCAGTCGCTGGGCCTGGGCGTGGTGCAAATGAACACCGGGCTGAGCCAGGTGTTTGACCTGCCCATCAGCCTgggg                                                                                                                           <  2:1080679/150‑1 (MQ=255)
                             gTGAGCGCATCTATGGGCCCATCGGACATGTGGTGGACATCCTTGCCGTGGTGATTACCGCGTTCGGTGTATCCCAGTCGCTGGGCCTGGGCGTGGTGCAAATGAACACCGGGCTGAGCCAGGTGTTTGACCTGCCCATCAGCCTggggg                                                                                                                          >  1:549923/1‑150 (MQ=255)
                                         aTGGGCCCATCGGACATGTGGTGGACATCCTTGCCGTGGTGATTACCGCGTTCGGTGTATCCCAGTCGCTGGGCCTGGGCGTGGTGCAAATGAACACCGGGCTGAGCCAGGTGTTTGACCTGCCCATCAGCCTGGGGGTGCAGATCACCt                                                                                                              <  1:2220432/150‑1 (MQ=255)
                                          tGGGCCCATCGGACATGTGGTGGACATCCTTGCCGTGGTGATTACCGCGTTCGGTGTATCCCAGTCGCTGGGCCTGGGCGTGGTGCAAATGAACACCGGGCTGAGCCAGGTGTTTGACCTGCCCATCAGCCTGGGGGTGCAGATCACCtt                                                                                                             >  2:1913040/1‑150 (MQ=255)
                                                 aTCGGACATGTGGTGGACATCCTTGCCGTGGTGATTACCGCGTTCGGTGTATCCCCGTCGCTGGGCCTGGGCGTGGTGCAAATGAACACCGGGCTGAGCCAGGTGTTTGACCTGc                                                                                                                                         <  2:271986/115‑1 (MQ=255)
                                                 aTCGGACATGTGGTGGACATCCTTGCCGTGGTGATTACCGCGTTCGGTGTATCCCAGTCGCTGGGCCTGGGCGTGGTGCAAATGAACACCGGGCTGAGCCAGGTGTTTGACCTGCCCATCAGcc                                                                                                                                >  2:732170/1‑124 (MQ=255)
                                                 aTCGGACATGTGGTGGACATCCTTGCCGTGGTGATTACCGCGTTCGGTGTATCCCAGTCGCTGGGCCTGGGCGTGGTGCAAATGAACACCGGGCTGAGCCAGGTGTTTGACCTGCCCATCAGCCTGGGGGTGCAGATCACCTTGATCGTg                                                                                                      <  1:260920/150‑1 (MQ=255)
                                                  tCGGACATGTGGTGGACATCCTTGCCGTGGTGATTACCGCGTTCGGTGTATCCCAGTCGCTGGGCCTGGGCGTGGTGCAAATGAACACCGGGCTGAGCCAGGTGTTTGACCTGCCCATCAGCCTGGGGGTGCAGATCACCTTGATCGTGc                                                                                                     >  2:2660173/1‑150 (MQ=255)
                                                        aTGTGGTGGACATCCTTGCCGTGGTGATTACCGCGTTCGGTGTATCCCAGTCGCTGGGCCTGGGCGTGGTGCAAATGAACACCGGGCTGAGCCAGGTGTTTGACCTGCCCATCAGCCTGGGGGTGCAGATCACCTTGATCGTGCTCATTa                                                                                               <  2:2232519/150‑1 (MQ=255)
                                                        aTGTGGTGGACATCCTTGCCGTGGTGATTACCGCGTTCGGTGTATCCCAGTCGCTGGGCCTGGGCGTGGTGCAAATGAACACCGGGCTGAGCCAGGTGTTTGACCTGCCCATCAGCCTGGGGGTGCAGATCACCTTGATCGTGCTCATTa                                                                                               <  1:2553380/150‑1 (MQ=255)
                                                         tgtgGTGGACATCCTTGCCGTGGTGATTACCGCGTTCGGTGTATCCCAGTCGCTGGGCCTGGGCGTGGTGCAAATGAACACCGGGCTGAGCCAGGTGTTTGACCTGCCCATCAGCCTGGGGGTGCAGATCACCTTGATCGTGCTCATTAc                                                                                              >  2:2561362/1‑150 (MQ=255)
                                                         tgtgGTGGACATCCTTGCCGTGGTGATTACCGCGTTCGGTGTATCCCAGTCGATGGGCTGGGTCGTGGTGCAAATGAACACCGGGCTGAGCCAGGTGTTTGACCTGCCCATCAGCCTGGGGGTGCAGATCACCTTGATCGTGCTCATTAc                                                                                              <  2:416327/150‑1 (MQ=255)
                                                           tggtggACATCCTTGCCGTGGTGATTACCGCGTTCGGTGTATCCCAGTCGCTGGGCCTGGGCGTGGTGCAAATGAACACCGGGCTGAGCCAGGTGTTTGACCTGCCCATCAGCCTGGGGGTGCAGATCACCTTGATCGTGCTCATTAccc                                                                                            >  2:2557415/1‑150 (MQ=255)
                                                              tggACATCCTTGCCGTGGTGATTACCGCGTTCGGTGTATCCCAGTCGCTGGGCCTGGGCGTGGTGCAAATGAACACCGGGCTGAGCCAGGTGTTTGACCTGCCCATCAGCCTGGGGGTGCAGATCACCTTGATCGTGCTCATTACCCTCg                                                                                         >  1:2805084/1‑150 (MQ=255)
                                                                 aCATCCTTGCCGTGGTGATTACCGCGTTCGGTGTATCCCAGTCGCTGGGCCTGGGCGTGGTGCAAATGAACACCGGGCTGAGCCAGGTGTTTGACCTGCCCATCAGCCTGGGGGTGCAGATCACCTTGATCGTGCTCATTACCCTCGTta                                                                                      <  1:2456594/150‑1 (MQ=255)
                                                                  cATCCTTGCCGTGGTGATTACCGCGTTCGGTGTATCCCAGTCGCTGGGCCTGGGCGTGGTGCAAATGAACACCGGGCTGAGCCAGGTGTTTGACCTGCCCATCAGCCTGGGGGTGCAGATCACCTTGATCGTGCTCATTACCCTCGTtac                                                                                     <  2:201044/150‑1 (MQ=255)
                                                                   aTCCTTGCCGTGGTGATTACCGCGTTCGGTGTATCCCAGTCGCTGGGCCTGGGCGTGGTGCAAATGAACACCGGGCTGAGCCAGGTGTTTGACCTGCCCATCAGCCTGGGGGTGCAGATCACCTTGATCGTGCTCATTACCCTCGTtact                                                                                    >  1:974277/1‑150 (MQ=255)
                                                                      cTTGCCGTGGTGATTACCGCGTTCGGTGTATCCCAGTCGCTGGGCCTGGGCGTGGTGCAAATGAACACCGGGCTGAGCCAGGTGTTTGACCTGCCCATCAGCCTGGGGGTGCAGATCACCTTGATCGTGCTCATTACCCTCGTTACTAcc                                                                                 >  2:1998912/1‑150 (MQ=255)
                                                                      cTTGCCGTGGTGATTACCGCGTTCGGTGTATCCCAGTCGCTGGGCCTGGGCGTGGTGCAAATGAACACCGGGCTGAGCCAGGTGTTTGACCTGCCCATCAGCCTGGGGGTGCAGATCACCTTGATCGTGCTCATTACCCTCGTTACTAcc                                                                                 >  1:2457885/1‑150 (MQ=255)
                                                                      cTTGCCGTGGTGATTACCGCGTTCGGTGTATCCCAGTCGCTGGGCCTGGGCGTGGTGCAAATGAACACCGGGCTGAGCCAGGTGTTTGACCTGCCCATCAGCCTGGGGGTGCAGATCACCTTGATCGTGCTCATTACCCTCGTTACTAcc                                                                                 <  1:2557415/150‑1 (MQ=255)
                                                                      cTTGCCGTGGTGATTACCGCGTTCGGTGTATCCCAGTAGCTGGGCCTGGGCGTGGTGCAAATGAACACCGGGCTGAGCCAGGTGTTTGACCTGCCCATCAGCCTGGGGGTGCAGATCACCTTGATCGTGCTCATTACCCTCGTTACTAcc                                                                                 >  2:1183546/1‑150 (MQ=255)
                                                                       ttGCCGTGGTGATTACCGCGTTCGGTGTATCCCAGTCGCTGGGCCTGGGCGTGGTGCAAATGAACACCGGGCTGAGCCAGGTGTTTGACCTGCCCATCAGCCTGGGGGTGCAGATCACCTTGATCGTGCTCATTACCCTCGTTACTACCg                                                                                >  1:2751140/1‑150 (MQ=255)
                                                                            gtggtgATTACCGCGTTCGGTGTATCCCAGTCGCTGGGCCTGGGCGTGGTGCAAATGAACACCGGGCTGAGCCAGGTGTTTGACCTGCCCATCAGCCTGGGGGTGCAGATCACCTTGATCGTGCTCATTACCCTCGTTACTACCGTATCg                                                                           <  2:1666520/150‑1 (MQ=255)
                                                                             tggtgATTACCGCGTTCGGTGTATCCCAGTCGCTGGGCCTGGGCGTGGTGCAAATGAACACCGGGCTGAGCCAGgt                                                                                                                                                    <  2:72913/76‑1 (MQ=255)
                                                                             tggtgATTACCGCGTTCGGTGTATCCCAGTCGCTGGGCCTGGGCGTGGTGCAAATGAACACCGGGCTGAGCCAGgt                                                                                                                                                    >  1:72913/1‑76 (MQ=255)
                                                                             tggtgATTACCGCGTTCGGTGTATCCCAGTCGCTGGGCCTGGGCGTGGTGCAAATGAACACCGGGCTGAGCCAGGTGTTTGACCTGCCCATCAGCCTGGGGGTGCAGATCACCTTGATCGTGCTCATTACCCTCGTTACTACCGTATCgg                                                                          >  1:1579377/1‑150 (MQ=255)
                                                                              ggtgATTACCGCGTTCGGTGTATCCCAGTCGCTGGGCCTGGGCGTGGTGCAAATGAACACCGGGCTGAGCCAGGTGTTTGACCTGCCCATCAGCCTGGGGGTGCAGATCACCTTGATCGTGCTCATTACCCTCGTTACTACCGTATCg                                                                           <  1:310609/148‑1 (MQ=255)
                                                                              ggtgATTACCGCGTTCGGTGTATCCCAGTCGCTGGGCCTGGGCGTGGTGCAAATGAACACCGGGCTGAGCCAGGTGTTTGACCTGCCCATCAGCCTGGGGGTGCAGATCACCTTGATCGTGCTCATTACCCTCGTTACTACCGTATCg                                                                           >  2:310609/1‑148 (MQ=255)
                                                                                tgATTACCGCGTTCGGTGTATCCCAGTCGCTGGGCCTGGGCGTGGTGCAAATGAACACCGGGCTGAGCCAGGTGTTTGACCTGCCCATCAGCCTGGGGGTGCAGATCATCTTGATCGTGCTCATTACCCTCGTTACTACCGTATCGGTGa                                                                       <  1:2669959/150‑1 (MQ=255)
                                                                                tgATTACCGCGTTCGGTGTATCCCAGTCGCTGGGCCTGGGCGTGGTGCAAATGAACACCGGGCTGAGCCAGGTGTTTGACCTGCCCATCAGCCTGGGGGTGCAGATCACCTTGATCGTGCTCATTACCCTCGTTACTACCGTATCGGTGa                                                                       <  1:3029541/150‑1 (MQ=255)
                                                                                 gATTACCGCGTTCGGTGTATCCCAGTCGCTGGGCCTGGGCGTGGTGCAAATGAACACCGGGCTGAGCCAGGTGTTTGACCTGCCCATCAGCCTGGGGGTGCAGATCACCTTGATCGTGCTCATTACCCTCGTTACTACCGTATCGGTGAt                                                                      <  2:875884/150‑1 (MQ=255)
                                                                                 gATTACCGCGTTCGGTGTATCACAGTCGCTGGGCCTGGGCGTGGTGCAAATGAACACCGGGCTGAGCCAGGTGTTTGACCTGCCCATCAGCCTGGGGGTGCAGATCACCTTGATCGTGCTCATTACCCTCGTTACTACCGTATCGGTGAt                                                                      >  1:1283773/1‑150 (MQ=255)
                                                                                   ttACCGCGTTCGGTGTATCCCAGTCGCTGGGCCTGGGCGTGGTGCAAATGAACACCGGGCTGAGCCAGGTGTTTGACCTGCCCATCAGCCTGGGGGTGCAGATCACCTTGATCGTGCTCATTAccc                                                                                            >  2:1174115/1‑126 (MQ=255)
                                                                                   ttACCGCGTTCGGTGTATCCCAGTCGCTGGGCCTGGGCGTGGTGCAAATGAACACCGGGCTGAGCCAGGTGTTTGACCTGCCCATCAGCCTGGGGGTGCAGATCACCTTGATCGTGCTCATTAccc                                                                                            <  1:1174115/126‑1 (MQ=255)
                                                                                      ccGCGTTCGGTGTATCCCAGTCGCTGGGCCTGGGCGTGGTGCAAATGAACACCGGGCTGAGCCAGGTGTTTGACCTGCCCATCAGCCTGGGGGTGCAGATCACCTTGATCGTGCTCATTACCCTCGTTACTACCGTATCGGTGATGGCCg                                                                 >  2:2245736/1‑150 (MQ=255)
                                                                                       cgcgTTCGGTGTATCCCAGTCGCTGGGCCTGGGCGTGGTGCAAATGAACACCGGGCTGAGCCAGGTGTTTGACGTGCCCATCAGCCTGGGGGTGCAGATCACTTTGATCGTGCTCATTACCCTCGTTACTACCGTATCGGTGATGGCCgg                                                                >  1:1456128/1‑150 (MQ=255)
                                                                                       cgcgTTCGGTGTATCCCAGTCGCTGGGCCTGGGCGTGGTGCAAATGAACACCGGGCTGAGCCAGGTGTTTGACCTGCCCATCAGCCTGGGGGTGCAGATCACCTTGATCGTGCTCATTACCCTCGTTACTACCGTAt                                                                             <  1:733671/137‑1 (MQ=255)
                                                                                       cgcgTTCGGTGTATCCCAGTCGCTGGGCCTGGGCGTGGTGCAAATGAACACCGGGCTGAGCCAGGTGTTTGACCTGCCCATCAGCCTGGGGGTGCAGATCACCTTGATCGTGCTCATTACCCTCGTTACTACCGTATCGGTGATGGCCgg                                                                >  2:1243500/1‑150 (MQ=255)
                                                                                           ttCGGTGTATCCCAGTCGCTGGGCCTGGGCGTGGTGCAAATGAACACCGGGCTGAGCCAGGTGTTTGACCTGCCCATCAGCCTggg                                                                                                                            >  1:1991582/1‑86 (MQ=255)
                                                                                             cGGTGTATCCCAGTCGCTGGGCCTGGGCGTGGTGCAAATGAACACCGGGCTGAGCCAGGTGTTTGACCTGCCCATCAGCCTGGGGGTGCAGATCACCTTGATCGTGCTCATTACCCTCGTTACTACCGTATCGGTGATGGCCGGCGTGTc                                                          <  2:1444977/150‑1 (MQ=255)
                                                                                               gtgtATCCCAGTCGCTGGGCCTGGGCGTGGTGCAAATGAACACCGGGCTGAGCCAGGTGTTTGACCTGCCCATCAGCCTGGGGGTGCAGATCACCTTGATCGTGCTCATTACCCTCGTTACTACCGTATCGGTGATGGCCGGCGTGTcgc                                                        >  2:951284/1‑150 (MQ=255)
                                                                                                 gtATCCCAGTCGCTGGGCCTGGGCGTGGTGCAAATGAACACCGGGCTGAGCCAGGTGTTTGACCTGCCCATCAGCCTGGGGGTGCAGATCACCTTGATc                                                                                                         <  2:2032479/99‑1 (MQ=255)
                                                                                                 gtATCCCAGTCGCTGGGCCTGGGCGTGGTGCAAATGAACACCGGGCTGAGCCAGGTGTTTGACCTGCCCATCAGCCTGGGGGTGCAGATCACCTTGATc                                                                                                         >  1:2032479/1‑99 (MQ=255)
                                                                                                     cccAGTCGCTGGGCCTGGGCGTGGTGCAAATGAACACCGGGCTGAGCCAGGTGTTTGACCTGCCCATCAGCCTGGGGGTGCAGATCACCTTGATCGTGCTCATTACCCTCGTTACTACCGTATCGGTGATGGCCGGCGTGTCGCGTGGCa                                                  >  2:422145/1‑150 (MQ=255)
                                                                                                      ccAGTCGCTGGGCCTGGGCGTGGTGCAAATGAACACCGGGCTGAGCCAGGTGTTTGACCTGCCCATCAGCCTGGGGGTGCAGATCACCTTGCTCGTGCTCATTACCCTCGTTACTACCGTATCGGTGATGACCGGCGTGTAGCGTGGCAt                                                 >  1:1075829/1‑150 (MQ=255)
                                                                                                      ccAGTCGCTGGGCCTGGGCGTGGTGCAAATGAACACCGGGCTGAGCCAGGTGTTTGACCTGCCCATCAGCCTGGGGGTGCAGATCACCTTGATCGTGCTCATTACCCTCGTTACTACCGTATCGGTGATGGCCGGCGTGTCGCGTGGCAt                                                 >  2:2355066/1‑150 (MQ=255)
                                                                                                       cAGTCGCTGGGCCTGGGCGTGGTGCAAATGAACACCGGGCTGAGCCAGGTGTTTGACCTGCCCAACAGCCTGGAGGTGCAGATCACCTTGATCGTGCTCACTACCCTCCTGACTACCTGATCGGTGAAGGCCGGCGTGTCGAGCGGCATg                                                >  1:2981584/1‑150 (MQ=255)
                                                                                                                  cctgggcGTGGTGCAAATGAACACCGGGCTGAGCCAGGTGTTTGACCTGCCCATCAGCCTGGGGGTGCAGATCACCTTGATCGTGCTCATTACCCTCGTTACTACCGTATCGGTGATGGCCGGCGTGTCGCGTGGCATGAAGCGCCTTTc                                     <  2:368755/150‑1 (MQ=255)
                                                                                                                  cctgggcGTGGTGCAAATGAACACCGGGCTGAGCCAGGTGTTTGACCTGCCCATCAGCCTGGGGGTGCAGATCACCTTGATCGTGCTCATTACCCTCGTTACTACCGTATCGGTGATGGCCGGCGTGTCGCGTGGCATGAAGCGCCTTTc                                     >  2:38313/1‑150 (MQ=255)
                                                                                                                  cctgggcGTGGTGCAAATGAACACCGGGCTGAGCCAGGTGTTTGACCTGCCCATCAGCCTGGGGGTGCAGATCACCTTGATCGTGCTCATTACCCTCGTTACTACCGTATCGGTGATGGCCGGCGTGTCGCGTGGCATGAAGCGCCTTTc                                     <  1:2825582/150‑1 (MQ=255)
                                                                                                                  cctgggcGTGGTGCAAATGAACACCGGGCTGAGCCAGGTGTTTGACCTGCCCATCAGCCTGGGGGTGCAGATCACCTTGATCGAACTCCTTACACACTTTATTACCTTGTCTGTGATGGCCGGCGTGTCGCGTGGCATGAAGCGCCTTTc                                     >  1:2218150/1‑150 (MQ=255)
                                                                                                                   ctgggcGTGGTGCAAATGAACACCGGGCTGAGCCAGGTGTTTGACCTGCCCATCAGCCTGGGGGTGCAGATCACCTTGATCGTGCTCATTACCCTCGTTACTACCGTATCGGTGATGGCCGGCGTGTCGCGTGGCATGAAGCGCCTTTCg                                    <  1:2245736/150‑1 (MQ=255)
                                                                                                                   ctgggcGTGGTGCAAATGAACACCGGGCTGAGCCAGGTGTTTGACCTGCCCATCAGCCTGGGGGTGCAGATCACCTTGATCGTGCTCATTACCCTCGTTACTACCGTATCGGTGATGGCCGGCGTGTCGCGTGGCATGAAGCGCCTTTCg                                    <  1:2657560/150‑1 (MQ=255)
                                                                                                                    tgggcGTGGTGCAAATGAACACCGGGCTGAGCCAGGTGTTTGACCTGCCCATCAGCCTGGGGGTGCAGATCACCTTGATCGTGCTCATTACCCTCGTTACTACCGTATCGGTGATGGCCGGCGTGTCGCGTGGCATGAAGCGCCTTTCg                                    >  2:1673020/1‑149 (MQ=255)
                                                                                                                    tgggcGTGGTGCAAATGAACACCGGGCTGAGCCAGGTGTTTGACCTGCCCATCAGCCTGGGGGTGCAGATCACCTTGATCGTGCTCATTACCCTCGTTACTACCGTATCGGTGATGGCCGGCGTGTCGCGTGGCATGAAGCGCCTTTCg                                    <  1:1673020/149‑1 (MQ=255)
                                                                                                                        cGTGGTGCAAATGAACACCGGGCTGAGCCAGGTGTTTGACCTGCCCATCAGCCTGGGGGTGCAGATCACCTTGATCGTGCTCATTACCCTCGTTACTACCGTATCGGTGATGGCCGGCGTGTCGCGTGGCATGAAGCGCCTTTCggaatg                               <  1:1913040/150‑1 (MQ=255)
                                                                                                                           ggtgCAAATGAACACCGGGCCGAGCCAGGTGTTTGACCTGCCCATCAGCCTGGGGGTGCAGATCACCTTGATCGTGCTCATTAcc                                                                                             >  1:2776044/1‑85 (MQ=255)
                                                                                                                               cAAATGAACACCGGGCTGAGCCAGGTGTTTGACCTGCCCATCAGCCTGGGGGTGCAGATCACCTTGATCGTGCTCATTACCCTCGTTACTACCGTATCGGTGATGGCCGGCgt                                                             >  1:145620/1‑113 (MQ=255)
                                                                                                                               cAAATGAACACCGGGCTGAGCCAGGTGTTTGACCTGCCCATCAGCCTGGGGGTGCAGATCACCTTGATCGTGCTCATTACCCTCGTTACTACCGTATCGGTGATGGCCGGCGTGTCGCGTGGCATGAAGCGCCTTTCGGAATGGAATATg                        >  1:584175/1‑150 (MQ=255)
                                                                                                                                  aTGAACACCGGGCTGAGCCAGGTGTTTGACCTGCCCATCAGCCTGGGGGTGCAGATCACCTTGATCGTGCTCATTACCCTCGTTACTACCGTATCGGTGATGGCCGGCGTGTCGCGTGGCATGAAgcg                                           <  2:1411176/128‑1 (MQ=255)
                                                                                                                                  aTGAACACCGGGCTGAGCCAGGTGTTTGACCTGCCCATCAGCCTGGGGGTGCAGATCACCTTGATCGTGCTCATTACCCTCGTTACTACCGTATCGGTGATGGCCGGCGTGTCGCGTGGCATGAAgcg                                           >  1:1411176/1‑128 (MQ=255)
                                                                                                                                  aTGAACACAGGGCTGAGCCAGGTGTTTGACCTGCACATCAGACTGGGGGTGAAGATCACCTTGATCGTGATCATTACCCTCGTTACTACCGTATCGGTGATGGCCGGCGTGTCGCGTGGCATGAAATGTCGTTGGGAATGGAATATgatt                     >  2:1475711/1‑147 (MQ=255)
                                                                                                                                    gAACACCGGGCTGAGCCAGGTGTTTGACCTGCCCATCAGCCTGGGGGTGCAGATCACCTTGATCGTGCTCATTACCCTCGTTACTACCGTATCGGTGATGGCCGGCGTGTCGCGTGGCATGAAGCGCCTTTCGGAATGGAATATGCtgtt                   <  2:3069749/150‑1 (MQ=255)
                                                                                                                                    aaaCACCTGGCTCAGCCCGGTGTTTGACCTGCCCATCAGCCTGGGGGTGCAGATCACCTTGATCGTGGTCATTACCCTCGTTACTACCGTATCGGGGATGGCCGGCGGGTGGCGTGGCGTGAATCGCCTTGCGGGGTGGGATATGGtgtt                   >  2:2942027/2‑150 (MQ=255)
                                                                                                                                       cacCGGGCTGAGCCAGGTGTTTGACCTGCCCATCAGCCTGGGGGTGCAGATCACCTTGATCGTGCTCATTACCCTCGTTACTACCGTATCGGTGATGGCCGGCGTGTCGCGTGGCATGAAGCGCCTTTCGGAATGGAATATGCTGTTGTc                >  1:432032/1‑150 (MQ=255)
                                                                                                                                        acCGGGCTGAGCCAGGTGTTTGACCTGCCCATCAGCCTGGGGGTGCAGATCACCTTGATCGt                                                                                                       >  1:584381/1‑62 (MQ=255)
                                                                                                                                        acCGGGCTGAGCCAGGTGTTTGACCTGCCCATCAGCCTGGGGGTGCAGATCACCTTGATCGt                                                                                                       <  2:584381/62‑1 (MQ=255)
                                                                                                                                        acCGGGCTGAGCCAGGTGTTTGACCTGCCCATCAGCCTGGGGGTGCAGATCACCTTGATCGTGCTCATTACCCTCGTTACTACCGTATCGGTGATGGCCGGCGTGTc                                                          >  1:1093220/1‑107 (MQ=255)
                                                                                                                                        acCGGGCTGAGCCAGGTGTTTGACCTGCCCATCAGCCTGGGGGTGCAGATCACCTTGATCGTGCTCATTACCCTCGTTACTACCGTATCGGTGATGGCCGGCGTGTCGCGTGGCATGAAGCGcc                                         <  2:2818631/124‑1 (MQ=255)
                                                                                                                                        acCGGGCTGAGCCAGGTGTTTGACCTGCCCATCAGCCTGGGGGTGCAGATCACCTTGATCGTGCTCATTACCCTCGTTACTACCGTATCGGTGATGGCCGGCGTGTCGCGTGGCATGAAGCGcc                                         >  1:2818631/1‑124 (MQ=255)
                                                                                                                                        acCGGGCTGAGCCAGGTGTTTGACCTGCCCATCAGCCTGGGGGTGCAGATCACCTTGATCGTGCTCATTACCCTCGTTACTACCGTATCGGTGATGGCCGGCGTGTCGCGTGGCATGAAGCGCCTTTCGGAATGGAATATGCTGTTGTCg               >  2:1303621/1‑150 (MQ=255)
                                                                                                                                          cGGGCTGAGCCAGGTGTTTGACCTGCCCATCAGCCTGGGGGTGCAGATCACCTTGATCGTGCTCATTACCCTCGTTACTACCGTATCGGTGATGGCCGGCGTGTCGCGTGGCATGAAGCGCCTTTCGGAATGGAATATGCTGTTGTCggt             <  2:584175/150‑1 (MQ=255)
                                                                                                                                           gggCTGAGCCAGGTGTTTGACCTGCCCATCAGCCTGGGGGTGCAGATCACCTTGATCGTGCTCATTACCCTCGTTACTACCGTATCGGTGATGGCCGGCGTGTCGCGTGGCATGAAGCGCCTTTCGGAATGGAATATGCTGTTGTCggtg            <  1:993700/150‑1 (MQ=255)
                                                                                                                                           gggCTGAGCCAGGTGTTTGACCTGCCCATCAGCCTGGGGGTGCAGATCACCTTGAGCGTGCTCATTACCCTCGTTACTACCGTATCGGTGATGGTCGGCGTGTCGCGTGGCATGAAGCGCCTTTCGGCATGGAATATGCTGTTGTCggcg            <  2:1283773/150‑3 (MQ=255)
                                                                                                                                              cTGAGCCAGGTGTTTGACCTGCACGTCAGGCCGGGGGTGCAGATCACCTTGATCGTGCTCATTACCCTCGTTACTACCGTATCGGTGATGGCCGGCGTGTCGCGTGGCATGAAGCGCCTTTCGGAATGAAAAATAATGTTGTCggggggg         <  2:2805084/150‑3 (MQ=255)
                                                                                                                                               tGAGCCAGGTGTTTGCCCTGCCCATCAGCGTGGGGGTGCAGATCACCTTGATCGTGCTCATTACCCTCGTTACTACCGTATCGGTGATGGCCGGCGTGTCGCGTGGCATGAAGCGCCTTTCGGAATGGAATATGCTCTTGTCGGTGGTGc        <  2:677612/150‑1 (MQ=255)
                                                                                                                                                 aGCCCGGTGTTTGACCTGCCCATCAGCCTGGGGGTGCAGATCACCTTGATCGTGCTCATTACCCTCGTTACTACCGTATCGGTGATGGCCGGCGTGTCGCGTGGCATGAAGCGCCTTTCGGAATGGAatat                         >  1:135932/1‑131 (MQ=255)
                                                                                                                                                   ccAGGTGTTTGACCTGCCCATCAGCCTGGGGGTGCAGATCACCTTGATCGTGCTCATTACCCTCGTTACTACCGTATCGGTGATGGCCGGCGTGTCGCGTGGCATGAAGCGCCTTTCGGAATGGAATATGCTGTTGTCGGTGGTGCtggt    >  2:3005695/1‑150 (MQ=255)
                                                                                                                                                   ccAGGTGTTTGACCTGCCCATCAGCCTGGGGGTGCAGATCACCTTGATCGTGCTCATTACCCTCGTTACTACCGTATCGGTGATGGCCGGCGTGTCGCGTGGCATCAAGCGCCTTTCGGAATGGAATATGCTGTTGTCGATGGTGCcggt    >  1:2186176/1‑150 (MQ=255)
                                                                                                                                                    cAGGTGTTTGACCTGCCCATCAGCCTGGGGGTGCAGATCACCTTGATCGTGCTCATTACCCTCGTTACTACCGTATCg                                                                           >  2:1482895/1‑78 (MQ=255)
                                                                                                                                                     aGGTGTTTGACCTGCCCATCAGCCTGGGGGTGCAGATCACCTTGATCGTGCTCATTACCCTCGTTACTACCGTATCGGTGATGGCCGGCGTGTCGCGTGGCATGAAGCGCCTTTCGGAATGGAATATGCTGTTGTCGGTGGTGCtggtgg  >  1:966833/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CATGCGCTCGATCCTTTACCCGCTGATTGGTGAGCGCATCTATGGGCCCATCGGACATGTGGTGGACATCCTTGCCGTGGTGATTACCGCGTTCGGTGTATCCCAGTCGCTGGGCCTGGGCGTGGTGCAAATGAACACCGGGCTGAGCCAGGTGTTTGACCTGCCCATCAGCCTGGGGGTGCAGATCACCTTGATCGTGCTCATTACCCTCGTTACTACCGTATCGGTGATGGCCGGCGTGTCGCGTGGCATGAAGCGCCTTTCGGAATGGAATATGCTGTTGTCGGTGGTGCTGGTGG  >  NC_002947/6125774‑6126072

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: