Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,124,287 G→A 5.6% A240T (GCC→ACC)  PP_3628 → choline/carnitine/betaine transporter family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,124,2870GA 5.6% 211.8 / 2.5 71A240T (GCC→ACC) PP_3628choline/carnitine/betaine transporter family protein
Reads supporting (aligned to +/- strand):  ref base G (31/36);  new base A (2/2);  total (33/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.85e-01

GCCGCTCGCAGCTGATTTCCGCCGCCTTCCGGCCGCTGATCGGCAAGCACGCAAACGGCCCTCTCGGCCGCCTGATCGACATGATGGCAATCTTCGCCACGCTGTTCGGCTCGGCCGCCTCGCTGGGCCTTGGGGCGCTGCAGATTGCCGGCGGCCTGGAATACAACGGCTGGATCGAAAGCCCCGGTAAACTGTTCTACATCTCGATCATCACTCTGCTGACCATTGCCTTCGTGACCTCGGCGGTGTCGGGCATCGGCAAGGGTATCCAGTGG  >  NC_002947/4124141‑4124415
                                                                                                                                                  |                                                                                                                                
gCCGCTCGCAGCTGATTTCCGCCGCCTTCCGGCCGCTGATCGGCAAGCACGCAAACGGCCCTCTCGGCCGCCTGATCGACATGATGGCAATCTTCGCCACGCTGTTCGGCTCGGCCGCCTCGCTGGGCCTTGGGGCGCTGCAGATTACcg                                                                                                                               >  2:1672030/1‑150 (MQ=255)
gCCGCTCGCAGCTGATTTCCGCCGCCTTCCGGCCGCTGATCGGCAAGCACGCAAACGGCCCTCTCGGCCGCCTGATCGACATGATGGCAATCTTCGCCACGCTGTTCGGCTCGGCCGCCTCGCTGGGCCTTGGGGCGCTGCAGATTACcg                                                                                                                               >  2:1671819/1‑150 (MQ=255)
  cGCTCGCAGCTGATTTCCGCCGCCTTCCGGCCGCTGATCGGCAAGCACGCAAACGGCCCTCTCGGCCGCCTGATCGACATGATGGCAATCTTCGCCACGCTGTTCGGCTCGGCCGCCTCGCTGGGCCTTGGGGCGCTGCAGATTGCcggc                                                                                                                             >  2:1623217/1‑150 (MQ=255)
            tGATTTCCGCCGCCTTCCGGCCGCTGATCGGCAAGCACGCAAACGGCCCTCTCGGCCGCCTGATCGACATGATGGCAATCTTCGCCACGCTGTTCGGCTCGGCCTCCTCGCTGGGCCTTGGGGCGCTGCAGATTGCCGGCGGCCt                                                                                                                        >  2:1581541/1‑145 (MQ=255)
            tGATTTCCGCCGCCTTCCGGCCGCTGATCGGCAAGCACGCAAACGGCCCTCTCGGCCGCCTGATCGACATGATGGCAATCTTCGCCACGCTGTTCGGCTCGGCCTCCTCGCTGGGCCTTGGGGCGCTGCAGATTGCCGGCGGCCt                                                                                                                        <  1:1581541/145‑1 (MQ=255)
            tGATTTCCGCCGCCTTCCGGCCGCTGATCGGCAAGCACGCAAACGGCCCTCTCGGCCGCCTGATCGACATGATGGCAATCTTCGCCACGCTGTTCGGCTCGGCCGCCTCGCTGGGCCTTGGGGCGCTGCAGATTg                                                                                                                                  <  2:1611519/135‑1 (MQ=255)
            tGATTTCCGCCGCCTTCCGGCCGCTGATCGGCAAGCACGCAAACGGCCCTCTCGGCCGCCTGATCGACATGATGGCAATCTTCGCCACGCTGTTCGGCTCGGCCGCCTCGCTGGGCCTTGGGGCGCTGCAGATTg                                                                                                                                  >  1:1611519/1‑135 (MQ=255)
            tGATTTCCGCCGCCTTCCGGCCGCTGATCGGCAAGCACGCAAACGGCCCTCTCGGCCGCCTGATCGACATGATGGCAATCTTCGCCACGCTGTTCGGCTCGGCCGCCTCGCTGGGCCTTGGGGCGCTGCAGATTGCCGGCGGCCTGGAAt                                                                                                                   >  2:1153423/1‑150 (MQ=255)
             gATTTCCGCCGCCTTCCGGCCGCTGATCGGCAAGCACGCAAACGGCCCTCTCGGCCGCCTGATCGACATGATGGCAATCTTCGCCACGCTGTTCGGCTCGGCCGCCTCGCTGGGCCTTGGGGCGCTGCAGATTGCCGGCGGCCTGGAATa                                                                                                                  >  2:515307/1‑150 (MQ=255)
                ttCCGCCGCCTTCCGGCCGCTGATCGGCAAGCACGCAAACGGCCCTCTCGGCCGCCTGATCGACATGATGGCAATCTTCGCCACGCTGTTCGGCTCGGCCGCCTCGCTGGGCCTTGGGGCGCTGCAGATTGCCGGCGGCCTGGAATACaa                                                                                                               <  2:1388833/150‑1 (MQ=255)
                 tCCGCCGCCTTCCGGCCGCTGATCGGCAAGCACGCAAACGGCCCTCTCGGCCGCCTGATCGACATGATGGCAATCTTCGCCACGCTGTTCGGCTCGGCCGCCTCGCTGGGCCTTGGGGCGCTGCAGATTGCCGGCGGCCTGGAATACAAc                                                                                                              <  1:1435208/150‑1 (MQ=255)
                     ccgccTTCCGGCCGCTGATCGGCAAGCACGCAAACGGCCCTCTCGGCCGCCTGATCGACATGATGGCAATCTTCGCCACGCTGTTCGGCTCGGCCGCCTCGCTGGGCCTTGGGGCGCTGCAGATTGCCGGCGGCCTGGAATACAACGGCt                                                                                                          <  1:1652820/150‑1 (MQ=255)
                         cTTCCGGCCGCTGATCGGCAAGCACGCAAACGGCCCTCTCGGCCGCCTGATCGACATGATGGCAATCTTCGCCACGCTGTTCGGCTCGGCCGCCTCGCTGGGCCTTGGGGCGCTGCAGATTGCCGGCGGCCTGGAATACAACGGCTGGAt                                                                                                      <  1:555168/150‑1 (MQ=255)
                          ttCCGGCCGCTGATCGGCAAGCACGCAAACGGCCCTCTCGGCCGCCTGATCGACATGATGGCAATCTTCGCCACGCTGTTCGGCTCGGCCGCCTCGCTGGGCCTTGGGGCGCTGCAGATTGCCGGCGGCCTGGAATACAACGGCTGGATc                                                                                                     <  2:498471/150‑1 (MQ=255)
                          ttCCGGCCGCTGATCGGCAAGCACGCAAACGGCCCTCTCGGCCGCCTGATCGACATGATGGCAATCTTCGCCACGCTGTTCGGCTCGGCCGCCTCGCTGGGCCTTGGGGCGCTGCAGATTGCCGGCGGCCTGGAATACAACGGCTGGATc                                                                                                     >  1:1159889/1‑150 (MQ=255)
                             cGGCCGCTGATCGGCAAGCACGCAAACGGCCCTCTCGGCCGCCTGATCGACATGATGGCAATCTTCGCCACGCTGTTCGGCTCGGCCGCCTCGCTGGGCCTTGGGGCGCTGCAGATTGCCGGCGGCCTGGAATACAACGGCTGGATCGaa                                                                                                  <  2:1604069/150‑1 (MQ=255)
                              ggCCGCTGATCGGCAAGCACGCAAACGGCCCTCTCGGCCGCCTGATCGACATGATGGCAATCTTCGCCACGCTGTTCGGCTCGGCCGCCTCGCTGGGCCTTGGGGCGCTGCAGATTGCCGGCGGCCTGGAATACAACGGCTGGATCGaaa                                                                                                 >  2:258549/1‑150 (MQ=255)
                              ggCCGCTGATCGGCAAGCACGCAAACGGCCCTCTCGGCCGCCTGATCGACATGATGGCAATCTTCGCCACGCTGTTCGGCTCGGCCGCCTCGCTGGGCCTTGGGGCGCTGCAGATTGCCGGAGGCATGGGATACAACGGCTGGATCGaaa                                                                                                 >  2:932202/1‑150 (MQ=255)
                               gCCGCTGATCGGCAAGCACGCAAACGGCCCTCTCGGCCGCCTGATCGACATGATGGCAATCTTCGCCACGCTGTTCGGCTCGGCCGCCTCGCTGGGCCTTGGGGCGCTGCAGATTGCCGGCGGCCTGGAATACAACGGCTGGATCGAAAg                                                                                                >  1:879635/1‑150 (MQ=255)
                                  gCTGATCGGCAAGCACGCAAACGGCCCTCTCGGCCGCCTGATCGACATGATGGCAATCTTCGCCACGCTGTTCGGCTCGGCCGCCTCGCTGGGCCTTGGGGCGCTGCAGATTGCCGGCGGCCTGGAATACAACg                                                                                                             >  2:938006/1‑134 (MQ=255)
                                  gCTGATCGGCAAGCACGCAAACGGCCCTCTCGGCCGCCTGATCGACATGATGGCAATCTTCGCCACGCTGTTCGGCTCGGCCGCCTCGCTGGGCCTTGGGGCGCTGCAGATTGCCGGCGGCCTGGAATACAACg                                                                                                             <  1:938006/134‑1 (MQ=255)
                                    tGATCGGCAAGCACGCCAGTCTCCCTCTCGGCCGCCTGATCGACATGATGGCAATCTTCGCCACGCTGTTCGGCTCGGCCGCCTCGCTGGGCCTTGGGGCGCTGCAGATTGCCGGCGGCCTGGAATACAACGGCTGGATCGAAAGCCCCg                                                                                           <  1:99571/150‑1 (MQ=255)
                                      aTCGGCAAGCACGCAAACGGCCCTCTCGGCCGCCTGATCGACATGATGGCAATCTTCGCCACGCTGTTCGGCTCGGCCGCCTCGCTGGGCCTTGGGGCGCTGCAGATTGCCGGCGGCCTGGAATACAACGGCTGGATCGAAAGCCCCGGt                                                                                         >  1:652205/1‑150 (MQ=255)
                                         ggCAAGCACGCAAACGGCCCTCTCGGCCGCCTGATCGACATGATGGCAATCTTCGCCACGCTGTTCGGCTCGGCCGCCTCGCTGGGCCTTGGGGCGCTGCAGATTGCCGGCGGCCTGGAATACAACGGCTGGATCGAAAGCCCCGGTaaa                                                                                      >  1:1250528/1‑150 (MQ=255)
                                          gCAAGCACGCAAACGGCCCTCTCGGCCGCCTGATCGACATGATGGCAATCTTCGCCACGCTGTTCGGCTCGGCCGCCTCGCTGGGCCTTGGGGCGCTGCAGATTGCCGGCGGCCTGGAATACAACg                                                                                                             <  1:1739087/126‑1 (MQ=255)
                                          gCAAGCACGCAAACGGCCCTCTCGGCCGCCTGATCGACATGATGGCAATCTTCGCCACGCTGTTCGGCTCGGCCGCCTCGCTGGGCCTTGGGGCGCTGCAGATTGCCGGCGGCCTGGAATACAACg                                                                                                             >  2:1739087/1‑126 (MQ=255)
                                            aaGCACGCAAACGGCCCTCTCGGCCGCCTGATCGACATGATGGCAATCTTCGCCACGCTGTTCGGCTCGGCCGCCTCGCTGGGCCTTGGGGCGCTGCAGATTGCcg                                                                                                                               >  1:314641/1‑106 (MQ=255)
                                            aaGCACGCAAACGGCCCTCTCGGCCGCCTGATCGACATGATGGCAATCTTCGCCACGCTGTTCGGCTCGGCCGCCTCGCTGGGCCTTGGGGCGCTGCAGATTGCcg                                                                                                                               <  2:314641/106‑1 (MQ=255)
                                              gCACGCAAACGGCCCTCTCGGCCGCCTGATCGACATGATGGCAATCTTCGCCACGCTGTTCGGCTCGGCCGCCTCGCTGGGCCTTGGGGCGCTGCAGATTGCCGGCGGCCTGGAATACAACGGCTGGATCGAAAGCCCCGGTAAACTGtt                                                                                 <  2:1159889/150‑1 (MQ=255)
                                               cACGCAAACGGCCCTCTCGGCCGCCTGATCGACATGATGGCAATCTTCGCCACGCTGTTCGGCTCGGCCGCCTCGCTGGGCCTTGGGGCGCTGCAGATTGCCGGCGGCCTGGAATACAACGGCTGGATCGAAAGCCCCGGTAAACTGTTc                                                                                <  1:1819929/150‑1 (MQ=255)
                                               cACGCAAACGGCCCTCTCGGCCGCCTGATCGACATGATGGCAATCTTCGCCACGCTGTTCGGCTCGGCCGCCTCGCTGGGCCTTGGGGCGCTGCAGATTGCCGGCGGCCTGGAATACAACGGCTGGATCGAAAGCCCCGGTAAACTGTTc                                                                                <  1:515307/150‑1 (MQ=255)
                                                    aaaCGGCCCTCTCGGCCGCCTGATCGACATGATGGCAATCTTCGCCACGCTGTTCGGCTCGGCCGCCTCGCTGGGCCTTGGGGCGCTGCAGATTGCCGGCGGCCTGGAATACAACGGCTGGATCGAAAGccc                                                                                             >  1:856708/1‑132 (MQ=255)
                                                    aaaCGGCCCTCTCGGCCGCCTGATCGACATGATGGCAATCTTCGCCACGCTGTTCGGCTCGGCCGCCTCGCTGGGCCTTGGGGCGCTGCAGATTGCCGGCGGCCTGGAATACAACGGCTGGATCGAAAGccc                                                                                             <  2:856708/132‑1 (MQ=255)
                                                         gCCCTCTCGGCCGCCTGATCGACATGATGGCAATCTTCGCCACGCTGTTCGGCTCGGCCGCCTCGCTGGGCCTTGGGGCGCTGCAGATTGCCGGCGGCCTGGAATACAACGGCTGGATCGAAAGCCCCGGTAAACTGTTCTACATCTCGa                                                                      <  1:901413/150‑1 (MQ=255)
                                                          cccTCTCGGCCGCCTGATCGACATGATGGCAATCTTCGCCACGCTGTTCGGCTCGGCCGCCTCGCTGGGCCTTGGGGCGCTGCAGATTGCCGGCGGCCTGGAATACAACGGCTGGATCGAAAGCCCCGGTAAACTGTTCTACATCTCGat                                                                     <  1:1623217/150‑1 (MQ=255)
                                                           ccTCTCGGCCGCCTGATCGACATGATGGCAATCTTCGCCACGCTGTTCGGCTCGGCCGCCTCGCTGGGCCTTGGGGCGCTGCAGATTGCCGGCGGCCTGGAATACAACGGCTGGATCGAAAGCCCCGGTAAACTGTTCTACATCTCGatc                                                                    <  2:640009/150‑1 (MQ=255)
                                                           ccTCTCGGCCGCCTGATCGACATGATGGCAATCTTCGCCACGCTGTTCGGCTCGGCCGCCTCGCTGGGCCTTGGGGCGCTGCAGATTGCCGGCGGCCTGGAATACAACGGCTGGATCGAAAGCCCCGGTAAACTGTTCTACATCTCGatc                                                                    >  1:242910/1‑150 (MQ=255)
                                                             tctcGGCCGCCTGATCGACATGATGGCAATCTTCGCCACGCTGTTCGGCTCGGCCGCCTCGCTGGGCCTTGGGGCGCTGCAGATTGCCGGCGGCCTGGAATACAACGGCTGGATCGAAAGCCCCGGTAAACTGTTCTACATCTCGatcat                                                                  <  2:923769/150‑1 (MQ=255)
                                                               tcGGCCGCCTGATCGACATGATGGCAATCTTCGCCACGCTGTTCGGCTCGGCCGCCTCGCTGGGCCTTGGGGCGCTGCAGATTGCCGGCGGCCTGGAATACAACGGCTGGATCGAAAGCCCCGGTAAACTGTTCTACATCTCGAtcatca                                                                >  2:625586/1‑150 (MQ=255)
                                                                    cgccTGATCGACATGATGGCAATCTTCGCCACGCTGTTCGGCTCGGCCGCCTCGCTGGGCCTTGGGGCGCTGCAGATTGCCGGCGGCCTGGAATACAACGGCTGGATCGAAAGCCCCGGTAAACTGTTCTACATCTCGATCATCACTctg                                                           >  2:1930727/1‑150 (MQ=255)
                                                                     gccTGATTGACATGATGGCAATCTTCGCCACGCTGTTCGGCTCGGCCGCCTCGCTGGGCCTTGGGGCGCTGCAGATTGCCGGCGGCCTGGAATACAACGGCTGGATCGAAAGCCCCGGTAAACTGTTCTACATCTCGATCATCACTctgc                                                          >  2:192416/1‑150 (MQ=255)
                                                                     gccTGATCGACATGATGGCAATCTTCGCCACGCTGTTCGGCTCGGCCGCCTCGCTGGGCCTTGGGGCGCTGCAGATTGCCGGCGGCCTGGAATACAACGGCTGGATCGAAAGCCCCGGTAAACTGTTCTACATCTCGATCATCACTctgc                                                          <  1:325357/150‑1 (MQ=255)
                                                                     gccTGATCGACATGATGGCAATCTTCGCCACGCTGTTCCGGTCGGCCGCCTCGCTGGGCCTTGGGGCGCTGCAGATTGCCGGCGGCCTGGAATACAACGGCTGGATCGAAAGCCCCGGTAAACTGTTCTACATCTCGATCATCACTctgc                                                          <  1:101397/150‑1 (MQ=255)
                                                                          aTCGACATGATGGCAATCTTCGCCACGCTGTTCGGCTCGGCCGCCTCGCTGGGCCTTGGGGCGCTGCAGATTGCCGGCGGCCTGGAATACAACGGCTGGATCGAAAGCCCCGGTAAACTGtt                                                                                 <  2:1541527/122‑1 (MQ=255)
                                                                          aTCGACATGATGGCAATCTTCGCCACGCTGTTCGGCTCGGCCGCCTCGCTGGGCCTTGGGGCGCTGCAGATTGCCGGCGGCCTGGAATACAACGGCTGGATCGAAAGCCCCGGTAAACTGtt                                                                                 >  1:1541527/1‑122 (MQ=255)
                                                                              aCATGATGGCAATCTTCGCCACGCTGTTCGGCTCGGCCGCCTCGCTGGGCCTTGGGGCGCTGCAGATTGCCGGCGGCCTGGAATACAACGGCTGGATc                                                                                                     <  2:1887698/98‑1 (MQ=255)
                                                                              aCATGATGGCAATCTTCGCCACGCTGTTCGGCTCGGCCGCCTCGCTGGGCCTTGGGGCGCTGCAGATTGCCGGCGGCCTGGAACACAACGGCTGGATc                                                                                                     >  1:1887698/1‑98 (MQ=255)
                                                                               cATGATGGCAATCTTCGCCACGCTGTTCGGCTCGGCCGCCTCGCTGGGCCTTGGGGCGCTGCAGATTGCCGGCGGCCTGGAATACAACGGCTGGATCGAAAGCCCCGGTAAACTGTTCTACATCTCGATCATCACTCTGCTGACCATTGc                                                <  2:652205/150‑1 (MQ=255)
                                                                               cATGATGGCAATCTTCGCCACGCTGTTCGGCTCGGCCGCCTCGCTGGGCCTTGGGGCGCTGCAGATTGCCGGCGGCCTGGAATACAACGGCTGGATCGAAAGCCCCGGTAAACTGTTCTACATCTCGATCATCACTCTGCTGACCATTGc                                                <  1:673471/150‑1 (MQ=255)
                                                                                atgatgGCAATCTTCGCCACGCTGTTCGGCTCGGCCGCCTCGCTGGGCCTTGGGGCGCTGCAGATTGCCGGCGGCCTGGAATa                                                                                                                  <  2:679894/83‑1 (MQ=255)
                                                                                atgatgGCAATCTTCGCCACGCTGTTCGGCTCGGCCGCCTCGCTGGGCCTTGGGGCGCTGCAGATTGCCGGCGGCCTGGAATa                                                                                                                  >  1:679894/1‑83 (MQ=255)
                                                                                     ggCAATCTTCGCCACGCTGTTCGGCTCGGCCGCCTCGCTGGGCCTTGGGGCGCTGCAGATTGCCGGCGGCCTGGAATACAACGGCTGGATCGAAAGCCCCGGTAAACTGTTCTACATCTCGatcat                                                                  <  2:15241/126‑1 (MQ=255)
                                                                                     ggCAATCTTCGCCACGCTGTTCGGCTCGGCCGCCTCGCTGGGCCTTGGGGCGCTGCAGATTGCCGGCGGCCTGGAATACAACGGCTGGATCGAAAGCCCCGGTAAACTGTTCTACATCTCGatcat                                                                  >  1:15241/1‑126 (MQ=255)
                                                                                     ggCAATCTTCGCCACGCTGTTCGGCTCGGCCGCCTCGCTGGGCCTTGGGGCGCTGCAGATTGCCGGCGGCCTGGAATACAACGGCTGGATCGAAAGCCCCGGTAAACTGTTCTACATCTCGATCATCACTCTGCTGACCATTGCCTTCGt                                          >  2:1403251/1‑150 (MQ=255)
                                                                                               gCCACGCTGTTCGGCTCGGCCGCCTCGCTGGGCCTTGGGGCGCTGCAGATTGCCGGCGGCCTGGAATACAACGGCTGGATCGAAAGCCCCGGTAAACTGTTCTACATCTCGATCATCACTCTGCTGACCATTGCCTTCGTGACCTcggcg                                >  2:1150845/1‑150 (MQ=255)
                                                                                                 cACGCTGTTCGGCTCGGCCGCCTCGCTGGGCCTTGGGGCGCTGCAGATTGCCGGCGGCCTGGAATACAACGGCTGGATCGAAAGCCCCGGTAAACTGTTCTACATCTCGATCATCACTCTGCTGACCATTGCCTTCGTGACCTcggcg                                >  2:21266/1‑148 (MQ=255)
                                                                                                 cACGCTGTTCGGCTCGGCCGCCTCGCTGGGCCTTGGGGCGCTGCAGATTGCCGGCGGCCTGGAATACAACGGCTGGATCGAAAGCCCCGGTAAACTGTTCTACATCTCGATCATCACTCTGCTGACCATTGCCTTCGTGACCTcggcg                                <  1:21266/148‑1 (MQ=255)
                                                                                                  aCGCTGTTCGGCTCGGCCGCCTCGCTGGGCCTTGGGGCGCTGCGGATTGCCGGCGGCCTGGAATACAACGGCTGGATCGAAAGCCCCGGTAAACTGTTCTACAtc                                                                          <  1:727550/105‑1 (MQ=255)
                                                                                                  aCGCTGTTCGGCTCGGCCGCCTCGCTGGGCCTTGGGGCGCTGCGGATTGCCGGCGGCCTGGAATACAACGGCTGGATCGAAAGCCCCGGTAAACTGTTCTACAtc                                                                          >  2:727550/1‑105 (MQ=255)
                                                                                                           ggctcggcCGCCTCGCTGGGCCTTGGGGCGCTGCAGATTGCCGGCGGCCTGGAATACAACGGCTGGATCGAAAGCCCCGGTAAACTGTTCTACATCTCGATCATCACTCTGCTGACCATTGCCTTCGTGACCTCGGCGGTGTCGGGCATc                    <  1:1153423/150‑1 (MQ=255)
                                                                                                            gctcggcCGCCTCGCTGGGCCTTGGGGCGCTGCAGATTGCCGGCGGCCTGGAATACAACGGCTGGATCGAAAGCCCCGGTAAACTGTTCTACATCTCGATCATCACTCTGCTGACCATTGCCTTCGTGACCTCGGCGGTGTCGGGCATCg                   >  2:1753614/1‑150 (MQ=255)
                                                                                                                 gccgccTCGCTGGGCCTTGGGGCGCTGCAGATTGCCGGCGGCCTGGAATACAACGGCTGGATCGAAAGCCCCGGTAAACTGTTCTACATCTCGATCATCACTCTGCTGACCATTGCCTTCGTGACCTCGGCGGTGTCgg                         <  2:375578/139‑1 (MQ=255)
                                                                                                                 gccgccTCGCTGGGCCTTGGGGCGCTGCAGATTGCCGGCGGCCTGGAATACAACGGCTGGATCGAAAGCCCCGGTAAACTGTTCTACATCTCGATCATCACTCTGCTGACCATTGCCTTCGTGACCTCGGCGGTGTCgg                         >  1:375578/1‑139 (MQ=255)
                                                                                                                 gccgccTCGCTGGGCCTTGGGGCGCTGCAGATTGCCGGCGGCCTGGAATACAACGGCTGGATCGAAAGCCCCGGTAAACTGTTCTACATCTCGATCATCACTCTGCTGACCATTGCCTTCGTGACCTCGGCGGTGTCGGGCATCGGCAAg              >  2:1846969/1‑150 (MQ=255)
                                                                                                                 gccgccTCGCTGGGCCTTGGGGCGCTGCAGATTGCCGGCGGCCTGGAATACAACGGCTGGATCGAAAGCCCCGGTAAACTGTTCTACATCTCGATCATCACTCTGCTGACCATTGCCTTCGTGACCTCGGCGGTGTCGGGCATCGGCAAg              <  1:258549/150‑1 (MQ=255)
                                                                                                                  ccgccTCGCTGGTGCTTGGGGCGCTGCAGATTGCCGGCGGCCTGGAATACAACGGCTGGATCGAAAGCCCCGGTAAACTGTTCTACATCTCGATCATCACTCTGCTGACCATTGCCTTCGTGACCTCGGCGGTGTCGGGCATCGGCAAgg             <  1:1753614/150‑1 (MQ=255)
                                                                                                                  ccgccTCGCTGGGCCTTGGGGCGCTGCAGATTACCGGCGGCCTGGAATACAACGGCTGGATCGAAAGCCCCGGTAAACTGTTCTACATCTCGATCATCACTCTGCTGACCATTGCCTTCGTGACCTCGGCGGTGTCGGGCATCGGCAAgg             <  1:1671819/150‑1 (MQ=255)
                                                                                                                  ccgccTCGCTGGGCCTTGGGGCGCTGCAGATTACCGGCGGCCTGGAATACAACGGCTGGATCGAAAGCCCCGGTAAACTGTTCTACATCTCGATCATCACTCTGCTGACCATTGCCTTCGTGACCTCGGCGGTGTCGGGCATCGGCAAgg             <  1:1672030/150‑1 (MQ=255)
                                                                                                                      cTCGCTGGGCCTTGGCGCGCTGCAGGTTGCCGGCGGCCTGGAATACAACGGCTGGAGCGAAAGCCCCGGTAAACTGTTCTACATCTCGATGATCACTCTGCTGACCATTGCCTTCGGGACCTCGGCGGTGTCGGGGATCGGCTAGGGTAt         >  1:1314423/1‑150 (MQ=255)
                                                                                                                       tCGCTGGGCCTTGGGGCGCTGCAGATTGCCGGCGGCCTGGAATACAACGGCTGGATCGAAAGCCCCGGTAAACTGTTCTACATCTCGATCATCACTCTGCTGACCATTGCCTTCGTGACCTCGGCGGTGTCGGGCATCGGCaa               >  1:1001502/1‑143 (MQ=255)
                                                                                                                       tCGCTGGGCCTTGGGGCGCTGCAGATTGCCGGCGGCCTGGAATACAACGGCTGGATCGAAAGCCCCGGTAAACTGTTCTACATCTCGATCATCACTCTGCTGACCATTGCCTTCGTGACCTCGGCGGTGTCGGGCATCGGCaa               <  2:1001502/143‑1 (MQ=255)
                                                                                                                           tGGGCCTTGGGGCGCTGCAGATTGCCGGCGGCCTGGAATACAACGGCTGGATCGAAAGCCCCGGTAAACTGTTCTACATCTCGATCATCACTCTGCTGACCATTGCCTTCGTGACCTCGGCGGTGTCGGGCATCGGCAAGGGTATCCAGt    <  1:932202/150‑1 (MQ=255)
                                                                                                                             ggCCTTGGGGCGCTGCAGATTGCCGGCGGCCTGGAATACAACGGCTGGATCGAAAGCCCCGGTAAACTGTTCTACATCTCGATCATCACTCTGCTGACCATTGCCTTCGTGACCTCGGCGGTGTCGGGCATCGGCAAGGGTATCCAGTgg  <  1:192416/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                
GCCGCTCGCAGCTGATTTCCGCCGCCTTCCGGCCGCTGATCGGCAAGCACGCAAACGGCCCTCTCGGCCGCCTGATCGACATGATGGCAATCTTCGCCACGCTGTTCGGCTCGGCCGCCTCGCTGGGCCTTGGGGCGCTGCAGATTGCCGGCGGCCTGGAATACAACGGCTGGATCGAAAGCCCCGGTAAACTGTTCTACATCTCGATCATCACTCTGCTGACCATTGCCTTCGTGACCTCGGCGGTGTCGGGCATCGGCAAGGGTATCCAGTGG  >  NC_002947/4124141‑4124415

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