Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,376,178 C→A 92.6% E1004* (GAA→TAA)  mscM ← mechanosensitive channel protein, miniconductance

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,376,1780CA92.6% 838.2 / 57.2 284E1004* (GAA→TAA) mscMmechanosensitive channel protein, miniconductance
Reads supporting (aligned to +/- strand):  ref base C (10/11);  new base A (135/128);  total (145/139)
Fisher's exact test for biased strand distribution p-value = 8.22e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

Only 100 of 289 total aligned reads displayed.

TACGGTGACCCATCTCAGCGGCGTAAATACGCAGCTCGAAAATCTGAATCCCCTGTTGCAGATCCACCAGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTCGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATAAACTGCTCGGTAATAAACGCCTTGTTCGGCACGATTATCTCTTTACGGTCCCAGTCGCTGATGGTGGTGGCGC  >  CP009273/4376031‑4376324
                                                                                                                                                   |                                                                                                                                                  
tACGGTGACCCATCTCAGCGGCGTAAATACGCAGCTCGAAAATCTGAATCCCCTGTTGCAGATCCACCAGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGc                                                                                                                                                  <  2:1387536/150‑1 (MQ=255)
 aCGGTGACCCATCTCAGCGGCGTAAATACGCAGCTCGAAAATCTGAATCCCCTGTTGCAGATCCACCAGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCt                                                                                                                                                 <  1:537022/150‑1 (MQ=255)
 aCGGTGACCCATCTCAGCGGCGTAAATACGCAGCTCGAAAATCTGAATCCCCTGTTGCAGATCCACCAGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCt                                                                                                                                                 <  2:3083396/150‑1 (MQ=255)
   ggTGACCCATCTCAGCGGCGTAAATACGCAGCTCGAAAATCTGAATCCCCTGTTGCAGATCCACCAGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTAt                                                                                                                                               >  2:1443320/1‑150 (MQ=255)
    gTGACCCATCTCAGCGGCGTAAATACGCAGCTCGAAAATCTGAATCCCCTGTTGCAGATCCACCAGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTAtt                                                                                                                                              >  1:1311103/1‑150 (MQ=255)
    gTGACCCATCTCAGCGGCGTAAATACGCAGCTCGAAAATCTGAATCCCCTGTTGCAGATCCACCAGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTAtt                                                                                                                                              <  2:2060661/150‑1 (MQ=255)
      gACCCATCTCAGCGGCGTAAATACGCAGCTCGAAAATCTGAATCCCCTGTTGCAGATCCACCAGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTgg                                                                                                                                            >  1:2552437/1‑150 (MQ=255)
       aCCCATCTCAGCGGCGTAAATACGCAGCTCGAAAATCTGAATCCCCTGTTGCAGATCCACCAGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGc                                                                                                                                           >  1:700925/1‑150 (MQ=255)
           aTCTCAGCGGCGTAAATACGCAGCTCGAAAATCTGAATCCCCTGTTGCAGATCCACCAGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCg                                                                                                                                       >  1:1672590/1‑150 (MQ=255)
             ctcAGCGGCGTAAATACGCAGCTCGAAAATCTGAATCCCCTGTTGCAGATCCACCAGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGc                                                                                                                                     >  2:1599935/1‑150 (MQ=255)
                    gcgTAAATACGCAGCTCGAAAATCTGAATCCCCTGTTGCAGATCCACCAGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTCGCTATTGGCATCGGCAGGGGcc                                                                                                                              >  1:2663131/1‑150 (MQ=255)
                      gTAAATACGCAGCTCGAAAATCTGAATCCCCTGTTGCAGATCCACCAGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCgg                                                                                                                            >  1:2377035/1‑150 (MQ=255)
                       tAAATACGCAGCTCGAAAATCTGAATCCCCTGTTGCAGATCCACCAGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGt                                                                                                                           <  2:3974450/150‑1 (MQ=255)
                         aaTACGCAGCTCGAAAATCTGAATCCCCTGTTGCAGATCCACCAGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTAt                                                                                                                         <  1:3002718/150‑1 (MQ=255)
                                aGCTCGAAAATCTGAATCCCCTGTTGCAGATCCACCAGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAac                                                                                                                  >  1:1518277/1‑150 (MQ=255)
                                 gCTCGAAAATCTGAATCCCCTGTTGCAGATCCACCAGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAaca                                                                                                                 <  1:1082368/150‑1 (MQ=255)
                                  cTCGAAAATCTGAATCCCCTGTTGCAGATCCACCAGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAacac                                                                                                                <  2:3938018/150‑1 (MQ=255)
                                   tCGAAAATCTGAATCCCCTGTTGCAGATCCACCAGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAAcacc                                                                                                               <  2:1878248/150‑1 (MQ=255)
                                       aaaTCTGAATCCCCTGTTGCAGATCCACCAGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCAcgc                                                                                                           >  1:4018259/1‑150 (MQ=255)
                                       aaaTCTGAATCCCCTGTTGCAGATCCACCAGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCAcgc                                                                                                           >  1:3580452/1‑150 (MQ=255)
                                         aTCTGAATCCCCTGTTGCAGATCCACCAGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCAcgcgc                                                                                                         >  2:1289153/1‑150 (MQ=255)
                                         aTCTGAATCCCCTGTTGCAGATCCACCAGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCAcgcgc                                                                                                         <  2:4575839/150‑1 (MQ=255)
                                            tGAATCCCCTGTTGCAGATCCACCAGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTg                                                                                                      <  1:2503115/150‑1 (MQ=255)
                                            tGAATCCCCTGGTGCAGATCCACCAGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTg                                                                                                      <  2:3495033/150‑1 (MQ=255)
                                                 ccccTGTTGCAGATCCACCAGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGa                                                                                                 <  2:202350/150‑1 (MQ=255)
                                                 ccccTGTTGCAGATCCACCAGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGa                                                                                                 >  2:4326882/1‑150 (MQ=255)
                                                  cccTGTTGCAGATCCACCAGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTCGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAg                                                                                                >  2:3708951/1‑150 (MQ=255)
                                                    cTGTTGCAGATCCACCAGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTc                                                                                              <  2:1745803/150‑1 (MQ=255)
                                                      gTTGCAGATCCACCAGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCag                                                                                            <  1:4266046/150‑1 (MQ=255)
                                                      gTTGCAGATCCACCAGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCag                                                                                            <  1:4187007/150‑1 (MQ=255)
                                                         gCAGATCCACCAGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCagaga                                                                                         <  1:3782114/150‑1 (MQ=255)
                                                            gATCCACCAGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCg                                                                                      <  2:2552437/150‑1 (MQ=255)
                                                              tCCACCAGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGAc                                                                                    >  2:1464508/1‑150 (MQ=255)
                                                               ccaccaGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGAcc                                                                                   <  2:570847/150‑1 (MQ=255)
                                                                 accaGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAg                                                                                 >  2:861735/1‑150 (MQ=255)
                                                                    aGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTg                                                                              <  2:4398149/150‑1 (MQ=255)
                                                                      gAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGAt                                                                            >  1:3761063/1‑150 (MQ=255)
                                                                      gAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGAt                                                                            >  1:116146/1‑150 (MQ=255)
                                                                      gAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGAt                                                                            <  2:2917364/150‑1 (MQ=255)
                                                                      gAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGAt                                                                            >  1:4500867/1‑150 (MQ=255)
                                                                       aaGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATa                                                                           >  2:2707133/1‑150 (MQ=255)
                                                                         gACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATaaa                                                                         <  1:797980/150‑1 (MQ=255)
                                                                         gACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATaaa                                                                         >  1:2039242/1‑150 (MQ=255)
                                                                          aCTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATAAAc                                                                        >  2:1496897/1‑150 (MQ=255)
                                                                             tCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATAAACTGc                                                                     >  2:2653269/1‑150 (MQ=255)
                                                                              ccGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATAAACTGCt                                                                    <  1:4580121/150‑1 (MQ=255)
                                                                              ccGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGGGACCGAGTCAGAGAGCGACCAGTTGATAAACTGCt                                                                    <  2:4247553/150‑1 (MQ=255)
                                                                                     cGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATAAACTGCTCGGtaat                                                             >  1:3459243/1‑150 (MQ=255)
                                                                                      gggTTGTCGATCACCAGTGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATAAACTGCTCGGtaata                                                            >  2:301008/1‑150 (MQ=255)
                                                                                       ggTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTCGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATAAACTGCTCGGtaataa                                                           <  1:2700891/150‑1 (MQ=255)
                                                                                       ggTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTCGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGACTCAGAGAGCGACCAGTTGATAAACTGCTCGGtaataa                                                           <  1:565659/150‑1 (MQ=255)
                                                                                          tGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATAAACTGCTCGGTAATAAACg                                                        <  2:3047111/150‑1 (MQ=255)
                                                                                          tGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATGGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGCCAGAGAGCGACCAGTTGATACACTGCTCGGTAATAAACg                                                        <  2:2301170/150‑1 (MQ=255)
                                                                                            tCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATAAACTGCTCGGTAATAAACGcc                                                      <  1:2567558/150‑1 (MQ=255)
                                                                                             cGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATAAACTGCTCGGTAATAAACGCCt                                                     <  2:1588511/150‑1 (MQ=255)
                                                                                              gATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATAAACTGCTCGGTAATAAACGCCtt                                                    >  1:218886/1‑150 (MQ=255)
                                                                                                tCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATAAACTGCTCGGTAATAAACGCCTTGt                                                  >  1:4559642/1‑150 (MQ=255)
                                                                                                 cACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATAAACTGCTCGGTAATAAACGCCTTGtt                                                 <  2:3109702/150‑1 (MQ=255)
                                                                                                  aCCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATAAACTGCTCGGTAATAAACGCCTTGTTc                                                <  2:4300775/150‑1 (MQ=255)
                                                                                                   ccAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTCGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATAAACTGCTCGGTAATAAACGCCTTGTTCg                                               >  1:958078/1‑150 (MQ=255)
                                                                                                   ccAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATAAACTGCTCGGTAATAAACGCCTTGTTCg                                               >  2:3388323/1‑150 (MQ=255)
                                                                                                     aGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATAAACTGCTCGGTAATAAACGCCTTGTTCGGc                                             <  1:3235208/150‑1 (MQ=255)
                                                                                                      gCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATAAACTGCTCGGTAATAAACGCCTTGTTCGGCa                                            <  2:4178633/150‑1 (MQ=255)
                                                                                                      gCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATAAACTGCTCGGTAATAAACGCCTTGTTCGGCa                                            <  1:1930285/150‑1 (MQ=255)
                                                                                                       cGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATAAACTGCTCGGTAATAAACGCCTTGTTCGGCAc                                           <  1:2859937/150‑1 (MQ=255)
                                                                                                       cGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATAAACTGCTCGGTAATAAACGCCTTGTTCGGCAc                                           <  1:33459/150‑1 (MQ=255)
                                                                                                        gAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATAAACTGCTCGGTAATAAACGCCTTGTTCGGCACg                                          <  1:278064/150‑1 (MQ=255)
                                                                                                         agcagcGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATAAACTGCTCGGTAATAAACGCCTTGTTCGGCACGa                                         >  2:2129880/1‑150 (MQ=255)
                                                                                                           cagcGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATAAACTGCTCGGTAATAAACGCCTTGTTCGGCACGAtt                                       <  2:2883607/150‑1 (MQ=255)
                                                                                                            agcGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTCGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATAAACTGCTCGGTAATAAACGCCTTGTTCGGCACGATTa                                      >  2:465517/1‑150 (MQ=255)
                                                                                                                      tGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATAAACTGCTCGGTAATAAACGCCTTGTTCGGCACGATTATCTCTTTACg                            >  2:1979772/1‑150 (MQ=255)
                                                                                                                       gCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATAAACTGCTCGGTAATAAANGCCTTGTTCGGCACGATTATCTCTTTACgg                           <  1:4281329/150‑1 (MQ=255)
                                                                                                                        cGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATAAACTGCTCGGTAATAAACGCCTTGTTCGGCACGATTATCTCTTTACGGt                          >  1:2287257/1‑150 (MQ=255)
                                                                                                                         ggTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATAAACTGCTCGGTAATAAACGCCTTGTTCGGCACGATTATCTCTTTACGGTc                         >  1:1425170/1‑150 (MQ=255)
                                                                                                                            gAGCAGGATTTCCGTCACTTCTTCGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATAAACTGCTCGGTAATAAACGCCTTGTTCGGCACGATTATCTCTTTACGGTCCCa                      <  2:3233568/150‑1 (MQ=255)
                                                                                                                            gAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATAAACTGCTCGGTAATAAACGCCTTGTTCGGCACGATTATCTCTTTACGGTCCCa                      <  1:4393699/150‑1 (MQ=255)
                                                                                                                            gAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATAAACTGCTCGGTAATAAACGCCTTGTTCGGCACGATTATCTCTTTACGGTCCCa                      >  2:4109602/1‑150 (MQ=255)
                                                                                                                            gAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATAAACTGCTCGGTAATAAACGCCTTGTTCGGCACGATTATCTCTTTACGGTCCCa                      <  1:3098363/150‑1 (MQ=255)
                                                                                                                            gAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATAAACTGCTCGGTAATAAACGCCTTGTTCGGCACGATTATCTCTTTACGGTCCCa                      >  2:2249750/1‑150 (MQ=255)
                                                                                                                             aGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATAAACTGCTCGGTAATAAACGCCTTGTTCGGCACGATTATCTCTTTACGGTCCCAg                     >  2:2343723/1‑150 (MQ=255)
                                                                                                                             aGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATAAACTGCTCGGTAATAAACGCCTTGTTCGGCACGATTATCTCTTTACGGTCCCAg                     >  1:2665188/1‑150 (MQ=255)
                                                                                                                             aGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATAAACTGCTCGGTAATAAACGCCTTGTTCGGCACGATTATCTCTTTACGGTCCCAg                     >  2:32657/1‑150 (MQ=255)
                                                                                                                             aGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATAAACTGCTCGGTAATAAACGCCTTGTTCGGCACGATTATCTCTTTACGGTCCCAg                     >  2:164584/1‑150 (MQ=255)
                                                                                                                              gCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATAAACTGCTCGGTAATAAACGCCTTGTTCGGCACGATTATCTCTTTACGGTCCCAGt                    >  1:913209/1‑150 (MQ=255)
                                                                                                                                aGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGTGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATAAACTGCTCGGTAATAAACGCCTTGTTCGGCACGATTATCTCTTTAAGGTCCCAGTCg                  >  1:3075302/1‑150 (MQ=255)
                                                                                                                                aGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATAAACTGCTCGGTAATAAACGCCTTGTTCGGCACGATTATCTCTTTACGGTCCCAGTCg                  >  1:1134001/1‑150 (MQ=255)
                                                                                                                                 ggATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATAAACTGCTCGGTAATAAACGCCTTGTTCGGCACGATTATCTCTTTACGGTCCCAGTCGc                 >  2:4429623/1‑150 (MQ=255)
                                                                                                                                   aTTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATAAACTGCTCGGTAATAAACGCCTTGTTCGGCACGATTATCTCTTTACGGTCCCAGTAGCTg               >  1:4077270/1‑150 (MQ=255)
                                                                                                                                    tttCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATAAACTGCTCGGTAATAAACGCCTTGTTCGGCACGATTATCTCTTTACGGTCCCAGTCGCTGa              <  1:3321730/150‑1 (MQ=255)
                                                                                                                                      tCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATAAACTGCTCGGTAATAAACGCCTTGTTCGGCACGATTATCTCTTTACGGTCCCAGTCGCTGAtg            <  1:1599935/150‑1 (MQ=255)
                                                                                                                                      tCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATAAACTGCTCGGTAATAAACGCCTTGTTCGGCACGATTATCTCTTTACGGTCCCAGTCGCTGAtg            >  2:1538596/1‑150 (MQ=255)
                                                                                                                                       ccGTCACTTCTTAGCTATTGGCATTGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATAAACTGCTCGGTAATAAACGCCTTGTTCGGCACGATTATCTCTTTACGGTCCCAGTCGCTGAtgg           <  1:1815549/150‑1 (MQ=255)
                                                                                                                                         gTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATAAACTGCTCGGTAATAAACGCCTTGTTCGGCACGATTATCTCTTTACGGTCCCAGTCGCTGAtggtg         >  1:2069641/1‑150 (MQ=255)
                                                                                                                                          tCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATAAACTGCTCGGTAATAAACGCCTTGTTCGGCACGATTATCTCTTTACGGTCCCAGTCGCTGAtggtgg        >  2:3837596/1‑150 (MQ=255)
                                                                                                                                          tCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATAAACTGCTCGGTAATAAACGCCTTGTTCGGCACGATTATCTCTTTACGGTCCCAGTCGCTGAtggtgg        <  1:3535167/150‑1 (MQ=255)
                                                                                                                                            aCTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATAAACTGCTCGGTAATAACCGCCTTGTTCGGCACGATTATCTCTTTACGGTCCCAGTCGCTGAtggtggtg      <  2:153208/150‑1 (MQ=255)
                                                                                                                                            aCTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATAAACTGCTCGGTAATAAACGCCTTGTTCGGCACGATTATCTCTTTACGGTCCCAGTCGCTGAtggtggtg      <  2:3638146/150‑1 (MQ=255)
                                                                                                                                             cttcttAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATAAACTGCTCGGTAATAAACGCCTTGTTCGGCACGATTATCTCTTTACGGTCCCAGTCGCTGAtggtggtgg     >  1:927593/1‑150 (MQ=255)
                                                                                                                                             cttcttAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATAAACTGCTCGGTAATAAACGCCTTGTTCGGCACGATTATCTCTTTACGGTCCCAGTCGCTGAtggtggtgg     <  2:2012327/150‑1 (MQ=255)
                                                                                                                                                cttaGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATAAACTGCTCGGTAATAAACGCCTTGTTCGGCACGATTATCTCTTTACGGTCCCAGTCGCTGATGGTGGTGgcgc  >  1:420941/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
TACGGTGACCCATCTCAGCGGCGTAAATACGCAGCTCGAAAATCTGAATCCCCTGTTGCAGATCCACCAGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTCGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATAAACTGCTCGGTAATAAACGCCTTGTTCGGCACGATTATCTCTTTACGGTCCCAGTCGCTGATGGTGGTGGCGC  >  CP009273/4376031‑4376324

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: ‑