Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,033,806 T→A 6.1% D141E (GAT→GAA PP_0891 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,033,8060TA 6.1% 323.6 / 9.1 99D141E (GAT→GAAPP_0891hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (48/45);  new base A (4/2);  total (52/47)
Fisher's exact test for biased strand distribution p-value = 6.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.20e-01

CCACAGTGAAGGCACAGATCGACCTGCAGCCGATCAATGACCTGGCCGACGGTGCACAGATCGAATTGCGCCTGCCGCTGACGGTCACCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGATGCGCGGCGCTTCCAGGTGGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCAAGTTCGCCAGGTTGAAGTCCATCAACCCGTCGGTACCGGTTAGCGCGGTGCTGATCTTCAC  >  NC_002947/1033664‑1033955
                                                                                                                                              |                                                                                                                                                     
ccACAGTGAAGGCACAGATCGACCTGCAGCCGATCAATGACCTGGCCGACGGTGCACAGATCGAATTGCGCCTGCCGCTGACGGTCACCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGATGCgcggc                                                                                                                                                <  2:1136236/150‑1 (MQ=255)
 cacaGTGAAGGCACAGATCGACCTGCAGCCGATCAATGACCTGGCCGACGGTGCACAGATCGAATTGCGCGTGCCGCTGACGGTCACCCTGCATGGCCAGTCGCACAGTTACAGTGTCCTGTTGCTGGCCACCCGCCAGGATGCgcggcg                                                                                                                                               >  2:541470/1‑150 (MQ=255)
 cacaGTGAAGGCACAGATCGACCTGCAGCCGATCAATGACCTGGCCGACGGTGCACAGATCGAATTGCGCCTGCCGCTGACGGTCACCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGATGCgcggcg                                                                                                                                               >  1:1395345/1‑150 (MQ=255)
   caGTGAAGGCACAGATCGACCTGCAGCCGATCAATGACCTGGCCGACGGTGCACAGATCGAATTGCGCCTGCCGATGAGGGTCACCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGATGCGCGGCGCt                                                                                                                                             >  2:805139/1‑150 (MQ=255)
          ggCACAGATCGACCTGCAGCCGATCAATGACCTGGCCGACGGTGCACAGATCGAATTGCGCCTGCCGCTGACGGTCACCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGCTGCGTTCCACCCGCCTGGATGCGCGGCGCTTCCAggt                                                                                                                                      <  1:2371895/150‑1 (MQ=255)
           gCACAGATCGACCTGCAGCCGATCAATGACCTGGCCGACGGTGCACAGATCGAATTGCGCCTGCCGCTGACGGTCACCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGATGCGCGGCGCTTCCAggtg                                                                                                                                     <  2:1395345/150‑1 (MQ=255)
            cacaGATCGACCTGCAGCCGATCAATGACCTGGCCGACGGTGCACAGATCGAATTGCGCCTGCCGCTGACGGTCACCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGATGCGCGGCGCTTCCAggtgg                                                                                                                                    <  1:1844985/150‑1 (MQ=255)
            cacaGATCGACCTGCAGCCGATCAATGACCTGGCCGACGGTGCACAGATCGAATTGCGCCTGCCGCTGACGGTCACCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGATGCGCGGCGCTTCCAggt                                                                                                                                      >  2:2172078/1‑148 (MQ=255)
            cacaGATCGACCTGCAGCCGATCAATGACCTGGCCGACGGTGCACAGATCGAATTGCGCCTGCCGCTGACGGTCACCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGATGCGCGGCGCTTCCAggt                                                                                                                                      <  1:2172078/148‑1 (MQ=255)
            cacaGATCGACCTGCAGCCGATCAATGACCTGGCCGACGGTGCACAGATCGAATTGCGCCTGCCGCTGACGGTCACCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGATGCGCGGCGCTTCCAggcgg                                                                                                                                    <  2:642570/150‑1 (MQ=255)
                gATCGACCTGCAGCCGATCAATGACCTGGCCGACGGTGCACAGATCGAATTGCGCCTGCCGCTGACGGTCACCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGATGCGCGGCGCTTCCAGGTGGTGAc                                                                                                                                >  2:2121351/1‑150 (MQ=255)
                 aTCGACCTGCAGCCGATCAATGACCTGGCCGACGGTGCACAGATCGAATTGCGCCTGCCGCTGACGGTCACCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGATGCGCGGCGCTTCCAGGTGGTGAcc                                                                                                                               >  1:1782429/1‑150 (MQ=255)
                  tCGACCTGCAGCCGATCAATGACCTGGCCGACGGTGCACAGATCGAATTGCGCCTGCCGCTGACGGTCACCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGATGCGCGGCGCTTCCAGGTGGTGAccc                                                                                                                              <  2:1277345/150‑1 (MQ=255)
                    gACCTGCAGCCGATCAATGACCTGGCCGACGGTGCACAGATCGAATTGCGCCTGCCGCTGACGGTCACCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCCCCTGGATGCGCTGCGCTTCCAGGTGGTGACCCtt                                                                                                                            <  2:1234678/150‑1 (MQ=255)
                     aCCTGCAGCCGATCAATGACCTGGCCGACGGTGCACAGATCGAATTGCGCCTGCCGCTGACGGTCACCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGATGCGCGGCGCTTCCAGGTGGTGACCCTTg                                                                                                                           >  1:713048/1‑150 (MQ=255)
                        tGCAGCCGATCAATGACCTGGCCGACGGTGCACAGATCGAATTGCGCCTGCCGCTGACGGTCACCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGATGCGCGGCGCTTCCAGGTGGTGACCCTTGAGc                                                                                                                        >  2:1784409/1‑150 (MQ=255)
                        tGCAGCCGATCAATGACCTGGCCGACGGTGCACAGATCGAATTGCGCCTGCCGCTGACGGTCACCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGATGCGCGGCGCTTCCAGGTGGTGACCCTTGAGc                                                                                                                        >  2:1089210/1‑150 (MQ=255)
                              cGATCAATGACCTGGCCGACGGTGCACAGATCGAATTGCGCCTGCCGCTGACGGTCACCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGATGCGCGGCGCTTCCAGGTGGTGACCCTTGAGCCTCTgg                                                                                                                  <  2:1356225/150‑1 (MQ=255)
                               gATCAATGACCTGGCCGACGGTGCACAGATCGAATTGCGCCTGCCGCTGACGGTCACCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGATGCGCGGCGCTTCCAGGTGGTGACCCTTGAGCCTCTGGt                                                                                                                 >  1:987213/1‑150 (MQ=255)
                               gATCAATGACCTGGCCGACGGTGCACAGATCGAATTGCGCCTGCCGCTGACGGTCACCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGATGCGCGGCGCTTCCAGGTGGTGACCCTTGAGCCTCTGGt                                                                                                                 >  1:1288312/1‑150 (MQ=255)
                                aTCAATGACCTGGCCGACGGTGCACAGATCGAATTGCGCCTGCCGCTGACGGTCACCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGATGCGCGGCGCTTCCAGGTGGTGACCCTTGAGCCTCTGGtg                                                                                                                >  1:1805055/1‑150 (MQ=255)
                                aTCAATGACCTGGCCGACGGTGCAACGCTCGAGTTGCGGCTGGCGCTCACGTTCGCCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCCCCCCCCTGGCTGCGCGGTGCTTTCTGGTGCTGAGTCGGGAGCctcttttg                                                                                                                <  2:529742/150‑5 (MQ=255)
                                   aaTGACCTGGCCGACGGTGCACAGATCGAATTGCGCCTGCCGCTGACGGTCACCCTGCATGGCCAGTCGCACAGCTAAAATGTCCTGTTGCTGGCCACCCGCCTGGATGCGCGGCGCTTCCAGGTGGTGACCCTTGAGCCTCTGGtgctg                                                                                                             >  2:1399694/1‑150 (MQ=255)
                                      gACCTGGCCGACGGTGCACAGATCGAATTGCGCCTGCCGCTGACGGTCACCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGATGCGCGGCGCTTCCAggtg                                                                                                                                     <  1:21732/123‑1 (MQ=255)
                                      gACCTGGCCGACGGTGCACAGATCGAATTGCGCCTGCCGCTGACGGTCACCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGATGCGCGGCGCTTCCAggtg                                                                                                                                     >  2:21732/1‑123 (MQ=255)
                                            gCCGACGGTGCACAGATCGAATTGCGCCTGCCGCTGACGGTCACCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGATGCGCGGCGCTTCCAGGTGGTGACCCTTGAGCCTCTgg                                                                                                                  >  1:1861630/1‑136 (MQ=255)
                                            gCCGACGGTGCACAGATCGAATTGCGCCTGCCGCTGACGGTCACCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGATGCGCGGCGCTTCCAGGTGGTGACCCTTGAGCCTCTgg                                                                                                                  <  2:1861630/136‑1 (MQ=255)
                                               gACGGTGCACAGATCGAATTGCGCCTGCCGCTGACGGTCACCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGATGCGCGGCGCTTCCAGGTGGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGAt                                                                                                 >  1:1615178/1‑150 (MQ=255)
                                                    tGCACAGATCTAATTGCGCCTGCCGCTGACGGTCACCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGATGCGCGGCGCTTCCAGGTGGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTgg                                                                                            <  1:715346/150‑1 (MQ=255)
                                                    tGCACAGATCGAATTGCGCCTGCCGCTGACGGTCACCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGATGCGCGGCGCTTCCAGGTGGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTgg                                                                                            <  1:905426/150‑1 (MQ=255)
                                                       acaGATCGAATTGCGCCTGCCGCTGACGGTCACCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGATGCGCGGCGCTTCCAGGTGGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCtt                                                                                         >  1:277196/1‑150 (MQ=255)
                                                       acaGATCGAATTGCGCCTGCCGCTGACGGTCACCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGATGCGCGGCGCTTCCAGGTGGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCtt                                                                                         >  2:289757/1‑150 (MQ=255)
                                                       acaGATCGAATTGCGCCTGCCGCTGACGGTCACCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGATGCGCGGCGCTTCCAGGTGGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCtt                                                                                         >  2:2342485/1‑150 (MQ=255)
                                                       acaGATCGAATTGCGCCTGCCGCTGACGGTCACCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGAAGCGCCGCGCATCCAggc                                                                                                                                      <  1:1051258/105‑2 (MQ=255)
                                                       acaGATCGAATTGCGCCTGCCGCTGACGGTCACCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGAAGCGCCGCGCATCCAggc                                                                                                                                      >  2:1051258/1‑104 (MQ=255)
                                                       acaGATCGAATTGCGCCTGCCGCTGACGGTCACACTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGATGCGCGGCGCTTCCAGGTGGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCtt                                                                                         <  2:570366/150‑1 (MQ=255)
                                                           aTCGAATTGCGCCTGCCGCTGACGGTCACCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGATGCGCGGCGCTTCCAGGTGGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTtgctg                                                                                     >  2:268841/1‑150 (MQ=255)
                                                           aTCGAATTGCGCCTGCCGCTGACGGTCACCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGATGCGCGGCGCTTCCAGGTGGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTtgctg                                                                                     >  2:222218/1‑150 (MQ=255)
                                                            tCGAATTGCGCCTGCCGCTGACGGTCACCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGATGCGCGGCGCTTCCAGGTGGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTg                                                                                             >  2:754930/1‑141 (MQ=255)
                                                            tCGAATTGCGCCTGCCGCTGACGGTCACCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGATGCGCGGCGCTTCCAGGTGGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTg                                                                                             <  1:754930/141‑1 (MQ=255)
                                                            tCGAATTGCGCCTGCCGCTGACGGTCACCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGATGCGCGGCGCTTCCAGGTGGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTtgctgc                                                                                    <  1:268841/150‑1 (MQ=255)
                                                               aaTTGCGCCTGCCGCTGACGGTCACCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGATGCGCGGCGCTTCCAGGTGGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCgg                                                                                 <  2:1261524/150‑1 (MQ=255)
                                                               aaTTGCGCCTGCCGCTGACGGTCACCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGATGCGCGGCGCTTCCAGGTGGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCgg                                                                                 <  2:987213/150‑1 (MQ=255)
                                                                aTTGCGCCTGCCGCTGACGGTCACCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGATGCGCGGCGCTTCCAGGTGGTGACCCTTGAGCCTTTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCggg                                                                                <  1:1949609/150‑1 (MQ=255)
                                                                    cgcCTGCCGCTGACGGTCACCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGATGCGCGGCGCTTCCAGGTGGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTg                                                                            <  1:222218/150‑1 (MQ=255)
                                                                    cgcCTGCCGCTGACGGTCACCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGATGCGCGGCGCTTCCAGGTGGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTg                                                                            >  2:1456765/1‑150 (MQ=255)
                                                                     gcCTGCCGCTGACGGTCACCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGATGCGCGGCGCTTCCAGGTGGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTgg                                                                                            <  1:361550/133‑1 (MQ=255)
                                                                     gcCTGCCGCTGACGGTCACCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGATGCGCGGCGCTTCCAGGTGGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTgg                                                                                            >  2:361550/1‑133 (MQ=255)
                                                                     gcCTGCCGCTGACGGTCACCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGATGCGCGGCGCTTCCAGGTGGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTgg                                                                           >  1:2181929/1‑150 (MQ=255)
                                                                     gcCTGCCGCTGACGGTCACCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGATGCGCGGCGCTTCCAGGTGGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTgg                                                                           >  1:2036692/1‑150 (MQ=255)
                                                                      ccTGCCGCTCACGGTCACCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGATGCGCGGCGCTTCCAGGTGGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGa                                                                          <  1:1309133/150‑1 (MQ=255)
                                                                       cTGCCGCTGACGGTCACCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGATGCGCGGCGCTTCCAGGTGGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGaa                                                                         >  1:1722395/1‑150 (MQ=255)
                                                                       cTGCCGCTGACGGTCACCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGATGCGCGGCGCTTCCAGGTGGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGaa                                                                         >  2:1962722/1‑150 (MQ=255)
                                                                              tGACGGTCACCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGATGCGCGGCGCTTCCAGGTGGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTgc                                                                  <  1:1089210/150‑1 (MQ=255)
                                                                                  ggTCACCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGATGCGCGGCGCTTCCAGGTGGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAATCCTGCGCaa                                                              <  2:713048/150‑1 (MQ=255)
                                                                                      aCCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGATGCGCGGCGCTTCCAGGTGGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCAAGTTc                                                          <  2:2181929/150‑1 (MQ=255)
                                                                                        ccTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGATGCGCGGCGCTTCCAGGTGGTGAccc                                                                                                                              >  1:1052253/1‑80 (MQ=255)
                                                                                        ccTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGATGCGCGGCGCTTCCAGGTGGTGAccc                                                                                                                              <  2:1052253/80‑1 (MQ=255)
                                                                                          tGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGATGCGCGGCGCTTCCAGGTGGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCAAGTTCGCCa                                                      <  2:374395/150‑1 (MQ=255)
                                                                                             aTGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGATGCGCGGCGCTTCCAGGTGGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCAAGTTCGCCAGGt                                                   <  1:2121351/150‑1 (MQ=255)
                                                                                              tGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGATGCGCGGCGCTTCCAggt                                                                                                                                      >  1:1806969/1‑66 (MQ=255)
                                                                                              tGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGATGCGCGGCGCTTCCAggt                                                                                                                                      <  2:1806969/66‑1 (MQ=255)
                                                                                                 ccAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGATGCGCGGCGCTTCCAGGTGGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGATGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCAAGTTCGCCAGGTTGaa                                               >  1:1144643/1‑150 (MQ=255)
                                                                                                 ccAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGATGCGCGGCGCTTCCAGGTGGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGaaa                                                                        >  2:143340/1‑125 (MQ=255)
                                                                                                 ccAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGATGCGCGGCGCTTCCAGGTGGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGaaa                                                                        <  1:143340/125‑1 (MQ=255)
                                                                                                   aGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGATGCGCGGCGCTTCCAGGTGGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCAAGTTCGCCAGGTTGAAGt                                             <  2:1805055/150‑1 (MQ=255)
                                                                                                       gCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGATGCGCGGCGCTTCCAGGTGGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCAAGTTCGCCAGGTTGAAGTCCAt                                         >  2:2289202/1‑150 (MQ=255)
                                                                                                        cacaGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGATGCGCGGCGCTTCCAGGTGGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCAAGTTCGCCAGGTTGAAGTCCATc                                        <  2:277196/150‑1 (MQ=255)
                                                                                                        cacaGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGATGCGCGGCGCTTCCAGGTGGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCAAGTTCGCCAGGTTGAAGTCCATc                                        <  1:2342485/150‑1 (MQ=255)
                                                                                                              tACAATGTCCTGTTGCTGGCCACCCGCCTGGATGCGCGGCGCTTCCAGGTGGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCAAGTTCGCCAGGTTGAAGTCCATCAACCCg                                  >  1:1003028/1‑150 (MQ=255)
                                                                                                                cAATGTCCTGTTGCTGGCCACCCGCCTGGATGCGCGGCGCTTCCAGGTGGTGACCCTTGa                                                                                                                          >  1:44319/1‑60 (MQ=255)
                                                                                                                cAATGTCCTGTTGCTGGCCACCCGCCTGGATGCGCGGCGCTTCCAGGTGGTGACCCTTGa                                                                                                                          <  2:44319/60‑1 (MQ=255)
                                                                                                                       cTGTTGCTGGCCACCCGCCTGGATGCGCGGCGCTTCCAGGTGGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCAAGTTCGCCAGGTTGAAGTCCATCAACCCGTCGGTACCg                         <  1:944350/150‑1 (MQ=255)
                                                                                                                       cTGTTGCTGGCCACCCGCCTGGATGCGCGGCGCTTCCAGGTGGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCAAGTTCGCCAGGTTGAAGTCCATCAACCCGTCGGTACCg                         >  1:1267308/1‑150 (MQ=255)
                                                                                                                        tGTTGCTGGCCACCCGCCTGGATGCGCGGCGCTTCCAGGTGGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTgg                                                                           >  1:472261/1‑99 (MQ=255)
                                                                                                                        tGTTGCTGGCCACCCGCCTGGATGCGCGGCGCTTCCAGGTGGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTgg                                                                           <  2:472261/99‑1 (MQ=255)
                                                                                                                         gTTGCTGGCCACCCGCCTGGATGCGCGGCGCTTCCAGGTGGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCAAGTTCGCCAGGTTGAAGTCCATCAACCCGTCGGTACCGGt                       <  2:2361559/150‑1 (MQ=255)
                                                                                                                           tGCTGGCCACCCGCCTGGATGCGCGGCGCTTCCAGGTGGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCAAGTTCGCCAGGTTGAAGTCCATCAACCCGTCGGTACCGGTTa                     >  2:45855/1‑150 (MQ=255)
                                                                                                                             cTGGCCACCCGCCTGGATGCGCGGCGCTTCCAGGTGGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGc                                                                              >  1:977596/1‑91 (MQ=255)
                                                                                                                             cTGGCCACCCGCCTGGATGCGCGGCGCTTCCAGGTGGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGc                                                                              <  2:977596/91‑1 (MQ=255)
                                                                                                                             cTGGCCACCCGCCTGGATGCGCGGCGCTTCCAGGTGGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAg                                                                       <  2:1625576/98‑1 (MQ=255)
                                                                                                                             cTGGCCACCCGCCTGGATGCGCGGCGCTTCCAGGTGGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAg                                                                       >  1:1625576/1‑98 (MQ=255)
                                                                                                                               ggCCACCCGCATGGATGCGCGGCTCTTCCAGGTGGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCAAGTTCGCCAGGTTGAAGTCCATCAACCCGTCGGTACCGGTTAgcgc                 >  2:595213/1‑150 (MQ=255)
                                                                                                                                gCCACCCGCCTGGATGCGCGGCGCTTCCAGGTGGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCAAGTTCGCCAGGTTGAAGTCCATCAACCCGTCGGTACCGGTTAGcgcg                >  2:503320/1‑150 (MQ=255)
                                                                                                                                  tcaccaccTGGAAGCGCCGCGCATCCAGGTGGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCAAGTTCGCCAGGTTGAAGTCCATCAACCCGTCGGTACCGGTTAGCGCGGt              >  2:322132/7‑150 (MQ=255)
                                                                                                                                  cACCCGCCTGGATGCGCGGCGCTTCCAGGTGGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGaaa                                                                        >  1:888404/1‑92 (MQ=255)
                                                                                                                                  cACCCGCCTGGATGCGCGGCGCTTCCAGGTGGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGaaa                                                                        <  2:888404/92‑1 (MQ=255)
                                                                                                                                  cACCCGCCTGGATGCGCGGCGCTTCCAGGTGGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCAAGTTCGCCAGGTTGAAGTcc                                           <  1:1687875/121‑1 (MQ=255)
                                                                                                                                  cACCCGCCTGGATGCGCGGCGCTTCCAGGTGGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCAAGTTCGCCAGGTTGAAGTcc                                           >  2:1687875/1‑121 (MQ=255)
                                                                                                                                      caccTGGAAGCGCCGCGCATCCAGGTGGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCAAGTTCGCCAGGTTGAAGTCCATCAACCCGTCGGTACCGGTTAGCGCGGTGCTg          >  2:2287344/3‑150 (MQ=255)
                                                                                                                                       gCCTGGATGCGCGGCGCTTCCAGGTGGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCAAGTTCGCCAGGTTGAAGTCCATCaa                                      >  2:1635407/1‑121 (MQ=255)
                                                                                                                                       gCCTGGATGCGCGGCGCTTCCAGGTGGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCAAGTTCGCCAGGTTGAAGTCCATCaa                                      <  1:1635407/121‑1 (MQ=255)
                                                                                                                                        ccTGGAAGCGCCGCGCATCCAGGTGGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCAAGTTCGCCAGGTTGAAGTCCATCAACCCGTCGGTACCGGTTAGCGCGGTGCTGAt        <  1:322132/150‑1 (MQ=255)
                                                                                                                                        ccTGGAAGCGCCGCGCATCCAGGTGGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCAAGTTCGCCAGGTTGAAGTCCATCAACCCGTCGGTACCGGTTAGCGCGGTGCTGAt        >  2:737947/1‑150 (MQ=255)
                                                                                                                                         cTGGATGCGCGGCGCTTCCAGGTGGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGc                                                                      >  1:1817876/1‑87 (MQ=255)
                                                                                                                                         cTGGATGCGCGGCGCTTCCAGGTGGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGc                                                                      <  2:1817876/87‑1 (MQ=255)
                                                                                                                                            gATGCGCGGCGCTTCCAGGTGGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCAAGTTCGCCAGGTTGAAGTCCATCAACCCGTCGGTACCGGTTAGCGCGGTGCTGATCTTc    >  2:445693/1‑150 (MQ=255)
                                                                                                                                            gATGCGCGGCGCTTCCAGGTGGTGACCCTTGAGCCTCTGGTCCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCAAGTTCGCCAGTTTGAAGTCCATCAACCCGTCGGTACCGGTTAGCGCGGTGCTGATCTTc    >  2:86114/1‑150 (MQ=255)
                                                                                                                                             aTGCGCGGCGCTTCCAGGTGGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCAAGTTCGCCAGGTTGAAGTCCATCAACCCGTCGGTACCGGTTAGCGCGGTGCTGATCTTCa   <  2:1288312/150‑1 (MQ=255)
                                                                                                                                              tGCGCGGCGCTTCCAGGTGGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCAAGTTCGCCAGGTTGAAGTCCATCAACCCGTCGGTACCGGTTAGCGCGGTGCTGATCTTCAc  >  1:910991/1‑150 (MQ=255)
                                                                                                                                              |                                                                                                                                                     
CCACAGTGAAGGCACAGATCGACCTGCAGCCGATCAATGACCTGGCCGACGGTGCACAGATCGAATTGCGCCTGCCGCTGACGGTCACCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGATGCGCGGCGCTTCCAGGTGGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCAAGTTCGCCAGGTTGAAGTCCATCAACCCGTCGGTACCGGTTAGCGCGGTGCTGATCTTCAC  >  NC_002947/1033664‑1033955

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