Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,070,302 C→G 7.4% G325R (GGT→CGT) ‡ PP_4465 ← porin

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,070,3020CG 7.4% 644.8 / 30.7 199G325R (GGT→CGT) ‡PP_4465porin
Reads supporting (aligned to +/- strand):  ref base C (85/98);  new base G (5/10);  total (90/109)
Fisher's exact test for biased strand distribution p-value = 4.22e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 224 total aligned reads displayed.

CCCCGCCTTTAGGCCTGGAACACCTAAGGCGGCGAAATCGTAGCTGTAAAGGCCAAAAGTAGAATTTTCGCCCGCGCGCGTGAAGCTATTCACCATAACGTCGGTATATAGATATAGGGATTGGACGCCCCCTTCACCTTCTGGTCGACCTCGACCATCAACCACGTTGCTACTATTGACGGTTGCAAAGCCAGACTTATCACTGACTTGTTGATGCCCTAAAGTCAGCGCATGACCACCCGCTTTATAGGTGAAAGCCAAAGACCAGGTGGTATTATCTATTTTCCCATTTTCATTGT  >  NC_002947/5070153‑5070451
                                                                                                                                                     |                                                                                                                                                     
ccccGCCTTTAGGCCTGGAACACCTAAGGCGGCGAAATCGTAGCTGTAAAGGCCAAAAGTAGAATTTTCGCCCGCGCGCGTGAAGCTATTCACCATAACGTCGGTATATAGATATAGGGATTGGACGCCCCCTTCACCTTCTGGtcgacc                                                                                                                                                       >  2:1447704/1‑150 (MQ=255)
 cccGCCTTTAGGCCTGGAACACCTAAGGCGGCGAAATCGTAGCTGTAAAGGCCAAAAGTAGAATTTTCGCCCGCGCGCGTGAAGCTATTCACCATAACGTCGGTATATAGATATAGGGATTGGACGCCCCCTTCACCTTCTGGtcgacct                                                                                                                                                      >  2:523305/1‑150 (MQ=255)
 cccGCCTTTAGGCCTGGAACACCTAAGGCGGCGAAATCGTAGCTGTAAAGGCCAAAAGTAGAATTTTCGCCCGCGCGCGTGAAGCTATTCACCATAACGTCGGTATATAGATATAGGGATTGGACGCCCCCTTCACCTTCTGGtcgacct                                                                                                                                                      >  2:3950924/1‑150 (MQ=255)
 cccGCCTTTAGGCCTGGAACACCTAAGGCGGCGAAATCGTAGCTGTAAAGGCCAAAAGTAGAATTTTCGCCCGCGCGCGTGAAGCTATTCACCATAACGTCGGTATATAGATATAGGGATTGGACGCCCCCTTCACCTTCTGGtcgacct                                                                                                                                                      >  1:8837494/1‑150 (MQ=255)
  ccGCCTTTAGGCCTGGAACACCTAAGGCGGCGAAATCGTAGCTGTAAAGGCCAAAAGTAGAATTTTCGCCCGCGCGCGTGAAGCTATTTACCATAACGTCGGTATATAGATATAGGGATTGGACGCCCCCTTCACCTTCTGGtcgacctc                                                                                                                                                     <  1:4524773/150‑1 (MQ=255)
  ccGCCTTTAGGCCTGGAACACCTAAGGCGGCGAAATCGTAGCTGTAAAGGCCAAAAGTAGAATTTTCGCCCGCGCGCGTGAAGCTATTCACCATAACGTCGGTATATAGATATAGGGATTGGACGCCCCCTTCACCTTCTGGtcgacctc                                                                                                                                                     >  1:4064817/1‑150 (MQ=255)
  ccGCCTTTAGGCCTGGAACACCTAAGGCGGCGAAATCGTAGCTGTAAAGGCCAAAAGTAGAATTTTCGCCCGCGCGCGTGAAGCTATTCACCATAACGTCGGTATATAGATATAGGGATTGGACGCCCCCTTCACCTTCTGGtcgacctc                                                                                                                                                     >  2:4024089/1‑150 (MQ=255)
    gCCTTTAGGCCTGGAACACCTAAGGCGGCGAAATCGTAGCTGTAAAGGCCAAAAGTAGAATTTTCGCCCGCGCGCGTGAAGCTATTCACCATAACGTCGGTATATAGATATAGGGATTGGACGCCCCCTTCACCTTCTGGtcgacctcga                                                                                                                                                   <  1:7837505/150‑1 (MQ=255)
    gCCTTTAGGCCTGGAACACCTAAGGCGGCGAAATCGTAGCTGTAAAGGCCAAAAGTAGAATTTTCGCCCGCGCGCGTGAAGCTATTCACCATAACGTCGGTATATAGATATAGGGATTGGACGCCCCCTTCACCTTCTGGtcgacctcga                                                                                                                                                   >  1:337525/1‑150 (MQ=255)
     ccTTTAGGCCTGGAACACCTAAGGCGGCGAAATCGTAGCTGTAAAGGCCAAAAGTAGAATTTTCGCCCGCGCGCGTGAAGCTATTCACCATAACGTCGGTATATAGATATAGGGATTGGACGCCCCCTTCACCTTCTGGtcgacctcgac                                                                                                                                                  <  2:6258031/150‑1 (MQ=255)
     ccTTTAGGCCTGGAACACCTAAGGCGGCGAAATCGTAGCTGTAAAGGCCAAAAGTAGAATTTTCGCCCGCGCGCGTGAAGCTATTCACCATAACGTCGGTATATAGATATAGGGATTGGACGCCCCCTTCACCTTCTGGtcgacctcgac                                                                                                                                                  <  2:7966049/150‑1 (MQ=255)
       tttAGGCCTGGAACACCTAAGGCGGCGAAATCGTAGCTGTAAAGGCCAAAAGTAGAATTTTCGCCCGCGCGCGTGAAGCTATTCACCATAACGTCGGTATATAGATATAGGGATTGGACGCCCCCTTCACCTTCTGGTCGACCTCGACCa                                                                                                                                                <  2:4673264/150‑1 (MQ=255)
       tttAGGCCTGGAACACCTAAGGCGGCGAAATCGTAGCTGTAAAGGCCAAAAGTAGAATTTTCGCCCGCGCGCGTGAAGCTATTCACCATAACGTCGGTATATAGATATAGGGATTGGACGCCCCCTTCACCTTCTGGTCGACCTCGACCa                                                                                                                                                >  1:2055029/1‑150 (MQ=255)
         tAGGCCTGGAACACCTAAGGCGGCGAAATCGTAGCTGTAAAGGCCAAAAGTAGAATTTTCGCCCGCGCGCGTGAAGCTATTCACCATAACGTCGGTATATAGATATAGGGATTGGACGCCCCCTTCACCTTCTGGTCGACCTCGACCATc                                                                                                                                              <  2:3908028/150‑1 (MQ=255)
          aGGCCTGGAACACCTACGGCGGCGAAATCGTAGCTGTAAAGGCCAAAAGTAGAATTTTCGCCCGCGCGCGTGAAGCTCTTCACCATAACGTCGGTATATAGATATAGGGATTGGACGCCCCCTTCACCTTCTGGTCGAGGTCGACCagaa                                                                                                                                             <  2:3006654/150‑4 (MQ=255)
          aGGCCTGGAACACCTAAGGCGGCGAAATCGTAGCTGTAAAGGCCAAAAGTAGAATTTTCGCCCGCGCGCGTGAAGCTATTCACCATAACGTCGGTATATAGATATAGGGATTGGACGCCCCCTTCACCTTCTGGTCGAGGTCGACCagaa                                                                                                                                             <  2:7588005/150‑4 (MQ=255)
          aGGCCTGGAACACCTAAGGCGGCGAAATCGTAGCTGTAAAGGCCAAAAGTAGAATTTTCGCCCGCGCGCGTGAAGCTATTCACCATAACGTCGGTATATAGATAGAGGGATTGGACGCCCCCTTCACCTTCTGGTCGAGGTCGACCagaa                                                                                                                                             <  1:8711157/150‑4 (MQ=255)
           ggCCTGGAACACCTAAGGCGGCGAAATCGTAGCTGTAAAGGCCAAAAGTAGAATTTTCGCCCGCGCGCGTGAAGCTATTCACCATAACGTCGGTATATAGATATAGGGATTGGACGCCCCCTTCACCTTCTGGTCGAGGTCGACCagaag                                                                                                                                            <  1:406516/150‑5 (MQ=255)
           ggCCTGGAACACCTAAGGCGGCGAAATCGTAGCTGTAAAGGCCAAAAGTAGAATTTTCGCCCGCGCGCGTGAAGCTATTCACCATAACGTCGGTATATAGATATAGGGATTGGACGCCCCCTTCACCTTCTGGTCGAGGTCGACCagaag                                                                                                                                            <  1:2697930/150‑5 (MQ=255)
            gCCTGGAACACCTAAGGCGGCGAAATCGTAGCTGTAAAGGCCAAAAGTAGAATTTTCGCCCGCGCGCGTGAAGCTATTCACCATAACGTCGGTATATAGATATAGGGATTGGACGCCCCCTTCACCTTCTGGTCGAGGTCGACCagaagg                                                                                                                                           <  1:2791747/150‑6 (MQ=255)
               tGGAACACCTAAGGCGGCGAAATCGTAGCTGTAAAGGCCAAAAGTAGAATTTTCGCCCGCGCGCGTGAAGCTATTCACCATAACGTCGGTATATAGATATAGGGATTGGACGCCCCCTTCACCTTCTGGTCGACCTCGACCATCAACCAc                                                                                                                                        >  1:1987007/1‑150 (MQ=255)
                 gAACACCTAAGGCGGCGAAATCGTAGCTGTAAAGGCCAAAAGTAGAATTTTCGCCCGCGCGCGTGAAGCTATTCACCATAACGTCGGTATATAGATATAGGGATTGGACGCCCCCTTCACCTTCTGGTCGACCTCGACCATCAACCACGt                                                                                                                                      <  2:2708236/150‑1 (MQ=255)
                  aaCACCTAAGGCGGCGAAATCGTAGCTGTAAAGGCCAAAAGTAGAATTTTCGCCCGCGCGCGTGAAGCTATTCACCATAACGTCGGTATATAGATATAGGGATTGGACGCCCCCTTCACCTTCTGGTCGACCTCGACCATCAACCACGtt                                                                                                                                     <  2:1036599/150‑1 (MQ=255)
                    cacCTAAGGCGGCGAAATCGTAGCTGTAAAGGCCAAAAGTAGAATTTTCGCCCGCGCGCGTGAAGCTATTCACCATAACGTCGGTATATAGATATAGGGATTGGACGCCCCCTTCACCTTCTGGTCGACCTCGACCATCAACCACGTTGc                                                                                                                                   <  2:7080368/150‑1 (MQ=255)
                      ccTAAGGCGGCGAAATCGTAGCTGTAAAGGCCAAAAGTAGAATTTTCGCCCGCGCGCGTGAAGCTATTCACCATAACGTCGGTATATAGATATAGGGATTGGACGCCCCCTTCACCTTCTGGTCGACCTCGACCATCAACCACGTTGcta                                                                                                                                 <  2:1812947/150‑1 (MQ=255)
                       cTAAGGCGGCGAAATCGTAGCTGTAAAGGCCAAAAGTAGAATTTTCGCCCGCGCGCGTGAAGCTATTCACCATAACGTCGGTATATAGATATAGGGATTGGACGCCCCCTTCACCTTCTGGTCGACCTCGACCATCAACCACGTTGctac                                                                                                                                >  1:6576953/1‑150 (MQ=255)
                       cTAAGGCGGCGAAATCGTAGCTGTAAAGGCCAAAAGTAGAATTTTCGCCCGCGCGCGTGAAGCTATTCACCATAACGTCGGTATATAGATATAGGGATTGGACGCCCCCTTCACCTTCTGGTCGACCTCGACCATCAACCACGTTGctac                                                                                                                                >  2:6584535/1‑150 (MQ=255)
                        tAAGGCGGCGAAATCGTAGCTGTAAAGGCCAAAAGTAGAATTTTCGCCCGCGCGCGTGAAGCTATTCACCATAACGTCGGTATATAGATATAGGGATTGGACGCCCCCTTCACCTTCTGGTCGACCTCGACCATCAACCACGTTGctact                                                                                                                               >  1:9132232/1‑150 (MQ=255)
                                 gAAATCGTAGCTGTAAAGGCCAAAAGTAGAATTTTCGCCCGCGCGCGTGAAGCTATTCACCATAACGTCGGTATATAGATATAGGGATTGGACGCCCCCTTCACCTTCTGGTCGACCTCGACCATCAACCACGTTGCTACTATTGACGGt                                                                                                                      <  1:1606454/150‑1 (MQ=255)
                                  aaaTCGTAGCTGTAAAGGCCAAAAGTAGAATTTTCGCCCGCGCGCGTGAAGCTATTCACCATAACGTCGGTATATAGGTATAGGGATTGGACGCCCCCTTCACCTTCTGGTCGACCTCGACCATCAACCACGTTGCTACTATTGACGGtt                                                                                                                     <  2:2827324/150‑1 (MQ=255)
                                  aaaTCGTAGCTGTAAAGGCCAAAAGTAGAATTTTCGCCCGCGCGCGTGAAGCTATTCAACATAACGTCGGTATATAGATATAGGGATTGGACGCCCCCTTCACCTTCTGGTCGACCTCGACCATCAACCACGTTGCTACTATTGACGGtt                                                                                                                     <  1:3786552/150‑1 (MQ=255)
                                    aTCGTAGCTGTAAAGGCCAAAAGTAGAATTTTCGCCCGCGCGCGTGAAGCTATTCACCATAACGTCGGTATATAGATATAGGGATTGGACGCCCCCTTCACCTTCTGGTCGACCTCGACCATCAACCACGTTGCTACTATTGACGGTTGc                                                                                                                   <  1:2915004/150‑1 (MQ=255)
                                       gTAGCTGTAAAGGCCAAAAGTAGAATTTTCGCCCGCGCGCGTGAAGCTATTCACCATAACGTCGGTATATAGATATAGGGATTGGACGCCCCCTTCACCTTCTGGTCGACCTCGACCATCAACCACGTTGCTACTATTGACGGTTGCaaa                                                                                                                <  2:8873419/150‑1 (MQ=255)
                                       gTAGCTGTAAAGGCCAAAAGTAGAATTTTCGCCCGCGCGCGTGAAGCTATTCACCATAACGTCGGTATATAGATATAGGGATTGGACGCCCCCTTCACCTTCTGGTCGACCTCGACCATCAACCACGTTGCTACTATTGACGGTTGCaaa                                                                                                                >  2:2946508/1‑150 (MQ=255)
                                       gTAGCTGTAAAGGCCAAAAGTAGAATTTTCGCCCGCGCGCGTGAAGCTATTCACCATAACGTCGGTATATAGATATAGGGATTGGACGCCCCCTTCACCTTCTGGTCGACCTCGACCATCAACCACGTTGCTACTATTGACGGTTGCaaa                                                                                                                <  2:2762139/150‑1 (MQ=255)
                                        tAGCTGTAAAGGCCAAAAGTAGAATTTTCGCCCGCGCGCGTGAAGCTATTCACCATAACGTCGGTATATAGATATAGGGATTGGACGCCCCCTTCACCTTCTGGTCGACCTCGACCATCAACCACGTTGCTACTATTGACGGTTGCAAAg                                                                                                               <  2:2524362/150‑1 (MQ=255)
                                             gTAAAGGCCAAAAGTAGAATTTTCGCCCGCGCGCGTGAAGCTATTCACCATAACGTCGGTATATAGATATAGGGATTGGACGCCCCCTTCACCTTCTGGTCGACTTCGACCATCAACCACGTTGCTACTATTGACGGTTGCAAAGCCAGa                                                                                                          <  2:2986516/150‑1 (MQ=255)
                                             gTAAAGGCCAAAAGTAGAATTTTCGCCCGCGCGCGTGAAGCTATTCACCATAACGTCGGTATATAGATATAGGGATTGGACGCCCCCTTCACCTTCTGGTCGACCTCGACCATCAACCACGTTGCTACTATTGACGGTTGCAAAGCCAGa                                                                                                          <  2:6265684/150‑1 (MQ=255)
                                             gTAAAGGCCAAAAGTAGAATTTTCGCCCGCGCGCGTGAAGCTATTCACCATAACGTCGGTATATAGATATAGGGATTGGACGCCCCCTTCACCTTCTGGTCGACCTCGACCATCAACCACGTTGCTACTATTGACGGTTGCAAAGCCAGa                                                                                                          <  2:9016500/150‑1 (MQ=255)
                                                     cAAAAGTAGAATTTTCGCCCGCGCGCGTGAAGCTATTCACCATAACGTCGGTATATAGATATAGGGATTGGACGCCCCCTTCACCTTCTGGTCGACCTCGACCATCAACCACGTTGCTACTATTGACGGTTGCAAAGCCAGACTTATCAc                                                                                                  >  2:7193856/1‑150 (MQ=255)
                                                      aaaaGTAGAATTTTCGCCCGCGCGCGTGAAGCTATTCACCATAACGTCGGTATATAGATATAGGGATTGGACGCCCCCTTCACCTTCTGGTCGACCTCGACCATCAACCACGTTGCTACTATTGACGGTTGCAAAGCCAGACTTATCACt                                                                                                 <  1:645842/150‑1 (MQ=255)
                                                           tAGAATTTTCGCCCGCGCGCGTGAAGCTATTCACCATAACGTCGGTATATAGATATAGGGATTGGACGCCCCCTTCACCTTCTGGTCGACCTCGACCATCAACCACGTTGCTACTATTGACGGTTGCAAAGCCAGACTTATCACTGACtt                                                                                            >  1:7673330/1‑150 (MQ=255)
                                                           tAGAATTTTCGCCCGCGCGCGTGAAGCTATTCACCATAACGTCGGTATATAGATATAGGGATTGGACGCCCCCTTCACCTTCTGGTCGACCTCGACCATCAACCACGTTGCTACTATTGACGGTTGCAAAGCCAGACTTATCACTGACtt                                                                                            <  2:5670686/150‑1 (MQ=255)
                                                           tAGAATTTTCGCCCGCGCGCGTGAAGCTATTCACCATAACGTCGGTATATAGATATAGGGATTGGACGCCCCCTTCACCTTCTGGTCGACCTCGACCATCAACCACGTTGCTACTATTGACGGTTGCAAAGCCAGACTTATCACTGACtg                                                                                            >  1:3053767/1‑149 (MQ=255)
                                                             gAATTTTCGCCCGCGCGCGTGAAGCTATTCACCATAACGTCGGTATATAGATATAGGGATTGGACGCCCCCTTCACCTTCTGGTCGACCTCGACCATCAACCACGTTGCTACTATTGACGGTTGCAAAGCCAGACTTATCACTGACttgt                                                                                          >  2:7344759/1‑150 (MQ=255)
                                                             gAATTTTCGCCCGCGCGCGTGAAGCTATTCACCATAACGTCGGTATATAGATATAGGGATTGGACGCCCCCTTCACCTTCTGGTCGACCTCGACCATCAACCACGTTGCTACTATTGACGGTTGCAAAGCCAGACTTATCACAGACttgt                                                                                          >  2:3317571/1‑150 (MQ=255)
                                                              aaTTTTCGCCCGCGCGCGTGAAGCTATTCACCATAACGTCGGTATATAGATATAGGGATTGGACGCCCCCTTCACCTTCTGGTCGACCTCGACCATCAACCACGTTGCTACTATTGACGGTTGCAAAGCCAGACTTATCACTGACttgtt                                                                                         <  2:4062591/150‑1 (MQ=255)
                                                               aTTTTCGCCCGCGCGCGTGAAGCTATTCACCATAACGTCGGTATATAGATATAGGGATTGGACGCCCCCTTCACCTTCTGGTCGACCTCGACCATCAACCACGTTGCTACTATTGACGGTTGCAAAGCCAGACTTATCACTGACttgttg                                                                                        <  1:5038346/150‑1 (MQ=255)
                                                                ttttCGCCCGCGCGCGTGAAGCTATTCACCATAACGTCGGTATATAGATATAGGGATTGGACGCCCCCTTCACCTTCTGGTCGACCTCGACCATCAACCACGTTGCTACTATTGACGGTTGCAAAGCCAGACTTATCACTGACTTGTTGa                                                                                       >  2:5980418/1‑150 (MQ=255)
                                                                ttttCGCCCGCGCGCGTGAAGCTATTCACCATAACGTCGGTATATAGATATAGGGATTGGACGCCCCCTTCACCTTCTGGTCGACCTCGACCATCAACCACGTTGCTACTATTGACGGTTGCAAAGCCAGACTTATCACTGACTTGTTGa                                                                                       <  1:6784229/150‑1 (MQ=255)
                                                                ttttCGCCCGCGCGCGTGAAGCTATTCACCATAACGTCGGTATATAGATATAGGGATTGGACGCCCCCTTCACCTTCTGGTCGACCTCGACCATCAACCACGTTGCTACTATTGACGGTTGCAAAGCCAGACTTATCACTGACTTGTTGa                                                                                       >  2:6978346/1‑150 (MQ=255)
                                                                        cgcgcgcgTGAAGCTATTCACCATAACGTCGGTATATAGATATAGGGATTGGACGCCCCCTTCACCTTCTGGTCGACCTCGACCATCAACCACGTTGCTACTATTGACGGTTGCAAAGCCAGACTTATCACTGACTTGTTGATGCCCTaa                                                                               >  2:5019299/1‑150 (MQ=255)
                                                                         gcgcgcgTGAAGCTATTCACCATAACGTCGGTATATAGATATAGGGATTGGACGCCCCCTTCACCTTCTGGTCGACCTCGACCATCAACCACGTTGCTACTATTGACGGTTGCAAAGCCAGACTTATCACTGACTTGTTGATGCCCTaaa                                                                              <  2:3773401/150‑1 (MQ=255)
                                                                         gcgcgcgTGAAGCTATTCACCATAACGTCGGTATATAGATATAGGGATTGGACGCCCCCTTCACCTTCTGGTCGACCTCGACCATCAACCACGTTGCTACTATTGACGGTTGCAAAGCCAGACTTATCACTGACTTGTTGATGCCCTaaa                                                                              >  1:1007289/1‑150 (MQ=255)
                                                                         gcgcgcgTGAAGCTATTCACCATAACGTCGGTATATAGATATAGGGATTGGACGCCCCCTTCACCTTCTGGTCGACCTCGACCATCAACCACGTTGCTACTATTGACGGTTGCAAAGCCAGACTTATCACTGACTTGTTGATGCCCTaaa                                                                              <  2:1987007/150‑1 (MQ=255)
                                                                         gcgcgcgTGAAGCTATTCACCATAACGTCGGTATATAGATATAGGGATTGGACGCCCCCTTCACCTTCTGGTCGACCTCGACCATCAACCACGTTGCTACTATTGACGGTTGCAAAGCCAGACTTATCACTGACTTGTTGATGCCCTaaa                                                                              <  2:1340002/150‑1 (MQ=255)
                                                                         gcgcgcgTGAAGCTATTCACCATAACGTCGGTATATAGATATAGGGATTGGACGCCCCCTTCACCTTCTGGTCGACCTCGACCATCAACCACGTTGCTACTATTGACGGTTGCAAAGCCAGACTTATCACTGACTTGTTGATGCCCTaaa                                                                              <  2:846294/150‑1 (MQ=255)
                                                                             gcgTGAAGCTATTCACCATAACGTCGGTATATAGATATAGGGATTGGACGCCCCCTTCACCTTCTGGTCGACCTCGACCATCAACCACGTTGCTACTATTGACGGTTGCAAAGCCAGACTTATCACTGACTTGTTGATGCCCTAAAGTCa                                                                          <  2:8971965/150‑1 (MQ=255)
                                                                             gcgTGAAGCTATTCACCATAACGTCGGTATATAGATATAGGGATTGGACGCCCCCTTCACCTTCTGGTCGACCTCGACCATCAACCACGTTGCTACTATTGACGGTTGCAAAGCCAGACTTATCACTGACTTGTTGATGCCCTAAAGTCa                                                                          <  1:6363776/150‑1 (MQ=255)
                                                                              cgTGAAGCTATTCACCATAACGTCGGTATATAGATATAGGGATTGGACGCCCCCTTCACCTTCTGGTCGACCTCGACCATCAACCACGTTGCTACTATTGACGGTTGCAAAGCCAGACTTATCACTGACTTGTTGATGCCCTAAAGTCAg                                                                         <  1:6902750/150‑1 (MQ=255)
                                                                               gTGAAGCTATTCACCATAACGTCGGTATATAGATATAGGGATTGGACGCCCCCTTCACCTTCTGGTCGACCTCGACCATCAACCACGTTGCTACTATTGACGGTTGCAAAGCCAGACTTATCACTGACTTGTTGATGCCCTAAAGTCAgc                                                                        <  1:2529677/150‑1 (MQ=255)
                                                                               gTGAAGCTATTCACCATAACGTCGGTATATAGATATAGGGATTGGACGCCCCCTTCACCTTCTGGTCGACCTCGACCATCAACCACGTTGCTACTATTGACGGTTGCAAAGCCAGACTTATCACTGACTTGTTGATGCCCTAAAGTCAgc                                                                        <  1:9103982/150‑1 (MQ=255)
                                                                                   aGCTATTCACCATAACGTCGGTATATAGATATAGGGATTGGACGCCCCCTTCACCTTCTGGTCGACCTCGACCATCAACCACGTTGCTACTATTGACGGTTGCAAAGCCAGACTTATCACTGACTTGTTGATGCCCTAAAGTCAGCGCAt                                                                    <  1:5410367/150‑1 (MQ=255)
                                                                                    gCTATTCACCATAACGTCGGTATATAGATATAGGGATTGGACGCCCCCTTCACCTTCTGGTCGACCTCGACCATCAACCACGTTGCTACTATTGACGGTTGCAAAGCCAGACTTATCACTGACTTGTTGATGCCCTAAAGTCAGCGCATg                                                                   >  1:8958584/1‑150 (MQ=255)
                                                                                            ccATAACGTCGGTATATAGATATAGGGATTGGACGCCCCCTTCACCTTCTGGTCGACCTCGACCATCAACCACGTTGCTACTATTGACGGTTGCAAAGCCAGACTTATCACTGACTTGTTGATGCCCTAAAGTCAGCGCATGACCACCCg                                                           >  1:1263259/1‑150 (MQ=255)
                                                                                            ccATAACGTCGGTATATAGATATAGGGATTGGACGCCCCCTTCACCTTCTGGTCGACCTCGACCATCAACCACGTTGCTACTATTGACGGTTGCAAAGCCAGACTTATCACTGACTTGTTGATGCCCTAAAGTCAGCGCATGACCACCCg                                                           >  1:5099999/1‑150 (MQ=255)
                                                                                               tAACGTCGGTATATAGATATAGGGATTGGACGCCCCCTTCACCTTCTGGTCGACCTCGACCATCAACCACGTTGCTACTATTGACGGTTGCAAAGCCAGACTTATCACTGACTTGTTGATGCCCTAAAGTCAGCGCATGACCACCCGCtt                                                        <  2:8482927/150‑1 (MQ=255)
                                                                                                aaCGTCGGTATATAGATATAGGGATTGGACGCCCCCTTCACCTTCTGGTCGACCTCGACCATCAACCACGTTGCTACTATTGACGGTTGCAAAGCCAGACTTATCACTGACTTGTTGATGCCCTAAAGTCCGCGCATGACCACCCGCttt                                                       >  2:8771302/1‑150 (MQ=255)
                                                                                                   gTCGGTATATAGATATAGGGATTGGACGCCCCCTTCACCTTCTGGTCGACCTCGACCATCAACCACGTTGCTACTATTGACGGTTGCAAAGCCAGACTTATCACTGACTTGTTGATGCCCTAAAGTCAGCGCATGACCACCCGCTTtata                                                    <  1:4782610/150‑1 (MQ=255)
                                                                                                   gTCGGTATATAGATATAGGGATTGGACGCCCCCTTCACCTTCTGGTCGACCTCGACCATCAACCACGTTGCTACTATTGACGGTTGCAAAGCCAGACTTATCACTGACTTGTTGATGCCCTAAAGTCAGCGCATGACCACCCGCTTtata                                                    <  1:8599656/150‑1 (MQ=255)
                                                                                                     cGGTATATAGATATAGGGATTGGACGCCCCCTTCACCTTCTGGTCGACCTCGACCATCAACCACGTTGCTACTATTGACGGTTGCAAAGCCAGACTTATCACTGACTTGTTGATGCCCTAAAGTCAGCGCATGACCACCCGCTTTATAgg                                                  <  1:1713554/150‑1 (MQ=255)
                                                                                                     cGGTATATAGATATAGGGATTGGACGCCCCCTTCACCTTCTGGTCGACCTCGACCATCAACCACGTTGCTACTATTGACGGTTGCAAAGCCAGACTTATCACTGACTTGTTGATGCCCTAAAGTCAGCGCATGACCACCCGCTTTATAgg                                                  <  2:488645/150‑1 (MQ=255)
                                                                                                            tagatatagGGATTGGACGCCCCCTTCACCTTCTGGTCGACCTCGACCATCAACCACGTTGCTACTATTGACGGTTGCAAAGCCAGACTTATCACTGACTTGTTGATGCCCTAAAGTCAGCGCATGACCACCCGCTTTATAGGTGAAAGc                                           <  1:806320/150‑1 (MQ=255)
                                                                                                             agatatagGGATTGGACGCCCCCTTCACCTTCTGGTCGACGTCGACCATCAACCACGTTGCTACTATTGACGGTTGCAAAGCCAGTCTTATCACTGTGTTGTTGTTGCCGTATAGTCAGCGCATGACCACCCGCTTTATAGGTGTATGcc                                          >  2:4630761/1‑150 (MQ=255)
                                                                                                             agatatagGGATTGGACGCCCCCTTCACCTTCTGGTCGACCTCGACCATCAACCACGTTGCTACTATTGACGGTTGCAAAGCCAGACTTATCACTGACTTGTTGATGCCCTAAAGTCAGCGCATGACCACCCGCTTTATAGGTGAAAGcc                                          <  2:4394952/150‑1 (MQ=255)
                                                                                                                 atagGGATTGGACGCCCACTTCACCTTCTGGTCGACCTCGACCATCAACCACGTTGCTACTATTGACGGTTGCAAAGCCAGACTTATCACTGACTTGTTGATGCCCTAAAGTCAGCGCATGACCACCCGCTTTATAGGTGAAAGCCAAAg                                      >  1:8139956/1‑150 (MQ=255)
                                                                                                                    gggATTGGACGCCCCCTTCACCTTCTGGTCGACCTCGACCATCAACCACGTTGCTACTATTGACGGTTGCAAAGCCAGACTTATCACTGACTTGTTGATGCCCTAAAGTCAGCGCATGACCACCCGCTTTATAGGTGAAAGCCAAAGAcc                                   >  1:754587/1‑150 (MQ=255)
                                                                                                                     ggATTGGACGCCCCCTTCACCTTCTGGTCGACCTCGACCATCAACCACGTTGCTACTATTGACGGTTGCAAAGCCAGACTTATCACTGACTTGTTGATGCCCTAAAGTCAGCGCATGACCACCCGCTTTATAGGTGAAAGCCAAAGACCa                                  >  1:5768845/1‑150 (MQ=255)
                                                                                                                     ggATTGGACGCCCCCTTCACCTTCTGGTCGACCTCGACCATCAACCACGTTGCTACTATTGACGGTTGCAAAGCCAGACTTATCACTGACTTGTTGATGCCCTAAAGTCAGCGCATGACCACCCGCTTTATAGGTGAAAGCCAAAGACCa                                  >  1:2358917/1‑150 (MQ=255)
                                                                                                                       aTTGGACGCCCCCTTCACCTTCTGGTCGACCTCGACAATCAACCACGTTGCTACTATTGACGGTTGCAAAGCCAGACTTATCACTGACTTGTTGATGCCCTAAAGTCAGCGCATGACAACCCGCTTTATAGGTGAAAGCCAAAGACCAgg                                >  1:145751/1‑150 (MQ=255)
                                                                                                                         tGGACGCCCCCTTCACCTTCTGGTCGACCTCGACCATCAACCACGTTGCTACTATTGACGGTTGCAAAGCCAGACTTATCACTGACTTGTTGATGCCCTAAAGTCAGCGCATGACCACCCGCTTTATAGGTGAAAGCCAAAGACCAggtg                              <  1:3590816/150‑1 (MQ=255)
                                                                                                                              gCCCCCTTCACCTTCTGGTCGACCTCGACCATCAACCACGTTGCTACTATTGACGGTTGCAAAGCCAGACTTATCACTGACTTGTTGATGCCCTAAAGTCAGCGCATGACCACCCGCTTTATAGGTGAAAGCCAAAGACCAGGTGGtatt                         >  2:1156461/1‑150 (MQ=255)
                                                                                                                              gCCCCCTTCACCTTCTGGTCGACCTCGACCATCAACCACGTTGCTACTATTGACGGTTGCAAAGCCAGACTTATCACTGACTTGTTGATGCCCTAAAGTCAGCGCATGACCACCCGCTTTATAGGTGAAAGCCAAAGACCAGGTGGtatt                         <  1:3950924/150‑1 (MQ=255)
                                                                                                                              gCCCCCTTCACCTTCTGGTCGACCTCGACCATCAACCACGTTGCTACTATTGACGGTTGCAAAGCCAGACTTATCACTGACTTGTTGATGCCCTAAAGTCAGCGCATGACCACCCGCTTTATAGGTGAAAGCCAAAGACCAGGTGGtatt                         >  2:2597267/1‑150 (MQ=255)
                                                                                                                              gCCCCCTTCACCTTCTGGTCGACCTCGACCATCAACCACGTTGCTACTATTGACGGTTGCAAAGCCAGACTTATCACTGACGTGTTGATGCCCTAAAGTCAGCGCATGACCACCCGCTTTATAGGTGAAAGCCAAAGACCAGGTGGtatt                         >  2:2384042/1‑150 (MQ=255)
                                                                                                                                ccccTTCACCTTCTGGTCGACCTCGACCATCAACCACGTTGCTACTATTGACGGTTGCAAAGCCAGACTTATCACTGACTTGTTGATGCCCTAAAGTCAGCGCATGACCACCCGCTTTATAGGTGAAAGCCAAAGACCAGGTGGtattat                       >  1:4376439/1‑150 (MQ=255)
                                                                                                                                ccccTTCACCTTCTGGTCGACCTCGACCATCAACCACGTTGCTACTATTGACGGTTGCAAAGCCAGACTTATCACTGACTTGTTGATGCCCTAAAGTCAGCGCATGACCACCCGCTTTATAGGTGAAAGCCAAAGACCAGGTGGtattat                       >  2:7197595/1‑150 (MQ=255)
                                                                                                                                   cTTCACCTTCTGGTCGACCTCGACCATCAACCACGTTGCTACTATTGACGGTTGCAAAGCCAGACTTATCACTGACTTGTTGATGCCCTAAAGTCAGCGCATGACCACCCGCTTTATAGGTGAAAGCCAAAGACCAGGTGGTATTATCTa                    <  2:2921122/150‑1 (MQ=255)
                                                                                                                                    ttCACCTTCTGGTCGACCTCGACCATCAACCACGTTGCTACTATTGACGGTTGCAAAGCCAGACTTATCACTGACTTGTTGATGCCCTAAAGTCAGCGCATGACCACCCGCTTTATAGGTGAAAGCCAAAGACCAGGTGGTATTATCTAt                   >  1:5658181/1‑150 (MQ=255)
                                                                                                                                        ccTTCTGGTCGACCTCGACCATCAACCACGTTGCTACTATTGACGGTTGCAAAGCCAGACTTATCACTGACTTGTTGATGCCCTAAAGTCAGCGCATGACCACCCGCTTTATAGGTGAAAGCCAAAGACCAGGTGGTATTATCTATTTTc               >  1:8283361/1‑150 (MQ=255)
                                                                                                                                         cTTCTGGTCGACCTCGACCATCAACCACGTTGCTACTATTGACGGTTGCAAAGCCAGACTTATCACTGACTTGTTGATGCCCTAAAGTCAGCGCATGACCACCCGCTTTATAGGTGAAAGCCAAAGACCAGGTGGTATTATCTATTTTcc              >  1:7440976/1‑150 (MQ=255)
                                                                                                                                          ttCTGGTCGACCTCGACCATCAACCACGTTGCTACTATTGACGGTTGCAAAGCCAGACTTATCACTGACTTGTTGATGCCCTAAAGTCAGCGCATGACCACCCGCTTTATAGGTGAAAGCCAAAGACCAGGTGGTATTATCTATTTTccc             <  1:4868638/150‑1 (MQ=255)
                                                                                                                                               gTCGACCTCGACCATCAACCACGTTGCTACTATTGACGGTTGCAAAGCCAGACTTATCATTGACTTGTTGATGCCCTAAAGTCAGCGCATGACCACCCGCTTTATAGGTGGAAGCCAAAGACCAGGTGGTATTATCTATTTTCCCAtttt        >  2:1554752/1‑150 (MQ=255)
                                                                                                                                               gTCGACCTCGACCATCAACCACGTTGCTACTATTGACGGTTGCAAAGCCAGACTTATCACTGACTTGTTGATGCCCTAAAGTCAGCGCATGACCACCCGCTTTATAGGTGAAAGCCAAAGACCAGGTGGTATTATCTATTTTCCCAtttt        >  1:5456698/1‑150 (MQ=255)
                                                                                                                                               gTCGACCTCGACCATCAACCACGTTGCTACTATTGACGGTTGCAAAGCCAGACTTATCACTGACTTGTTGATGCCCTAAAGTCAGCGCATGACCACCCGCTTTATAGGTGAAAGCCAAAGACCAGGTGGTATTATCTATTTTCCCAtttt        >  2:3241874/1‑150 (MQ=255)
                                                                                                                                               gTCGACCTCGACCATCAACCACGTTGCTACTATTGACGGTTGCAAAGCCAGACTTATCACTGACTTGTTGATGCCCTAAAGTCAGCGCATGACCACCCGCTTTATAGGTGAAAGCCAAAGACCAGGTGGTATTATCTATTTTCCCAtttt        >  2:913492/1‑150 (MQ=255)
                                                                                                                                                tcgacctcgaccATCAACCACGTTGCTACTATTGACGGTTGCAAAGCCAGACTTATCACTGACTTGTTGATGCCCTAAAGTCAGCGCATGACCACCCGCTTTATAGGTGAAAGCCAAAGACCAGGTGGTATTATCTATTTTCCCATTTTc       <  2:8282203/150‑1 (MQ=255)
                                                                                                                                                  gacctcgaccATCAACCACGTTGCTACTATTGACGGTTGCAAAGCCAGACTTATCACTGACTTGTTGATGCCCTAAAGTCAGCGCATGACCACCCGCTTTATAGGTGAAAGCCAAAGACCAGGTGGTATTATCTATTTTCCCATTTTCAt     >  2:5452418/1‑150 (MQ=255)
                                                                                                                                                   acctcgaccATCAACCACGTTGCTACTATTGACGGTTGCAAAGCCAGACTTATCACTGACTTGTTGATGCCCTAAAGTCAGCGCATGACCACCCGCTTTATAGGTGAAAGCCAAAGACCAGGTGGTATTATCTATTTTCCCATTTTCAtt    <  2:7788346/150‑1 (MQ=255)
                                                                                                                                                   acctcgaccATCAACCACGTTGCTACTATTGACGGTTGCAAAGCCAGACTTATCACTGACTTGTTGATGCCCTAAAGTCAGCGCATGACCACCCGCTTTATAGGTGAAAGCCAAAGACCAGGTGGTATTATCTATTTTCCCATTTTCAtt    <  1:6000724/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CCCCGCCTTTAGGCCTGGAACACCTAAGGCGGCGAAATCGTAGCTGTAAAGGCCAAAAGTAGAATTTTCGCCCGCGCGCGTGAAGCTATTCACCATAACGTCGGTATATAGATATAGGGATTGGACGCCCCCTTCACCTTCTGGTCGACCTCGACCATCAACCACGTTGCTACTATTGACGGTTGCAAAGCCAGACTTATCACTGACTTGTTGATGCCCTAAAGTCAGCGCATGACCACCCGCTTTATAGGTGAAAGCCAAAGACCAGGTGGTATTATCTATTTTCCCATTTTCATTGT  >  NC_002947/5070153‑5070451

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