Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 233,931 G→A 28.9% S37N (AGT→AAT)  yafT → lipoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273233,9310GA28.9% 410.6 / 203.0 237S37N (AGT→AAT) yafTlipoprotein
Reads supporting (aligned to +/- strand):  ref base G (100/67);  new base A (41/27);  total (141/96)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.52e-01

Only 100 of 243 total aligned reads displayed.

AATATTCTAATCATTATGACAGGCGAGGGAGTGTCCAATTATGAATTCAAAAAAGCTTTGTTGCATATGTGTGTTATTCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAGTAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGATATTCAAGCTGCGCTTGATGGAGATGAGTTCCGCGTTCCCCTTAATTCTCCTGTAATACTTGTTC  >  CP009273/233782‑234080
                                                                                                                                                     |                                                                                                                                                     
 atagTCTAATCATTATGACAGGCGAGGGAGTGTCCAATTATGAATTCAAAAAAGCTTTGTTGCATATGTGTGTTATTCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAGt                                                                                                                                                      >  1:1483564/1‑150 (MQ=255)
       tAATCATTATGACAGGCGAGGGAGTGTCCAATTATGAATTCAAAAAAGCTTTGTTGCATATGTGTGTTATTCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAGTAATAtt                                                                                                                                                <  1:3812793/150‑1 (MQ=255)
          tCATTATGACAGGCGAGGGAGTGTCCAATTATGAATTCAAAAAAGCTTTGTTGCATATGTGTGTTATTCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAGTAATATtaat                                                                                                                                             <  1:1218821/150‑1 (MQ=255)
          tCATTATGACAGGCGAGGGAGTGTCCAATTATGAATTCAAAAAAGCTTTGTTGCATATGTGTGTTATTCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAGTAATATtaat                                                                                                                                             <  2:1599534/150‑1 (MQ=255)
          tCATTATGACAGGCGAGGGAGTGTCCAATTATGAATTCAAAAAAGCTTTGTTGCATATGTGTGTTATTCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAATAATATtaat                                                                                                                                             >  2:3357477/1‑150 (MQ=255)
               aTGACAGGCGAGGGAGTGTCCAATTATGAATTCAAAAAAGCTTTGTTGCATATGTGTGTTATTCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAGTAATATTAATaaaaa                                                                                                                                        >  2:1218621/1‑150 (MQ=255)
                tGACAGGCGAGGGAGTGTCCAATTATGAATTCAAAAAAGCTTTGTTGCATATGTGTGTTATTCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAGTAATATTAATAAAAac                                                                                                                                       >  1:3655471/1‑150 (MQ=255)
                  aCAGGCGAGGGAGTGTCCAATTATGAATTCAAAAAAGCTTTGTTGCATATGTGTGTTATTCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAGTAATATTAATAAAAacac                                                                                                                                     <  1:1755687/150‑1 (MQ=255)
                   cAGGCGAGGGAGTGTCCAATTATGAATTCAAAAAAGCTTTGTTGCATATGTGTGTTATTCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAGTAATATTAATAAAAACACt                                                                                                                                    >  1:4057955/1‑150 (MQ=255)
                    aGGCGAGGGAGTGTCCAATTATGAATTCAAAAAAGCTTTGTTGCATATGTGTGTTATTCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAGTAATATTAATAAAAACACTc                                                                                                                                   >  1:4033397/1‑150 (MQ=255)
                    aGGCGAGGGAGTGTCCAATTATGAATTCAAAAAAGCTTTGTTGCATATGTGTGTTATTCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAATAATATTAATAAAAACACTc                                                                                                                                   >  1:247242/1‑150 (MQ=255)
                     ggCGAGGGAGTGTCCAATTATGAATTCAAAAAAGCTTTGTTGCATATGTGTGTTATTCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAGTAATATTAATAAAAACACTcc                                                                                                                                  >  2:2368482/1‑150 (MQ=255)
                     ggCGACGGAGTGTCCAATTATGAATTCAAAAAAGCTTTGTTGCATATGTGTGTTATTCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATTGATGAAAATAAGAAAAAAGCTCAAGTCACACAAAGAAATATTAATAAAAAAACTcc                                                                                                                                  >  2:2944921/1‑150 (MQ=255)
                       cGAGGGAGTGTCCAATTATGAATTCAAAAAAGCTTTGTTGCATATGTGTGTTATTCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAGTAATATTAATAAAAACACTcccc                                                                                                                                <  1:1722407/150‑1 (MQ=255)
                         aGGGAGTGTCCAATTATGAATTCAAAAAAGCTTTGTTGCATATGTGTGTTATTCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAATAATATTAATAAAAACACTCCCCAg                                                                                                                              <  2:3286426/150‑1 (MQ=255)
                          gggAGTGTCCAATTATGAATTCAAAAAAGCTTTGTTGCATATGTGTGTTATTCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAGTAATATTAATAAAAACACTCCCCAGc                                                                                                                             >  1:1627266/1‑150 (MQ=255)
                           ggAGTGTCCAATTATGAATTCAAAAAAGCTTTGTTGCATATGTGTGTTATTCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAATAATATTAATAAAAACACTCCCCAGCa                                                                                                                            >  2:3043683/1‑150 (MQ=255)
                              gtgtCCAATTATGAATTCAAAAAAGCTTTGTTGCATATGTGTGTTATTCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAGTAATATTAATAAAAACACTCCCCAGCAACt                                                                                                                         <  1:224033/150‑1 (MQ=255)
                               tgtCCAATTATGAATTCAAAAAAGCTTTGTTGCATATGTGTGTTATTCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAGTAATATTAATAAAAACACTCCCCAGCAACTg                                                                                                                        <  2:3444648/150‑1 (MQ=255)
                                gtCCAATTATGAATTCAAAAAAGCTTTGTTGCATATGTGTGTTATTCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAGTAATATTAATAAAAACACTCCCCAGCAACTga                                                                                                                       >  1:346116/1‑150 (MQ=255)
                                 tCCAATTATGAATTCAAAAAAGCTTTGTTGCATATGTGTGTTATTCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAATAATATTAATAAAAACACTCCCCAGCAACTgac                                                                                                                      <  2:4122087/150‑1 (MQ=255)
                                   cAATTATGAATTCAAAAAAGCTTTGTTGCATATGTGTGTTATTCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAGTAATATTAATAAAAACACTCCCCAGCAACTgacag                                                                                                                    <  1:3608223/150‑1 (MQ=255)
                                   cAATTATGAATTCAAAAAAGCTTTGTTGCATATGTGTGTTATTCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAGTAATATTAATAAAAACACTCCCCAGCAACTgacag                                                                                                                    <  1:3506099/150‑1 (MQ=255)
                                     aTTATGAATTCAAAAAAGCTTTGTTGCATATGTGTGTTATTCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAATAATATTAATAAAAACACTCCCCAGCAAATgacagac                                                                                                                  >  2:463878/1‑150 (MQ=255)
                                      ttATGAATTCAAAAAAGCTTTGTTGCATATGTGTGTTATTCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAGTAATATTAATAAAAACACTCCCCAGCAACTgacagaca                                                                                                                 <  2:1362863/150‑1 (MQ=255)
                                           aaTTCAAAAAAGCTTTGTTGCATATGTGTGTTATTCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAGTAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGAt                                                                                                            >  2:1987353/1‑150 (MQ=255)
                                           aaTTCAAAAAAGCTTTGTTGCATATGTGTGTTATTCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAGTAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGAt                                                                                                            <  2:1642011/150‑1 (MQ=255)
                                            aTTCAAAAAAGCTTTGTTGCATATGTGTGTTATTCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAGTAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGAtt                                                                                                           >  2:2978895/1‑150 (MQ=255)
                                              tCAAAAAAGCTTTGTTGCATATGTGTGTTATTCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAGTAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTa                                                                                                         <  2:147574/150‑1 (MQ=255)
                                              tCAAAAAAGCTTTGTTGCATATGTGTGTTATTCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAATAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTa                                                                                                         <  1:741260/150‑1 (MQ=255)
                                                aaaaaaGCTTTGTTGCATATGTGTGTTATTCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAGTAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTAtt                                                                                                       <  2:2536860/150‑1 (MQ=255)
                                                aaaaaaGCTTTGTTGCATATGTGTGTTATTCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAGTAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTAtt                                                                                                       <  1:743690/150‑1 (MQ=255)
                                                    aaGCTTTGTTGCATATGTGTGTTATTCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAGTAATAGTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCTGt                                                                                                   >  1:4033476/1‑150 (MQ=255)
                                                       cTTTGTTGCATATGTGTGTTATTCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGAAAGATGCTCAAGTCACACAAAGTAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAAt                                                                                                <  1:1009532/150‑1 (MQ=255)
                                                       cTTTGTTGCATATGTGTGTTATTCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAATAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAAt                                                                                                <  1:2705772/150‑1 (MQ=255)
                                                         ttgttgCATATGTGTGTTATTCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAGTAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGa                                                                                              >  2:2493264/1‑150 (MQ=255)
                                                          tgttgCATATGTGTGTTATTCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAGTAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGaa                                                                                             >  1:838753/1‑150 (MQ=255)
                                                          tgttgCATATGTGTGTTATTCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAGTAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGaa                                                                                             >  1:1289395/1‑150 (MQ=255)
                                                           gttgCATATGTGTGTTATTCTCAATTATTTAAATATTTTAATATTAATATTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAGTAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGaaa                                                                                            >  2:1584715/1‑150 (MQ=255)
                                                             tgCATATGTGTGTTATTCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAGTAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGAAAcc                                                                                          >  1:1194550/1‑150 (MQ=255)
                                                             tgCATATGTGTGTTATTCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAGTAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGAAAcc                                                                                          >  1:609999/1‑150 (MQ=255)
                                                               cATATGTGTGTTATTCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAGTAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGAAACCAc                                                                                        >  2:1094079/1‑150 (MQ=255)
                                                                atatGTGTGTTATTCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAATAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTCATGAAACCAct                                                                                       >  2:3156521/1‑150 (MQ=255)
                                                                  atGTGTGTTATTCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAGTAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGAAACCActct                                                                                     >  2:1179188/1‑150 (MQ=255)
                                                                      gtgtTATTCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAGTAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGAAACCACTCTGGcc                                                                                 <  2:2137866/150‑1 (MQ=255)
                                                                       tgtTATTCTCGCTGCTTGCAGTATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAGTAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGAAACCACTCTGGCCg                                                                                <  2:782359/150‑1 (MQ=255)
                                                                       tgtTATTCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAGTAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGAAACCACTCTGGCCg                                                                                <  1:1851982/150‑1 (MQ=255)
                                                                           aTTCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAATAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGAAACCACTCTGGCCGTATc                                                                            <  2:1879481/150‑1 (MQ=255)
                                                                            ttCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAGTAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGAAACCACTCTGGCCGTATcc                                                                           <  1:318152/150‑1 (MQ=255)
                                                                                 gctgctTGCAGGATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAATAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGAAACCACTCTGGCCGTATCCGAaga                                                                      >  1:4083083/1‑150 (MQ=255)
                                                                                  ctgctTGCAGGATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAGTAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTCATGAAACCACTCTGGCCGTATCCGAagag                                                                     >  2:2548630/1‑150 (MQ=255)
                                                                                  ctgctTGCAGGATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAATAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGAAACCACTCTGGCCGTATCCGAagag                                                                     >  1:3991078/1‑150 (MQ=255)
                                                                                       tGCAGGATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAGTAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGatat                                                                >  2:1243443/1‑150 (MQ=255)
                                                                                        gCAGGATGTGCCTCTTAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAGTAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGATAtt                                                               <  2:3304685/150‑1 (MQ=255)
                                                                                        gCAGGATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAGTAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGATAtt                                                               >  2:482737/1‑150 (MQ=255)
                                                                                         cAGGATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAATAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGATATTc                                                              >  1:854631/1‑150 (MQ=255)
                                                                                          aGGATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAGTAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGATATTCa                                                             >  2:798601/1‑150 (MQ=255)
                                                                                          aGGATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAGTAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGATATTCa                                                             <  1:2385163/150‑1 (MQ=255)
                                                                                           ggATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAGTAATATTAATAAAAACACTACCCAGCAACTGACAGACAAAGATTTATTCGGTAATGAAACCACTCTGGCCGTATACGAAGAGGATATTCaa                                                            >  2:2306932/1‑150 (MQ=255)
                                                                                            gATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAATAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGATATTCAAg                                                           >  1:671881/1‑150 (MQ=255)
                                                                                             aTGTGCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAGTAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGATATTCAAGc                                                          >  2:1615954/1‑150 (MQ=255)
                                                                                                    tctGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAATAATATTAATAAAAACACTCCCCCGCAACTGACAGACAAAGATTTATTCGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGATATTCAAGCTGCGCtt                                                   <  2:373554/150‑1 (MQ=255)
                                                                                                            ttctATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAGTAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGATATTCAAGCTGCGCTTGATGgaga                                           >  2:3179742/1‑150 (MQ=255)
                                                                                                             tctATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAATAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGATATTCAAGCTGCGCTTGATGGAGAt                                          <  1:3247692/150‑1 (MQ=255)
                                                                                                              ctATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAATAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGATATTCAAGCTGCGCTTGATGGAGATg                                         <  1:944871/150‑1 (MQ=255)
                                                                                                               tATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAGTAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGATATTCAAGCTGCGCTTGATGGAGATGa                                        <  1:2994829/150‑1 (MQ=255)
                                                                                                                aTTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAGTAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGATATTCAAGCTGCGCTTGATGGAGATGAg                                       >  1:4141420/1‑150 (MQ=255)
                                                                                                                  tgatgaAAAGAAGAAAAAAGCTCAAGTCACACAAAGTAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGATATTCAAGCTGCGCTTGATGGAGATGAGtt                                     <  1:2935531/150‑1 (MQ=255)
                                                                                                                    atgaAAAGAAGAAAAAAGCTCAAGTCACACAAAGTAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGATATTCAAGCTGCGCTTGATGGAGATGAGTTcc                                   <  2:131978/150‑1 (MQ=255)
                                                                                                                    atgaAAAGAAGAAAAAAGCTCAAGTCACACAAAGTAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGATATTCAAGCTGCGCTTGATGGAGATGAGTTcc                                   <  2:2898447/150‑1 (MQ=255)
                                                                                                                        aaaGAAGAAAAAAGCTCAAGTCACACAAAGTAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGATATTCAAGCTGCGCTTGATGGAGATGAGTTCCGCGt                               >  2:2245809/1‑150 (MQ=255)
                                                                                                                        aaaGAAGAAAAAAGCTCAAGTCACACAAAATAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGATATTCAAGCTGCGCTTGATGGAGATGAGTTCCGCGt                               >  2:4162036/1‑150 (MQ=255)
                                                                                                                        aaaGAAGAAAAAAGCTCAAGTCACACAAAATAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGATATTCAAGCTGCGCTTGATGGAGATGAGTTCCGCGt                               >  2:2630348/1‑150 (MQ=255)
                                                                                                                           gaagaaAAAAGCTCAAGTCACACAAAGTAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGATATTCAAGCTGCGCTTGATGGAGATGAGTTCCGCGTTcc                            <  2:1794795/150‑1 (MQ=255)
                                                                                                                             agaaAAAAGCTCAAGTCACACAAAATAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGATATTCAAGCTGCGCTTGATGGAGATGAGTTCCGCGTTcccc                          >  2:2642935/1‑150 (MQ=255)
                                                                                                                              gaaAAAAGCTCAAGTCACACAAAATAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGATATTCAAGCTGCGCTTGATGGAGATGAGTTCCGCGTTCCCCt                         >  1:1861338/1‑150 (MQ=255)
                                                                                                                               aaaaaaGCTCAAGTCACACAAAGTAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGATATTCAAGCTGCGCTTGATGGAGATGAGTTCCGCGTTCCCCtt                        <  1:1557323/150‑1 (MQ=255)
                                                                                                                               aaaaaaGCTCAAGTCACACAAAGTAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGATATTCAAGCTGCGCTTGATGGAGATGAGTTCCGCGTTCCCCtt                        >  2:3914236/1‑150 (MQ=255)
                                                                                                                                  aaaGCTCAAGTCACACAAAGTAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGATATTCAAGCTGCGCTTGATGGAGATGAGTTCCGCGTTCCCCTTAAt                     <  2:3354705/150‑1 (MQ=255)
                                                                                                                                  aaaGCTCAAGTCACACAAAATAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGATATTCAAGCTGCGCTTGATGGAGATGAGTTCCGCGTTCCCCTTAAt                     <  2:2092002/150‑1 (MQ=255)
                                                                                                                                    aGCTCAAGTCACACAAAATAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGATATTCAAGCTGCGCTTGATGGAGATGAGTTCCGCGTTCCCCTTAATtc                   >  1:517359/1‑150 (MQ=255)
                                                                                                                                      cTCAAGTCACACAAAATAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGATATTCAAGCTGCGCTTGATGGAGATGAGTTCCGCGTTCCCCTTAATtctc                 <  2:2820013/150‑1 (MQ=255)
                                                                                                                                        caggTCACACAAAGTAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGATATTCAAGCTGCGCTTGATGGAGATGAGTTCCGCGTTCCCCTTAATTCTCCt               >  1:1971785/4‑150 (MQ=255)
                                                                                                                                        cAAGTCACACAAAGTAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGATATTCAAGCTGCGCTTGATGGAGATGAGTTCCGCGTTCCCCTTAATTCTCCt               <  1:443922/150‑1 (MQ=255)
                                                                                                                                        cAAGTCACACAAAATAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGATATTCAAGCTGCGCTTGATGGAGATGAGTTCCGCGTTCCCCTTAATTCTCCt               >  2:168893/1‑150 (MQ=255)
                                                                                                                                         aaGTCACACAAAGTAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGATATTCAAGCTGCGCTTGATGGAGATGAGTTCCGCGTTCCCCTTAATTCTCCTg              <  2:976273/150‑1 (MQ=255)
                                                                                                                                            tCACACAAAATAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGCTATTCAAGCTGCGCTTGATGGAGATGAGTTCCGCGTTCCCCTTAATTCTCCTGTaa           >  2:933229/1‑150 (MQ=255)
                                                                                                                                            tCACACAAAATAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGATATTCAAGCTGCGCTTGATGGAGATGAGTTCCGCGTTCCCCTGAATTCTCCTGTaa           >  2:6211/1‑150 (MQ=255)
                                                                                                                                             cacacaAAGTAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGATATTCAAGCTGCGCTTGATGGAGATGAGTTCCGCGTTCCCCTTAATTCTCCTGTAAt          <  2:2240474/150‑1 (MQ=255)
                                                                                                                                             cacacaAAATAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGATATTCAAGCTGCGCTTGATGGAGATGAGTTCCGCGTTCCCCTTAATTCTCCTGTAAt          <  1:3778729/150‑1 (MQ=255)
                                                                                                                                              acacaAAGTAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGATATTCAAGCTGCGCTTGATGGAGATGAGTTCCGCGTTCCCCTTAATTCTCCTGTAATa         <  2:1502321/150‑1 (MQ=255)
                                                                                                                                              acacaAAGTAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGATATTCAAGCTGCGCTTGATGGAGATGAGTTCCGCGTTCCCCTTAATTCTCCTGTAATa         >  1:3649644/1‑150 (MQ=255)
                                                                                                                                              acacaAAATAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGATATTCAAGCTGCGCTTGATGGAGATGAGTTCCGCGTTCCCCTTAATTCTCCTGTAATa         >  2:2679607/1‑150 (MQ=255)
                                                                                                                                               cacaAAGTAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGATATTCAAGCTGCGCTTGATGGAGATGAGTTCCGCGTTCCCCTTAATTCTCCTGTAATAc        >  1:2619831/1‑150 (MQ=255)
                                                                                                                                                 caAAGTAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGATATTCAAGCTGCGCTTGATGGAGATGAGTTCCGCGTTCCCCTTAATTCTCCTGTAATACtt      >  1:2254162/1‑150 (MQ=255)
                                                                                                                                                 caAAGTAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGATATTCAAGCTGCGCTTGATGGAGATGAGTTCCGCGTTCCCCTTAATTCTCCTGTAATACtt      >  1:1831843/1‑150 (MQ=255)
                                                                                                                                                  aaaaTAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGATATTCAAGCTGCGCTTGATGGAGATGAGTTCCGCGTTCCCCTTAATTCTCCTGTAATACTTg     >  1:2402987/1‑150 (MQ=255)
                                                                                                                                                  aaaGTAATATTAATAAAAACACTCCCCAGCAACTTACAGACAAAGATTTATTCGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGATATTCAAGCTGCGCTTGATGGAGATGAGTTCCGCGTTCCCCTTAATTCTCCTGTAATACTTg     >  2:629856/1‑150 (MQ=255)
                                                                                                                                                  aaaGTAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGATATTCAAGCTGCGCTTGATGGAGATGAGTTCCGCGTTCCCCTTAATTCTCCTGTAATACTTg     <  1:2011917/150‑1 (MQ=255)
                                                                                                                                                    aGTAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGATATTCAAGCTGCGCTTGATGGAGATGAGTTCCGCGTTCCCCTTAATTCTCCTGTAATACTTGtt   >  1:2633598/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
AATATTCTAATCATTATGACAGGCGAGGGAGTGTCCAATTATGAATTCAAAAAAGCTTTGTTGCATATGTGTGTTATTCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAGTAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGATATTCAAGCTGCGCTTGATGGAGATGAGTTCCGCGTTCCCCTTAATTCTCCTGTAATACTTGTTC  >  CP009273/233782‑234080

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