Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,088,734 G→T 37.3% V135F (GTT→TTT)  ycdT → diguanylate cyclase, membrane‑anchored

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,088,7340GT37.3% 348.3 / 280.1 271V135F (GTT→TTT) ycdTdiguanylate cyclase, membrane‑anchored
Reads supporting (aligned to +/- strand):  ref base G (82/88);  new base T (50/51);  total (132/139)
Fisher's exact test for biased strand distribution p-value = 9.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.24e-01

Only 100 of 271 total aligned reads displayed.

ATTTACTTTGTTATAACCGTCATTATCATTCAGCAAATTATTGAGGAGCGTTTGACAAGCAGTGTTGTCCAAAATGACATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTTTGTGCATAAGTTTGTTTTTTGTTTTTGGGGGGCCGATAGTTGCTCACATACTTTCTAGTCATTATGAGAGCTACAATTTACACATTGCCGAGTTAACC  >  CP009273/1088587‑1088880
                                                                                                                                                   |                                                                                                                                                   
 tttACTTTGTTATAACCGTCATTATCATTCAGCAAATTATTGAGGAGCGTTTGACAAGCAGTGTTGTCCAAAATGACATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTa                                                                                                                                                  >  2:4067334/1‑150 (MQ=255)
   tACTTTGTTATAACCGTCATTATCATTCAGCAAATTATTGAGGAGCGTTTGACAAGCAGTGTTGTCCAAAATGACATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGt                                                                                                                                                >  2:3937141/1‑150 (MQ=255)
   tACTTTGTTATAACCGTCATTATCATTCAGCAAATTATTGAGGAGCGTTTGACAAGCAGTGTTGTCCAAAATGACATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGt                                                                                                                                                >  2:2466166/1‑150 (MQ=255)
   tACTTTGTTATAACCGTCATTATCATTCAGCAAATTATTGAGGAGCGTTTGACAAGCAGTGTTGTCCAAAATGACATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGt                                                                                                                                                >  2:3937320/1‑150 (MQ=255)
     cTTTGTTATAACCGTCATTATCATTCAGCAAATTATTGAGGAGCGTTTGACAAGCAGTGTTGTCCAAAATGACATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGa                                                                                                                                              <  1:2606498/150‑1 (MQ=255)
        tGTTATAACCGTCATTATCATTCAGCAAATTATTGAGGAGCGTTTGACAAGCAGTGTTGTCCAAAATGACATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAATTTAGTGaaaa                                                                                                                                           >  2:3321619/1‑150 (MQ=255)
           tataACCGTCATTATCATTCAGCAAATTATTGAGGAGCGTTTGACAAGCAGTGTTGTCCAAAATGACATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAacac                                                                                                                                        >  1:3413786/1‑150 (MQ=255)
              aaCCGTCATTATCATTCAGCAAATTATTGAGGAGCGTTTGACAAGCAGTGTTGTCCAAAATGACATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAATTTAGTGAAAACACaaa                                                                                                                                     >  1:1354422/1‑150 (MQ=255)
              aaCCGTCATTATCATTCAGCAAATTATTGAGGAGCGTTTGACAAGCAGTGTTGTCCAAAATGACATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACaaa                                                                                                                                     >  2:3461380/1‑150 (MQ=255)
                ccGTCATTATCATTCAGCAAATTATTGAGGAGCGTTTGACAAGCAGTGTTGTCCAAAATGACATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAAc                                                                                                                                   <  1:2528875/150‑1 (MQ=255)
                 cGTCATTATCATTCAGCAAATTATTGAGGAGCGTTTGACAAGCAGTGTTGTCCAAAATGACATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAATTTAGTGAAAACACAAAACa                                                                                                                                  <  1:1143419/150‑1 (MQ=255)
                   tCATTATCATTCAGAAAATTATTGAGGAGCGTTTGACAAGCAGTGTTGTCCAAAATGACATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAATTTAGTGAAAACACAAAACAgc                                                                                                                                <  2:757047/150‑1 (MQ=255)
                         tCATTCAGCAAATTATTGAGGAGCGTTTGACAAGCAGTGTTGTCCAAAATGACATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAatt                                                                                                                          <  1:3484401/150‑1 (MQ=255)
                          cATTCAGCAAATTATTGAGGAGCGTTTGACAAGCAGTGTTGTCCAAAATGACATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAattt                                                                                                                         >  1:2421285/1‑150 (MQ=255)
                          cATTCAGCAAATTATTGAGGAGCGTTTGACAAGCAGTGTTGTCCAAAATGACATTGCAATCTATTATTTGTTTCGTCAGAGGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAattt                                                                                                                         >  1:3400173/1‑150 (MQ=255)
                              cAGCAAATTATTGAGGAGCGTTTGACAAGCAGTGTTGTCCAAAATGACATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAATTTAGTGAAAACACAAAACAGCGCAatttattt                                                                                                                     >  2:290482/1‑150 (MQ=255)
                               aGCAAATTATTGAGGAGCGTTTGACAAGCAGTGTTGTCCAAAATGACATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAATTTAGTGAAAACACAAAACAGCGCAATTTAtttt                                                                                                                    >  1:1872411/1‑150 (MQ=255)
                               aGCAAATTATTGAGGAGCGTTTGACAAGCAGTGTTGTCCAAAATGACATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTAtttt                                                                                                                    >  1:1689316/1‑150 (MQ=255)
                                 cAAATTATTGAGGAGCGTTTGACAAGCAGTGTTGTCCAAAATGACATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTcc                                                                                                                  >  1:3605781/1‑150 (MQ=255)
                                  aaaTTATTGAGGAGCGTTTGACAAGCAGTGTTGTCCAAAATGACATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCa                                                                                                                 <  2:200123/150‑1 (MQ=255)
                                    attattGAGGAGCGTTTGACAAGCAGTGTTGTCCAAAATGACATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCaaa                                                                                                               >  1:286652/1‑150 (MQ=255)
                                     ttattGAGGAGCGTTTGACAAGCAGTGTTGTCCAAAATGACATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCaaaa                                                                                                              >  1:2746700/1‑150 (MQ=255)
                                              agCGTTTGACAAGCAGTGTTGTCCAAAATGACATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAATTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACtt                                                                                                     <  2:1738119/150‑1 (MQ=255)
                                                cGTTTGACAAGCAGTGTTGTCCAAAATGACATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAATTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTTtg                                                                                                   >  1:2057818/1‑150 (MQ=255)
                                                cGTTTGACAAGCAGTGTTGTCCAAAATGACATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTTtg                                                                                                   <  1:1158188/150‑1 (MQ=255)
                                                 gTTTGACAAGCAGTGTTGTCCAAAATGACATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAATTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTTtgt                                                                                                  >  1:4487239/1‑150 (MQ=255)
                                                 gTTTGACAAGCAGTGTTGTCCAAAATGACATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTTtgt                                                                                                  >  1:3784837/1‑150 (MQ=255)
                                                 gTTTGACAAGCAGTGTTGTCCAAAATGACATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTTtgt                                                                                                  >  1:3783411/1‑150 (MQ=255)
                                                  tttGACAAGCAGTGTTGTCCAAAATGACATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTTtgtg                                                                                                 <  1:352957/150‑1 (MQ=255)
                                                   ttGACAAGCAGTGTTGTCCAAAATGACATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTTTGTGc                                                                                                <  1:1383407/150‑1 (MQ=255)
                                                    tGACAAGCAGTGTTGTCCAAAATGACATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAATTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTTTGTGCa                                                                                               <  2:1110731/150‑1 (MQ=255)
                                                     gACAAGCAGTGTTGTCCAAAATGACATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATTACTTTGTGCAt                                                                                              >  2:533857/1‑150 (MQ=255)
                                                     gACAAGCAGTGTTGTCCAAAATGACATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTTTGTGCAt                                                                                              >  2:533789/1‑150 (MQ=255)
                                                      aCAAGCAGTGTTGTCCAAAATGACATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTTTGTGCATa                                                                                             <  2:3414067/150‑1 (MQ=255)
                                                          gCAGTGTTGTCCAAAATGACATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAATTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTTTGTGCATAAgtt                                                                                         <  1:2131113/150‑1 (MQ=255)
                                                          gCAGTGTTGTCCAAAATGACATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAATTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTTTGTGCATAAgtt                                                                                         <  1:2131331/150‑1 (MQ=255)
                                                          gCAGTGTTGTCCAAAATGACAATGCAATCTATTATATGTTTCGTCAGAAGTGTTTGTGCATATTAATATTTCTGGCATGGGTGAATAAAGTTAGTGAAAACACAAAACAGAGCACTTTATTTTCCAAAAAAATGACGATGGGCCTAAgat                                                                                         >  2:1102601/1‑148 (MQ=255)
                                                            aGTGTTGTCCAAAATGACATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTTTGTGCATAAgtttg                                                                                       >  1:91892/1‑150 (MQ=255)
                                                             gtgtTGTCCAAAATGACATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAATTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTTTGTGCATAAgtttgt                                                                                      <  2:1381539/150‑1 (MQ=255)
                                                             gtgtTGTCCAAAATGACATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTGTGTGCATAAgtttgt                                                                                      >  2:1219908/1‑150 (MQ=255)
                                                             gtgtTGTACAAAATGACATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTTTGTGCATAAgtttgt                                                                                      <  2:1496695/150‑1 (MQ=255)
                                                               gttgtCCAAAATGACATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAATTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTTTGTGCATAAgtttgttt                                                                                    <  2:2529366/150‑1 (MQ=255)
                                                                ttgtCCAAAATGACATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAATTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTTTGTGCATAAGtttgtttt                                                                                   <  1:4296167/150‑1 (MQ=255)
                                                                ttgtCCAAAATGACATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTTTGTGCATAAGtttgtttt                                                                                   <  2:1702076/150‑1 (MQ=255)
                                                                 tgtCCAAAATGACATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTTTGTGCATAAGtttgttttt                                                                                  <  2:2272699/150‑1 (MQ=255)
                                                                    ccAAAATGACATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTTTGTGCATAAGtttgttttttgt                                                                              <  2:106020/150‑1 (MQ=255)
                                                                    ccAAAATGACACTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTTTGTGCATAAGtttgttgtttgt                                                                              >  2:4498549/1‑150 (MQ=255)
                                                                      aaaaTGACATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAATTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTTTGTGCATAAGtttgttttttgttt                                                                            <  1:2731451/150‑1 (MQ=255)
                                                                       aaaTGACATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAATTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTTTGTGCATAAGTttgttttttgtttt                                                                           >  2:3384924/1‑150 (MQ=255)
                                                                       aaaTGACATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTTTGTGCATAAGTttgttttttgtttt                                                                           >  2:868863/1‑150 (MQ=255)
                                                                          tGACATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTTTGTGCATAAGTTTGTTTTTTGTTTTTgg                                                                        <  1:3341120/150‑1 (MQ=255)
                                                                            aCATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTTTGTGCATAAGTTTGTTTTTTGTTTTTgggg                                                                      <  1:2826074/150‑1 (MQ=255)
                                                                            aCATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTTTGTGCATAAGTTTGTTTTTTGTTTTTgggg                                                                      <  2:1453785/150‑1 (MQ=255)
                                                                               ttGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTTATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTTTGTGCATAAGTTTGTTTTTTGTTTTTGGGGGGc                                                                   <  1:3461348/150‑1 (MQ=255)
                                                                               ttGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAATTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTTTGTGCATAAGTTTGTTTTTTGTTTTTGGGGGGc                                                                   >  2:329877/1‑150 (MQ=255)
                                                                               ttGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAATTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTTTGTGCATAAGTTTGTTTTTTGTTTTTGGGGGGc                                                                   >  2:329652/1‑150 (MQ=255)
                                                                                tGCAATCTATCATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTTTGTGCATAAGTTTGTTTTTTGTTTTTGGGGGGcc                                                                  >  2:3614400/1‑150 (MQ=255)
                                                                                 gCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAATTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTTTGTGCATAAGTTTGTTTTTTGTTTTTGGGGGGCCg                                                                 >  1:3408210/1‑150 (MQ=255)
                                                                                  cAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCGGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTTTGTGCATAAGTTTGTTTTTTGTTTTTGGGGGGCCGa                                                                <  1:2558936/150‑1 (MQ=255)
                                                                                    aTCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTTTGTGCATAAGTTTGTTTTTTGTTTTTGGGGGGCCGATa                                                              >  2:649039/1‑150 (MQ=255)
                                                                                    aTCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTTTGTGCATAAGTTTGTTTTTTGTTTTTGGGGGGCCGATa                                                              >  2:4127783/1‑150 (MQ=255)
                                                                                     tCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTTTGTGCATAAGTTTGTTTTTTGTTTTTGGGGGGCCGATAg                                                             <  2:807597/150‑1 (MQ=255)
                                                                                       tattatTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTTTGTGCATAAGTTTGTTTTTTGTTTTTGGGGGGCCGATAGtt                                                           <  2:2937615/150‑1 (MQ=255)
                                                                                          tattTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAATTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTTTGTGCATAAGTTTGTTTTTTGTTTTTGGGGGGCCGATAGTTGCt                                                        <  1:138471/150‑1 (MQ=255)
                                                                                          tattTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTTTGTGCATACGTTTGTGTTTTGTTTTTGGGGGGCCGATAGTTGCt                                                        >  2:1570266/1‑150 (MQ=255)
                                                                                            tttGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTTTGTGCATAAGTTTGTTTTTTGTTTTTGGGGGGCCGATAGTTGCTca                                                      <  1:741086/150‑1 (MQ=255)
                                                                                             ttGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAATTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTTTGTGCATAAGTTTGTTTTTTGTTTTTGGGGGGCCGATAGTTGCTcac                                                     >  1:173949/1‑150 (MQ=255)
                                                                                              tGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACATCGCAATTAATTTTCCAAAAAAATGAATATGTGCATAAGTTTGTTTTTTGTTTTTGGGGGGCCGATAGATGCTcaca                                                    >  1:2903605/1‑150 (MQ=255)
                                                                                                      cAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAATTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTTTGTGCATAAGTTTGTTTTTTGTTTTTGGGGGGCCGATAGTTGCTCACATACTTTCt                                            <  2:1297640/150‑1 (MQ=255)
                                                                                                      cAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTTTGTGCATAAGTTTGTTTTTTGTTTTTGGGGGGCCGATAGTTGCTCACATACTTTCt                                            <  1:2242434/150‑1 (MQ=255)
                                                                                                       aGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAATTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTTTGTGCATAAGTTTGTTTTTTGTTTTTGGGGGGCCGATAGTTGCTCACATACTTTCTa                                           >  2:3750062/1‑150 (MQ=255)
                                                                                                            aGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAATTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTTTGTGCATAAGTTTGTTTTTTGTTTTTGGGGGGCCGATAGTTGCTCACATACTTTCTAGTCAt                                      >  2:678734/1‑150 (MQ=255)
                                                                                                             gTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTTTGTGCATAAGTTTGTTTTTTGTTTTTGGGGGGCCGATAGTTGCTCACATACTTTCTAGTCAtt                                     >  2:1558020/1‑150 (MQ=255)
                                                                                                               ttGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTTTGTGCATAAGTTTGTTTTTTGTTTTTGGGGGGCCGATAGTTGCTCACATACTTTCTAGTCATTAt                                   <  1:3313544/150‑1 (MQ=255)
                                                                                                               ttGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTTTGTGCATAAGTTTGTTTTTTGTTTTTGGGGGGCCGATAGTTGCTCACATACTTTCTAGTCATTAt                                   >  1:4002357/1‑150 (MQ=255)
                                                                                                               ttGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTTTGTGCATAAGTTTGTTTTTTGTTTTTGGGGGGCCGATAGTTGCTCACATACTTTCTAGTCATTAt                                   >  1:4002291/1‑150 (MQ=255)
                                                                                                                   gCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTTTGTGCATAAGTTTGTTTTTTGTTTTTGGGGGGCCGATAGTTGCTCACATACTTTCTAGTCATTATgaga                               <  2:1862450/150‑1 (MQ=255)
                                                                                                                      tatTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTTTGTGCATAAGTTTGTTTTTTGTTTTTGGGGGGCCGATAGTTGCTCACATACTTTCTAGTCATTATGAGAGCt                            <  2:1385228/150‑1 (MQ=255)
                                                                                                                      tatTAATAATTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTTTGTGCATAAGTTTGTTTTTTGTTTTTGGGGGGCCGATAGTTGCTCACATACTTTCTAGTCATTATGAGAGCt                            <  2:4371182/150‑1 (MQ=255)
                                                                                                                             atTTCTGGCATTGGTGAATAAATTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTTTGTGCATAAGTTTGTTTTTTGTTTTTGGGGGGCCGATAGTTGCTCACATACTTTCTAGTCATTATGAGAGCTACAAttt                     <  1:1133901/150‑1 (MQ=255)
                                                                                                                             atTTCTGGCATTGGTGAATAAATTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTTTGTGCATAAGTTTGTTTTTTGTTTTTGGGGGGCCGATAGTTGCTCACATACTTTCTAGTCATTATGAGAGCTACAAttt                     <  1:2951716/150‑1 (MQ=255)
                                                                                                                              tttCTGGCATTGGTGAATAAATTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTTTGTGCATAAGTTTGTTTTTTGTTTTTGGGGGGCCGATAGTTGCTCACATACTTTCTAGTCATTATGAGAGCTACAATTTa                    <  2:2234661/150‑1 (MQ=255)
                                                                                                                               ttCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTTTGTGCATAAGTTTGTTTTTTGTTTTTGGGGGGCCGATAGTTGCTCACATACTTTCTAGTCATTATGAGAGCTACAATTTac                   <  2:1902313/150‑1 (MQ=255)
                                                                                                                                tCTGGCATTGGTGAATAAATTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTTTGTGCATAAGTTTGTTTTTTGTTTTTGGGGGGCCGATAGTTGCTCACATACTTTCTAGTCATTATGAGAGCTACAATTTaca                  >  1:1666248/1‑150 (MQ=255)
                                                                                                                                tCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTTTGTGCATAAGTTTGTTTTTTGTTTTTGGGGGGCCGATAGTTGCTCACATACTTTCTAGTCATTATGAGAGCTACAATTTaca                  <  2:3303786/150‑1 (MQ=255)
                                                                                                                                  tGGCATTGGTGAATAAATTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTTTGTGCATAAGTTTGTTTTTTGTTTTTGGGGGGCCGATAGTTGCTCACATACTTTCTAGTCATTATGAGAGCTACAATTTAcaca                >  2:1382113/1‑150 (MQ=255)
                                                                                                                                    gCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTTTGTGCATAAGTTTGTTTTTTGTTTTTGGGGGGCCGATAGTTGCTCACATACTTTCTAGTCATTATGAGAGCTACAATTTACACAtt              >  1:4461183/1‑150 (MQ=255)
                                                                                                                                         ggTGAATAAATTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTTTGTGCATAAGTTTGTTTTTTGTTTTTGGGGGGCCGATAGTTGCTCACATACTTTCTAGTCATTATGAGAGCTACAATTTACACATTGCCGa         <  1:3948992/150‑1 (MQ=255)
                                                                                                                                          gTGAATAAATTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTTTGTGCATAAGTTTGTTTTTTGTTTTTGGGGGGCCGATAGTTGCTCACATACTTTCTAGTCATTATGAGAGCTACAATTTACACATTGCCGAg        >  1:4087531/1‑150 (MQ=255)
                                                                                                                                          gTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTTTGTGCATAAGTTTGTTTTTTGTTTTTGGGGGGCCGATAGTTGCTCACATACTTTCTAGTCATTATGAGAGCTACAATTTACACATTGCCGAg        >  1:3042955/1‑150 (MQ=255)
                                                                                                                                          gTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTTTGTGCATAAGTTTGTTTTTTGTTTTTGGGGGGCCGATAGTTGCTCACATACTTTCTAGTCATTATGAGAGCTACAATTTACACATTGCCGAg        >  2:4319895/1‑150 (MQ=255)
                                                                                                                                            gAATAAATTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTTTGTGCATAAGTTTGTTTTTTGTTTTTGGGGGGCCGATAGTTGCTCACATACTTTCTAGTCATTATGAGAGCTACAATTTACACATTGCCGAGtt      >  2:273287/1‑150 (MQ=255)
                                                                                                                                            gAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTTTGTGCATAAGTTTGTTTTTTGTTTTTGGGGGGCCGATAGTTGCTCACATACTTTCTAGTCATTATGAGAGCTACAATTTACACATTGCCGAGtt      >  1:3299797/1‑150 (MQ=255)
                                                                                                                                              aTAAATTTAGTTAAAACACAAAACAGCGCAATTTATTTTCCAAAAATATGACTTTGTGCATAAGTTTGTTTTTTGTTTTTGGGGGGCCGATAGTTGCTCAAATACTTTCTAGTCATTATGAGAGCTACAATTTACACATTGCCGAGTTaa    <  1:132272/150‑1 (MQ=255)
                                                                                                                                              aTAAATTTAGTGAAAACACAAAACAGCGCAATTTATTTTCAAAAAAAATGACTTTGTGCATAAGTTTGTTTTTTGTTTTTGGGGGGCCGATAGTTGCTCACATACTTTCTAGTCATTATGAGAGCTACAATTTACACATTGCCGAGTTaa    <  1:132986/150‑1 (MQ=255)
                                                                                                                                              aTAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTTTGTGCATAAGTTTGTTTTTTGTTTTTTGGGGGGCCGATAGTTGCTCACATACTTTCTAGTCATTATGAGAGCTACAATTTACACATTGCCGAGTTa     >  2:1581121/1‑150 (MQ=255)
                                                                                                                                              aTAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTTTGTGCATAAGTTTGTTTTTTGTTTTTGGGGGGCCGATAGTTGCTCACATACTTTCTAGTCATTATGAGAGCTACAATTTACACATTGCCGAGTTaa    <  2:1001392/150‑1 (MQ=255)
                                                                                                                                               tAAATTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTTTGTGCATAAGTTTGTTTTTTGTTTTTGGGGGGCCGATAGTTGCTCACATACTTTCTAGTCATTATGAGAGCTACAATTTACACATTGCCGAGTTAAc   >  2:2383166/1‑150 (MQ=255)
                                                                                                                                               tAAATTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTTTGTGCATAAGTTTGTTTTTTGTTTTTGGGGGGCCGATAGTTGCTCACATACTTTCTAGTCATTATGAGAGCTACAATTTACACATTGCCGAGTTAAc   >  1:651368/1‑150 (MQ=255)
                                                                                                                                                aaaTTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTTTGTGCATAAGTTTGTTTTTTGTTTTTGGGGGGCCGATAGTTGCTCACATACTTTCTAGTCATTATGAGAGCTACAATTTACACATTGCCGAGTTAAcc  >  2:2749541/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
ATTTACTTTGTTATAACCGTCATTATCATTCAGCAAATTATTGAGGAGCGTTTGACAAGCAGTGTTGTCCAAAATGACATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTTTGTGCATAAGTTTGTTTTTTGTTTTTGGGGGGCCGATAGTTGCTCACATACTTTCTAGTCATTATGAGAGCTACAATTTACACATTGCCGAGTTAACC  >  CP009273/1088587‑1088880

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: