Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,396,403 T→A 5.4% F278Y (TTC→TAC)  hsdM → type I restriction modification system methylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,396,4030TA 5.4% 219.6 / 3.5 74F278Y (TTC→TAC) hsdMtype I restriction modification system methylase
Reads supporting (aligned to +/- strand):  ref base T (30/40);  new base A (2/2);  total (32/42)
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.70e-01

ACCGCTGCCGCACCTGCTGGCCACCACCAACATGATCCTGCATGGCATCGAAGTGCCCAGCCAGATCCGTCACGACAACACCCTGAGCAAACCGCTGATCAGCTGGGGCCCAAGCGAGCGCGTGCATTGTATCGTCGCCAACCCGCCGTTCGGCGGCATGGAAGAAGACGGTATCGAAACAAATTTTCCTGCCGCTTTCCGCACCCGGGAAACCGCCGATTTGTTCTTGGTATTGATTATGCAGCTGCTCAAAGATGGTGGCCGCGCCGCAGTGGTACTGCCCGATGGCTTCCTTTTT  >  NC_002947/5396254‑5396551
                                                                                                                                                     |                                                                                                                                                    
aCCGCTGCCGCACCTGCTGGCCACCACCAACATGATCCTGCATGGCATCGAAGTGCCCAGCCAGATCCGTCACGACAACACCCTGAGCAAACCGCTGATCAGCTGGGGCCCAAGCGAGCGCGTGCATTGTATCGTCGCCAACCCGCCGtt                                                                                                                                                      >  2:1741217/1‑150 (MQ=255)
aCCGCTGCCGCACCTGCTGGCCACCACCAACATGATCCTGCATGGCATCGAAGTGCCCAGCCAGATCCGTCACGACAACACCCTGAGCAAACCGCTGATCAGCTGGGGCCCAAGCGAGCGCGTGCATTGTATCGTCGCCAACCCGCCGtt                                                                                                                                                      >  1:1224734/1‑150 (MQ=255)
 ccGCTGCCGCACCTGCTGGCCACCACCAACATGATCCTGCATGGCATCGAAGTGCCCAGCCAGATCCGTCACGACAACACCCTGAGCAAACCGCTGATCAGCTGGGGCCCAAGCGAGCGCGTGCATTGTATCGTCGCCAACCCGCCGTTc                                                                                                                                                     <  2:1284978/150‑1 (MQ=255)
   gCTGCCGCACCTGCTGGCCACCACCAACATGATCCTGCATGGCATCGAAGTGCCCAGCCAGATCCGTCACGACAACACCCTGAGCAAACCGCTGATCAGCTGGGGCCCAAGCGAGCGCGTGCATTGTATCGTCGCCAACCCGCCGTTcgg                                                                                                                                                   >  2:1025849/1‑150 (MQ=255)
    cTGCCGCACCTGCTGGCCACCACCAACATGATCCTGCATGGCATCGAAGTGCCCAGCCAGATCCGTCACGACAACACCCTGAGCAAACCGCTGATCAGCTGGTGCCCAAGCGAGCGCGTGCATTGTATCGTCGCCAACCCGCCGTTcggc                                                                                                                                                  <  2:572421/150‑1 (MQ=255)
     tGCCGCACCTGCTGGCCACCACCAACATGATCCTGCATGGCATCGAAGTGCCCAGCCAGATCCGTCACGACAACACCCTGAGCAAACCGCTGATCAGCTGGGGCCCAAGCGAGCGCGTGCATTGTATCGTCGCCAACCCGCCGTTcggcg                                                                                                                                                 <  1:811874/150‑1 (MQ=255)
              tgctgGCCACCACCAACATGATCCTGCATGGCATCGAAGTGCCCAGCCAGATCCGTCACGACAACACCCTGAGCAAACCGCTGATCAGCTGGGGCCCAAGCGAGCGCGTGCATTGTATCGTCGCCAACCCGCCGTTCggcggc                                                                                                                                               <  2:940723/143‑1 (MQ=255)
              tgctgGCCACCACCAACATGATCCTGCATGGCATCGAAGTGCCCAGCCAGATCCGTCACGACAACACCCTGAGCAAACCGCTGATCAGCTGGGGCCCAAGCGAGCGCGTGCATTGTATCGTCGCCAACCCGCCGTTCggcggc                                                                                                                                               >  1:940723/1‑143 (MQ=255)
              tgctgGCCACCACCAACATGATCCTGCATGGCATCGAAGTGCCCAGCCAGATCCGTCACGACAACACCCTGAGCAAACCGCTGATCAGCTGGGGCCCAAGCGAGCGCGTGCATTGTATCGTCGCCAACCCGCCGTTCGGCGGCATGgaag                                                                                                                                        >  1:1763668/1‑150 (MQ=255)
               gctgGCCACCACCAACATGATCCTGCATGGCATCGAAGTGCCCAGCCAGATCCGTCACGACAACACCCTGAGCAAACCGCTGATCAGCTGGGGCCCAAGCGAGCGCGTGCATTGTATCGTCGCCAACCCGCCGTTCGGCGGCATGgaaga                                                                                                                                       >  1:893781/1‑150 (MQ=255)
                   gCCACCACCAACATGATCCTGCATGGCATCGAAGTGCCCAGCCAGATCCGTCACGACAACACCCTGAGCAAACCGCTGATCAGCTGGGGCCCAAGCGAGCGCGTGCATTGTATCGTCGCCAACCCGCCGTTCGGCGGCATGGAAGAAGAc                                                                                                                                   >  1:1312017/1‑150 (MQ=255)
                   gCCACCACCAACATGATCCTGCATGGCATCGAAGTGCCCAGCCAGATCCGTCACGACAACACCCTGAGCAAACCGCTGATCAGCTGGGGCCCAAGCGAGCGCGTGCATTGTATCGTCGCCAACCCGCCGTTCGGCGGCATGGAAGAAGAc                                                                                                                                   >  1:1311664/1‑150 (MQ=255)
                   gCCACCACCAACATGATCCTGCATGGCATCGAAGTGCCCAGCCAGATCCGTCACGACAACACCCTGAGCAAACCGCTGATCAGCTGGGGCCCAAGCGAGCGCGTGCATTGTATCGTCGCCAACCCGCCGAAcggcgggt                                                                                                                                              >  2:911756/1‑137 (MQ=255)
                   gCCACCACCAACATGATCCTGCATGGCATCGAAGTGCCCAGCCAGATCCGTCACGACAACACCCTGAGCAAACCGCTGATCAGCTGGGGCCCAAGCGAGCGCGTGCATTGTATCGTCGCCAACCCGCCGAAcggcgggt                                                                                                                                              <  1:911756/139‑3 (MQ=255)
                   gCCACCACCAACATGATCCTGCATGGCAACGAAGTGCCCAGCCAGATCCGTCACGACAACACCCTGAGCAAACCGCTGATCAGCTGGGGCCCAAGCGAGCGCGTGCATTGTATCGTCGCCAACCCGCCGTTCGGCGGCATGGAAGAAGAc                                                                                                                                   <  2:1361452/150‑1 (MQ=255)
                    ccaccaccaACATGATCCTGCATGGCATCGAAGTGCCCTGCCAGATCCGTCACGACAACACCCTGAGCAAACCGCTGATCAGCTGGGGCCCAAGCGAGCGCGTGCATTGTATCGTCGCCAACCCGCCGTTCGGCGGCATGGAAGAAGACg                                                                                                                                  <  1:1535042/150‑1 (MQ=255)
                    ccaccaccaACATGATCCTGCATGGCATCGAAGTGCCCAGCCAGATCCGTCACGACAACACCCTGAGCAAGCAGCCGATCAGCGTGGTGGCGGCCGTGCGGGTGCGTTGTATCGTCGCCAACCCGCCGTTCGGCGGCATGGAAGAAGACg                                                                                                                                  >  1:1382728/1‑150 (MQ=255)
                     caccaccaACATGATCCTGCATGGCATCGAAGTGCCCAGCCAGATCCGTCTCGACAACACCCTGCGCAAACCGCTGATCAGTGAGGGCCCAAGCGAGCGCGTGCATTGTATCGTCGCCAACCAGCATTTTGGCGCCCTTGAGGAAGTCgg                                                                                                                                 >  1:33773/1‑150 (MQ=255)
                     caccaccaACATGATCCTGCATGGCATCGAAGTGCCCAGCCAGATCCGTCACGACAACACCCTGAGCAAACCGCTGATCAGCTGGGGCCCAAGCGAGCGCGTGCATTGTATCGTCGCCAACCCGCCGTTCGGCGGCATGGAAGAAGACgg                                                                                                                                 >  1:1815034/1‑150 (MQ=255)
                         accaACATGATCCTGCATGGCATCGAAGTGCCCAGCCAGATCCGTCACGACAACACCCTGAGCAAACCGCTGATCAGCTGGGGCCCAAGCGAGCGCGTGCATTGTATCGTCGCCAACCCGCCGAAcggcgggtt                                                                                                                                             >  1:1540870/1‑131 (MQ=255)
                         accaACATGATCCTGCATGGCATCGAAGTGCCCAGCCAGATCCGTCACGACAACACCCTGAGCAAACCGCTGATCAGCTGGGGCCCAAGCGAGCGCGTGCATTGTATCGTCGCCAACCCGCCGAAcggcgggtt                                                                                                                                             <  2:1540870/134‑4 (MQ=255)
                             aCATGATCCTGCATGGCATCGAAGTGCCCAGCCAGATCCGTCACGACAACACCCTGAGCAAACCGCTGATCAGCTGGGGCCCAAGCGAGCGCGTGCATTGTATCGTCGCCAACCCGCCGTTCGGCGGCATGGAAGAAGACGGTATCGaaa                                                                                                                         <  2:1516567/150‑1 (MQ=255)
                             aCATGATCCTGCATGGCATCGAAGTGCCCAGCCAGATCCGTCACGACAACACCCTGAGCAAACCGCTGATCAGCTGGGGCCCAAGCGAGCGCGTGCATTGTATCGTCGCCAACCCGCCGTTCGGCGGCATGGAAGAAGACGGTATCGaaa                                                                                                                         <  1:1491979/150‑1 (MQ=255)
                                   tCCTGCATGGCATCGAAGTGCCCAGCCAGATCCGTCACGACAACACCCTGAGCAAACCGCTGATCAGCTGGGGCCCAAGCGAGCGCGTGCATTGTATCGTCGCCAACCCGCCGTTCGGCGGCATGGAAGAAGACGGTATCGAAACAAAtt                                                                                                                   >  2:395206/1‑150 (MQ=255)
                                      tGCATGGCATCGAAGTGCCCAGCCAGATCCGTCACGACAACACCCTGAGCAAACCGCTGATCAGCTGGGGCCCAAGCGAGCGCGTGCATTGTATCGTCGCCAACCCGCCGTTCGGCGGCATGGAAGAAGACGGTATCGAAACAAATTTTc                                                                                                                >  1:1168688/1‑150 (MQ=255)
                                      tGCATGGCATCGAAGTGCCCAGCCAGATCCGTCACGACAACACCCTGAGCAAACCGCTGATCAGCTGGGGCCCAAGCGAGCGCGTGCATTGTATCGTCGCCAACCCGCCGTTCGGCGGCATGGAAGAAGACGGTATCGAAACAAATTTTc                                                                                                                >  1:1171778/1‑150 (MQ=255)
                                            gCATCGAAGTGCCCAGCCAGATCCGTCACGACAACACCCTGAGCAAACCGCTGATCAGCTGGGGCCCAAGCGAGCGCGTGCATTGTATCGTCGCCAACCCGCCGTTCGGCGGCATGGAAGAAGACGGTATCGAAACAAATTTTCCTGCCg                                                                                                          >  2:936282/1‑150 (MQ=255)
                                             cATCGAAGTGCCCAGCCAGATCCGTCACGACAACACCCTGAGCAAACCGCTGATCAGCTGGGGCCCAAGCGAGCGCGTGCATTGTATCGTCGCCAACCCGCCGTTc                                                                                                                                                     >  2:213983/1‑106 (MQ=255)
                                             cATCGAAGTGCCCAGCCAGATCCGTCACGACAACACCCTGAGCAAACCGCTGATCAGCTGGGGCCCAAGCGAGCGCGTGCATTGTATCGTCGCCAACCCGCCGTTc                                                                                                                                                     <  1:213983/106‑1 (MQ=255)
                                              aTCGAAGTGCCCAGCCAGATCCGTCACGACAACACCCTGAGCAAACCGCTGATCAGCTGGGGCCCAAGCGAGCGCGTGCATTGTATCGTCGCCAACCCGCCGTTCGGCGGCATGGAAGAAGACGGTATCGAAACAAATTTTCCTGCCGCt                                                                                                        >  1:690834/1‑150 (MQ=255)
                                                         cagccagATCCGTCACGACAACACCCTGAGCAAACCGCTGATCAGCTGGGGCCCAAGCGAGCGCGTGCATTGTATCGTCGCCAACCCGCCGTTcgg                                                                                                                                                   >  2:1872731/1‑96 (MQ=255)
                                                         cagccagATCCGTCACGACAACACCCTGAGCAAACCGCTGATCAGCTGGGGCCCAAGCGAGCGCGTGCATTGTATCGTCGCCAACCCGCCGTTcgg                                                                                                                                                   <  1:1872731/96‑1 (MQ=255)
                                                            ccagATCCGTCACGACAACACCCTGAGCAAACCGCTGCTTAGCTGGGGCCCAAGCGAGTGCGTGCATTGTATCGTCGCCAACCCGCCGTTCGGCGGCATGGAAGAAGACGGTATCGAAACAAATTTTCCTGCCGCTTTCCGCACCCGGGa                                                                                          <  1:958380/150‑1 (MQ=255)
                                                            ccagATCCGTCACGACAACACCCTGAGCAAACCGCTGATCAGCTGGGGCCCAAGCGAGAGCGTGCATTGTATCGTCGCCAACCCGCCGTTCGGCGGCATGGAAGAAGACGGTATCGAAACAAATTTTCCTGCCGCTTTCCGCACCCGGGa                                                                                          >  1:1854829/1‑150 (MQ=255)
                                                                    gTCACGACAACACCCTGAGCAAACCGCTGATCAGCTGGGGCCCAAGCGAGCGCGTGCATTGTATCGTCGCCAACCCGCCGTTCGGCGGCATGGAAGAAGACGGTATCGAAACAAATTTTCCTGCCGCTTTCCGCACCCGGGAAAccgccg                                                                                  >  2:660699/1‑150 (MQ=255)
                                                                    gTCACGACAACACCCTGAGCAAACCGCTGATCAGCTGGGGCCCAAGCGAGCGCGTGCATTGTATCGTCGCCAACCCGCCGTTCGGCGGCACAGAACAGGACGGTATCGAAACAAATTTTCCTGCCGCTTTCCGCACCCGGGAAAccgccg                                                                                  >  2:52349/1‑150 (MQ=255)
                                                                     tCACGACAACACCCTGAGCAAACCGCTGATCAGCTGGGGCCCAAGCGAGCGCGTGCATTGTATCGTCGCCAACCCGCCGTTCGGCGGCATGGAAGAAGACGGTATCGAAACAAATTTTCCTGCCGCTTTCCGCACCCGGGAAACCGCCGa                                                                                 <  1:660699/150‑1 (MQ=255)
                                                                       aCGACAACACCCTGAGCAAACCGCTGATCAGCTGGGGCCCAAGCGAGCGCGTGCATTGTATCGTCGCCAACCCGCCGTTCGGCGGCATGGAAGAAGACGGTATCGAAACAAATTTTCCTGCCGCTTTCCGCACCCGGGAAACCGCCGAtt                                                                               <  1:236158/150‑1 (MQ=255)
                                                                       aCGACAACACCCTGAGCAAACCGCTGATCAGCTGGGGCCCAAGCGAGCGCGTGCATTGTATCGTCGCCAACCCGCCGTTCGGCGGCATGGAAGAAGACGGTATCGAAACAAATTTTCCTGCCGCTTTCCGCACCCGGGAAACCGCCGAtt                                                                               <  2:649159/150‑1 (MQ=255)
                                                                          acaacaCCCTGAGCAAACCGCTGATCAGCTGGGGCCCAAGCGAGCGCGTGCATTGTATCGTCGCCAACCCCCCGTTCGGCGGCATGGACGAAGACGGTATCGAAACAAATTTTCCTGCCGCTTTCCGCACCCGGGAAACCGCCGATTTGt                                                                            >  1:1077004/1‑150 (MQ=255)
                                                                              cacCCTGAGCAAACCGCTGATCAGCTGGGGCCCAAGCGAGCGCGTGCATTGTATCGTCGCCAACCCGCCGTTCGGCGGCATAGAAGAAGACGGTATCGAAACAAATTTTCCTGCCGCTTTCCGCACCCGGGAAACCGCCGATTTGTTCtt                                                                        <  1:63636/150‑1 (MQ=255)
                                                                               acCCTGAGCAAACCGCTGATCAGCTGGGGCCCAAGCGAGCGCGTGCATTGTATCGTCGCCAACCCGCCGTTCGGCGGCATGGAAGAAGACGGTATCGAAACAAATTTTCCTGCCGCTTTCCGCACCCGGGAAACCGCCGATTTGTTCTTg                                                                       >  1:262393/1‑150 (MQ=255)
                                                                               acCCTGAGCAAACCGCTGATCAGCTGGGGCCCAAGCGAGCGCGTGCATTGTATCGTCGCCAACCCGCCGTTCGGCGGCATGGAAGAAGACGGTATCGAAACAAATTTTCCTGCCGCTTTCCGCACCCGGGAAACCGCCGATTTGTTCTTg                                                                       <  1:395206/150‑1 (MQ=255)
                                                                                 ccTGAGCAAACCGCTGATCAGCTGGGGCCCAAGCGAGCGCGTGCATTGTATCGTCGCCAACCCGCCGTTCGGCGGCATGGAAGAAGACGGTATCGAAACaaa                                                                                                                     <  2:1206798/102‑1 (MQ=255)
                                                                                 ccTGAGCAAACCGCTGATCAGCTGGGGCCCAAGCGAGCGCGTGCATTGTATCGTCGCCAACCCGCCGTTCGGCGGCATGGAAGAAGACGGTATCGAAACaaa                                                                                                                     >  1:1206798/1‑102 (MQ=255)
                                                                                 ccTGAGCAAACCGCTGATCAGCTGGGGCCCAAGCGAGCGCGTGCATTGTATCGTCGCCAACCCGCCGTTCGGCGGCATGGAAGAAGACGGTATCGAAACAAATTTTCCTGCCGCTTTCCGCACCCGGGAAACCGCCGATTTGTTCTTGGt                                                                     <  2:1224734/150‑1 (MQ=255)
                                                                                  gtGAGCAAACCGTTGAGCAGCTGGGGCCCAAGCGGGCGCGTGCATTGTATCGTCGCCAACCCGCCGTTCGGCGACATGGAAGAAGACGGTATCGTAAGTATTATTCCTCGCGCTTTCCGCACCCGGGACATCGCCTATTTGTTCTAGGTa                                                                    <  2:126826/149‑1 (MQ=255)
                                                                                  cTGAGCAAACGGCTGATCAGCTGGGGCCCAAGCGAGCGCGTGCATTGTATCGTCGCCAACCCGCCGTTCGGCGGCATGGAAGAAGACGGTATCGAAACAAATTTTCCTGCCGCTTTCCGCACCCGGGAAACCGCCGATTTGTTCTTGGTa                                                                    <  1:1854315/150‑1 (MQ=255)
                                                                                  cTGAGCAAACCGCTGATCAGCTGGGGCCCAAGCGAGCGCGTGCATTGTATCGTCGCCAACCCGCCGTTCGGCGGCATGGAAGAAGACGGTATCGAAACAAATTTTCCTGCCGCTTTCCGCACCCGGGAAACCGCCGATTTGTTCTTGGTa                                                                    <  2:1311664/150‑1 (MQ=255)
                                                                                  cTGAGCAAACCGCTGATCAGCTGGGGCCCAAGCGAGCGCGTGCATTGTATCGTCGCCAACCCGCCGTTCGGCGGCATGGAAGAAGACGGTATCGAAACAAATTTTCCTGCCGCTTTCCGCACCCGGGAAACCGCCGATTTGTTCTTGGTa                                                                    <  2:1312017/150‑1 (MQ=255)
                                                                                  cTGAGCAAACCGCTGATCAGCTGGGGCCCAAGCGAGCGCGTGCATTGTATCGTCGCCAACCCGCCGTTCGGCGGCATGGAAGAAGACGGTATCGAAACAAATTTTCCTGCCGCTTTCCGCACCCGGGAAACCGCCGATTTGTTCTTGGTa                                                                    >  1:428882/1‑150 (MQ=255)
                                                                                   tGAGCAAACCGCTGATCAGCTGGGGCCCAAGCGAGCGCGTGCATTGTATCGTCGCCAACCCGCCGTTCGGCGGCATGGAAGAAGACGGTATCGAAACAAATTTTCCTGCCGCTTTCCGCACCCGGGAAACCGCCGATTTGTTCTTGGTAt                                                                   <  2:1171778/150‑1 (MQ=255)
                                                                                      gaaatccGCTGATCAGCTGGGGCCCAAGCGAGCGCGTGCATTGTATCGTCGCCAACCCGCCGTTCGGCGGCATGGAAGAAGACGGTATCGAAACAAATTTTCCTGCCGCTTTCCGCACCCGGGAAACCGCCGATTTGTTCTTGGTATTGa                                                                <  1:928309/145‑1 (MQ=255)
                                                                                      gCAAACCGCTGATCAGCTGGGGCCCAAGCGAGCGCGTGCATTGTATCGTCGCCAACCCGCCGTTCGGCGGCATGGAAGAAGACGGTATCGAAACAAATTTTCCTGCCGCTTTCCGCACCCGGGAAACCGCCGATTTGTTCTTGGTATTGa                                                                <  2:1763668/150‑1 (MQ=255)
                                                                                      ccAAACCGCTGATCAGCTGGGGCCCAAGCGAGCGCGTGCATTGTATCGTCGCCAACCCGCCGTTCGGCGGCATGGAAGAAGACGGTATCGAAACAAATTTTCCTGCCGCTTTCCGCACCCGGGAAACCGCCGATTTGTTCTTGGTATTGa                                                                <  1:928313/149‑1 (MQ=255)
                                                                                       cAAACCGCTGATCAGCTGGGGCCCAAGCGAGCGCGTGCATTGTATCGTCGCCAACCCGCCGTTCGGCGGCATGGAAGAAGACGGTATCGAAACAAATTTTCCTGCCGCTTTCCGCACCCGGGAAACCGCCGATTTGTTCTTGGTATTGAt                                                               <  1:1741217/150‑1 (MQ=255)
                                                                                            cGCTGATCAGCTGGGGCCCAAGCGAGCGCGTGCATTGTATCGTCGCCAACCCGCCGTTCGGCGGCATGGAAGAAGACTGTATCGAAAGAAATTTTCCTGCCTCTTTCCGCACCCGGTAAACCGCCGATGTTGGCTTGGTATTGATTATGc                                                          >  2:1841879/1‑150 (MQ=255)
                                                                                             gCTGATCAGCTGGGGCCCAAGCGAGCGCGTGCATTGTATCGTCGCCAACCCGCCGTTCGGCGGCATGGAAGAAGACGGTATCGAAACAAATTTTCCTGCCGCTTTCCGCACCCGGGAAACCGCCGATTTGTTCTTGGTATTGATTATGCa                                                         <  2:1601491/150‑1 (MQ=255)
                                                                                             gCTGATCAGCTGGGGCCCAAGCGAGCGCGTGCATTGTATCGTCGCCAACCCGCCGTTCGGCGGCATGGAAGAAGACGGTATCGAAACAAATTTTCCTGCCGCTTTCCGCACCCGGGAAACCGCCGATTTGTTCTTGGTATTGATTATGCa                                                         <  2:690834/150‑1 (MQ=255)
                                                                                               tGATCAGCTGGGGCCCAAGCGAGCGCGTGCATTGTATCGTCGCCAACCCGCCGTTCGGCGGCATGGAAGAAGACGGTATCGAAACAAATTTTCCTGCCGCTTTCCGCACCCGGGAAACCGCCGATTTGTTCTTGGTATTGATTATGCAgc                                                       >  1:785716/1‑150 (MQ=255)
                                                                                                       tGGGGCCCAAGCGAGCGCGTGCATTGTATCGTCGCCAACCCGCCGTTCGGCGGCATGGAAGAAGACGGTATCGAAACAAATTTTCCTGCCGCTTTCCGCACCCGGGAAACCGCCGATTTGTTCTTGGTATTGATTATGCAGCTGCTCaaa                                               >  1:195697/1‑150 (MQ=255)
                                                                                                           gCCCAAGCGAGCGCGTGCCTTTTATCGTCGCCAACCCGCCGTTCGGCGGCATGGAAGAAGACGGTATCGAAACAAATTTTCCTGCCGCTTTCCGCACCCGGGAAACCGCCGATTTGTTCTTGGTATTGATTATGCAGCTGCTCAAAGAtg                                           >  2:465603/1‑150 (MQ=255)
                                                                                                            cccAAGCGAGCGCGTGCATTGTATCGTCGCCAACCCGCCGTTCGGCGGCATGGAAGAAGACGGTATCGAAACAAATTTTCCTGCCGCTTTCCGCACCCGGGAAACCGCCGATTTGTTCTTGGTATTGATTATGCAGCTGCTCAAAGAtgg                                          <  1:445733/150‑1 (MQ=255)
                                                                                                             ccAAGCGAGCGCGTGCATTGTATCGTCGCCAACCCGCCGTTCGGCGGCATGGAAGAAGACGGTATCGAAACAAATTTTCCTGCCGCTTTCCGCACCCGGGAAACCGCCGATTTGTTCTTGGTATTGATTATGCAGCTGCTCAAAGAtggt                                         <  2:195697/150‑1 (MQ=255)
                                                                                                                  cgagagCGTGCATTGTATCGTCGCCAACCCGCCGTTCGGCGGCATGGAAGAAGACGGTATCGAAACAAATTTTCCTGCCGCTTTCCGCACCCGGGAAACCGCCGATTTGTTCTTGGTATTGATTATGCAGCTGCTCAAAGATGGTGgccg                                    <  2:1854829/150‑1 (MQ=255)
                                                                                                                     gcgcGTGCATTGTATCGTCGCCAACCCGCCGTTCGGCGGCATGGAAGAAGACGGTATCGAAACAAATTTTCCTGCCGCTTTCCGCACCCGGGAAACCGCCGATTTGTTCTTGGTATTGATTATGCAGCTGCTCAAAGATGGTGgccgcgc                                 <  1:1706488/150‑1 (MQ=255)
                                                                                                                          ggCATTGCATCGCCGCCAACCCGCCGTTCGGCGGCATGGAAGAAGACGGTATCCAACCAAATTATCCTGCCGCTCTCCGCCCCCGGGAAACCGCCGATTTGTTCTTCGTACTGATTATGCAGCTGCGCAAAGATGGTGGGCGCGCCGCAg                            <  2:1923058/149‑1 (MQ=255)
                                                                                                                                   tcgtcgCCAACCCGCCGTTCGGCGGCATGGAAGAAGACGGTATCGAAACAAATTTTCCTGCCGCTTTCCGCACCCGGGAAACCGCCGATTTGTTCTTGGTATTGATTATGCAGCTGCTCAAAGATGGTGGCCGCGCCGCAGTGGTACTGc                   >  1:1192739/1‑150 (MQ=255)
                                                                                                                                       cgCCAACCCGCCGTTCGGCGGCATGGAAGAAGACGGTATCGAAACAAATTTTCCTGCCGCTTTCCGCACCCGGGAAACCGCCGATTTGTTCTTGGTATTGATTATGCAGCTGCTCAAAGATGGTGGCCGCGCCGCAGTGGTACTGCCCGa               <  2:1815034/150‑1 (MQ=255)
                                                                                                                                       cgCCAACCCGCCGTTCGGCGGCATGGAAGAAGACGGTATCGAAACAAATTTTCCTGCCGCTTTCCGCACCCGGGAAACCGCCGATTTGTTCTTGGTATTGATTATGCAGCTGCTCAAAGATGGTGGCCGCGCCGCAGTGGTACGGcgctg               <  2:1168688/150‑5 (MQ=255)
                                                                                                                                        gCCAACCCGCCGTTCGGCGGCATGGAAGAAGACGGTATCGAAACAAATTTTCCTGCCGCTTTCCGCACCCGGGAAACCGCCGATTTGTTCTTGGTATTGATTATGCAGCTGCTCAAAGATGGTGGCCGCGCCGCAGTGGTACTGCCCGAt              <  2:428882/150‑1 (MQ=255)
                                                                                                                                        gCCAACCCGCCGTTCGGCGGCATGGAAGAAGACGGTATCGAAACAAATTTTCCTGCCGCTTTCCGCACCCGGGAAACCGCCGATTTGTTCTTGGTATTGATTATGCAGCTGCTCAAAGATGGTGGCCGCGCCGCAGTGGTACTGCCCGAt              >  1:1488238/1‑150 (MQ=255)
                                                                                                                                         ccAACCCGCCGTTCGGCGGCATGGAAGAAGACGGTATCGAAACAAATTTTCCTGCCGCTTTCCGCACCCGGGAAACCGCCGATTTGTTCTTGGTATTGATTATGCAGCTGCTCAAAGATGGTGGCCGCGCCGCAGTGGTACTGCCCGATg             <  1:1025849/150‑1 (MQ=255)
                                                                                                                                                gccgTTCGGCGGCATGGAAGAAGACGGTATCGAAACAAATTTGCCTGCCGCTTTCCGCACCCGGCAAACCGCCGATTTGTTCTTGGTATTGATTATGCAGCTGCTCAAAGATGGTGGCCGCGCCGCAGTGGTACTGCCCGATGGCTTCCt      <  1:782035/150‑1 (MQ=255)
                                                                                                                                                 ccgTTCGGCGGCATGGAAGAAGACGGTATCGAAACAAATTTTCCTGCCGCTTTCCGCACCCGGGAAACCGCCGATTTGTTCTTGGTATTGATTATGCAGCTGCTCAAAGAtg                                           >  1:725021/1‑112 (MQ=255)
                                                                                                                                                 ccgTTCGGCGGCATGGAAGAAGACGGTATCGAAACAAATTTTCCTGCCGCTTTCCGCACCCGGGAAACCGCCGATTTGTTCTTGGTATTGATTATGCAGCTGCTCAAAGAtg                                           <  2:725021/112‑1 (MQ=255)
                                                                                                                                                    ttCGGCGGCATGGAAGAAGACGGTATCGAAACAAATTTTCCTGCCGCTTTCCGCACCCGGGAAACCGCCGATTTGTTCTTGGTATTGATTATGCAGCTGCTCAAAGATGGTGGCCGCGCCGCAGTGGTACTGCCCGATGGCTTCCttttt  <  2:1488238/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
ACCGCTGCCGCACCTGCTGGCCACCACCAACATGATCCTGCATGGCATCGAAGTGCCCAGCCAGATCCGTCACGACAACACCCTGAGCAAACCGCTGATCAGCTGGGGCCCAAGCGAGCGCGTGCATTGTATCGTCGCCAACCCGCCGTTCGGCGGCATGGAAGAAGACGGTATCGAAACAAATTTTCCTGCCGCTTTCCGCACCCGGGAAACCGCCGATTTGTTCTTGGTATTGATTATGCAGCTGCTCAAAGATGGTGGCCGCGCCGCAGTGGTACTGCCCGATGGCTTCCTTTTT  >  NC_002947/5396254‑5396551

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: