Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,576,745 G→T 20.4% pseudogene (348/381 nt) yrhA → pseudogene, interrupted by IS1E

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,576,7450GT20.4% 1106.3 / 189.6 421pseudogene (348/381 nt)yrhApseudogene, interrupted by IS1E
Reads supporting (aligned to +/- strand):  ref base G (165/170);  new base T (40/46);  total (205/216)
Fisher's exact test for biased strand distribution p-value = 7.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.81e-01

Only 100 of 431 total aligned reads displayed.

TTAAGAGCAATGGATGGATTCGAAGTTAATGGATTACGATTATTTAGCCTCTCGATTCCAGAACCTTCAGTTAAAAACCTTTTTGCCGTAAATGAATTTTATAGAAATAATGATGATTTCATAAACCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGACGGGGTGG  >  CP009273/3576596‑3576894
                                                                                                                                                     |                                                                                                                                                     
ttAAGAGCAATGGATGGATTCGAAGTTAATGGATTACGATTATTTAGCCTCTCGATTCCAGAACCTTCAGTTAAAAACCTTTTTGCCGTAAATGAATTTTATAGAAATAATGATGATTTCATAAACCCTGATCTACAAGAACGGTTAGTg                                                                                                                                                       <  2:5295551/150‑1 (MQ=255)
ttAAGAGCAATGGATGGATTCGAAGTTAATGGATTACGATTATTTAGCCTCTCGATTCCAGAACCTTCAGTTAAAAACCTTTTTGCCGTAAATGAATTTTATAGAAATAATGATGATTTCATAAACCCTGATCTACAAGAACGGTTAGTg                                                                                                                                                       <  2:4006081/150‑1 (MQ=255)
 tAAGAGCAATGGATGGATTCGAAGTTAATGGATTACGATTATTTAGCCTCTCGATTCCAGAACCTTCAGTTAAAAACCTTTTTGCCGTAAATGAATTTTATAGAAATAATGATGATTTCATAAACCCTGATCTACAAGAACGGTTAGTGa                                                                                                                                                      <  2:5068581/150‑1 (MQ=255)
   agagCAATGGATGGATTCGAAGTTAATGGATTACGATTATTTAGCCTCTCGATTCCAGAACCTTCAGTTAAAAACCTTTTTGCCGTAAATGAATTTTATAGAAATAATGATGATTTCATAAACCCTGATCTACAAGAACGGTTAGTTATc                                                                                                                                                    >  1:2295265/1‑150 (MQ=255)
   agagCAATGGATGGATTCGAAGTTAATGGATTACGATTATTTAGCCTCTCGATTCCAGAACCTTCAGTTAAAAACCTTTTTGCCGTAAATGAATTTTATAGAAATAATGATGATTTCATAAACCCTGATCTACAAGAACGGTTAGTGATc                                                                                                                                                    <  2:3912173/150‑1 (MQ=255)
   agagCAATGGATGGATTCGAAGTTAATGGATTACGATTATTTAGCCTCTCGATTCCAGAACCTTCAGTTAAAAACCTTTTTGCCGTAAATGAATTTTATAGAAATAATGATGATTTCATAAACCCTGATCTACAAGAACGGTTAGTGATc                                                                                                                                                    <  1:1978909/150‑1 (MQ=255)
       cAATGGATGGATTCGAAGTTAATGGATTACGATTATTTAGCCTCTCGATTCCAGAACCTTCAGTTAAAAACCTTTTTGCCGTAAATGAATTTTATAGAAATAATGATGATTTCATAAACCCTGATCTACAAGAACGGTTAGTGATCgggg                                                                                                                                                >  1:2205043/1‑150 (MQ=255)
         atggatggATTCGAAGTTAATGGATTACGATTATTTAGCCTCTCGATTCCAGAACCTTCAGTTAAAAACCTTTTTGCCGTAAATGAATTTTATAGAAATAATGATGATTTCATAAACCCTGATCTACAAGAACGGTTAGTGATCGGGGAt                                                                                                                                              <  1:709206/150‑1 (MQ=255)
           ggatggatTCGAAGTTAATGGATTACGATTATTTAGCCTCTCGATTCCAGAACCTTCAGTTAAAAACCTTTTTGCCGTAAATGAATTTTATAGAAATAATGATGATTTCATAAACCCTGATCTACAAGAACGGTTAGTGATCGGGGATta                                                                                                                                            >  1:4255417/1‑150 (MQ=255)
             atggatTCGAAGTTAATGGATTACGATTATTTAGCCTCTCGAGTCCAGAACCTTCAGTTAAAAACCTTTTTGCCGTAAATGAATTTTATAGAAATAATGATGATTTCATAAACCCTGATCTACAAGAACGGTTAGTTATCGGGGATtata                                                                                                                                          <  2:3392630/150‑1 (MQ=255)
              tggatTCGAAGTTAATGGATTACGATTATTTAGCCTCTCGATTCCAGAACCTTCAGTTAAAAACCTTTTTGCCGTAAATGAATTTTATAGAAATAATGATGATTTCATAAACCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAg                                                                                                                                         <  2:2579749/150‑1 (MQ=255)
              tggatTCGAAGTTAATGGATTACGATTATTTAGCCTCTCGATTCCAGAACCTTCAGTTAAAAACCTTTTTGCCGTAAATGAATTTTATAGAAATAATGATGATTTCATAAACCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAg                                                                                                                                         <  2:1575964/150‑1 (MQ=255)
              tggatTCGAAGTTAATGGATTACGATTATTTAGCCTCTCGATTCCAGAACCTTCAGTTAAAAACCTTTTTGCCGTAAATGAATTTTATAGAAATAATGATGATTTCATAAACCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAg                                                                                                                                         <  1:5251641/150‑1 (MQ=255)
                 atTCGAAGTTAATGGATTACGATTATTTAGCCTCTCGATTCCAGAACCTTCAGTTAAAAACCTTTTTGCCGTAAATGAATTTTATAGAAATAATGATGATTTCATAAACCCTGATCTACAAGAACGGTTAGTTATCGGGGATTATAGCAt                                                                                                                                      <  1:122543/150‑1 (MQ=255)
                 atTCGAAGTTAATGGATTACGATTATTTAGCCTCTCGATTCCAGAACCTTCAGTTAAAAACCTTTTTGCCGTAAATGAATTTTATAGAAATAATGATGATTTCATAAACCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCAt                                                                                                                                      <  1:1515907/150‑1 (MQ=255)
                      aaGTTAATGGATTACGATTATTTAGCCTCTCGATTCCAGAACCTTCAGTTAAAAACCTTTTTGCCGTAAATGAATTTTATAGAAATAATGATGATTTCATAAACCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCaa                                                                                                                                 >  2:5160398/1‑150 (MQ=255)
                         ttAATGGATTACGATTATTTAGCCTCTCGATTCCAGAACCTTCAGTTAAAAACCTTTTTGCCGTAAATGAATTTTATAGAAATAATGATGATTTCATAAACCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAatat                                                                                                                              <  2:1701957/150‑1 (MQ=255)
                          tAATGGATTACGATTATTTAGCCTCTCGATTCCAGAACCTTCAGTTAAAAACCTTTTTGCCGTAAATGAATTTTATAGAAATAATGATGATTTCATAAACCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATAtt                                                                                                                             >  2:1298913/1‑150 (MQ=255)
                            aTGGATTACGATTATTTAGCCTCTCGATTCCAGAACCTTCAGTTAAAAACCTTTTTGCCGTAAATGAATTTTATAGAAATAATGATGATTTCATAAACCCTGATCTACAAGAACGGTTAGTTATCGGGGATTATAGCATTTCAATATTTa                                                                                                                           >  1:1173122/1‑150 (MQ=255)
                             tGGATTACGATTATTTAGCCTCTCGATTCCAGAACCTTCAGTTACAAACCTTTTTGCCGTAAATGAATTTTATAGAAATAATGATGATTTCATAAACCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTAc                                                                                                                          <  2:4042517/150‑1 (MQ=255)
                              ggATTACGATTATTTAGCCTCTCGATTCCAGAACCTTCAGTTAAAAACCTTTTTGCCGTAAATGAATTTTATAGAAATAATGATGATTTCATAAACCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACt                                                                                                                         <  2:441468/150‑1 (MQ=255)
                              ggATTACGATTATTTAGCCTCTCGATTCCAGAACCTTCAGTTAAAAACCTTTTTGCCGTAAATGAATTTTATAGAAATAATGATGATTTCATAAACCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACt                                                                                                                         <  2:1158414/150‑1 (MQ=255)
                                 ttACGATTATTTAGCCTCTCGATTCCAGAACCTTCAGTTAAAAACCTTTTTGCCGTAAATGAATTTTATAGAAATAATGATGATTTCATAAACCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTAt                                                                                                                      <  2:5024454/150‑1 (MQ=255)
                                     gATTATTTAGCCTCTCGATTCCAGAACCTTCAGTTAAAAACCTTTTTGCCGTAAATGAATTTTATAGAAATAATGATGATTTCATAAACCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACa                                                                                                                  >  2:1812226/1‑150 (MQ=255)
                                     gATTATTTAGCCTCTCGATTCCAGAACCTTCAGTTAAAAACCTTTTTGCCGTAAATGAATTTTATAGAAATAATGATGATTTCATAAACCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACa                                                                                                                  <  1:4669468/150‑1 (MQ=255)
                                     gATTATTTAGCCTCTCGATTCCAGAACCTTCAGTTAAAAACCTTTTTGCCGTAAATGAATTTTATAGAAATAATGATGATTTCATAAACCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACa                                                                                                                  >  2:4994486/1‑150 (MQ=255)
                                       ttattTAGAGTCTCGAGTCCAGAACCTTCAGTTAAAAATGTTTTTGCAGTTAATGTATTTTATAGTTATATTGATGATTTCATAAACCCTGTTCTACAAGATCGGTTAGTGATCGGGGATTATAGCATTTCAATATTTTCTTATGTCAtt                                                                                                                <  1:2442112/150‑1 (MQ=255)
                                            tAGCCTCTCGATTCCAGAACCTTCAGTTAAAAACCTTTTTGCCGTAAATGAATTTTATAGAAATAATGATGATTTCATAAACCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAgg                                                                                                           >  1:548093/1‑150 (MQ=255)
                                              gCCTCTCGATTCCAGAACCTTCAGTTAAAAACCTTTTTGCCGTAAATGAATTTTATAGAAATAATGATGATTTCATAAACCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTg                                                                                                         <  2:1602628/150‑1 (MQ=255)
                                                   tcGATTCCAGAACCTTCAGTTAAAAACCTTTTTGCCGTAAATGAATTTTATAGAAATAATGATGATTTCATAAACCCTGATCTACAAGAACGGTTAGTTATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGAtgct                                                                                                    <  1:5155761/150‑1 (MQ=255)
                                                      aTTCCAGAACCTTCAGTTAAAAACCTTTTTGCCGTAAATGAATTTTATAGAAATAATGATGATTTCATAAACCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGcc                                                                                                 <  1:4343022/150‑1 (MQ=255)
                                                      aTTCCAGAACCTTCAGTTAAAAACCTTTTTGCCGTAAATGAATTTTATAGAAATAATGATGATTTCATAAACCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGcc                                                                                                 >  1:4654412/1‑150 (MQ=255)
                                                      aTTCCAGAACCTTCAGTTAAAAACCTTTTTGCCGTAAATGAATTTTATAGAAATAATGATGATTTCATAAACCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGATGcc                                                                                                 >  2:536417/1‑150 (MQ=255)
                                                        tCCAGAACCTTCAGTTAAAAACCTTTTTGCCGTAAATGAATTTTATAGAAATAATGATGATTTCATAAACCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCaa                                                                                               >  2:4439669/1‑150 (MQ=255)
                                                           aGAACCTTCAGTTAAAAACCTTTTTGCCGTAAATGAATTATATAGAAATAATGATGATTTCATAAACCCTGATCTACAAGAACGGTTAGAGATCGGGGATTATAGCATTTCAATATTTACTTATGACATAAAATGTGATTCTGACAACtt                                                                                            >  2:2834411/1‑150 (MQ=255)
                                                             aaCCTTCAGTTAAAAACCTTTTTGCCGTAAATGAATTTTATAGAAATAATGATGATTTCATAAACCCTGATCTACAAGAACGGTTAGTTATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTAc                                                                                          >  2:375251/1‑150 (MQ=255)
                                                               ccTTCAGTTAAAAACCTTTTTGCCGTAAATGAATTTTATAGAAATAATGATGATTTCATAAACCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTg                                                                                        <  1:4584197/150‑1 (MQ=255)
                                                                 ttCAGTTAAAAACCTTTTTGCCGTAAATGAATTTTATAGAAATAATGATGATTTCATAAACCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGAt                                                                                      <  1:1696092/150‑1 (MQ=255)
                                                                  tCAGTTAAAAACCTTTTTGCCGTAAATGAATTTTATAGAAATAATGATGATTTCATAAACCCTGATCTACAAGAACGGTTAGTTATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGAtt                                                                                     <  2:207136/150‑1 (MQ=255)
                                                                   cAGTTAAAAACCTTTTTGCCGTAAATGAATTTTATAGAAATAATGATGATTTCATAAACCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGAttt                                                                                    <  1:4742394/150‑1 (MQ=255)
                                                                     gTTAAAAACCTTTTTGCCGTAAATTAATTTTATAGAAATAATGATGATTTCATAAACCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAg                                                                                  <  2:3270596/150‑1 (MQ=255)
                                                                     gTTAAAAACCTTTTTGCCGTAAATGAATTTTATAGAAATAATGATGATTTCATAAACCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAg                                                                                  >  2:1465492/1‑150 (MQ=255)
                                                                        aaaaaCCTTTTTGCCGTAAATGAATTTTATAGAAATAATGATGATTTCATAAACCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAgtgt                                                                               <  2:2477713/150‑1 (MQ=255)
                                                                        aaaaaCCTTTTTGCCGTAAATGAATTTTATAGAAATAATGATGATTTCATAAACCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAgtgt                                                                               >  1:4041958/1‑150 (MQ=255)
                                                                          aaaCCTTTTTGCCGTAAATGAATTTTATAGAAATAATGATGATTTCATAAACCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTat                                                                             >  2:4336678/1‑150 (MQ=255)
                                                                           aaCCTTTTTGCCGTAAATGAATTTTATAGAAATAATGATGATTTCATAAACCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTatg                                                                            >  1:4455376/1‑150 (MQ=255)
                                                                           aaCCTTTTTGCCGTAAATGAATTTTATAGAAATAATGATGATTTCATAAACCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTatg                                                                            >  1:3897252/1‑150 (MQ=255)
                                                                             catttttGCCGTAAATGAATTTTATAGAAATAATGATGATTTCATAAACCCTGATCTACAAGAACGGTTAGTTATCGGGGATTATAGCTTTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTatgat                                                                          <  1:1397830/148‑1 (MQ=255)
                                                                              cTTTTTGCCGTAAATGAATTTTATAGAAATAATGATGATTTCATAAACCCTGATCTACAAGAACGGTTAGTTATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTatgatg                                                                         >  1:1015249/1‑150 (MQ=255)
                                                                              cTTTTTGCCGTAAATGAATTTTATAGAAATAATGATGATTTCATAAACCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGAGTTAGTGTatgatg                                                                         >  1:4526943/1‑150 (MQ=255)
                                                                                    gCCGTAAATGAATTTTATAGAAATAATGATGATTTCATAAACCCTGATCTACAAGAACGGTTAGTTATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGttt                                                                   >  1:873049/1‑150 (MQ=255)
                                                                                      cGTAAATGAATTTTATAGAAATAATGATGATTTCATAAACCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGttttt                                                                 >  2:4197376/1‑150 (MQ=255)
                                                                                      cGTAAATGAATTTTATAGAAATAATGATGATTTCATAAACCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGttttt                                                                 <  1:306515/150‑1 (MQ=255)
                                                                                       gTAAATGAATTTTATAGAAATAATGATGATTTCATAAACCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTg                                                                >  2:899708/1‑150 (MQ=255)
                                                                                         aaaTGAATTTTATAGAAATAATGATGATTTCATAAACCCTGATCTACAAGAACGGGTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAg                                                              >  2:2667496/1‑150 (MQ=255)
                                                                                          aaTGAATTTTATAGAAATAATGATGATTTCATAAACCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAgg                                                             >  2:3666315/1‑150 (MQ=255)
                                                                                           aTGAATTTTATAGAAATAATGATGATTTCATAAACCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGt                                                            >  1:3336399/1‑150 (MQ=255)
                                                                                           aTGAATTTTATAGAAATAATGATGATTTCATAAACCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGt                                                            <  2:1877935/150‑1 (MQ=255)
                                                                                           aTGAATTTTATAGAAATAATGATGATTTCATAAACCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGt                                                            >  1:937240/1‑150 (MQ=255)
                                                                                           aTGAATTTTATAGAAATAATGATGATTTCATAAACCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTGAGTGTATGATGGTGTTTTTGAGGt                                                            >  2:3537140/1‑150 (MQ=255)
                                                                                              aaTTTTATAGAAATAATGATGATTTCATAAACCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCt                                                         >  2:4512536/1‑150 (MQ=255)
                                                                                               aTTTTATAGAAATAATGATGATTTCATAAACCCTGATCTACAAGAACGGTTAGTTATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTc                                                        >  1:4969659/1‑150 (MQ=255)
                                                                                                 tttATAGAAATAATGATGATTTCATAAACCCTGATCTACAAGAACGGTTAGTTATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCa                                                      <  1:2535267/150‑1 (MQ=255)
                                                                                                 tttATAGAAATAATGATGATTTCATAAACCCTGATCTACAAGAACGGTTAGTTATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCa                                                      <  1:2495089/150‑1 (MQ=255)
                                                                                                 tttATAGAAATAATGATGATTTCATAAACCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCa                                                      <  1:3445960/150‑1 (MQ=255)
                                                                                                  ttATAGAAATAATGATGATTTCATAAACCCTGATCTACAAGAACGGTTAGTTATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAg                                                     >  1:4407863/1‑150 (MQ=255)
                                                                                                    ataGAAATAATGATGATTTCATAAACCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTg                                                   <  1:2729393/150‑1 (MQ=255)
                                                                                                       gAAATAATGATGATTTCATAAACCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCt                                                >  2:690152/1‑150 (MQ=255)
                                                                                                         aataatGATGATTTCATAAACCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTc                                              >  2:1089246/1‑150 (MQ=255)
                                                                                                          ataatGATGATTTCATAAACCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCt                                             >  1:1616269/1‑150 (MQ=255)
                                                                                                           taatGATGATTTCATAAACCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTg                                            >  2:2946883/1‑150 (MQ=255)
                                                                                                              tgatgaTTTCATAAACCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGttt                                         <  1:3720313/150‑1 (MQ=255)
                                                                                                              tgatgaTTTCATAAACCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGttt                                         >  2:2122993/1‑150 (MQ=255)
                                                                                                                   atTTCATAAACCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATc                                    <  1:2348972/150‑1 (MQ=255)
                                                                                                                   atTTCATAAACCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCa                                   >  2:3467402/1‑150 (MQ=255)
                                                                                                                     ttCATAAACCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAg                                  <  1:1893156/150‑1 (MQ=255)
                                                                                                                      tCATAAACCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCatc                                 >  2:1893162/1‑148 (MQ=255)
                                                                                                                        aTAAACCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTTTTTCTATAAGCTg                               >  1:2939213/1‑150 (MQ=255)
                                                                                                                        aTAAACCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTg                               >  2:4912760/1‑150 (MQ=255)
                                                                                                                          aaaCCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTc                             >  1:2987351/1‑150 (MQ=255)
                                                                                                                          aaaCCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTc                             >  1:2642727/1‑150 (MQ=255)
                                                                                                                           aaCCCTGATCTACAAGAACGGTTAGTTATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTcc                            >  1:1786057/1‑150 (MQ=255)
                                                                                                                             cccTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCcct                          >  2:4752231/1‑150 (MQ=255)
                                                                                                                             cccTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCcct                          <  1:4006199/150‑1 (MQ=255)
                                                                                                                              ccTGATCTACAAGAACGGTTAGTTATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCcctc                         >  2:2274374/1‑150 (MQ=255)
                                                                                                                               cTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCcctcc                        >  2:3107996/1‑150 (MQ=255)
                                                                                                                                 gATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTg                      <  2:3311063/150‑1 (MQ=255)
                                                                                                                                  aTCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGt                     >  1:3804291/1‑150 (MQ=255)
                                                                                                                                        aaGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCt               >  2:4646805/1‑150 (MQ=255)
                                                                                                                                        aaGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGAGTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCt               >  2:846334/1‑150 (MQ=255)
                                                                                                                                         aGAACTGTTAGTGATCGGTGATTATAGAATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTa              <  1:3214608/150‑1 (MQ=255)
                                                                                                                                           aaCGGTTAGTTATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACt            <  1:3437989/150‑1 (MQ=255)
                                                                                                                                            aCGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTg           <  2:571759/150‑1 (MQ=255)
                                                                                                                                            aCGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGAAATTAAAGGTGATGCTGCCAACTCACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTg           <  2:2762796/150‑1 (MQ=255)
                                                                                                                                              ggTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGAc         <  2:4287088/150‑1 (MQ=255)
                                                                                                                                              ggTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCAACTGAc         >  1:58585/1‑150 (MQ=255)
                                                                                                                                               gTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGACg        <  2:2648712/150‑1 (MQ=255)
                                                                                                                                               gTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGACg        <  2:2008535/150‑1 (MQ=255)
                                                                                                                                                    tGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGACGGggtg   <  1:4941568/150‑1 (MQ=255)
                                                                                                                                                     gATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGACGGggtgg  >  1:2387713/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TTAAGAGCAATGGATGGATTCGAAGTTAATGGATTACGATTATTTAGCCTCTCGATTCCAGAACCTTCAGTTAAAAACCTTTTTGCCGTAAATGAATTTTATAGAAATAATGATGATTTCATAAACCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGACGGGGTGG  >  CP009273/3576596‑3576894

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: