Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,127,845 A→C 6.4% 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,127,8450AC 6.4% 557.1 / 21.0 173Q221P (CAG→CCG) PP_5374choline/carnitine/betaine transporter family protein
Reads supporting (aligned to +/- strand):  ref base A (82/80);  new base C (5/6);  total (87/86)
Fisher's exact test for biased strand distribution p-value = 7.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.84e-01

Only 100 of 174 total aligned reads displayed.

ATGCGCTCGATCCTTTACCCGCTGATTGGTGAGCGCATCTATGGGCCCATCGGACATGTGGTGGACATCCTTGCCGTGGTGATTACCGCGTTCGGTGTATCCCAGTCGCTGGGCCTGGGCGTGGTGCAAATGAACACCGGGCTGAGCCAGGTGTTTGACCTGCCCATCAGCCTGGGGGTGCAGATCACCTTGATCGTGCTCATTACCCTCGTTACTACCGTATCGGTGATGGCCGGCGTGTCGCGTGGCATGAAGCGCCTTTCGGAATGGAATATGCTGTTGTCGGTGGTGCTGGTGG  >  NC_002947/6127697‑6127994
                                                                                                                                                    |                                                                                                                                                     
aTGCGCTCGATCCTTTACCCGCTGATTGGTGAGCGCATCTATGGGCCCATCGGACATGTGGTGGACATCCTTGCCGTGGTGATTACCGCGTTCGGTGTATCCCAGTCGCTGGGCCTGGGCGTGGTGCAAATGAACACCGGGCTGAGCCAg                                                                                                                                                      <  2:3389264/150‑1 (MQ=255)
    gcTCGATCCTTTACCCGCTGATTGGTGAGCGCATCTATGGGCCCATCGGACATGTGGTGGACATCCTTGCCGTGGTGATTACCGCGTTCGGTGTATCCCAGTCGCTGGGCCTGGGCGTGGTGCAAATGAACACCTGGCTCAGCCCGgtgt                                                                                                                                                  <  1:2773391/150‑1 (MQ=255)
      tCGATCCTTTACCCGCTGATTGGTGAGCGCATCTATGGGCCCATCGGACATGTGGTGGACATCCTTGCCGTGGTGATTACCGCGTTCGGTGTATCCCAGTCGCTGGGCCTGGGCGTGGTGCAAATGAACACCGGGCTGAGCCAGGTGttt                                                                                                                                                <  2:734929/150‑1 (MQ=255)
      tCGATCCTTTACCCGCTGATTGGTGAGCGCATCTATGGGCCCATCGGACATGTGGTGGACATCCTTGCCGTGGTGATTACCGCGTTCGGTGTATCCCAGTCGCTGGGCCTGGGCGTGGTGCAAATGAACACCGGGCTGAGCCAGGTGttt                                                                                                                                                <  2:592268/150‑1 (MQ=255)
      tCGATCCTTTACCCGCTGATTGGTGAGCGCATCTATGGGCCCATCGGACATGTGGTGGACATCCTTGCCGTGGTGATTACCGCGTTCGGTGTATCCCAGTCGCTGGGCCTGGGCGTGGTGCAAATGAACACCGGGCTGAGCCAGGTGttt                                                                                                                                                >  2:4531483/1‑150 (MQ=255)
      tCGATCCTTTACCCGCTGATTGGTGAGCGCATCTATGGGCCCATCGGACATGTGGTGGACATCCTTGCCGTGGTGATTACCGCGTTCGGTGTATCCCAGTCGCTGGGCCTGGGCGTGGTGCAAATGAACACCGGGCTGAGCCAGGTGttt                                                                                                                                                <  2:1318809/150‑1 (MQ=255)
        gATCCTTTACCCGCTGATTGGTGAGCGCATCTATGGGCCCATCGGACATGTGGTGGACATCCTTGCCGTGGTGATTACCGCGTTCGGTGTATCCCAGTCGCTGGGCCTGGGCGTGGTGCAAATGAACACCGGGCTGAGCCAGGTGTTTGa                                                                                                                                              <  2:4158842/150‑1 (MQ=255)
        gATCCTTTACCCGCTGATTGGTGAGCGCATCTATGGGCCCATCGGACATGTGGTGGACATCCTTGCCGTGGTGATTACCGCGTTCGGTGTATCCCAGTCGCTGGGCCTGGGCGTGGTGCAAATGAACACCGGGCTGAGCCAGGTGTTTGa                                                                                                                                              <  2:1354167/150‑1 (MQ=255)
        gATCCTTTACCCGCTGATTGGTGAGCGCATCTATGGGCCCATCGGACATGTGGTGGACATCCTTGCCGTGGTGATTACCGCGTTCGGTGTATCCCAGTCGCTGGGCCTGGGCGTGGTGCAAATGAACACCGGGCTGAGCCAGGTGTTTGa                                                                                                                                              <  1:846509/150‑1 (MQ=255)
        gATCCTTTACCCGCTGATTGGTGAGCGCATCTATGGGCCCATCGGACATGTGGTGGACATCCTTGCCGTGGTGATTACCGCGTTCGGTGTATCCCAGTCGCTGGGCCTGGGCGTGGTGCAAATGAACACCGGGCTGAGCCAGGTGTTTGa                                                                                                                                              <  1:4687455/150‑1 (MQ=255)
        gATCCTTTACCCGCTGATTGGTGAGCGCATCTATGGGCCCATCGGACATGTGGTGGACATCCTTGCCGTGGTGATTACCGCGTTCGGTGTATCCCAGTCGCTGGGCCTGGGCGTGGTGCAAATGAACACCGGGCTGAGCCAGGTGTTTGa                                                                                                                                              >  1:4266908/1‑150 (MQ=255)
        gATCCTTTACCCGCTGATTGGTGAGCGCATCTATGGGCCCATCGGACATGTGGTGGACATCCTTCCCGTGGTGATTACCGCGTTCGGTGTATCCCAGTCGCTGGGCCTGGGCGTGGTGCAAATGAACACCGGGCTGAGCCAGGTGTTTGa                                                                                                                                              <  2:941424/150‑1 (MQ=255)
          tCCTTTACCCGCTGATTGGTGAGCGCATCTATGGGCCCATCGGACATGTGGTGGACATCCTTGCCGTGGTGATTACCGCGTTCGGTGTATCCCAGTCGCTGGGCCTGGGCGTGGTGCAAATGAACACCTGGCTCAGCCCGGTGttcattt                                                                                                                                            <  2:2021228/150‑6 (MQ=255)
             tttACCCGCTGATTGGTGAGCGCATCTATGGGCCCATCGGACATGTGGTGGACATCCTTGCCGTGGTGATTACCGCGTTCGGTGTATCCCAGTCGCTGGGCCTGGGCGTGGTGCAAATGAACACCGGGCTGAGCCAGGTGTTTGACCTGc                                                                                                                                         >  1:1409452/1‑150 (MQ=255)
               tACCCGCTGATTGGTGAGCGCATCTATGGGCCCATCGGACATGTGGTGGACATCATTGCCGTGGTGATTACCGCGTTCGGTGTATCCCAGTCGCTGGGCCTGGGCGTGGTGCAAATGAACACCGGGCTGAGCCAGGTGTTTGACCTGccc                                                                                                                                       <  2:1881591/150‑1 (MQ=255)
                aCCCGCTGATTGGTGAGCGCATCTATGGGCCCATCGGACATGTGGTGGACATCCTTGCCGTGGTGATTACCGCGTTCGGTGTATCCCAGTCGCTGGGCCTGGGCGTGGTGCAAATGAACACCGGGCTGAGCCAGGTGTTTGACCTGCCCa                                                                                                                                      <  1:5012995/150‑1 (MQ=255)
                 cccGCTGATTGGTGAGCGCATCTATGGGCCCATCGGACATGTGGTGGACATCCTTGCCGTGGTGATTACCGCGTTCGGTGTATCCCAGTCGCTGGGCCTGGGCGTGGTGCAAATGAACACCGGGCTGAGCCAGGTGTTTGACCTGCCCAt                                                                                                                                     >  1:2112737/1‑150 (MQ=255)
                   cGCTGATTGGTGAGCGCATCTATGGGCCCATCGGACATGTGGTGGACATCCTTGCCGTGGTGATTACCGCGTTCGGTGTATCCCAGTCGCTGGGCCTGGGCGTGGTGCAAATGAACACCGGGCTGAGCCAGGTGTTTGACCTGCCCATCa                                                                                                                                   >  1:2010009/1‑150 (MQ=255)
                    gCTGATTGGTGAGCGCATCTATGGGCCCATCGGACATGTGGTGGACATCCTTGCCGTGGTGATTACCGCGTTCGGTGTATCCCAGTCGCTGGGCCTGGGCGTGGTGCAAATGAACACCGGGCTGAGCCAGGTGTTTGACCTGCCCATCAg                                                                                                                                  >  1:4151902/1‑150 (MQ=255)
                        aTTGGTGAGCGCATCTATGGGCCCATCGGACATGTGGTGGACATCCTTGCCGTGGTGATTACCGCGTTCGGTGTATCCCAGTCGCTGGGCCTGGGCGTGGTGCAAATGAACACCGGGCTGAGCCAGGTGTTTGACCTGCCCATCAGCCTg                                                                                                                              <  2:408566/150‑1 (MQ=255)
                         ttGGTGAGCGCATCTATGGGCCCATCGGACATGTGGTGGACATCCTTGCCGTGGTGATTACCGCGTTCGGTGTATCCCAGTCGCTGGGCCTGGGCGTGGTGCAAATGAACACCGGGCTGAGCCAGGTGTTTGACCTGCCCATCAGCCTgg                                                                                                                             <  2:2951920/150‑1 (MQ=255)
                          tGGTGAGCGCATCTATGGGCCCATCGGACATGTGGTGGACATCCTTGCCGTGGTGATTACCGCGTTCGGTGTATCCCAGTCGCTGGGCCTGGGCGTGGTGCAAATGAACACCTGGCTCAGCCCGGTGttcat                                                                                                                                              <  2:899277/132‑4 (MQ=255)
                          tGGTGAGCGCATCTATGGGCCCATCGGACATGTGGTGGACATCCTTGCCGTGGTGATTACCGCGTTCGGTGTATCCCAGTCGCTGGGCCTGGGCGTGGTGCAAATGAACACCGGGCTGAGCCAGGTGTTTGACCTGCCCATCAGCCTggg                                                                                                                            <  1:4594530/150‑1 (MQ=255)
                          tGGTGAGCGCATCTATGGGCCCATCGGACATGTGGTGGACATCCTTGCCGTGGTGATTACCGCGTTCGGTGTATCCCAGTCGCTGGGCCTGGGCGTGGTGCAAATGAACACCGGGCTGAGCCAGGTGTTTGACCTGCCCATCAGCCTggg                                                                                                                            <  1:5297132/150‑1 (MQ=255)
                           ggTGAGCGCATCTATGGGCCCATCGGACATGTGGTGGACATCCTTGCCGTGGTGATTACCGCGTTCGGTGTATCCCAGTCGCTGGGCCTGGGCGTGGTGCAAATGAACACCGGGCTGAGCCAGGTGTTTGACCTGCCCATCAGCCTgggg                                                                                                                           <  1:1134011/150‑1 (MQ=255)
                             tGAGCGCATCTATGGGCCCATCGGACATGTGGTGGACATCCTTGCCGTGGTGATTACCGCGTTCGGTGTATCCCAGTCGCTGGGCCTGGGCGTGGTGCAAATGAACACCGGGCTGAGCCAGGTGTTTGAc                                                                                                                                             <  1:4988549/130‑1 (MQ=255)
                               aGCGCATCTATGGGCCCATCGGACATGTGGTGGACATCCTTGCCGTGGTGATTACCGCGATCGGTGTATCCCAGTCGCTGGGCCTGGGCGTGGTGCAAATGAACACCGGGCTGAGCCAGGTGTTTGACCTGCCCATCAGCCTGGGGGTGc                                                                                                                       >  2:1529423/1‑150 (MQ=255)
                               aGCGCATCTATGGGCCAATCGGACATGTGGTGGACATCATTGCCGTGGTGATTACCGCGTTCGGTGTATCCCAGTCGCTGGGCCTGGGCGTGGTGCAAATGAACAACGGGCTGAGCCAGGTGTTTGACCTGCCCATCAGCCTGGGGGTGc                                                                                                                       >  2:1529303/1‑150 (MQ=255)
                                  gcATCTATGGGCCCATCGGACATGTGGTGGACATCCTTGCCGTGGTGATTACCGCGTTCGGTGTATCCCAGTCGCTGGGCCTGGGCGTGGTGCAAATGAACACCGGGCTGAGCCAGGTGTTTGACCTGCCCATCAGCCTGGGGGTGCAGa                                                                                                                    <  1:3531377/150‑1 (MQ=255)
                                    aTCTATGGGCCCATCGGACATGTGGTGGACATCCTTGCCGTGGTGATTACCGCGTTCGGTGTATCCCAGTCGCTGGGCCTGGGCGTGGTGCAAATGAACACCGGGCTGAGCCAGGTGTTTGACCTGCCCATCAGCCTGGGGGTGCAGATc                                                                                                                  >  1:319890/1‑150 (MQ=255)
                                     tCTATGGGCCCATCGGACATGTGGTGGACATCCTTGCCGTGGTGATTACCGCGTTCGGTGTATCCCAGTCGCTGGGCCTGGGCGTGGTGCAAATGACCACCGGGCTGAGCCAGGTGTTTGACCGGCCCATCAGCCTGGGGGTGCAGATCa                                                                                                                 <  1:1681049/150‑1 (MQ=255)
                                       tATGGGCCCATCGGACATGTGGTGGACATCCTTGCCGTGGTGATTACCGCGTTCGGTGTATCCCAGTCGCTGGGCCTGGGCGTGGTGCAAATGAACACCGGGCTGAGCCAGGTGTTTGACCTGCCCATCAGCCTGGGGGTGCAGATCAcc                                                                                                               >  2:2495879/1‑150 (MQ=255)
                                       tATGGGCCCATCGGACATGTGGTGGACATCCTTGCCGTGGTGATTACCGCGTTCGGTGTATCCCAGTCGCTGGGCCTGGGCGTGGTGCAAATGAACACCGGGCTGAGCCAGGTGTTTGACCTGCCCATCAGCCTGGGGGTGCAGATCAcc                                                                                                               >  2:2495885/1‑150 (MQ=255)
                                        aTGGGCCCATCGGACATGTGGTGGACATCCTTGCCGTGGTGATTACCGCGTTCGGTGTATCCCAGTCGCTGGGCCTGGGCGTGGTGCAAATGAACACCGGGCTGAGCCAGGTGTTTGACCTGCCCATCAGCCTGGGGGTGCAGATCACCt                                                                                                              <  1:2396232/150‑1 (MQ=255)
                                        aTGGGCCCATCGGACATGTGGTGGACATCCTTGCCGTGGTGATTACCGCGTTCGGTGTATCCCAGTCGCTGGGCCTGGGCGTGGTGCAAATGAACACCGGGCTGAGCCAGGTGTTTGACCTGCCCATCAGCCTGGGGGTGCAGATCACCt                                                                                                              <  2:2010009/150‑1 (MQ=255)
                                        aTGGGCCCATCGGACATGTGGTGGACATCCTTGCCGTGGTGATTACCGCGTTCGGTGTATCCCAGTCGCTGGGCCTGGGCGTGGTGCAAATGAACACCGGGCTGAGCCAGGTGTTTGACCTGCCCATCAGCCTGGGGGTGCAGATCACCt                                                                                                              >  1:4626642/1‑150 (MQ=255)
                                        aTGGGCCCATCGGACATGTGGTGGACATCCTTGCCGTGGTGATTACCGCGTTCGGTGTATCCCAGTCGCTGGGCCTGGGCGTGGTGCAAATGAACACCGGGCTGAGCCAGGTGTTTGACCTGCCCATCAGCCTGGGGGTGCAGATCACCt                                                                                                              >  1:3279096/1‑150 (MQ=255)
                                            gCCCATCGGACATGTGGTGGACATCCTTGCCGTGGTGATTACCGCGTTCGGTGTATCCCAGTCGCTGGGCCTGGGCGTGGTGCAAATGAACACCGGGCTGAGCCAGGTGTTTGACCTGCCCATCAGCCTGGGGGTGCAGATCACCTTGAt                                                                                                          >  2:3281426/1‑150 (MQ=255)
                                                aTCGGACATGTGGTGGACATCCTTGCCGTGGTGATTACCGCGTTCGGTGTATCCCAGTCGCTGGGCCTGGGCGTGGTGCAAATGAACACCGGGCTGAGCCAGGTGTTTGACCTGCCCATCAGCCTGGGGGTGCAGATCACCTTGATCGTg                                                                                                      >  1:381296/1‑150 (MQ=255)
                                                 tCGGACATGTGGTGGACATCCTTGCCGTGGTGATTACCGCGTTCGGTGTATCCCAGTCGCTGGGCCTGGGCGTGGTGCAAATGAACACCGGGCTGAGCCAGGTGTTTGACCTGCCCATCAGCCTGGGGGTGCAGATCACCTTGATCGTGc                                                                                                     <  1:2289295/150‑1 (MQ=255)
                                                 tCGGACATGTGGTGGACATCCTTGCCGTGGTGATTACCGCGTTCGGTGTATCCCAGTCGCTGGGCCTGGGCGTGGTGCAAATGAACACCGGGCTGAGCCAGGTGTTTGACCTGCCCATCAGCCTGGGGGTGCAGATCACCTTGATCGTGc                                                                                                     <  1:2954451/150‑1 (MQ=255)
                                                       aTGTGGTGGACATCCTTGCCGTGGTGATTACCGCGTTCGGTGTATCCCAGTCGCTGGGCCTGGGCGTGGTGCAAATGAACACCGGGCTGAGCCAGGTGTTTGACCTGCCCATCAGCCTGGGGGTGCAGATCACCTTGATCGTGCTCATTa                                                                                               <  2:2878972/150‑1 (MQ=255)
                                                       aTGTGGTGGACATCCTTGCCGTGGTGATTACCGCGTTCGGTGTATCCCAGTCGCTGGGCCTGGGCGTGGTGCAAATGAACACCGGGCTGAGCCAGGTGTTTGACCTGCCCATCAGCCTGGGGGTGCAGATCACCTTGATCGTGCTCATTa                                                                                               >  2:4844328/1‑150 (MQ=255)
                                                          tggtggACATCCTTGCCGTGGTGATTACCGCGTTCGGTGTATCCCAGTCGCTGGGCCTGGGCGTGGTGCAAATGAACACCGGGCTGAGCCAGGTGTTTGACCTGCCCATCAGCCTGGGGGTGCAGATCACCTTGATCGTGCTCATTAccc                                                                                            >  2:3497114/1‑150 (MQ=255)
                                                             tggACATCCTTGCCGTGGTGATTACCGCGTTCGGTGTATCCCAGTCGCTGGGCCTGGGCGTGGTGCAAATGAACACCGGGCTGAGCCAGGTGTTTGACCTGCCCATCAGCCTGGGGGTGCAGATCACCTTGATCGTGCTCATTACCCTCg                                                                                         >  1:1820994/1‑150 (MQ=255)
                                                              ggACATCCTTGCCGTGGTGATTACCGCGTTCGGTGTATCCCAGTCGCTGGGCCTGGGCGTGGTGCAAATGAACACCGGGCTGAGCCAGGTGTTTGACCTGCCCATCAGCCTGGGGGTGCAGATCACCTTGATCGTGCTCATTACCCTCGt                                                                                        >  1:5233262/1‑150 (MQ=255)
                                                              ggACATCCTTGCCGTGGTGATTACCGCGTTCGGTGTATCCCAGTCGCTGGGCCTGGGCGTGGTGCAAATGAACACCGGGCTGAGCCAGGTGTTTGACCTGCCCATCAGCCTGGGGGTGCAGATCACCTTGATCGTGCTCATTACCCTCGt                                                                                        <  1:4216060/150‑1 (MQ=255)
                                                              ggACATCCTTGCCGTGGTGATTACCGCGTTCGGTGTATCCCAGTCGCTGGGCCTGGGCGTGGTGCAAATGAACACCGGGCTGAGCCAGGTGTTTGACCTGCCCATCAGCCTGGGGGTGCAGATCACCTTGATCGTGCTCATTACCCTCGt                                                                                        >  1:1001126/1‑150 (MQ=255)
                                                                aCATCCTTGCCGTGGTGATTACCGCGTTCGGTGTATCCCAGTCGCTGGGCCTGGGCGTGGTGCAAATGAACACCGGGCTGAGCCAGGTGTTTGACCTGCCCATCAGCCTGGGGGTGCAGATCACCTTGATCGTGCTCATTACCCTCGTta                                                                                      <  1:2500808/150‑1 (MQ=255)
                                                                 cATCCTTGCCGTGGTGATTACCGCGTTCGGTGTATCCCAGTCGCTGGGCCTGGGCGTGGTGCAAATGAACACCGGGCTGAGCCAGGTGTTTGACCTGCCCATCAGCCTGGGGGTGCAGATCACCTTGATCGTGCTCATTACCCTCGTtac                                                                                     <  2:4549815/150‑1 (MQ=255)
                                                                 cATCCTTGCCGTGGTGATTACCGCGTTCGGTGTATCCCAGTCGCTGGGCCTGGGCGTGGTGCAAATGAACACCGGGCTGAGCCAGGTGTTTGACCTGCCCATCAGCCTGGGGGTGCAGATCACCTTGATCGTGCTCATTACCCTCGTtac                                                                                     >  2:2518380/1‑150 (MQ=255)
                                                                  aTCCTTGCCGTGGTGATTACCGCGTTCGGTGTATCCCAGTCGCTGGGCCTGGGCGTGGTGCAAATGAACACCGGGCTGAGCCAGGTGTTTGACCTGCCCATCAGCCTGGGGGTGCAGATCACCTTGATCGTGCTCATTACCCTCGTtact                                                                                    >  2:4396160/1‑150 (MQ=255)
                                                                      ttGCCGTGGTGATTACCGCGTTCGGTGTATCCCAGTCGCTGGGCCTGGGCGTGGTGCAAATGAACACCGGGCTGAGCCAGGTGTTTGACCTGCCCATCAGCCTGGGGGTGCAGATCACCTTGATCGTGc                                                                                                     <  1:629035/129‑1 (MQ=255)
                                                                           gtggtgATTACCGCGTTCGGTGTATCCCAGTCGCTGGGCCTGGGCGTGGTGCAAATGAACACCGGGCTGAGCCAGGTGTTTGACCTGCCCATCAGCCTGGGGGTGCAGATCACCTTGATCGTGCTCATTACCCTCGTTACTACCGTATCg                                                                           >  2:5646/1‑150 (MQ=255)
                                                                              gtgATTACCGCGTTCGGTGTATCCCAGTCGCTGGGCCTGGGCGTGGTGCAAATGAACACCGGGCTGAGCCAGGTGTTTGACCTGCCCATCAGCCTGGGGGTGCAGATCACCTTGATCGTGCTCATTACCCTCGTTACTACCGTATCGGTg                                                                        >  2:330056/1‑150 (MQ=255)
                                                                               tgATTACCGCGTTCGGTGTATCCCAGTCGCTGGGCCTGGGCGTGGTGCAAATGAACACCGGGCTGAGCCAGGTGTTTGACCTGCCCATCAGCCTGGGGGTGCAGATCACCTTGATCGTGCTCATTACCCTCGTTACTACCGTATCGGTGa                                                                       <  2:638162/150‑1 (MQ=255)
                                                                               tgATTACCGCGTTCGGTGTATCCCAGTCGCTGGGCCTGGGCGTGGTGCAAATGAACACCGGGCTGAGCCAGGTGTTTGACCTGCCCATCAGCCTGGGGGTGCAGATCACCTTGATCGTGCTCATTACCCTCGTTACTACCGTATCGGTGa                                                                       <  1:5366974/150‑1 (MQ=255)
                                                                               tgATTACCGCGTTCGGTGTATCCCAGTCGCTGGGCCTGGGCGTGGTGCAAATGAACACCGGGCTGAGCCAGGTGTTTGACCTGCCCATCAGCCTGGGGGTGCAGATCACCTTGATCGTGCTCATTACCCTCGTTACTACCGTATCGGTGa                                                                       <  2:2621571/150‑1 (MQ=255)
                                                                                 aTTACCGCGTTCGGTGTATCCCAGTCGCTGGGCCTGGGCGTGGTGCAAATGAACACCGGGCTGAGCCAGGTGTTTGACCTGCCCATCAGCCTGGGGGTGCAGATCACCTTGATCGTGCTCATTACCCTCGTTACTACCGTATCGGTGATg                                                                     <  2:4244979/150‑1 (MQ=255)
                                                                                  ttACCGCGTTCGGTGTATCCCAGTCGCTGGGCCTGGGCGTGGTGCAAATGAACACCGGGCTGAGCCAGGTGTTTGACCTGCCCATCAGCCTGGGGGTGCAGATCACCTTGATCGTGCTCATTACCCTCGTTACTACCGTATCGGTGATgg                                                                    >  2:3978625/1‑150 (MQ=255)
                                                                                       gcgTTCGGTGTATCCCAGTCGCTGGGCCTGGGCGTGGTGCAAATGAACACCGGGCTGAGCCAGGTGTTTGACCTGCCCATCAGCCTgg                                                                                                                             >  1:4488841/1‑88 (MQ=255)
                                                                                       gcgTTCGGTGTATCCCAGTCGCTGGGCCTGGGCGTGGTGCAAATGAACACCGGGCTGAGCCAGGTGTTTGACCTGCCCATCAGCCTgg                                                                                                                             <  2:4488841/88‑1 (MQ=255)
                                                                                                 tATCCCAGTCGCTGGGCCTGGGCGTGGTGCAAATGAACACCGGGCTGAGCCAGGTGTTTGACCTGCCCATCAGCCTGGGGGTGCAGATCACCTTGATCGTGCTCATTACCCTCGTTACTACCGTATCGGTGATGGCCGGCGTGTCGCGTg                                                     >  2:541820/1‑150 (MQ=255)
                                                                                                   tCCCAGTCGCTGGGCCTGGGCGTGGTGCAAATGAACACCGGGCTGAGCCAGGTGTTTGACCTGCCCATCAGCCTGGGGGTGCAGATCACCTTGATCGTGCTCATTACCCTCGTTACTACCGTATCGGTGATGGCCGGCGTGTCGCGTGGc                                                   <  1:4828740/150‑1 (MQ=255)
                                                                                                   tCCCAGTCGCTGGGCCTGGGCGTGGTGCAAATGAACACCGGGCTGAGCCAGGTGTTTGACCTGCCCATCAGCCTGGGGGTGCAGATCACCTTGATCGTGCTCATTACCCTCGTTACTACCGTATCGGTGATGGCCGGCGTGTCGCGTGGc                                                   >  1:196798/1‑150 (MQ=255)
                                                                                                    cccAGTCGCTGGGCCTGGGCGTGGTGCAAATGAACACCGGGCTGAGCCAGGTGTATGACCTGCCCATCAGCCTGGGGGTGCAGATCACCTTGATCGTGCTCATTACCCTCGTCACTACCGTATCGGTGATGGCCGGCGTGTCGCGTGGCa                                                  >  2:706034/1‑150 (MQ=255)
                                                                                                     ccAGTCGCTGGGCCTGGGCGTGGTGCAAATGAACACCGGGCTGAGCCAGGTGTTTGACCTGCCCATCAGCCTGGGGGTGCAGATCACCTTGATCGTGCTCATTACCCTCGTTACTACCGTATCGGTGATGGCCGGCGTGTCGCGTGGCAt                                                 >  1:4774551/1‑150 (MQ=255)
                                                                                                     ccAGTCGCTGGGCCTGGGCGTGGTGCAAATGAACACCGGGCTGAGCCAGGTGTTTGACCTGCCCATCAGCCTGGGGGTGCAGATCACCTTGATCGTGCTCATTACCCTCGTTACTACCGTATCGGTGATGGCCGGCGTGTCGCGTGGCAt                                                 >  2:1396454/1‑150 (MQ=255)
                                                                                                      cAGTCGCTGGGCCTGGGCGTGGTGCAAATGAACACCGGGCTGAGCCAGGTGTTTGACCTGCCCATCAGCCTGGGGGTGCAGATCACCTTGATCGTGCTCATTACCCTCGTTACTACCGTATCGGTGATGGCCGGCGTGTCGCGTGGCATg                                                >  1:3197689/1‑150 (MQ=255)
                                                                                                        gTCGCTGGGCCTGGGCGTGGTGCAAATGAACACCGGGCTGAGCCAGGTGTTTGACCTGCCCATCAGCCTGGGGGTGCAGATCACCTTGATCGTGCTCATTACCCTCGTTACTACCGTATCGGTGATGGCCGGCGTGTCGCGTGGCATGaa                                              >  2:2814136/1‑150 (MQ=255)
                                                                                                          cGCTGGGCCTGGGCGTGGTGCAAATGAACACCGGGCTGAGCCAGGTGTTTGACCTGCCCATCAGCCTGGGGGTGCAGATCACCTTGATCGTGCTCATTACCCTCGTTACTACCGTATCGGTGATGGCCGGCGTGTCGCGTGGCATGAAgc                                            >  1:4993927/1‑150 (MQ=255)
                                                                                                           gCTGGGCCTGGGCGTGGTGCAAATGAACACCGGGCTGAGCCAGGTGTTTGACCTGCCCATCAGCCTGGGGGTGCAGATCACCTTGATCGTGCTCATTACCCTCGTTACTACCGTATCGGTGATGGCCGGCGTGTCGCGTGGCATg                                                >  1:4766015/1‑145 (MQ=255)
                                                                                                           gCTGGGCCTGGGCGTGGTGCAAATGAACACCGGGCTGAGCCAGGTGTTTGACCTGCCCATCAGCCTGGGGGTGCAGATCACCTTGATCGTGCTCATTACCCTCGTTACTACCGTATCGGTGATGGCCGGCGTGTCGCGTGGCATGAAgcg                                           >  1:790688/1‑150 (MQ=255)
                                                                                                            ctgggcctgggcGTGGTGCAAATGAACACCGGGCTGAGCCAGGTGTTTGACCTGCCCATCAGCCTGGGGGTGCAGATCACCTTGATCGTGCTCATTACCCTCGTTACTACCGTATCGGTGATGGCCGGCGTGTCGCGTGGCATGAAgcgc                                          >  2:3724789/1‑150 (MQ=255)
                                                                                                                      gcGTGGTGCAAATGAACACCGGGCTGAGCCAGGTGTTTGACCTGCCCATCAGCCTGGGGGTGCAGATCACCTTGATCGTGCTCATTACCCTCGTTACTACCGTATCGGTGATGGCCGGCGTGTCGCGTGGCATGAAGCGCCTTTCggaat                                <  1:2397687/150‑1 (MQ=255)
                                                                                                                      gcGTGGTGCAAATGAACACCGGGCTGAGCCAGGTGTTTGACCTGCCCATCAGCCTGGGGGTGCAGATCACCTTGATCGTGCTCATTACCCTCGTTACTACCGTATCGGTGATGGCCGGCGTGTCGCGTGGCATGAAGCGCCTTTCggaat                                <  2:4575335/150‑1 (MQ=255)
                                                                                                                      gcGTGGTGCAAATGAACACCGGGCTGAGCCAGGTGTTTGACCTGCCCATCAGCCTGGGGGTGCAGATCACCTTGATCGTGCTCATTACCCTCGTTACTACCGTATCGGTGATGGCCGGCGTGTCGCGTGGCATGAAGCGCCTTTCggaat                                <  2:2941745/150‑1 (MQ=255)
                                                                                                                         tggtgCAAATGAACACCGGGCTGAGCCAGGTGTTTGACCTGCCCATCAGCCTGGGGGTGCAGATCACCTTGATCGTGCTCATTACCCTCGTTACTACCGTATCGGTGATGGCCGGCGTGTCGCGTGGCATGAAGCGCCTTTCggaatgga                             <  1:3978625/150‑1 (MQ=255)
                                                                                                                               aaaTGAACACCGGGCTGAGCCAGGTGTTTGACCTGCCCATCAGCCTGGGGGTGCAGATCACCTTGATCGTGCTCATTACCCTCGTTACTACCGTATCGGTGATGGCCGGCGTGTCGCGTGGCATGAAGCGCCTTTCGGAATGGAATATGc                       <  1:4285924/150‑1 (MQ=255)
                                                                                                                               aaaTGAACACCGGGCTGAGCCAGGTGTTTGACCTGCCCATCAGCCTGGGGGTGCAGATCACCTTGATCGTGCTCATTACCCTCGTTACTACCGTATCGGTGATGGCCGGCGTGTCGCGTGGCATGAAGCGCCTTTCGGAATGGAATATGc                       <  2:1133574/150‑1 (MQ=255)
                                                                                                                               aaaTGAACACCGGGCTGAGCCAGGTGTTTGACCTGCCCATCAGCCTGGGGGTGCAGATCACCTTGATCGTGCTCATTACCCTCGTTACTACCGTATCGGTGATGGCCGGCGTGTAGCGTGGCATGAAGCGCCTTTCGGAATGGAATATGc                       >  1:91248/1‑150 (MQ=255)
                                                                                                                                aaTGAACACCGGGCTGAGCCAGGTGTTTGACCTGCCCATCAGCCTGGGGGTGCAGATCACCTTGATCGTGCTCATTACCCTCGTTACTACCGTATCGGTGATGGCCGGCGTGTCGCGTGGCATGAAg                                             >  2:759680/1‑127 (MQ=255)
                                                                                                                                aaTGAACACCGGGCTGAGCCAGGTGTTTGACCTGCCCATCAGCCTGGGGGTGCAGATCACCTTGATCGTGCTCATTACCCTCGTTACTACCGTATCGGTGATGGCCGGCGTGTCGCGTGGCATGAAg                                             <  1:759680/127‑1 (MQ=255)
                                                                                                                                aaTGAACACCGGGCTGAGCCAGGTGTTTGACCTGCCCATCAGCCTGGGGGTGCAGATCACCTTGATCGTGCTCATTACCCTCGTTACGACCGTATCGGTGATGGCCGGCGTGTCGCGTGGCATGAAGCGCCTTTCGGAATGGAATATGCt                      <  2:190167/150‑1 (MQ=255)
                                                                                                                                  caaaCACCTGGCTCAGCCCGGTGTTTGACCTGCCCATCAGCCTGGGGGTGCAGATCACCTTGATCGTGCTCATTACCCTCGTTACTACCGTAt                                                                             >  2:309394/3‑93 (MQ=39)
                                                                                                                                  caaaCACCTGGCTCAGCCCGGTGTTTGACCTGCCCATCAGCCTGGGGGTGCAGATCACCTTGATCGTGCTCATTACCCTCGTTACTACCGTAt                                                                             <  1:309394/91‑1 (MQ=39)
                                                                                                                                   aaaCACCTGGCTCAGCCCGGTGTTTGACCTGCCCATCAGCCTGGGGGTGCAGATCACCTTGATCGTGCTCATTACCCTCGTTACTACCGTCTCGGTGATGGCCGGCGTGTCGCGGGGCATGAAGCGCCTTTCGGAATGGACTctgctggg                   >  2:5038957/2‑148 (MQ=255)
                                                                                                                                   aaaCACCTGGCTCAGCCCGGTGTTTGACCTGCCCATCAGCCTGGGGGTGCAGATCACCTTGATCGTGCCCATTACCCTCGTTACTACCGTATCGGTGATGGCCGGCGTGTCGCGTGGCATGAAGCGCCGTTCGGAATGGAATAGGCtgtt                   >  2:236442/2‑150 (MQ=255)
                                                                                                                                     acacCGGGCTGAGCCAGGTGTTTGACCTGCCCATCAGCCTGGGGGTGCAGATCACCTTGATCGTGCTCATTACCCTCGTTACTACCGTATCGGTGATGGCCGGCGTGTCGCGTGGCATGAAGCGCCTTTCGGAATGGAATATGCtgttgt                 <  1:263627/150‑1 (MQ=255)
                                                                                                                                      cacCGGGCTGAGCCAGGTGTTTGACCTGCCCATCAGCCTGGGGGTGCAGATCACCTTGATCGTGCTCATTACCCTCGTTACTACCGTATCGGTGATGGCCGGCGTGTCGCGTGGCATGAAGCGCCTTTCGGAATGGAATATGCTGTTGTc                >  1:4542289/1‑150 (MQ=255)
                                                                                                                                      cacCGGGCTGAGCCAGGTGTTTGACCTGCCCATCAGCCTGGGGGTGCAGATCACCTTGATCGTGCTCATTACCCTCGTTACTACCGTATCGGTGATGGCCGGCGTGTCGCGTGGCATGAAGCGCCTTTCGGAATGGAATATGCTGTTGTc                >  1:5213450/1‑150 (MQ=255)
                                                                                                                                       acCGGGCTGAGCCAGGTGTTTGACCTGCCCATCAGCCTGGGGGTGCAGATCACCTTGATCGTGCTCATTACCCTCGTTACTACCGTATCGGTGATGGCCGGCGTGTCGCGTGGCATGAAGCGCCTTTCGGAATGGAATATGCTGTTGTCg               >  2:2940234/1‑150 (MQ=255)
                                                                                                                                       acCGGGCTGAGCCAGGTGCTTGACCTGCCCATCAGCCTGGGGGTGCAGATCACCTTGATCGTGCTCATTACCCTCGTTACTACCGTATCGGTGATGGCCGGCGTGTCGCGTGGCATGAAGCGCCTTTCGGAATGGAATATGCTGTTGTCg               <  2:271721/150‑1 (MQ=255)
                                                                                                                                        ccGGGCTGAGCCAGGTGTTTGACCTGCCCATCAGCCTGGGGGTGCAGATCACCTTGATCGTGCTCATTACCCTCGTTACTACCGTATCGGTGATGGCCGGCGTGTCGCGTGGCATGAAGCGCCTTTCGGAATGGAATATGCTGTTGTCgg              <  1:1089638/150‑1 (MQ=255)
                                                                                                                                         cGGGCTGAGCCAGGTGTTTGACCTGCCCATCAGCCTGGGGGTGCAGATCACCTTGATCGTGCTCATTACCCTCGTTACTACCGTATCGGTGATGGCCGGCGTGTCGCGTGGCATGAAGCGCCTTTCGGAATGGAATATGCTGTTGTCggt             <  2:3143123/150‑1 (MQ=255)
                                                                                                                                           ggCTGAGCCAGGTGTTTGACCTGCCCATCAGCCTGGGGGTGCAGATCACCTTGATCGTGCTCATTACCCTCGTTACTACCGTATCGGTGATGGCCGGCGTGTCGCGTGGCATGAAGCGCCTTTCGGAATGGAATATGCTGTTGTCggtgg           <  2:2129256/150‑1 (MQ=255)
                                                                                                                                                 gCCAGGTGTTTGACCTGCCCATCAGCCTGGGGGTGCAGATCACCTTGATCGTGCTCATTACCCTCGTTACTACCGTATCGGTGATGGCCGGCGTGTCGCGTGGCATGAAGCGCCTTTCGGAATGGAATATGCTGTTGTCGGTGGTGCtgg     <  2:4774551/150‑1 (MQ=255)
                                                                                                                                                  ccAGGTGTTTGACCTGCCCATCAGCCTGGGGGTGCAGATCACCTTGATCGTGCTCATTACCCTCGTTACTACCGTATCGGTGATGGCCGGCGTGTCGCGTGGCATGAAGCGCCTTTCGGAATGGAATATGCTGTTGTCGGTGGTGCtggt    >  2:4710735/1‑150 (MQ=255)
                                                                                                                                                   cAGGTGTTTGACCTGCCCATCAGCCTGGGGGTGCAGATCACCTTGATCGTGCTCATTACCCTCGTTACTACCGTATCGGTGATGGCCGGCGTGTCGCGTGGCATGAAGCGCCTTTCGGAATGGAata                          >  1:2342154/1‑127 (MQ=255)
                                                                                                                                                    aGGTGTTTGACCTGCCCATCAGCCTGGGGGTGCAGATCACCTTGATCGTGCTCATTACCCTCGTTACTACCGTATCGGTGATGGCCGGCGTGTCGCGTGGCATGAAGCGCCTTTCGGAATGGAATATGCTGTTGTCGGTGGTGCtggtgg  <  1:2857915/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
ATGCGCTCGATCCTTTACCCGCTGATTGGTGAGCGCATCTATGGGCCCATCGGACATGTGGTGGACATCCTTGCCGTGGTGATTACCGCGTTCGGTGTATCCCAGTCGCTGGGCCTGGGCGTGGTGCAAATGAACACCGGGCTGAGCCAGGTGTTTGACCTGCCCATCAGCCTGGGGGTGCAGATCACCTTGATCGTGCTCATTACCCTCGTTACTACCGTATCGGTGATGGCCGGCGTGTCGCGTGGCATGAAGCGCCTTTCGGAATGGAATATGCTGTTGTCGGTGGTGCTGGTGG  >  NC_002947/6127697‑6127994

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: