Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,109,914 C→T 7.7% W1472* (TGG→TAG) ‡ PP_1882 ← membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,109,9140CT 7.7% 215.7 / 7.1 78W1472* (TGG→TAG) ‡PP_1882membrane protein
Reads supporting (aligned to +/- strand):  ref base C (40/32);  new base T (3/3);  total (43/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.90e-01

TTTGGAGCGCCGTACCGGTGACTGACAACGTGGCTGTAGCCACCGCGGCCGCCAACTTGATGCCGACGGCGAGCGATACGCTGGCAATTGCCCCGACACCTGCCGAAAGTGCCGCACCCGCGCCACTGATGCCGCCACCGGCGGCTGCCCAGGCAGCCGCGATTGCCGGCGCAGCTACGCCACCTGTCACCACGACGGCGGCCACGGAGGCAACTGCGGCGAAGAGTATGCCCAGGCCTATCCCCGTCCAGCCGATCGCACCTAAGGCCTTGCCCCGCCAGATGCTGTAGCTATTCAG  >  NC_002947/2109765‑2110062
                                                                                                                                                     |                                                                                                                                                    
tttGGAGCGCCGTACCGGTGACTGACAACGTGGCTGTAGCCACCGCGGCCGCCAACTTGATGCCGACGGCGAGCGATACGCTGGCAATTGCCCCGACACCTGCCGAAAGTGCCGCACCCGCGCCACTGATGCCGCCACCGGCGGCTGccc                                                                                                                                                      >  2:1705729/1‑150 (MQ=255)
  tGGAGCGCCGTACCGGTGACTGACAACGTGGCTGTAGCCACCGCGGCCGCCAACTTGATGCCGACGGCGAGCGATACGCTGGCAATTGCCCCGACACCTGCCGAAAGTGCCGCACCCGCGCCACTGATGCCGCCACCGGCGGCTGCCCAg                                                                                                                                                    <  1:657178/150‑1 (MQ=255)
     aGCGCCGTACCGGTGACTGACAACGTGGCTGTAGCCACCGCGGCCGCCAACTTGATGCCGACGGCGAGCGATACGCTGGCAATTGCCCCGACACCTGCCGAAAGTGCCGCACCCGCGCCACTGATGCCGCCACCGGCGGCTGCCCAGGCa                                                                                                                                                 <  2:365421/150‑1 (MQ=255)
      gcgcCGTACCGGTGACTGACAACGTGGCTGTAGCCACCGCGGCCGCCAACTTGATGCCGACGGCGAGCGATACGCTGGCAATTGCCCCGACACCTGCCGAAAGTGCCGCACCCGCGCCACTGATGCCGCCACCGGCGGCTGCCCAGGCAg                                                                                                                                                <  2:1210973/150‑1 (MQ=255)
      gcgcCGTACCGGTGACTGACAACGTGGCTGTAGCCACCGCGGCCGCCAACTTGATGCCGACGGCGAGCGATACGCTGGCAATTGCCCCGACACCTGCCGAAAGTGCCGCACCCGCGCCACTGATGCCGCCACCGGCGGCTGCCCAGGCAg                                                                                                                                                <  1:1705729/150‑1 (MQ=255)
         ccGTACCGGTGACTGACAACGTGGCTGTAGCCACCGCGGCCGCCAACTTGATGCCGACGGCGAGCGATACGCTGGCAATTGCCCCGACACCTGCCGAAAGTGCCGCACCCGCGCCACTGATGCCGCCACCGGCGGCTGCCCAGGCAGCcg                                                                                                                                             >  2:234230/1‑150 (MQ=255)
         ccGTACCGGTGACTGACAACGTGGCTGTAGCCACCGCGGCCGCCAACTTGATGCCGACGGCGAGCGATACGCTGGCAATTGCCCCGACACCTGCCGAAAGTGCCGCACCCGCGCCACTGATGCCGCCACCGGCGGCTGCCCAGGCAGCcg                                                                                                                                             >  1:1188228/1‑150 (MQ=255)
            tACCGGTGACTGACAACGTGGCTGTAGCCACCGCGGCCGCCAACTTGATGCCGACGGCGAGCGATACGCTGGCAATTGCCCCGACACCTGCTGAAAGTGCCGCACCCGCGCCACTGATGCCGCCACCGGCGGCTGCCCAGGCAGCCGCGa                                                                                                                                          <  2:984929/150‑1 (MQ=255)
            tACCGGTGACTGACAACGTGGCTGTAGCCACCGCGGCCGCCAACTTGATGCCGACGGCGAGCGATACGCTGGCAATTGCCCCGACACCTGCCGAAAGTGCCGCGCCCGCGCCACTGATGCCGCCACCGGCGGCTGCCCAGGCAGCCGCGa                                                                                                                                          <  2:1706593/150‑1 (MQ=255)
                      tgacAACGTGGCTGTAGCCACCGCGGCCGCCAACTTGATGCCGACGGCGAGCGATACGCTGGCAATTGCCCCGACACCTGCCGAAAGTGCCGCACCCGCGCCACTGATGCCGCCACCGGCGGCTGCCCAGGCAGCCGCGATTGCCGgcgc                                                                                                                                >  2:1126161/1‑150 (MQ=255)
                        acAACGTGGCTGTAGCCACCGCGGCCGCCAACTTGATGCCGACGGCGAGCGATACGCTGGCAATTGCCCCGACACCTGCCGAAAGTGCCGCACCCGCGCCACTGATGCCGCCACCGGCGGCTGCCCAGGCAGCCGCGATTGCCGgcgccg                                                                                                                              >  2:924458/1‑148 (MQ=255)
                        acAACGTGGCTGTAGCCACCGCGGCCGCCAACTTGATGCCGACGGCGAGCGATACGCTGGCAATTGCCCCGACACCTGCCGAAAGTGCCGCACCCGCGCCACTGATGCCGCCACCGGCGGCTGCCCAGGCAGCCGCGATTGCCGGCGCAg                                                                                                                              >  1:203836/1‑150 (MQ=255)
                            cGTGGCTGTAGCCACCGCGGCCGCCAACTTGATGCCGACGGCGAGCGATTCTCTGGCAATTGCCCCGACACCTGCCGAAAGTGCCGCACCCGCGCCACTGATGCCGCCACCGGCGGCTGCCCAGGCAGCCGCGATTGCCGGCGCAGCTAc                                                                                                                          <  1:593674/150‑1 (MQ=255)
                            cGTGGCTGTAGCCACAGCGGAATACCAGTTGTTGCTGATTGCGAGCGATACGCTGGCAATTGCCCCGACACCTGCCGAAAGTGCCGCACCCGCGCCACTGATGCCGCCACCGGCGGCTGCCCACGCAGCCGCGATTGCCGGCGCAGCTAc                                                                                                                          <  2:959063/150‑1 (MQ=255)
                              tGGCTGTAGCCACCGCGGCCGCCAACTTGATGCCGACGGCGAGCGATACGCTGGCAATTGCCCCGACACCTGCCGAAAGTGCCGCACCCGCGCCACTGATGCCGCCACCGGCGGCTGCCCAGGCAGCCGCGATTGCCGGCGCAGCTACGc                                                                                                                        >  1:394447/1‑150 (MQ=255)
                                   gTAGCCACCGCGGCCGCCAACTTGATGCCGACGGCGAGCGATACGCTGGCAATTGCCCCGACACCTGCCGAAAGTGCCGCACCCGCGCCACTGATGCCGCCACCGGCGGCTGCCCAGGCAGCCGCGATTGCCGGCCCAGCTACGCCACCt                                                                                                                   <  1:1004976/150‑1 (MQ=255)
                                    tAGCCACCGCGGCCGCCAACTTGATGCCGACGGCGAGCGATACGCTGGCAATTGCCCCGACACCTGCCGAAAGTGCCGCACCCGCGCCACTGATGCCGCCACCGGCGGCTGCCCAGGCAGCCGCGATTGCCGGCGCAGCTACGCCACCTg                                                                                                                  >  1:1103027/1‑150 (MQ=255)
                                     ggcagccGCGGCCGCCAATTTGATGCCGACGGCGAGCGAGACGCTGGCAATTGCCCCGACACCTGCCGAAAGTGCCGCACCCGCGCCACTGATGCCGCCACCGGCGGCTGCCCAGGCAGCCGCGATTGCCGGCGCAGCTACGCCACCTGt                                                                                                                 <  2:447067/145‑1 (MQ=255)
                                      gCCACCGCGGCCGCCACCTTGATGCCGACGGCGAGCGATACGCTGGCAATTGCCCCGACACCTGCCGAAAGTGCCGCACCCGCGCCACTGATGCCGCCACCGGCGGCTGCCCAGGCAGCCGCGATTGCCGGCGCAGCTACGCCACCTGTc                                                                                                                <  1:1265734/150‑1 (MQ=255)
                                       ccACCGCGGCCGCCAACTTGATGCCGACGGCGAGCGATACGCTGGCAATTGCCCCGACACCTGCCGAAAGTGCCGCACCCGCGCCACTGATGCCGCCACCGGCGGCTGCCCAGGCAGCCGCGATTGCCGGCGCAGCTACGCCACCTGTca                                                                                                               >  1:950571/1‑150 (MQ=255)
                                       ccACCGCGGCCGCCAACTTGATGCCGACGGCGAGCCATACGCTGGCAATTGCCCCGCCACCTGCCCAAAGTGCCGCACCCGCGCCACTGATGCCACCACCCGCCGCACCCCAGCCAGCCGAGATTAACGGCGCAAATACGACATATGTca                                                                                                               >  2:871470/1‑150 (MQ=255)
                                         aCCGCGGCCGCCAACTTGATGCCGACGGCGAGCGATACGCTGGCAATTGCCCCGACACCTGCCGAAAGTGCCGCACCCGCGCCACTGATGCCGCCACCGGCGGCTGCCCAGGCAGCCGCGATTGCCGGCGCAGCTACGCCACCTGTcacc                                                                                                             >  2:482701/1‑150 (MQ=255)
                                                ccgccAACTTGATGCCGACGGCGAGCGATACGCTGGCAATTGCCCCGACACCTGCCGAAAGTGCCGCACCCGCGCCACTGATGCCGCCACCGGCGGCTGCCCAGGCAGCCGCGATTGCCGGCGCAGCTACGCCACCTGTCAACAcggcgg                                                                                                      >  2:764249/1‑150 (MQ=255)
                                                 cgccAACTTGATGCCGACGGCGAGCGATACGCTGGCAATTGCCCCGACACCTGCCGAAAGTGCCGCACCCGCGCCACTGATGCCGCCACCGGCGGCTGCCCAGGCAGCCGCGATTGCCGGCGCAGCTACGCCACCTGTCACCACGAcggg                                                                                                     >  1:1843803/1‑149 (MQ=255)
                                                   ccAACTTGATGCCGACGGCGAGCGATACGCTGGCAATTGCCCCGACACCTGCCGAAAGTGCCGCACCCGCGCCACTGATGCCGCCACCGGCGGCTGCCCAGGCAGCCGCGATTGCCGGCGCAGCTACGCCACCTGTCACCACGAcggcgg                                                                                                   >  2:1344446/1‑150 (MQ=255)
                                                       cTTGATGCCGACGGCGAGCGATACGCTGGCAATTGCCCCGACACCTGCCGAAAGTGCCGCACCCGCGCCACTGATGCCGCCACCGGCGGCTGCCCAGGCAGCCGCGATTGCCGGCGCAGCTACGCCACCTGTCACCACGACGGCGGCCAc                                                                                               <  1:1126161/150‑1 (MQ=255)
                                                       cTTGATGCCGACGGCGAGCGATACGCTGGCAATTGCCCCGACACCTGCCGAAAGTGCCGCACCCGCGCCACTGATGCCGCCAACGGCGGCTGCCCAGGCAGCCGCGATTGCCGGCGCACCTACGCCACCTGTCACCACGACCGCGGCCAc                                                                                               >  2:1365429/1‑150 (MQ=255)
                                                        ttGATGCCGACGGCGAGCGATACGCTGGCAATTGCCCCGACACCTGCCGAAAGTGCCGCACCCGCGCCACTGATGCCGCCACCGGCGGCTGCCCAGGCAGCCGCGATTGCCGGCGCAGCTACGCCACCTGTCACCACGACGGCGGCCACg                                                                                              >  1:1423228/1‑150 (MQ=255)
                                                            tGCCGACGGCGAGCGATACGCTGGCAATTGCCCCGACACCTGCCGAAAGTGCCGCACCCGCGCCACTGATGCCGACACCGGCGGCTGCCCAGGCAGCCGCGATTGCCGGCGCAGCTACGCCACCTGTCACCACGACGGCGGCCACGGAgg                                                                                          >  1:124253/1‑150 (MQ=255)
                                                               cGACGGCGAGCGATACGCTGGCAATTGCCCCGACACCTGCCGAAAGTGCCGCACCCGCGCCACTGATGCCGCCACCGGCGGCTGCCTGGGCAGCcgc                                                                                                                                            >  2:1665215/1‑97 (MQ=255)
                                                               cGACGGCGAGCGATACGCTGGCAATTGCCCCGACACCTGCCGAAAGTGCCGCACCCGCGCCACTGATGCCGCCACCGGCGGCTGCCTGGGCAGCcgc                                                                                                                                            <  1:1665215/97‑1 (MQ=255)
                                                               cGACGGCGAGCGATACGCTGGCAATTGCCCCGACACCTGCCGAAAGTGCCGCACCCGCGCCACTGATGCCGCCACCGGCGGCTGCCCAGGCAGCCGCGATTGCCGGCGCAGCTACGCCACCTGTCACCACGACGGCGGCCACGGAGGCaa                                                                                       >  1:1113931/1‑150 (MQ=255)
                                                                      gagcgaTACGCTGGCAATTGCCCCGACACCTGCCGAAAGTGCCGCACCCGCGCCACTGATGCCGCCACCGGCGGCTGCCTGGGCAGCcgc                                                                                                                                            >  2:1629850/1‑90 (MQ=255)
                                                                      gagcgaTACGCTGGCAATTGCCCCGACACCTGCCGAAAGTGCCGCACCCGCGCCACTGATGCCGCCACCGGCGGCTGCCTGGGCAGCcgc                                                                                                                                            <  1:1629850/90‑1 (MQ=255)
                                                                      gagcgaTACGCTGGCAATTGCCCCGACACCTGCCGAAAGTGCCGCACCCGCGCCACTGATGCCGCCACCGGCGGCTGCCCAGGCAGCCGCGATTGCCGGCGCAGCTACGCCACCTGTCACCACGACGGCGGCCACGGAGGCAACTgcggc                                                                                <  2:1062809/150‑1 (MQ=255)
                                                                         cgaTACGCTGGCAATTGCCCCGACACCTGCCGAAAGTGCCGCACCCGCGCCACTGATGCCGCCACCGGCGGCTGCCCAGGCAGCCGCGATTGCCGGCGCAGCTACGCCACCTGTCACCACGACGGCGGCCACGGAGGCAACTGCGGCGaa                                                                             <  1:234230/150‑1 (MQ=255)
                                                                           aTACGCTGGCAATTGCCCCGACACCTGCCGAAAGTGCCGCACCCGCGCCACTGATGCCGCCACCGGCGGCTGCCCAGGCAGCCGCGATTGCCGGCGCAGCTACGCCACCTGTCACCACGACGGCGGCCACGGAGGCAACTGCGGCGAaga                                                                           >  2:1430652/1‑150 (MQ=255)
                                                                             aCGCTGGCAATTGCCCCGACACCTGCCGAAAGTGCCGCACCCGCGCCACTGATGCCGCCACCGGCGGCTGCCCAGGCAGCCGCGATTGCCGGCGCAGCTACGCCACCTGTCACCACGACGGCGGCCCCGGAGGCAACTGCGGCGAAGAGt                                                                         <  2:950571/150‑1 (MQ=255)
                                                                                     aaTTGCCCCGACACCTGCCGAAAGTGCCGCACCCGCGCCACTGATGCCGCCACCGGCGGCTGCCCAGGCAGCCGCGATTGCCGGCGCAGCTACGCCACCTGTCACCACGACGGCGGCCACGGAGGCAACTGCGGCGAAGAGTATGCCCAg                                                                 >  2:1652082/1‑150 (MQ=255)
                                                                                     aaTTGCCCCGACACCTGCCGAAAGTGCCGCACCCGCGCCACTGATGCCGCCACCGGCGGCTGCCCAGGCAGCCGCGATTGCCGGCGCAGCTACGCCACCTGTCACCACGACGGCGGCCACGGAGGCAACTGCGGCGAAGAGTATGCCCAg                                                                 >  2:1077916/1‑150 (MQ=255)
                                                                                       ttGCCCCGACACCTGCCGAAAGTGCCGCACCCGCGCCACTGATGCCGCCACCGGCGGCTGCCCAGGCAGCCGCGATTGCCGGCGCAGCTACGCCACCTGTCACCACGACGGCGGCCACGGAGGCAACTGCGGCGAAGAGTATGCCCAGGc                                                               <  2:1423228/150‑1 (MQ=255)
                                                                                           cccGACACCTGCCGAAAGTGCCGCACCCGCGCCACTGATGCCGCCACCGGCGGCTGCCCAGGCAg                                                                                                                                                >  1:859709/1‑65 (MQ=255)
                                                                                           cccGACACCTGCCGAAAGTGCCGCACCCGCGCCACTGATGCCGCCACCGGCGGCTGCCCAGGCAg                                                                                                                                                <  2:859709/65‑1 (MQ=255)
                                                                                           cccGACACCTGCCGAAAGTGCCGCACCCGCGCCACTGATGCCGCCACCGGCGGCTGCCCAGGCAGCCGCGATTGCCGGCGCAGCTACGCCACCTGTCACCACGACGGCGGCCACGGAGGCAACTGCGGCGAAGAGTATg                                                                      <  2:1688422/139‑1 (MQ=255)
                                                                                           cccGACACCTGCCGAAAGTGCCGCACCCGCGCCACTGATGCCGCCACCGGCGGCTGCCCAGGCAGCCGCGATTGCCGGCGCAGCTACGCCACCTGTCACCACGACGGCGGCCACGGAGGCAACTGCGGCGAAGAGTATg                                                                      >  1:1688422/1‑139 (MQ=255)
                                                                                           cccGACACCTGCCGAAAGTGCCGCACCCGCGCCACTGATGCCGCCACCGGCGGCTGCCCAGGCAGCCGCGATTGCCGGCGCAGCTACGCCACCTGTCACCACGACGGCGGCCACGGAGGCAACTGCGGCGAAGAGTATGCCCAGGCCTAt                                                           >  1:310995/1‑150 (MQ=255)
                                                                                           cccGACACCTGCCGAAAGTGCCGCACCCGCGCCACTGATGCCGCCACCGGCGGCTGCCCAGGCAGCCGCGATTGCCGGCGCAGCTACGCCACCTGTCACCACGACGGCGGCCACGGAGGCAACTGCGGCGAAGAGTATGCCCAGGCCTAt                                                           >  1:336640/1‑150 (MQ=255)
                                                                                             cGACACCTGCCGAAAGTGCCGCACCCGCGCCACTGATGCCGCCACCGGCGGCTGCCCAGGCAGCCGCGATTGCCGGCGCAGCTACGCCACCTGTCACCACGACGGCGGCCACGGAGGCAACTGCGGCGAAGAGTATGCCCAGGCCTATcc                                                         >  1:1761118/1‑150 (MQ=255)
                                                                                                 acCTGCCGAAAGTGCCGCACCCGCGCCACTGATGCCGCCACCGGCGGCTGCCCa                                                                                                                                                     >  2:565467/1‑54 (MQ=255)
                                                                                                 acCTGCCGAAAGTGCCGCACCCGCGCCACTGATGCCGCCACCGGCGGCTGCCCa                                                                                                                                                     <  1:565467/54‑1 (MQ=255)
                                                                                                 acCTGCCGAAAGTGCCGCACCCGCGCCACTGATGCCGCCACCGGCGGCTGCCCAGGCAGCCGCGATTGCCGGCGCAGCTACGCCACCTGTCACCACGACGGCGGCCACGGAGGCAACTGCGGCGAAGAGTATGCCCAGGCCTATCCCCGt                                                     >  2:195760/1‑150 (MQ=255)
                                                                                                  ccTGCCGAAAGTGCCGCACCCGCGCCACTGATGCCGCCACCGGCGGCTGCCCAGGCAGCCGCGATTGCCGGCGCAGCTACGCCACCTGTCACCACGACGGCGGCCACGGAGGCAACTGCGGCGAAGAGTATGCCCAGGCCTATCCCCGTc                                                    <  2:1103027/150‑1 (MQ=255)
                                                                                                    tGCCGAAAGTGCCGCACCCGCGCCACTGATGCCGCCACCGGCGGCGGCCCAGGCAGCCGCGATTGCCGGCGCAGCTACGCCACCTGTCACCACGACGGCGGCCACGGAGGCACCTGCGGCGAAGAGTATGCCCAGGCCTATCCCCGTCCa                                                  >  2:851609/1‑150 (MQ=255)
                                                                                                        gAAAGTGCCGCACCCGCGCCACTGATGCCGCCACCGGCGGCTGCCTGGGCAGCcgc                                                                                                                                            >  2:1321525/1‑56 (MQ=38)
                                                                                                        gAAAGTGCCGCACCCGCGCCACTGATGCCGCCACCGGCGGCTGCCTGGGCAGCcgc                                                                                                                                            <  1:1321525/56‑1 (MQ=38)
                                                                                                             ctgggCACCCGCGCCACTGATGCCGCCACCGGCGGCTGCCCAGGCAGCCGCGATTGCCGGCGCAGCTACGCCACCTGTCACCACGACGGCGGCCACGGAGGCAACTGCGGCGAAGAGTATGCCCAGGCCTATCCCCGTCCAGCCGATCGc                                         >  2:1085346/5‑150 (MQ=255)
                                                                                                                cGCACCCGCGCCACTGATGCCGCCACCGGCGGCTGCCCAGGCAGCCGCGATTGCCGGCGCAGCTACGCCACCTGTCACCACGACGGCGGCCACGGAGGCAACTGCGGCGAAGAGTATGCCCAGGCCTATCCCCGTCCAGCCGATCGCAcc                                      >  2:350924/1‑150 (MQ=255)
                                                                                                                 gCACCCGCGCCACTGATGCCGCCACCGGCGGCTGCCCAGGCAGCCGCGATTGCCGGCGCAGCTACGCCACCTGTCACCACGACGGCGGCCACGGAGGCAACTGCGGCGAAGAGTATGCCCAGGCCTATCCCCGTCCAGCCGATCGCACCt                                     >  1:1073482/1‑150 (MQ=255)
                                                                                                                    cccGCGCCACTGATGCCGCCACCGGCGGCTGCCCAGGCAGCCGCGATTGCCGGCGCAGCTACGCCACCTGTCACCACGACGGCGGCCACGGAGGCAACTGCGGCGAAGAGTATGCCCAGGCCTATCCCCGTCCAGCCGATCGCACCTAAg                                  <  1:1085346/150‑1 (MQ=255)
                                                                                                                     ccGCGCCACTGATGCCGCCACCGGCGGCTGCCCAGGCAGCCGCGATTGCCGGCGCAGCTACGCCACCTGTCACCACGACGGCGGACACGGAGGACACTACGGGGAATAGTAGGCCCAGGCCTATCCCCGTCCAGCCGATCGCACCTAAgg                                 >  1:1870447/1‑150 (MQ=255)
                                                                                                                        cgcCACTGATGCCGCCACCGGCGGCTGCCCAGGCAGCCGCGATTGCCGGCGCAGCTACGCCACCTGTCACCACGACGGCGGCCACGGAGGCAACTGCGGCGAAGAGTATGCCCAGGCCTATCCCCGTCCAACCGATCGCACCTAAGGCCt                              <  2:1352541/150‑1 (MQ=255)
                                                                                                                            aCTGATGCCGCCACCGGCGGCTGCCCAGGCAGCCGCGATTGCCGGCGCAGCTACGCCACCTGTCACCACGACGGCGGCCACGGAGGCAACTGCGGCGAAGAGTATGCCCAGGCCTATCCCCGTCCAGCCGATCGCACCTAAGGCCTTGcc                          >  2:1440441/1‑150 (MQ=255)
                                                                                                                            aCTGATGCCGCCACCGGCGGCTGCCCAGGCAGCCGCGATTGCCGGCGCAGCTACGCCACCTGTCACCACGACGGCGGCCACGGAGGCAACTGCGGCGAAGAGTATGCCCAGGCCTATCCCCGTCCAGCCGATCGCACCTAAGGCCTTGcc                          >  1:1630668/1‑150 (MQ=255)
                                                                                                                             cTGATGCCGCCACCGGCGGCTGCCCAGGCAGCCGCGATTGCCGGCGCAGCTACGCCACCTGTCACCACGACGGCGGCCACGGAGGCAACTgcg                                                                                  >  1:116966/1‑93 (MQ=255)
                                                                                                                             cTGATGCCGCCACCGGCGGCTGCCCAGGCAGCCGCGATTGCCGGCGCAGCTACGCCACCTGTCACCACGACGGCGGCCACGGAGGCAACTgcg                                                                                  <  2:116966/93‑1 (MQ=255)
                                                                                                                              tGATGCCGCCACCGGCGGCTGCCCAGGCAGCCGCGATTGCCGGCGCAGCTACGCCACCTGTCACCACGACGGCGGCCACGGAGGCAACTGCGGCGAAGAGTATGCCCAGGCCTATCCCCGTCCAGCCGATCGCACCTAAGGCCTTGcccc                        <  1:337180/150‑1 (MQ=255)
                                                                                                                                   ccgccACCGGCGGCTGCCCAGGCAGCCGCGATTGCCGGCGCAGCTACGCCACCTGTCACCACGACGGCGGCCACGGAGGCAACTGCGGCGAAGAGTATGCCCAGGCCTATCCCCGTCCAGCCGATCGCACCTAAGGCCTTGc                           <  1:170967/142‑1 (MQ=255)
                                                                                                                                   ccgccACCGGCGGCTGCCCAGGCAGCCGCGATTGCCGGCGCAGCTACGCCACCTGTCACCACGACGGCGGCCACGGAGGCAACTGCGGCGAAGAGTATGCCCAGGCCTATCCCCGTCCAGCCGATCGCACCTAAGGCCTTGc                           >  2:170967/1‑142 (MQ=255)
                                                                                                                                    cgccACCGGCGGCTGCCCAGGCAGCCGCGATTGCCGGCGCAGCTACGCCACCTGTCACCACGACGGCGGCCACGGAGGCAACTGCGGCGAAGAGTATGCCCAGGCCTATCCCCGTCCAGCCGATCGCACCTAAGGCCTTGCCCCGCCAGa                  >  2:480572/1‑150 (MQ=255)
                                                                                                                                     gccACCGGCGGCTGCCCAGGCAGCCGCGATTGCCGGCGCAGCTACGCCACCTGTCACCACGACGGCGGCCACGGAGGCAACTGCGGCGAAGAGTATGCCCAGGCCTATCCCCGTCCAGCCGATCGCACCTAAGGCCTTGCCCCGCCAGAt                 <  2:1761118/150‑1 (MQ=255)
                                                                                                                                       cACCGGCGGCTGCCCAGGCAGCCGCGATTGCCGGCGCAGCTACGCCACCTGTCACCACGACGGCGGCCACGGAGGCAACTGCGGCGAAGAGTATGCCCAGGCCTATCCCCGTCCAGCCGATCGCACCTAAGGCCTTGCCCCGCCAGa                  >  1:1396929/1‑147 (MQ=255)
                                                                                                                                       cACCGGCGGCTGCCCAGGCAGCCGCGATTGCCGGCGCAGCTACGCCACCTGTCACCACGACGGCGGCCACGGAGGCAACTGCGGCGAAGAGTATGCCCAGGCCTATCCCCGTCCAGCCGATCGCACCTAAGGCCTTGCCCCGCCAGa                  <  2:1396929/147‑1 (MQ=255)
                                                                                                                                         ccGGCGGCTGCCCAGGCAGCCGCGATTGCCGGCGCAGCTACGCCACCTGTCACCACGACGGCGGCCACGGAGGCAACTGCGGCGAAGAGTATGCCCAGGCCTATCCCCGTCCAGCCGATCGCACCTAAGGCCTTGCCCCGCCAGATGCTg             <  1:481203/150‑1 (MQ=255)
                                                                                                                                           ggcggcTGCCCAGGCAGCCGCGATTGCCGGCGCAGCTACGCCACCTGTCACCACGACGGCGGCCACGGAGGCAACTGCGGCGAAGAGTATGCCCAGGCCTATCCCCGTCCAGCCGATCGCACCTAAGGCCTTGCCCCGCCAGATGCTGTa           <  1:1430652/150‑1 (MQ=255)
                                                                                                                                           ggcggcTGCCCAGGCAGCCGCGATTGCCGGCGCAGCTACGCCACCTGTCACCACGACGGCGGCCACGGAGGCAACTGCGGCGAAGAGTATGCCCAGGCCTATCCCCGTCCAGCCGATCGCACCTAAGGCCTTGCCCCGCCAGATGCTGTa           <  1:1656757/150‑1 (MQ=255)
                                                                                                                                            gcggcTGCCCAGGCAGCCGCGATTGCCGGCGCAGCTACGCCACCTGTCACCACGACGGCGGCCACGGAGGCAACTGCGGCGAAGAGTATGCCCAGGCCTATCCCCGTCCAGCCGATCGCACCTAAGGCCTTGCCCCGCCAGATGCTGTAg          >  2:1408412/1‑150 (MQ=255)
                                                                                                                                              ggcTGCCCAGGCAGCCGCGATTGCCGGCGCAGCTACGCCACCTGTCACCACGACGGCGGCCACGGAGGCAACTGCGGCGAAGAGTATGCCCAGGCCTATCCCCGTCCAGcc                                               <  2:1252006/111‑1 (MQ=255)
                                                                                                                                              ggcTGCCCAGGCAGCCGCGATTGCCGGCGCAGCTACGCCACCTGTCACCACGACGGCGGCCACGGAGGCAACTGCGGCGAAGAGTATGCCCAGGCCTATCCCCGTCCAGcc                                               >  1:1252006/1‑111 (MQ=255)
                                                                                                                                                cTGCCCAGGCAGCCGCGATTGCCGGCGCAGCTACGCCACCTGTCACCACGACGGCGGCCACGGAGGCAACTGCGGCGAAGAGTATGCCCAGGCCTATCCCCGTCCAGCCGATCGCACCTAAGGCCTTGCCCCGCCAGATGCTGTAGCTAt      <  2:1073482/150‑1 (MQ=255)
                                                                                                                                                    ccAGGCAGCCGCGATTGCCGGCGCAGCTACGCCACCTGTCACCACGACGGCGGCCACGGAGGCAACTGCGGCGAAGAGTATGCCCAGGCCTATCGCCGTCCAGCCGATCGCACCTAAGGCCTTGCCCCGCCAGATGCTGTAGCTATTCAg  >  1:1530187/1‑150 (MQ=255)
                                                                                                                                                    ccAGGCAGCCGCGATTGCCGGCGCAGCTACGCCACCTGTCACCACGACGGCGGCCACGGAGGCAACTGCGGCGAAGAGTATGCCCAGGCCTATCCCCGTCCAGCCGATCGCACCTAAGGCCTTGCCCCGCCAGATGCTGTAGCTATTCAg  >  2:77183/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
TTTGGAGCGCCGTACCGGTGACTGACAACGTGGCTGTAGCCACCGCGGCCGCCAACTTGATGCCGACGGCGAGCGATACGCTGGCAATTGCCCCGACACCTGCCGAAAGTGCCGCACCCGCGCCACTGATGCCGCCACCGGCGGCTGCCCAGGCAGCCGCGATTGCCGGCGCAGCTACGCCACCTGTCACCACGACGGCGGCCACGGAGGCAACTGCGGCGAAGAGTATGCCCAGGCCTATCCCCGTCCAGCCGATCGCACCTAAGGCCTTGCCCCGCCAGATGCTGTAGCTATTCAG  >  NC_002947/2109765‑2110062

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