Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 733,047 G→A 22.9% pseudogene (767/1137 nt) ybfL → pseudogene, DDE domain transposase family;putative factor; Not classified; putative receptor protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273733,0470GA22.9% 322.1 / 96.8 153pseudogene (767/1137 nt)ybfLpseudogene, DDE domain transposase family;putative factor; Not classified; putative receptor protein
Reads supporting (aligned to +/- strand):  ref base G (59/57);  new base A (17/18);  total (78/75)
Fisher's exact test for biased strand distribution p-value = 8.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

Only 100 of 158 total aligned reads displayed.

TGAGGAAAAATTTCCGCTGAAAGAATTAAATAATCCAGAGCATGACAGTTACGCAATGAGTGAAAAGAGTCACGGCAGAGAAGAAATCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGACTGAAGAAATTATGCGTGGCAGTCTCCTTTCGGTCAATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCAGAAAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTG  >  CP009273/732898‑733195
                                                                                                                                                     |                                                                                                                                                    
   tgAAAAATTTCCGCTGAAAGAATTAAATAATCCAGAGCATGACAGTTACGCAATGAGTGAAAAGAGTCACGGCAGAGAAGAAATCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGACt                                                                                                                                                   <  1:4046868/149‑1 (MQ=21)
   gtaaaaaTTTCCGCTGAAAGAATTAAATAATCCAGAGCATGACAGTTACGCAATGAGTGAAAAGAGTCACGGCAGAGAAGAAATCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGACt                                                                                                                                                   <  1:2802737/148‑1 (MQ=21)
   ggAAAAATTTCCGCTGAAAGAATTAAATAATCCAGAGCATGACAGTTACGCAATGAGTGAAAAGAGTCACGGCAGAGAAGAAATCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGACt                                                                                                                                                   <  1:4046891/150‑1 (MQ=34)
    gAAAAATTTCCGCTGAAAGAATTAAATAATCCAGAGCATGACAGTTACGCAATGAGTGAAAAGAGTCACGGCAGAGAAGAAATCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGTACTg                                                                                                                                                  >  2:4272832/1‑150 (MQ=18)
     aaaaaTTTCCGCTGAAAGAATTAAATAATCCAGAGCATGACAGTTACGCAATGAGTGAAAAGAGTCACGGCAGAGAAGAAATCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGACTga                                                                                                                                                 >  1:1743758/1‑150 (MQ=34)
      aaaaTTTCCGCTGAAAGAATTAAATAATCCAGAGCATGACAGTTACGCAATGAGTGAAAAGAGTCACGGCAGAGAAGAAATCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGACTgaa                                                                                                                                                <  1:3057362/150‑1 (MQ=34)
      aaaaTTTCCGCTGAAAGAATTAAATAATCCAGAGCATGACAGTTACGCAATGAGTGAAAAGAGTCACGGCAGAGAAGAAATCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGAACTgaa                                                                                                                                                <  2:4371785/150‑1 (MQ=21)
      aaaaTTTCCGCTGAAAGAATTAAATAATCCAGAGCATGACAGTTACGCAATGAGTGAAAAGAGTCACGGCAGAGAAGAAATCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGAACTgaa                                                                                                                                                <  2:1920050/150‑1 (MQ=21)
       aaaTTTCCGCTGAAAGAATTAAATAATCCAGAGCATGACAGTTACGCAATGAGTGAAAAGAGTCACGGCAGAGAAGAAATCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATAGATTTCACGTTTGAATGGAAAGTACTgaag                                                                                                                                               >  1:4156586/1‑150 (MQ=21)
       aaaTTTCCGCTGAAAGAATTAAATAATCCAGAGCATGACAGTTACGCAATGAGTGAAAAGAGTCACGGCAGAGAAGAAATCCGTCTTCATATTGTTTGCGATGTCACTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGACTgaag                                                                                                                                               >  1:1660084/1‑150 (MQ=21)
          tttCCGCTGAAAGAATTAAATAATCCAGAGCATGACAGTTACGCAATGAGTGAAAAGAGTCACGGCAGAGAAGAAATCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGTAAGGACTGAAGaaa                                                                                                                                            >  1:4397435/1‑150 (MQ=21)
           ttCCGCTGAAAGAATTAAATAATCCAGAGCATGACAGTTACGCAATGAGTGAAAAGAGTCACGGCAGAGAAGAAATCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGACTGAAGAAAt                                                                                                                                           >  2:757012/1‑150 (MQ=34)
           ttCCGCTGAAAGAATTAAATAATCCAGAGCATGACAGTTACGCAATGAGTGAAAAGAGTCACGGCAGAGAAGAAATCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGACTGAAGAAAt                                                                                                                                           >  2:3608225/1‑150 (MQ=34)
           ttCCGCTGAAAGAATTAAATAATCCAGAGCATGACAGTTACGCAATGAGTGAAAAGAGTCACGGCAGAGAAGAAATCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGACTGAAGAAAt                                                                                                                                           >  2:103229/1‑150 (MQ=34)
              cGCTGAAAGAATTAAATAATCCAGAGCATGACAGTTACGCAATGAGTGAAAAGAGTCACGGCAGAGAAGAAATCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGACTGAAGAAATTAt                                                                                                                                        <  2:968860/150‑1 (MQ=34)
              cGCTGAAAGAATTAAATAATCCAGAGCATGACAGTTACGCAATGAGTGAAAAGAGTCACGGCAGAGAAGAAATCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGACTGAAGAAATTAt                                                                                                                                        >  1:3579481/1‑150 (MQ=34)
                cTGAAAGAATTAAATAATCCAGAGCATGACAGTTACGCAATGAGTGAAAAGAGTCACGGCAGAGAAGAAATCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGACTGAAGAAATTATGc                                                                                                                                      >  1:2348805/1‑150 (MQ=34)
                cTGAAAGAATTAAATAATCCAGAGCATGACAGTTACGCAATGAGTGAAAAGAGTCACGGCAGAGAAGAAATCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGACTGAAGAAATTATGc                                                                                                                                      >  1:2348702/1‑150 (MQ=34)
                cTGAAAGAATTAAATAATCCAGAGCATGACAGTTACGCAATGAGTGAAAAGAGTCACGGCAGAGAAGAAATCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGACTGAAGAAATTATGc                                                                                                                                      >  2:428110/1‑150 (MQ=34)
                 tGAAAGAATTAAATAATCCAGAGCATGACAGTTACGCAATGAGTGAAAAGAGTCACGGCAGAGAAGAAATCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGACTGAAGAAATTATGCg                                                                                                                                     <  2:1194360/150‑1 (MQ=34)
                  gAAAGAATTAAATAATCCAGAGCATGACAGTTACGCAATGAGTGAAAAGAGTCACGGCAGAGAAGAAATCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGACTGAAGAAATTATGCGt                                                                                                                                    <  1:3561913/150‑1 (MQ=34)
                  gAAAGAATTAAATAATCCAGAGCATGACAGTTACGCAATGAGTGAAAAGAGTCACGGCAGAGAAGAAATCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGACTGAAGAAATTATGCGt                                                                                                                                    >  1:4665528/1‑150 (MQ=34)
                  gAAAGAATTAAATAATCCAGAGCATGACAGTTACGCAATGAGTGAAAAGAGTCACGGCAGAGAAGAAATCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGACTGAAGAAATTATGCGt                                                                                                                                    >  2:3946988/1‑150 (MQ=34)
                  gAAAGAATTAAATAATCCAGAGCATGACAGTTACGCAATGAGTGAAAAGAGTCACGGCAGAGAAGAAATCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGACTGAAGAAATTATGCGt                                                                                                                                    >  2:3947461/1‑150 (MQ=34)
                  gAAAGAATTAAATAATCCAGAGCATGACAGTTACGCAATGAGTGAAAAGAGTCACGGCAGAGAAGAAATCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGACTGAAGAAATTATGCGt                                                                                                                                    <  1:1105761/150‑1 (MQ=34)
                  gAAAGAATTAAATAATCCAGAGCATGACAGTTACGCAATGAGTGAAAAGAGTCACGGCAGAGAAGAAATCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGAACTGAAGAAATTATGCGt                                                                                                                                    >  2:1743702/1‑150 (MQ=21)
                   aaaGAATTAAATAATCCAGAGCATGACAGTTACGCAATGAGTGAAAAGAGTCACGGCAGAGAAGAAATCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGACTGAAGAAATTATGCGTg                                                                                                                                   >  1:4712489/1‑150 (MQ=34)
                        aTTAAATAATCCAGAGCATGACAGTTACGCAATGAGTGAAAAGAGTCACGGCAGAGAAGAAATCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGACTGAAGAAATTATGCGTGGCAGt                                                                                                                              <  2:2066361/150‑1 (MQ=34)
                         ttAAATAATCCAGAGCATGACAGTTACGCAATGAGTGAAAAGAGTCACGGCAGAGAAGAAATCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGACTGAAGAAATTATGCGTGGCAGtc                                                                                                                             >  1:3394721/1‑150 (MQ=34)
                             ataatCCAGAGCATGACAGTTACGCAATGAGTGAAAAGAGTCACGGCAGAGAAGAAATCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGAACTGAAGAAATTATGCGTGGCAGTCTCCt                                                                                                                         <  2:2958231/150‑1 (MQ=21)
                               aatCCAGAGCATGACAGTTACGCAATGAGTGAAAAGAGTCACGGCAGAGAAGAAATCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGACTGAAGAAATTATGCGTGGCAGTCTCCttt                                                                                                                       >  1:4769299/1‑150 (MQ=34)
                                atCCAGAGCATGACAGTTACGCAATGAGTGAAAAGAGTCACGGCAGAGAAGAAATCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGACTGAAGAAATTATGCGTGGCAGTCTCCTTTc                                                                                                                      >  1:1119428/1‑150 (MQ=34)
                                    agagCATGACAGTTACGCAATGAGTGAAAAGAGTCACGGCAGAGAAGAAATCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGAACTGAAGAAATTATGCGTGGCAGTCTCCTTTCGGTc                                                                                                                  <  2:4542209/150‑1 (MQ=21)
                                    agagCATGACAGTTACGCAATGAGTGAAAAGAGTCACGGCAGAGAAGAAATCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGAACTGAAGAAATTATGCGTGGCAGTCTCCTTTCGGTc                                                                                                                  <  2:4542285/150‑1 (MQ=21)
                                    agagCATGACAGTTACGCAATGAGTGAAAAGAGTCACGGCAGAGAAGAAATCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGAACTGAAGAAATTATGCGTGGCAGTCTCCTTTCGGTc                                                                                                                  >  2:502191/1‑150 (MQ=21)
                                         aTGACAGTTACGCAATGAGTGAAAAGAGTCACGGCAGAGAAGAAATCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGACTGAAGAAATTATGCGTGGCAGTCTCCTTTCGGTCaataa                                                                                                             <  1:2376529/150‑1 (MQ=34)
                                         aTGACAGTTACGCAATGAGTGAAAAGAGTCACGGCAGAGAAGAAATCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGACTGAAGAAATTATGCGTGGCAGTCTCCTTTCGGTCaataa                                                                                                             <  1:2374795/150‑1 (MQ=34)
                                          tGACAGTTACGCAATGAGTGAAAAGAGTCACGGCAGAGAAGAGATCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGACTGAAGAAATTATGCGTGGCAGTCTCCTTTCGGTCaataat                                                                                                            >  1:4279940/1‑150 (MQ=21)
                                          tGACAGTTACGCAATGAGTGAAAAGAGTCACGGCAGAGAAGAAATCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGACTGAAGAAATTATGCGTGGCAGTCTCCTTTCGGTCaataat                                                                                                            >  1:4279932/1‑150 (MQ=34)
                                          tGACAGTTACGCAATGAGTGAAAAGAGTCACGGCAGAGAAGAAATCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTGTGAATGGAAAGGACTGAAGAAATTATGCGTGGCAGTCTCCTTTCGGTCaataat                                                                                                            <  2:1351844/150‑1 (MQ=21)
                                             cAGTTACGCAATGAGTGAAAAGAGTCACGGCAGAGAAGAAATCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGAACTGAAGAAATTATGCGTGGCAGTCTCCTTTCGGTCAATAATAGc                                                                                                         >  2:4127831/1‑150 (MQ=21)
                                                 tACGCAATGAGTGAAAAGAGTCACGGCAGAGAAGAAATCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGACTGAAGAAATTATGCGTGGCAGTCTCCTTTCGGTCAATAATAGCAGaa                                                                                                     >  1:1253273/1‑150 (MQ=34)
                                                 tACGCAATGAGTGAAAAGAGTCACGGCAGAGAAGAAATCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGACTGAAGAAATTATGCGTGGCAGTCTCCTTTCGGTCAATAATAGCAGaa                                                                                                     >  1:1255524/1‑150 (MQ=34)
                                                   cGCAATGAGTGAAAAGAGTCACGGCAGAGAAGAAATCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGACTGAAGAAATTATGCGTGGCAGTCTCCTTTCGGTCAATAATAGCAGAACa                                                                                                   >  1:1913489/1‑150 (MQ=34)
                                                     cAATGAGTGAAAAGAGTCACGGCAGAGAAGAAATCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGACTGTAGAAATTATGCGTGGCAGTCTCCTTTCGGTCAATCTTAGCCGAACaca                                                                                                 >  2:4658664/1‑148 (MQ=255)
                                                     cAATGAGTGAAAAGAGTCACGGCAGAGAAGAAATCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGACTGAAGAAATTATGCGTGGCAGTCTCCTTTCGGTCAATAATAGCAGAACaaa                                                                                                 >  2:332437/1‑150 (MQ=34)
                                                      aaTGAGTGAAAAGAGTCACGGCAGAGAAGAAATCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGACTGAAGAAATTATGCGTGGCAGTCTCCTTTCGGTCAATAATCGCAGAACaaaa                                                                                                >  2:4552545/1‑150 (MQ=21)
                                                       aTGAGTGAAAAGAGTCACGGCAGAGAAGAAATCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGACTGAAGAAATTATGCGTGGCAGTCTCCTTTCGGTCAATAATAGCAGAACAaaag                                                                                               <  2:882826/150‑1 (MQ=34)
                                                          aGTGAAAAGAGTCACGGCAGAGAAGAAATCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGACTGAAGAAATTATGCGTGGCAGTCTCCTTTCGGTCAATAATAGCAGAACAaaagaaa                                                                                            <  2:3788904/150‑1 (MQ=34)
                                                            tGAAAAGAGTCACGGCAGAGAAGAAATCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGAACTGAAGAAATTATGCGTGGCAGTCTCCTTTCGGTCAATAATAGCAGAACAAaagaaaga                                                                                          <  2:40815/150‑1 (MQ=21)
                                                              aaaaGAGTCACGGCAGAGAAGAAATCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGAACTGAAGAAATTATGCGTGGCAGTCTCCTTTCGGTCAATAATAGCAGAACAAAAGAAAGAGc                                                                                        >  1:3013980/1‑150 (MQ=21)
                                                                aaGAGTCACGGCAGAGAAGAAATCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGACTGAAGAAATTATGCGTGGCAGTCTCCTTTCGGTCAATAATAGCAGAACAAAAGAAAGAGCCa                                                                                      >  2:1151606/1‑150 (MQ=34)
                                                                 agagTCACGGCAGAGAAGAAATCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGAACTGAAGAAATTATGCGTGGCAGTCTCCTTTCGGTCAATAATAGCAGAACAAAAGAAAGAGCCAg                                                                                     <  1:3506441/150‑1 (MQ=21)
                                                                     tCACGGCAGAGAAGAAATCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGAACTGAAGAAATTATGCGTGGCAGTCTCCTTTCGGTCAATAATAGCAGAACAAAAGAAAGAGCCAGAAAt                                                                                 <  2:4246834/150‑1 (MQ=21)
                                                                        cGGCAGAGAAGAAATCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGACTGAAGAAATTATGCGTGGCAGTCTCCTTTCGGTCAATAATAGCAGAACAAAAGAAAGAGCCAGAAATGAc                                                                              <  2:1913489/150‑1 (MQ=34)
                                                                         ggCAGAGAAGAAATCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGACTGAAGAAATTATGCGTGGCAGTCTCCTTTCGGTCAATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACg                                                                             >  2:3799200/1‑150 (MQ=34)
                                                                            agagAAGAAATCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGACTGAAGAAATTATGCGTGGCAGTCTCCTTTCGGTCAATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTc                                                                          >  2:4736122/1‑150 (MQ=34)
                                                                                   aaaTCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGACTGAAGAAATTATGCGTGGCAGTCTCCTTTCGGTCAATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGAtatt                                                                   <  1:2608641/150‑1 (MQ=34)
                                                                                     aTCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGACTGAAGAAATTATGCGTGGCAGTCTCCTTTCGGTCAATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGAtattat                                                                 <  1:3983875/150‑1 (MQ=34)
                                                                                      tCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGACTGAAGAAATTATGCGTGGCAGTCTCCTTTCGGTCAATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATtata                                                                >  1:1945111/1‑150 (MQ=34)
                                                                                        cGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGACTGAAGAAATTATGCGTGGCAGTCTCCTTTCGGTCAATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATc                                                              >  1:3697295/1‑150 (MQ=34)
                                                                                        cGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGACTGAAGAAATTATGCGTGGCAGTCTCCTTTCGGTCAATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATc                                                              >  1:3698791/1‑150 (MQ=34)
                                                                                           cTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGAACTGAAGAAATTATGCGTGGCAGTCTCCTTTCGGTCAATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGt                                                           <  2:3938070/150‑1 (MQ=21)
                                                                                            ttCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGACTGAAGAAATTATGCGTGGCAGTCTCCTTTCGGTCAATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGtt                                                          >  1:3892239/1‑150 (MQ=34)
                                                                                            ttCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGAACTGAAGAAATTATGCGTGGCAGTCTCCTTTCGGTCAATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGtt                                                          >  2:3160505/1‑150 (MQ=21)
                                                                                             tCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGACTGAAGAAATTATGCGTGGCAGTCTCCTTTCGGTCAATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTc                                                         <  1:828482/150‑1 (MQ=34)
                                                                                             tCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGAACTGAAGAAATTATGCGTGGCAGTCTCCTTTCGGTCAATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTc                                                         <  1:52212/150‑1 (MQ=21)
                                                                                              cATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGACTGAAGAAATTATGCGTGGCAGTCTCCTTTCGGTCAATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTct                                                        <  2:4569111/150‑1 (MQ=34)
                                                                                              cATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGACTGAAGAAATTATGCGTGGCAGTCTCCTTTCGGTCAATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTct                                                        <  2:4569201/150‑1 (MQ=34)
                                                                                                  ttGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGACTGAAGAAATTATGCGTGGCAGTCTCCTTTCGGTCAATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTctgctg                                                    <  2:1119428/150‑1 (MQ=34)
                                                                                                     tttGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGACTGAAGAAATTATGCGTGGCAGTCTCCTTTCGGTCAATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGAtt                                                 <  2:82050/150‑1 (MQ=34)
                                                                                                      ttGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGACTGAAGAAATTATGCGTGGCAGTCTCCTTTCGGTCAATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGAttt                                                <  1:1857682/150‑1 (MQ=34)
                                                                                                      ttGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGACTGAAGAAATTATGCGTGGCAGTCTCCTTTCGGTCAATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGAttt                                                <  1:4851385/150‑1 (MQ=34)
                                                                                                      ttGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGACTGAAGAAATTATGCGTGGCAGTCTCCTTTCGGTCAATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGAttt                                                <  1:4851132/150‑1 (MQ=34)
                                                                                                      ttGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGACTGAAGAAATTATGCGTGGCAGTCTCCTTTCGGTCAATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGAttt                                                <  1:757012/150‑1 (MQ=34)
                                                                                                      ttGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGACTGAAGAAATTATGCGTGGCAGTCTCCTTTCGGTCAATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGAttt                                                <  1:1857661/150‑1 (MQ=34)
                                                                                                      ttGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGAACTGAAGAAATTATGCGTGGCAGTCTCCTTTCGGTCAATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGAttt                                                <  1:1917105/150‑1 (MQ=21)
                                                                                                         cGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGACTGAAGAAATTATGCGTGGCAGTCTCCTTTCGGTCAATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAAc                                             <  2:4040663/150‑1 (MQ=34)
                                                                                                         cGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGACTGAAGAAATTATGCGTGGCAGTCTCCTTTCGGTCAATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAAc                                             <  2:4040691/150‑1 (MQ=34)
                                                                                                               cccTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGACTGAAGAAATTATGCGTGGCAGTCTCCTTTCGGTCAATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCAGa                                       <  2:232061/150‑1 (MQ=37)
                                                                                                               cccTGATGAACTTATTGATTTCACGTTTGAATGGAAAGAACTGAAGAAATTATGCGTGGCAGTCTCCTTTCGGTCAATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCAGa                                       >  2:2104366/1‑150 (MQ=37)
                                                                                                                ccTGATGAACTTATTGATTTCACGTTTGAATGGAAAGAACTGAAGAAATTATGCGTGGCAGTCTCCTTTCGGTCAATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCAGaa                                      <  2:2499052/150‑1 (MQ=38)
                                                                                                                 cTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGACTGAAGAAATTATGCGTGGCAGTCTCCTTTCGGTCAATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCAGaaa                                     <  1:1603581/150‑1 (MQ=38)
                                                                                                                 cTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGACTGAAGAAATTATGCGTGGCAGTCTCCTTTCGGTCAATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCAGaaa                                     >  1:2934314/1‑150 (MQ=38)
                                                                                                                 cTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGACTGAAGAAATTATGCGTGGCAGTCTCCTTTCGGTCAATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCAGaaa                                     <  1:4363281/150‑1 (MQ=38)
                                                                                                                  tgatgaACTTATTGATTTCACGTTTGAATGGAAAGGACTGAAGAAATTATGCGTGGCAGTCTCCTTTCGGTCAATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCAGaaaa                                    >  2:2237053/1‑150 (MQ=39)
                                                                                                                   gatgaACTTATTGATTTCACGTTTGAATGGAAAGGACTGAAGAAATTATGCGTGGCAGTCTCCTTTCGGTCAATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCAGAAAAg                                   >  2:2585299/1‑150 (MQ=40)
                                                                                                                           tATTGATTTCACGTTTGAATGGAAAGGACTGAAGAAATTATGCGTGGCAGTCTCCTTTCGGTCAATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCAGAAAAGTTCGCcac                           >  2:2877752/1‑150 (MQ=255)
                                                                                                                            aTTGATTTCACGTTTGAATGGAAAGGACTGAAGAAATTATGCGTGGCAGTCTCCTTTCGGTCAATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCAGAAAAGTTCGCcaca                          <  2:1410953/150‑1 (MQ=255)
                                                                                                                            aTTGATTTCACGTTTGAATGGAAAGGACTGAAGAAATTATGCGTGGCAGTCTCCTTTCGGTCAATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCAGAAAAGTTCGCcaca                          <  2:3628742/150‑1 (MQ=255)
                                                                                                                            aTTGATTTCACGTTTGAATGGAAAGGACTGAAGAAATTATGCGTGGCAGTCTCCTTTCGGTCAATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCAGAAAAGTTCGCcaca                          <  2:3292181/150‑1 (MQ=255)
                                                                                                                                     aCGTTTGAATGGAAAGGACTGAAGAAATTATGCGTGGCAGTCTCCTTTCGGTCAATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCAGAAAAGTTCGCCACAGCAATCCGa                 <  1:1483411/150‑1 (MQ=255)
                                                                                                                                         ttGAATGGAAAGGACTGAAGAAATTATGCGTGGCAGTCTCCTTTCGGTCAATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCAGAAAAGTTCGCCACAGCAATCCGAAAcc             <  2:4656163/150‑1 (MQ=255)
                                                                                                                                          tGAATGGAAAGGACTGAAGAAATTATGCGTGGCAGTCTCCTTTCGGTCAATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCAGAAAAGTTCGCCACAGCAATCCGAAACCa            >  1:3204576/1‑150 (MQ=255)
                                                                                                                                          tGAATGGAAAGGACTGAAGAAATTATGCGTGGCAGTCTCCTTTCGGTCAATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCAGAAAAGTTCGCCACAGCAATCCGAAACCa            >  2:2704301/1‑150 (MQ=255)
                                                                                                                                                 aaaGGACTGAAGAAATTATGCGTGGCAGTCTCCTTTCGGTCAATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCAGAAAAGTTCGCCACAGCAATCCGAAACCACTGGCAc     <  2:4397435/150‑1 (MQ=255)
                                                                                                                                                  aaGGACTGAAGAAATTATGCGTGGCAGTCTCCTTTCGGTCAATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCAGAAAAGTTCGCCACAGCAATCCGAAACCACTGGCACg    <  1:3790838/150‑1 (MQ=255)
                                                                                                                                                    ggACTGAAGAAATTATGCGTGGCAGTCTCCTTTCGGTCAATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCAGAAAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTg  <  1:1639573/150‑1 (MQ=255)
                                                                                                                                                    ggACTGAAGAAATTATGCGTGGCAGTCTCCTTTCGGTCAATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCAGAAAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTg  <  1:4096213/150‑1 (MQ=255)
                                                                                                                                                    ggACTGAAGAAATTATGCGTGGCAGTCTCCTTTCGGTCAATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCAGAAAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTg  <  1:4096227/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
TGAGGAAAAATTTCCGCTGAAAGAATTAAATAATCCAGAGCATGACAGTTACGCAATGAGTGAAAAGAGTCACGGCAGAGAAGAAATCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGACTGAAGAAATTATGCGTGGCAGTCTCCTTTCGGTCAATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCAGAAAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTG  >  CP009273/732898‑733195

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: