Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 39,490 C→A 5.2% V73V (GTC→GTA PP_0036 → transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00294739,4900CA 5.2% 235.3 / 2.0 76V73V (GTC→GTAPP_0036transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base C (35/37);  new base A (2/2);  total (37/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.88e-01

CAATTGATGATCAGCTGGCCGTAGAGTTTTTCGCAACAGCACGTTGTGCCAGTTTCAAACAGGCTGCCCGAGGGCTCAATGTGCCAGTGGTAGGCCTGCGGAAACGCTTGGAGAAACTGGAAGAGCACATTGGTGCTCCCGTGTTTGTCTACAAACACAACAAGCTTGCCCTGACCCGCACAGGCGAGAGGGTCAGCCATTACCTGTGCCAGCTCTTCGGCCCAGACGGCCTTGGCAAGTCCGGTGAAAAGGAAGTCCCCAGGCTTACCATCGCGATCACAGAGCCAGTGCTGAACG  >  NC_002947/39342‑39638
                                                                                                                                                      |                                                                                                                                                     
cAATTGATGATCAGCTGGCCGTAGAGTTTTTCGCAACAGCACGTTGTGCCAGTTTCAAACAGGCTGCCCGAGGGCTCAATGTGCCAGTGGTAGGCCTGCGGAAACGCTTGGAGAAACTGGAAGAGCACATTGGTGCTCCCGTGTTTGTc                                                                                                                                                       <  2:343857/149‑1 (MQ=255)
  aTTGATGATCAGCTGGCCGTAGAGTTTTTCGCAACAGCACGTTGTGCCAGTTTCAAACAGGCTGCCCGAGGGCTCAATGTGCCAGTGGTAGGCCTGCGGAAACGCTTGGAGAAACTGGAAGAGCACATTGGTGCTCCCGTGTTTGTCTa                                                                                                                                                     >  2:54042/1‑149 (MQ=255)
          tCAGCTGGCCGTAGAGTTTTTCGCAACAGCACGTTGTGCCAGTTTCAAACAGGCTGCCCGAGGGCTCAATGTGCCAGTGGTAGGCCTGCGGAAACGCTTGGAGAAACTGGAAGAGCACATTGGTGCTCCCGTGTTTGTCTACAAacac                                                                                                                                              <  4:312442/148‑1 (MQ=255)
          tCAGCTGGCCGTAGAGTTTTTCGCAACAGCACGTTGTGCCAGTTTCAAACAGGCTGCCCGAGGGCTCAATGTGCCAGTGGTAGGCCTGCGGAAACGCTTGGAGAAACTGGAAGAGCACATTGGTGCTCCCGTGTTTGTCTACAAacac                                                                                                                                              >  8:361891/1‑148 (MQ=255)
             gCTGGCCGTAGAGTTTTTCGCAACAGCACGTTGTGCCAGTTTCAAACAGGCTGCCCGAGGGCTCAATGTGCCAGTGGTAGGCCTGCGGAAACGCTTGGAGAAACTGGAAGAGCACATTGGTGCTCCCGTGTTTGTAGACAAacacggga                                                                                                                                          <  1:380762/149‑5 (MQ=255)
             gCTGGCCGTAGAGTTTTTCGCAACAGCACGTTGTGCCAGTTTCAAACAGGCTGCCCGAGGGCTCAATGTGCCAGTGGTAGGCCTGCGGAAACGCTCGGAGAAACTGGAAGAGCACATTGGTGCTCCCGTGTTTGTAGACAAacacggga                                                                                                                                          <  5:378942/149‑5 (MQ=255)
              cTGGCCGTAGAGTTTTTCGCAACAGCACGTTGTGCCAGTTTCAAACAGGCTGCCCGAGGGCTCAATGTGCCAGTGGTAGGCCTGCGGAAACGCTTGGAGAAACTGGAAGAGCACATTGGTGCTCCCGTGTTTGTCTACAAACAcaacaa                                                                                                                                         <  5:38419/149‑1 (MQ=255)
              cTGGCCGTAGAGTTTTTCGCAACAGCACGTTGTGCCAGTTTCAAACAGGCTGCCCGAGGGCTCAATGTGCCAGTGGTAGGCCTGCGGAAACGCTTGGAGAAACTGGAAGAGCACATTGGTGCTCCCGTGTTTGTCTACAAACAcaacaa                                                                                                                                         <  4:197300/149‑1 (MQ=255)
                  ccGTAGAGTTTTTCGCAACAGCACGTTGTGCCAGTTTCAAACAGGCTGCCCGAGGGCTCAATGTGCCAGTGGTAGGCCTGCGGAAACGCTTGGTGAAACTGGAAGAGCACATTGGTGCTCCCGTGTTTGTCTACAAACACAACAAGCtt                                                                                                                                     <  1:305750/149‑1 (MQ=255)
                   cGTAGAGTTTTTCGCAACAGCACGTTGTGCCAGTTTCAAACAGGCTGCCCGAGGGCTCAATGTGCCAGTGGTAGGCCTGCGGAAACGCTTGGAGAAACTGGAAGAGCACATTGGTGCTCCCGTGTTTGTCTACAAACACAACAAGCTTg                                                                                                                                    >  2:90958/1‑149 (MQ=255)
                     tAGAGTTTTTCGCAACAGCACGTTGTGCCAGTTTCAAACAGGGCTGCCCGAGGGCTCAATGTGCCAGTGGTAGGCCTGCGGAAACGCTTGGAGAAACTGGAAGAGCACATTGGTGCTCCCGTGTTTGTCTACAAACACAACAAGCTTGc                                                                                                                                   >  1:273884/1‑149 (MQ=255)
                     tAGAGTTTTTCGCAACAGCACGTTGTGCCAGTTTCAAACAGGCTGCCCGAGGGCTCAATGTGCCAGTGGTAGGCCTGCGGAAACGCTTGGAGAAACTGGAAGAGCACATTGGTGCTCCCGTGTTTGTCTACAAACACAACAAGCTTGcc                                                                                                                                  >  6:384243/1‑149 (MQ=255)
                      agagTGTCTCGCCACAGCACGATCTGCCAGTTTCAAACAGGCTGCCCGAGGGCTCAATGTGCCAGTGGTAGGCCTGCGGAAACGCTTGGAGAAACTGGAAGAGCACATTGGTGCTCCCGTGTTTGTCTACAAACACAACAAGCTTGccc                                                                                                                                 <  4:246798/149‑1 (MQ=255)
                        agTTTTTCGCAACAGCACGTTGTGCCAGTTTCAAACAGGCTGCCCGAGGGCTCAATGTGCCAGTGGTAGGCCTGCGGAAACGCTTGGAGAAACTGGAAGAGCACATTGGTGCTCCCGTGTTTGTCTACAAACACAACAAGCTTGCCCt                                                                                                                                >  7:87081/1‑148 (MQ=255)
                             ttCGCAACAGCACGTTGTGCCAGTTTCAAACAGGCTGCCCGAGGGCTCAATGTGCCAGTGGTAGGCCTGCGGAAACGCTTGGAGAAACTGGAAGAGCACATTGGTGCTCCCGTGTTTGTCTACAAACACAACAATCTTGCCCTGAccc                                                                                                                           <  5:372965/148‑1 (MQ=255)
                                gCAACAGCACGTTGTGCCAGTTTCAAACAGGCTGCCCGAGGGCTCAATGTGCCAGTGGTAGGCCTGCGGAAACGCTTGGAGAAACTGGAAGAGCACATTGGTGCTCCCGTGTTTGTCTACAAACACAACAAGCTTGCCCTGACCCGca                                                                                                                        >  6:189931/1‑148 (MQ=255)
                                gCAACAGCACGTTGTGCCAGTTTCAAACAGGCTGCCCGAGGGCTCAATGTGCCAGTGGTAGGCCTGCGGAAACGCTTGGAGAAACTGGAAGAGCACATTGGTGCTCCCGTGTTTGTCTACAAACACAACAAGCTTGCCCTGACCCGca                                                                                                                        >  8:263333/1‑148 (MQ=255)
                                      gCACGTTGTGCCAGTTTCAAACAGGCTGCCCGAGGGCTCAATGTGCCAGTGGTAGGCCTGCGGAAACGCTTGGAGAAACTGGAAGAGCACATTGGTGCTCCCGTGTTTGTCTACAAACACAACAAGCTTGCCCTGACCCGCACAGGc                                                                                                                   >  5:353767/1‑147 (MQ=255)
                                       cACGTTGTGCCAGTTTCAAACAGGCTGCCCGAGGGCTCAATGTGCCAGTGGTAGGCCTGCGGAAACGCTTGGAGAAACTGGAAGAGCACATTGGTGCTCCCGTGTTTGTCTACAAACACAACAAGCTTGCCCTGACCCGCACAGGc                                                                                                                   >  4:312362/1‑146 (MQ=255)
                                        aCGTTGTGCCAGTTTCAAACAGGCTGCCCGAGGGCTCAATGTGCCAGTGGTAGGCCTGCGGAAACGCTTGGAGAAACTGGAAGAGCACATTGGTGCTCCCGTGTTTGTCTACAAACACAACAAGCTTGCCCTGACCCGCACAGGCgaga                                                                                                               >  4:198705/1‑149 (MQ=255)
                                        aCGTTGTGCCAGTTTCAAACAGGCTGCCCGAGGGCTCAATGTGCCAGTGGTAGGCCTGCGGAAACGCTTGGAGAAACTGGAAGAGCACATTGGTGCTACCGTGTTTTTCTACAAACACAACAAGCTTGCCCTGACCCGCACAGGCga                                                                                                                 >  2:52138/1‑147 (MQ=255)
                                                ccAGTTTCAAACAGGCTGCCCGAGGGCTCAATGTGCCAGTGGTAGGCCTGCGGAAACGCTTGGAGAAACTGGAAGAGCACATTGGTGCTCCCGTGTTTGTCTACAAACACAACAAGCTTGCCCTGACCCGCACAGGCGAGAGGGTCAg                                                                                                        <  1:54042/148‑1 (MQ=255)
                                                ccAGTTTCAAACAGGCTGCCCGAGGGCTCAATGTGCCAGTGGTAGGCCTGCGGAAACGCTTGGAGAAACTGGAAGAGCACATTGGTGCTCCCGTGTTTGTCTACAAACACAACAAGCTTGCCCTGACCCGCACAGGCGAGAGGGTCAGc                                                                                                       <  4:387481/149‑1 (MQ=255)
                                                  aGTTTCAAACAGGCTGCCCGAGGGCTCAATGTGCCAGTGGTAGGCCTGCGGAAACGCTTGGAGAAACTGGAAGAGCACATTGGTGCTCCCGTGTTTGTCTACAAACACAACAAGCTTGCCCTGACCCGCACAGGCGAGAGGGTCAGCCa                                                                                                     >  1:12184/1‑149 (MQ=255)
                                                       cAAACAGGCTGCCCGAGGGCTCAATGTGCCAGTGGTAGGCCTGCGGAAACGCTTGGAGAAACTGGAAGAGCACATTGGTGCTCCCGTGTTTGTCTACAAACACAACAAGCTTGCCCTGACCCGCACAGGCGAGAGGGTCAGCCATTAcc                                                                                                >  7:278513/1‑149 (MQ=255)
                                                       cAAACAGGCTGCCCGAGGGCTCAATGTGCCAGTGGTAGGCCTGCGGAAACGCTTGGAGAAACTGGAAGAGCACATTGGTGCTCCCGTGTTTGTCTACAAACACAACAAGCTTGCCATGACCCGCACAGGCGAGAGGGTCAGCCATTAcc                                                                                                >  5:257190/1‑149 (MQ=255)
                                                        aaaCAGGCTGCCCGAGGGCTCAATGTGCCAGTGGTAGGCCTGCGGAAACGCTTGGAGAAACTGGAAGAGCACATTGGTGCTCCCGTGTTTGTCTACAAACACAACAAGCTTGCCCTGACCCGCACAGGCGAGAGGGTCAGCCATTACCt                                                                                               <  5:271118/149‑1 (MQ=255)
                                                         aaCAGGCTGCCCGAGGGCTCAATGTGCCAGTGGTAGGCCTGCGGAAACGCTTGGAGAAACTGGAAGAGCACATTGGTGCTACCGTGTTTGTCTACAAACACAACAAGCTTGCCCTGACCCCCACAGCCgaga                                                                                                               >  2:294612/1‑132 (MQ=255)
                                                          aCAGGCTGCCCGAGGGCTCAATGTGCCAGTGGTAGGCCTGCGGAAACGCTTGGAGAAACTGGAAGAGCACATTGGTGCTCCCGTGTTTGTCTACAAACACAACAAGCTTGCCCTGACCCGCACAGGCGAGAGGGTCAGCCATTACCtgt                                                                                             <  3:312362/149‑1 (MQ=255)
                                                          aCAGGCTGCCCGAGGGCTCAATGTGCCAGTGGTAGGCCTGCGGAAACGCTTGGAGAAACTGGAAGAGCACATTGGTGCTCCCGTGTTTGTCTACAAACACAACAAGCTTGCCCTGACCCGCACAGGCGAGAGGGTCAGCCATTACCtgt                                                                                             <  8:222333/149‑1 (MQ=255)
                                                            aGGCTGCCCGAGGGCTCAATGTGCCAGTGGTAGGCCTGCGGAAACGCTTGGAGAAACTGGAAGAGCACATTGGTGCTCCCGTGTTTGTCTACAAACACAACAAGCTTGCCCTGACCCGCACAGGCGAGAGGGTCAGCCATTACCTGTGc                                                                                           <  8:87081/149‑1 (MQ=255)
                                                            aGGCTGCCCGAGGGCTCAATGTGCCAGTGGTAGGCCTGCGGAAACGCTTGGAGAAACTGGAAGAGCACATTGGTGCTCCCGTGTTTGTCTACAAACACAACAAGCTTGCCCTGACCCGCACAGGCGAGAGGGTCAGCCATTACCTGTGc                                                                                           <  6:353767/149‑1 (MQ=255)
                                                              ggCTGCCCGAGGGCTCAATGTGCCAGTGGTAGGCCTGCGGAAACGCTTGGAGAAACTGGAAGAGCACATTGGTGCTCCCGTGTTTGTCTACAAACACAACAAGCTTGCCCTGACCCGCACAGGCGAGAGGGTCAGCCATTACCTGTGcc                                                                                          <  4:136252/149‑1 (MQ=255)
                                                               gCTGCCCGAGGGCTCAATGTGCCAGTGGTAGGCCTGCGGAAACGCTTGGAGAAACTGGAAGAGCACATTGGTGCTCCCGTGTTTGTCTACAAACACAACAAGCTTGCCCTGACCCGCACAGGCGAGAGGGTCAGCCATTACCTGTGCCa                                                                                         <  5:384243/149‑1 (MQ=255)
                                                               gCTGCCCGAGGGCTCAATGTGCCAGTGGTAGGCCTGCGGAAACGCTTGGAGAAACTGGAAGAGCACATTGGTGCTCCCGTGTTTGTCTACAAACACAACAAGCTTGCCCTGACCCGCACAGGCGAGAGGGTCAGCCATTACCTGTGCCa                                                                                         <  1:352688/149‑1 (MQ=255)
                                                                  gCCCGAGGGCTCAATGTGCCAGTGGTAGGCCTGCGGAAACGCTTGGAGAAACTGGAAGAGCACATTGGTGCTCCCATGTTTGTCTACAAACACAACAAGCTTGCCCTGACCCGCACAGGCGAGAGGGTCAGCCATTACCTGTGc                                                                                           >  6:48985/1‑144 (MQ=255)
                                                                  gCCCGAGGGCTCAATGTGCCAGTGGTAGGCCTGCGGAAACGCTTGGAGAAACTGGAAGAGCACATTGGTGCTCCCATGTTTGTCTACAAACACAACAAGCTTGCCCTGACCCGCACAGGCGAGAGGGTCAGCCATTACCTGTGc                                                                                           <  5:48985/144‑1 (MQ=255)
                                                                   cccGAGGGCTCAATGTGCCAGTGGTAGGCCTGCGGAAACGCTTGGAGAAACTGGAAGAGCACATTGGTGCTCCCGTGTTTGTCTACAAACACAACAAGCTTGCCCTGACCCGCACAGGCGAGAGGGTCAGCCATTACCTGTGCCAGctc                                                                                     >  2:55492/1‑149 (MQ=255)
                                                                    ccGAGGGCTCAATGTGCCAGTGGTAGGCCTGCGGAAACGCTTGGAGAAACTGGAAGAGCACATTGGTGCTCCCGTGTTTGTCTACAAACACAACAAGCTTGCCCTGACCCGCACAGGCGAGAGGGTCAGCCATTACCTGTGCCAGctct                                                                                    >  5:198261/1‑149 (MQ=255)
                                                                     cGAGGGCTCAATGTGCCAGTGGTAGGCCTGCGGAAACGCTTGGAGAAACTGGAAGAGCACATTGGTGCTCCCGTGTTTGTCTACAAACACAACAAGCTTGCCCTGACCCGCACAGGCGAGAGGGTCAGCCATTACCTGTGCCAGCTCtt                                                                                   >  1:326891/1‑149 (MQ=255)
                                                                     cGAGGGCTCAATGTGCCAGTGGTAGGCCTGCGGAAAAGCTTGGAGATACTGTAAGAGCACATTGGTGCTCCCGTGTTTGGCTACAAACACAACAAGCTTGCCCTGACCCGCACAGGCGAGAGTGTCAGCCATTACCTGTGCCAGCTCtt                                                                                   <  4:334282/149‑1 (MQ=255)
                                                                           cTCAATGTGCGAGTGGTAGGCCTGCGGAAACTCTTGGAGAAACTGGAAGAGCACATTGGTGCTCCCGTGTTTGTCTACAAACACAACAAGCTTGCCCTGACCCGCACAGGCGAGAGGGTCAGcc                                                                                                      <  1:345806/124‑1 (MQ=255)
                                                                           cTCAATGTGCCAGTGGTAGGCCTGCGGAAACGCTTGGAGAAACTGGAAGAGCACATTGGTGCTCCCGTGTTTGTCTACAAACACAACAAGCTTGCCCTGACCCGCACAGGCGAGAGGGTCAGCCATTACCTGTGCCAGCTCTTCGGcc                                                                              >  7:49340/1‑148 (MQ=255)
                                                                              aaTGTGCCAGTGGTAGGCCTGCGGAAACGCTTGGAGAAAAATGGAAGAGCACATTGGTGCTCCCGTGTTTGTCTCCAAACACAACAAGCTTGCCCTGACCCGCACAGGCGAGAGGGTCAGCCATTACCTGTGCCAGCTCTTCGGcc                                                                              >  1:267435/1‑146 (MQ=255)
                                                                                tgtgCCAGTGGTAGGCCTGCGGAAACGCTTGGAGAAACTGGAAGAGCACATTGGTGCTCCCGTGTTTGTCTACAAACACAACAAGCTTGCCCTGACCCGCACAGGCGAGAGGGTCAGCCATTACCTGTGCCAGCTCTTCGGCCCAGACg                                                                        >  6:33909/1‑149 (MQ=255)
                                                                                    ccAGTGGTAGGCCTGCGGAAACGCTTGGAGAAACTGGAAGAGCACATTGGTGCTCCCGTGTTTGTCTACAAACACAACAAGCTTGCCCTGACCCGCACAGGCGAGAGGGTCAGCCATTACCTGTGCCAGCTCTTCGGCCCAGACGGCCt                                                                    <  2:218228/149‑1 (MQ=255)
                                                                                     cAGTGGTAGGCCTGCGGAAACGCTTGGAGAAACTGGAAGAGCACATTGGTGCTCCCGTGTTTGTCTACAAACACAACAAGCTTGCCCTGACCCGCACAGGCGAGAGGGTCAGCCATTACCTGTGCCAGCTCTTCGGCCCAGACGGCCtt                                                                   <  5:189931/149‑1 (MQ=255)
                                                                                           tAGGCCTGCGTAAACGCTTGTAGAAACTGTAAGAGCAAATTGGTGCTCCCGTGTTTGTCTACAAACAAAACAAGCTTGCCCTGACCCGCACAGGCGATAGTGTCAGCCATTACCTGTGCCAGCTCTTCGGCCCAGACGGCCTTGGCAAg                                                             <  1:90958/149‑1 (MQ=255)
                                                                                             ggCCTGCGGAAACGCTTGGAGAAACTGGAAGAGCACATTGGTGCTCCCGTGTTTGTCTACAAACACAACAAGCTTGCCCTGACCCGCACAGGCGAGAGGGTCAGCCATTACCTGTGCCAGCTCTTCGGCCCAGACGGCCTTGGCAAGt                                                            <  7:263333/148‑1 (MQ=255)
                                                                                              gCCTGCGGAAACGCTTGGAGAAACTGGAAGAGCACATTGGTGCTCCCGTGTTTGTCTACAAACACAACAAGCTTGCCCTGACCCGCACAGGCGAGAGGGTCAGCCATTACCTGTGCCAGCTCTTCGGCCCAGACGGCCTTGGCAAGTcc                                                          >  3:226265/1‑149 (MQ=255)
                                                                                              gCCTGCGGAAACGCTTGGAGAAACTGGAAGAGCACATTGGTGCTCCCGTGTTTGTCTACAAACACAACAAGCTTGCCCTGACCCGCACAGGCGAGAGGGTCAGCCATTACCTGTGCCAGCTCTTCGGCCCAGACGGCCTTGGCAAGTcc                                                          >  3:305837/1‑149 (MQ=255)
                                                                                                cTGCGGAAACGCTTGGAGAAACTGGAAGAGCACATTGGTGCTCCCGTGTTTGTCTACAAACACAACAAGCTTGCCCTGACCCGCACAGGCga                                                                                                                 <  4:91313/92‑1 (MQ=255)
                                                                                                cTGCGGAAACGCTTGGAGAAACTGGAAGAGCACATTGGTGCTCCCGTGTTTGTCTACAAACACAACAAGCTTGCCCTGACCCGCACAGGCga                                                                                                                 >  3:91313/1‑92 (MQ=255)
                                                                                                cTGCGGAAACGCTTGGAGAAACTGGAAGAGCACATTGGTGCTCCCGTGTTTGTCTACAAACACAACAAGCTTGCCCTGACCCGCACAGGCGAGAGGGTCAGCCATTACCTGTGCCAg                                                                                        >  3:364705/1‑117 (MQ=255)
                                                                                                cTGCGGAAACGCTTGGAGAAACTGGAAGAGCACATTGGTGCTCCCGTGTTTGTCTACAAACACAACAAGCTTGCCCTGACCCGCACAGGCGAGAGGGTCAGCCATTACCTGTGCCAg                                                                                        <  4:364705/117‑1 (MQ=255)
                                                                                                cTGCGGAAACGCTTGGAGAAACTGGAAGAGCACATTGGTGCTCCCGTGTTTGTCTACAAACACAACAAGCTTGCCCTGACCCGCACAGGCGAGAGGGTCAGCCATTACCTGTGCCAGCTCTTCGGCCCAGACGGcc                                                                     <  5:240781/136‑1 (MQ=255)
                                                                                                cTGCGGAAACGCTTGGAGAAACTGGAAGAGCACATTGGTGCTCCCGTGTTTGTCTACAAACACAACAAGCTTGCCCTGACCCGCACAGGCGAGAGGGTCAGCCATTACCTGTGCCAGCTCTTCGGCCCAGACGGcc                                                                     >  6:240781/1‑136 (MQ=255)
                                                                                                cTGCGGAAACGCTTGGAGAAACTGGAAGAGCACATTGGTGCTCCCGTGTTTGTCTACAAACACAACAAGCTTGCCCTGACCCGCACAGGCGAGAGGGTCAGCCATTACCTGTGCCAGCTCTTCGGCCCAGACGGCCTTGGCAAGTCCgg                                                        >  1:233645/1‑149 (MQ=255)
                                                                                                      cgattcTTGGAGAAACTGGAAGAGCACATTGGTGCTCCCGTGTTTGTCTACAAACACAACAAGCTTGCCCTGACCCGCACAGGCGAGAGGGTCAGCCATTACCTGTGCCAGCTCTTCGGCCCAGACGGCCTTGGCAAGTCCGGTGaaaa                                                  <  7:284595/144‑1 (MQ=255)
                                                                                                      aaaCGCTTGGAGAAACTGGAAGAGCACATTGGTGCTCCCGTGTTTGTCTACAAACACAACAAGCTTGCCCTGACCCGCACAGGCGAGAGGGTCAGCCATTACCTGTGCCAGCTCTTCGGCCCAGACGGCCTTGGCAAGTCCGGTGaaaa                                                  <  1:55492/149‑1 (MQ=255)
                                                                                                       aaCGCTTGGAGAAACTGGAAGAGCACATTGGTGCTCCCGTGTTTGTCTACAAACACAACAAGCTTGCCCTGACCCGCACAGGCGAGAGGGTCAGCCATTACCTGTGCCAGCTCTTCGGCCCAGACGGCCTTGGCAAGTCCGGTGAAAAg                                                 <  8:49340/149‑1 (MQ=255)
                                                                                                              ggAGAAACTGGAAGAGCACATTGGTGCTCCCGTGTTTGTCTACAAACACAACAAGCTTGCCCTGACCCGCACAGGCGAGAGGGTCAGCCATTACCTGTGCCAGCTCTTCGGCCCAGACGGCCTTGGCAAGTCCGGTGAAAAGGAAGt                                           <  2:12184/147‑1 (MQ=255)
                                                                                                                 gaAACTGGAAGAGCACATTGGTGCTCCCGTGTTTGTCTACAAACACAACAAGCTTGCCCTGACCCGCACAGGCGAGAGGGTCAGCCATTACCTGTGCCAGCTCTTCGGCCCAGACGGCCTTGGCAAGTCCGGTGAAAAGGAAGTCCCCa                                      <  4:305837/149‑1 (MQ=255)
                                                                                                                    aaCTGGAAGAGCACATTGGTGCTCCCGTGTTTGTCTACAAACACAACAAGCTTGCCCTGACCCGCACAGGCGAGAGGGTCAGCCATTACCTGTGCCAGCTCTTCGGCCCAGACGGCCTTGGCAAGTCCGGTGAAAAGGAAGTCCCCAgg                                    <  2:326891/149‑1 (MQ=255)
                                                                                                                           agagCACATTGGTGCTCCCGTGTTTGTCTACAAACACAACAAGCTTGCCCTGACCCGCACAGGCGAGAGGGTCAGCCATTACCTGTGCCAGCTCTTCGGCCCAGACGGCCTTGGCAAGTCCGGTGAAAAGGAAGTCCCCAGGCTTACCa                             >  5:126997/1‑149 (MQ=255)
                                                                                                                            gagCACATTGGTGCTCCCGTGTTTGTCTACAAACACAACAAGCTTGCCCTGACCCGCACAGGCGAGAGGGTCAGCCATTACCTGTGCCAGCTCTTCGGCCCAGACGGCCTTGGCAAGTCCGGTGAAAAGGAAGTCCCCAGGCTTACCAt                            >  3:171504/1‑149 (MQ=255)
                                                                                                                             agCACATTGGTGCTCCCGTGTTTGTCTACAAACACAACAAGCTTGCCCTGACCCGCACAGGCGAGAGGGTCAGCCATTACCTGTGCCAGCTCTTCGGCCCAGACGGCCTTGGCAAGTCCGGTGAAAAGGAAGTCCCCAGGCTTACCATc                           <  2:273884/149‑1 (MQ=255)
                                                                                                                              gCACATTGGTGCTCCCGTGTTTGTCTACAAACACAACAAGCTTGCCCTGACCCGCACAGGCGAGAGGGTCAGCCATTACCTGTGCCAGCTCTTCGGCCCAGACGGCCTTGGCAAGTCCGGTGAAAAGGAAGTCCCCAGGCTTACCATcg                          >  8:28538/1‑149 (MQ=255)
                                                                                                                               cacaTTGGTGCTCCCGTGTTTGTCTACAAACACAACAAGCTTGCCCTGACCCGCACAGGCGAGAGGGTCAGCCATTACCTGTGCCAGCTCTTCGGCCCAGACGGCCTTGGCAAGTCCGGTGAAAAGGAAGTCCCCAGGCTTACCATcgc                         >  7:138299/1‑149 (MQ=255)
                                                                                                                                   ttGGTGCTCCCGTGTTTGTCTACAAACACAACAAGCTTGCCCTGACCCGCACAGGCGAGAGGGTCAGCCATTACCTGTGCCAGCTCTTCGGCCCAGACGGCCTTGGCAAGTCCGGTGAAAAGGAAGTCCCCAGGCTTACCATCGCGATc                     <  7:28538/149‑1 (MQ=255)
                                                                                                                                     ggTGCTCCCGTGTTTGTCTACAAACACAACAAGCTTGCACTGACCCGCACAGGCGAGAGGTTCAGCCATTACCTGTGCCAGATCTTCGGCCCAGACGGACTAGGCAAGACCGGAGAAAAGGAAGTCCCCAGGCGTACAATcgataaca                    >  5:32537/1‑142 (MQ=255)
                                                                                                                                         cTCCCGTGTTTGTCTACAAACACAACAAGCTTGCCCTGACCCGCACAGGCGAGAGGGTCAGCCATTACCTGTGCCAGCTCTTCGGGCCAGACGGCCTTGGCAAGTCCGGTGAAATGGAAGTCCCCAGGCTTACCATCGCGATCACaga                <  4:226265/148‑1 (MQ=255)
                                                                                                                                           gttgtgtTTGTAGACAAACACAACAAGCTTGCCCTGACCCGCACAGGCGAGAGGGTCAGCCATTACCTGTGCCAGCTCTTCGGCCCAGACGGCCTTGGCAAGTCCGGTGAAAAGGAAGTCCCCAGGCTTACCATCGCGATCACAGAGc              >  1:49516/4‑148 (MQ=255)
                                                                                                                                           gttgtgtTTGTAGACAAACACAACAAGCTTGCCCTGACCCGCACACGCGAGAGGGTCAGCAATTAACAGTGACACATCTTCGGCCCAGACGGCCTTGCCAAGGAAGAAGAAAAGGAAGTCACCGGGCATACGATCGCGATCAgagaca              >  8:315628/4‑146 (MQ=255)
                                                                                                                                              gtgtTTGTCTACAAACACAACAAGCTTGCCCTGACCCGCACAGGCGAGAGGGTCAGCCATTACCTGTGCCAGCTCTTCGGCCCAGACGGCCTTGGCAAGTCCGGTGAAAAGGAATGTCCCCAGGCTTACCATCGCGATCACAGAGCCAg           <  1:52138/149‑1 (MQ=255)
                                                                                                                                                 tttGTCTACAAACACAACAAGCTTGCCCTGACCCGCACAGGCGAGAGGGTCAGCCATTACCTGTGCCAGCTCTTCGGCCCAGACGGCCTTGGCAAGTCCGGTGAAAAGGAAGTCCCCAGGCTTACCATCGCGATCACAGAGCCAGTGc        <  8:278513/148‑1 (MQ=255)
                                                                                                                                                     tCTACAAACACAACAAGCTTGCCCTGACCCGCACAGGCGAGAGGGTCAGCAATTACCTGTGCCAGCTCTTCGGCCCAGACGGACTTGGCAAGTACGGTGAAAAGGAAGTCCCCAGGCTTACAATCGAGAt                      >  3:227087/1‑130 (MQ=255)
                                                                                                                                                      cTACAAACACAACAAGCTTGCCCTGACCCGCACAGGCGAGAGGGTCAGCCATTACCTGTGCCAGCTCTTCGGCCCAGACGGCCTTGGCAAGTCCGGTGAAAAGGAAGTCCCCAGGCTTACCATCGCGATCACAGAGCCAGTGCTGAACg  >  8:193001/1‑149 (MQ=255)
                                                                                                                                                      |                                                                                                                                                     
CAATTGATGATCAGCTGGCCGTAGAGTTTTTCGCAACAGCACGTTGTGCCAGTTTCAAACAGGCTGCCCGAGGGCTCAATGTGCCAGTGGTAGGCCTGCGGAAACGCTTGGAGAAACTGGAAGAGCACATTGGTGCTCCCGTGTTTGTCTACAAACACAACAAGCTTGCCCTGACCCGCACAGGCGAGAGGGTCAGCCATTACCTGTGCCAGCTCTTCGGCCCAGACGGCCTTGGCAAGTCCGGTGAAAAGGAAGTCCCCAGGCTTACCATCGCGATCACAGAGCCAGTGCTGAACG  >  NC_002947/39342‑39638

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