Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,577,827 A→G 5.6% intergenic (+158/‑292) yrhA → / → yrhB pseudogene, interrupted by IS1E/stable heat shock chaperone

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,577,8270AG 5.6% 618.5 / 20.1 195intergenic (+158/‑292)yrhA/yrhBpseudogene, interrupted by IS1E/stable heat shock chaperone
Reads supporting (aligned to +/- strand):  ref base A (88/96);  new base G (8/3);  total (96/99)
Fisher's exact test for biased strand distribution p-value = 1.30e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.01e-01

Only 100 of 199 total aligned reads displayed.

AGGAAACAAAAATCATTGTAAACCATGCCATCTTTTTAGAGATTACCAGCAGAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGATTTATTTAGAATTTTACCCGAATCATAAGATTCTGATATCTTAAATAAGAAAATCCCGTTCGACGTTCTGAAAAAGAAATCATAAAATATTAGCAGATACT  >  CP009273/3577678‑3577976
                                                                                                                                                     |                                                                                                                                                     
 ggAAACAAAAATCATTGTAAACCATGCCATCTTTTTAGAGATTACCAGCAGAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTAt                                                                                                                                                      >  2:1745586/1‑150 (MQ=255)
  gAAACAAAAATCATTGTAAACCATGCCATCTTTTTAGAGATTACCAGCAGAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTGtt                                                                                                                                                     >  1:14038/1‑150 (MQ=255)
    aaCAAAAATCATTGTAAACCATGCCATCTTTTTAGAGATTACCAGCAGAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTAtttt                                                                                                                                                   >  2:4765720/1‑150 (MQ=255)
    aaCAAAAATCATTGTAAACCATGCCATCTTTTTAGAGATTACCAGCAGAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTAtttt                                                                                                                                                   <  2:1125197/150‑1 (MQ=255)
    aaCAAAAATCATTGTAAACCATGCCATATTTTTATAGATTACCAGCAGAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTAACAACGGAACTACACTCCTTAAGTGAGATGGACAGAAGACATGTAGACTCTAtttt                                                                                                                                                   >  1:3727572/1‑150 (MQ=255)
       aaaaaTCATTGTAAACCATGCCATCTTTTTAGAGATTACCAGCAGAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTGTTTTAGc                                                                                                                                                >  2:2978190/1‑150 (MQ=255)
          aaTCATTGTAAACCATGCCATCTTTTTAGAGATTACCAGCAGAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGa                                                                                                                                             >  1:4139982/1‑150 (MQ=255)
          aaTCATTGTAAACCATGCCATCTTTTTAGAGATTACCAGCAGAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGa                                                                                                                                             <  1:4047509/150‑1 (MQ=255)
           aTCATTGTAAACCATGCCATCTTTTTAGAGATTACCAGCAGAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGaa                                                                                                                                            <  1:4721924/150‑1 (MQ=255)
             cATTGTAAACCATGCCATCTTTTTAGAGATTACCAGCAGAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAAt                                                                                                                                          >  1:1892196/1‑150 (MQ=255)
                  tAAACCATGCCATCTTTTTAGAGATTACCAGCAGAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTa                                                                                                                                     >  2:5104392/1‑150 (MQ=255)
                   aaaCCATGCCATCTTTTTAGAGATTACCAGCAGAAAATTCCTTGTTCGAAATTTATCAGTTCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTaa                                                                                                                                    <  1:1897085/150‑1 (MQ=255)
                   aaaCCATGCCATCTTTTTAGAGATTACCAGCAGAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTaa                                                                                                                                    >  1:5081521/1‑150 (MQ=255)
                   aaaCCATGCCATCTTTTTAGAGATTACCAGCAGAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTaa                                                                                                                                    <  1:1831650/150‑1 (MQ=255)
                   aaaCCATGCCATCTTTTTAGAGATTACCAGCAGAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACAGGTAGTCTCTGTTTTAGCCGAAATAGCTaa                                                                                                                                    >  2:3812048/1‑150 (MQ=255)
                       cATGCCATCTTTTTAGAGATTACCAGCAGAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGc                                                                                                                                >  1:2349970/1‑150 (MQ=255)
                        aTGCCATCTTTTTAGAGATTACCAGCAGAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCt                                                                                                                               >  1:2623404/1‑150 (MQ=255)
                        aTGCCATCTTTTTAGAGATTACCAGCAGAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCt                                                                                                                               <  1:4618870/150‑1 (MQ=255)
                        aTGCCATCTTTTTAGAGATTACCAGCAGAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCt                                                                                                                               >  2:3329628/1‑150 (MQ=255)
                         tGCCATCTTTTTAGAGATTACCAGCAGAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTa                                                                                                                              >  1:1597437/1‑150 (MQ=255)
                         tGCCATCTTTTTAGAGATTACCAGCAGAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTa                                                                                                                              >  1:3582075/1‑150 (MQ=255)
                              tCTTTTTAGAGATTACCAGCAGAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTatgatg                                                                                                                         <  2:4458815/150‑1 (MQ=255)
                                tttttAGAGATTACCAGCAGAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTATgatgat                                                                                                                       <  1:4946044/150‑1 (MQ=255)
                                 ttttAGATATTACCAGCAGAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTATGATGAtt                                                                                                                      >  1:2246323/1‑150 (MQ=255)
                                  tttAGAGATTACCAGCAGAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTGTTTTAGCCGAAATAGCTAATGGCTATGATGAttt                                                                                                                     >  2:4048835/1‑150 (MQ=255)
                                  tttAGAGATTACCAGCAGAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTATGATGAttt                                                                                                                     <  1:545609/150‑1 (MQ=255)
                                  tttAGAGATTACCAGCAGAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTATGATGAttt                                                                                                                     >  1:2189513/1‑150 (MQ=255)
                                  tttAGAGATTACCAGCAGAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTATGATGAttt                                                                                                                     <  2:4563044/150‑1 (MQ=255)
                                    ttgcgaTTACCCGCTGAAAATTCCTTGTTCGATATTTAGCAGTGCAAGCTTTGGAAAAACTTACATCGGAACTACACTCCTTAAGTGAGATGGACCGAATACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTATGATGATTTTa                                                                                                                   <  2:564592/146‑1 (MQ=255)
                                     agagATTACCAGCAGAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTATGATGATTTTaa                                                                                                                  >  2:4241722/1‑150 (MQ=255)
                                     agagATTACCAGCAGAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTATGATGATTTTaa                                                                                                                  >  1:1802355/1‑150 (MQ=255)
                                        gaTTACCAGCAGAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGa                                                                                                               >  2:1347852/1‑150 (MQ=255)
                                          ttACCAGCAGAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACa                                                                                                             <  1:4332362/150‑1 (MQ=255)
                                              cagcagAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGtt                                                                                                         <  1:3933791/150‑1 (MQ=255)
                                              cagcagAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGtt                                                                                                         <  1:3933814/150‑1 (MQ=255)
                                              cagcagAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGtt                                                                                                         <  1:1378930/150‑1 (MQ=255)
                                                gcagAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGa                                                                                                       <  1:118865/150‑1 (MQ=255)
                                                 cagACAAGTCCTTGCTCGAAAGGGAGCACTCCAAGCTATGGAAACACTTACCACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGat                                                                                                      <  2:2667933/150‑1 (MQ=255)
                                                  agAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGatt                                                                                                     <  1:3659829/150‑1 (MQ=255)
                                                   gAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGattt                                                                                                    >  2:4491500/1‑150 (MQ=255)
                                                     aaaTTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGatttat                                                                                                  >  2:417323/1‑150 (MQ=255)
                                                       aTTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTGTTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGatttattt                                                                                                <  2:2488269/150‑1 (MQ=255)
                                                       aTTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGatttattt                                                                                                <  1:2522252/150‑1 (MQ=255)
                                                        ttCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGAtttattta                                                                                               >  2:3681560/1‑150 (MQ=255)
                                                             tGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGATTTATTTAGAAtt                                                                                          >  2:4972167/1‑150 (MQ=255)
                                                              gTTCGAAATTTAGCAGTTAAATCTATTGAAATACTTAAAACGGATATACACTCCTTAAGTGAGATGGACCGAAGACATGTTGTCTCTATTTGAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGATTTATTTAGAAttt                                                                                         <  2:3758846/150‑1 (MQ=255)
                                                              gTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGATTTATTTAGAAttt                                                                                         <  2:2242935/150‑1 (MQ=255)
                                                                 cGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGATTTATTTAGAATTTTAc                                                                                      >  2:607499/1‑150 (MQ=255)
                                                                 cGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGATTTATTTAGAATTTTAc                                                                                      <  1:4571550/150‑1 (MQ=255)
                                                                   aaaTTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGATTTATTTAGAATTTTAccc                                                                                    <  2:4017755/150‑1 (MQ=255)
                                                                   aaaTTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGATTTATTTAGAATTTTAccc                                                                                    >  1:775104/1‑150 (MQ=255)
                                                                   aaaTTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGATTTATTTAGAATTTTAccc                                                                                    >  2:1197320/1‑150 (MQ=255)
                                                                      tttAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGATTTATTTAGAATTTTACCCGaa                                                                                 <  2:2419909/150‑1 (MQ=255)
                                                                      tttAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGATTTATTTAGAATTTTACCCGaa                                                                                 >  2:1574206/1‑150 (MQ=255)
                                                                           cAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATATTGATTTATTTAGAATTTTACCCGAATCATa                                                                            >  1:3265225/1‑150 (MQ=255)
                                                                               gCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGATTTATTTAGAATTTTACCCGAATCATAAGAt                                                                        >  1:4698999/1‑150 (MQ=255)
                                                                                 aaGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGATTTATTTAGAATTTTACCCGAATCATAAGATTc                                                                      >  2:3390479/1‑150 (MQ=255)
                                                                                     tATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGATTTATTTAGAATTTTACCCGAATCATAAGATTCTGat                                                                  >  2:738670/1‑150 (MQ=255)
                                                                                      aTGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGATTTATTTAGAATTTTACCCGAATCATAAGATTCTGata                                                                 <  1:662680/150‑1 (MQ=255)
                                                                                       tGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGATTTATTTAGAATTTTACCCGAATCATAAGATTCTGatat                                                                <  1:276352/150‑1 (MQ=255)
                                                                                              aCTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGATTTATTTAGAATTTTACCCGAATCATAAGATTCTGATATCTTAAAt                                                         <  2:2793953/150‑1 (MQ=255)
                                                                                              aCTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGATTTATTTAGAATTTTACCCGAATCATAAGATTCTGATATCTTAAAt                                                         >  1:4886029/1‑150 (MQ=255)
                                                                                                  aCAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGATTTATTTAGAATTTTACCCGAATCATAAGATTCTGATATCTTAAATAAGa                                                     >  2:4873564/1‑150 (MQ=255)
                                                                                                       ggAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGATTTATTTAGAATTTTACCCGAATCATAAGATTCTGATATCTTAAATAAGAAAATc                                                >  2:1442087/1‑150 (MQ=255)
                                                                                                       ggAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGATTTATTTAGAATTTTACCCGAATCATAAGATTCTGATATCTTAAATAAGAAAATc                                                >  2:1442051/1‑150 (MQ=255)
                                                                                                           cTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGATTTATTTAGAATTTTACCCGAATCATAAGATTCTGATATCTTAAATAAGAAAATCCCGt                                            <  2:2928017/150‑1 (MQ=255)
                                                                                                             acacTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGATTTATTTAGAATTTTACCCGAATCATAAGATTCTGATATCTTAAATAAGAAAATCCCGTTc                                          <  1:3368500/150‑1 (MQ=255)
                                                                                                               acTCCTTAAGTGAGAAGGACCGAAGACATGTAGTCTCTATTTTATCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGATTTATTTAGAATTTTACCCGAATCATAAGATTCTGATATCTTAAATAATAAAATCCCGTTCGa                                        >  2:2124017/1‑150 (MQ=255)
                                                                                                                cTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGATTTATTTAGAATTTTACCCGAATCATAAGATTCTGATATCTTAAATAAGAAAATCCCGTTCGAc                                       <  2:2740077/150‑1 (MQ=255)
                                                                                                                  ccTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGATTTATTTAGAATTTTACCCGAATCATAAGATTCTGATATCTTAAATAAGAAAATCCCGTTCGACGt                                     >  2:811002/1‑150 (MQ=255)
                                                                                                                   cTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGATTTATTTAGAATTTTACCCGAATCATAAGATTCTGATATCTTAAATAAGAAAATCCCGTTCGACGtt                                    <  2:3582075/150‑1 (MQ=255)
                                                                                                                   cTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGATTTATTTAGAATTTTACCCGAATCATAAGATTCTGATATCTTAAATAAGAAAATCCCGTTCGACGtt                                    <  2:2349970/150‑1 (MQ=255)
                                                                                                                   cTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGATTTATTTAGAATTTTACCCGAATCATAAGATTCTGATATCTTAAATAAGAAAATCCCGTTCGACGtt                                    >  2:2893152/1‑150 (MQ=255)
                                                                                                                    ttAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGATTTATTTAGAATTTTACCCGAATCATAAGATTCTGATATCTTAAATAAGAAAATCCCGTTCGACGTTc                                   >  2:112686/1‑150 (MQ=255)
                                                                                                                      aaGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGATTTATTTAGAATTTTACCCGAATCATAAGATTCTGATATCTTAAATAAGAAAATCCCGTTCGACGTTCTg                                 >  1:3157091/1‑150 (MQ=255)
                                                                                                                      aaGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGATTTATTTAGAATTTTACCCGAATCATAAGATTCTGATATCTTAAATAAGAAAATCCCGTTCGACGTTCTg                                 <  2:458971/150‑1 (MQ=255)
                                                                                                                        gTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGATTTATTTAGAATTTTACCCGAATCATAAGATTCTGATATCTTAAATAAGAAAATCCCGTTCGACGTTCTGaa                               <  1:4242678/150‑1 (MQ=255)
                                                                                                                         tGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGATTTATTTAGAATTTTACCCGAATCATAAGATTCTGATATCTTAAATAAGAAAATCCCGTTCGACGTTCTGaaa                              <  1:1269132/150‑1 (MQ=255)
                                                                                                                            gaTGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGATTTATTTAGAATTTTACCCGAATCATAAGATTCTGATATCTTAAATAAGAAAATCCCGTTCGACGTTCTGAAAAAg                           <  1:2418631/150‑1 (MQ=255)
                                                                                                                            gaTGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGATTTATTTAGAATTTTACCCGAATCATAAGATTCTGATATCTTAAATAAGAAAATCCCGATCGACGTTCTGAAAAAg                           >  1:3687592/1‑150 (MQ=255)
                                                                                                                              tGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGATTTATTTAGAATTTTACCCGAATCATAAGATTCTGATATCTTAAATAAGAAAATCCCGTTCGACGTTCTGAAAAAGaa                         >  1:3419855/1‑150 (MQ=255)
                                                                                                                               ggACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGATTTATTTAGAATTTTACCCGAATCATAAGATTCTGATATCTTAAATAAGAAAATCCCGTTCGACGTTCTGAAAAAGaaa                        >  2:1662920/1‑150 (MQ=255)
                                                                                                                                 aCCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGATTTATTTAGAATTTTACCCGAATCATAAGATTCTGATATCTTAAATAAGAAAATCCCGTTCGACGTTCTGAAAAAGAAATc                      >  1:1754199/1‑150 (MQ=255)
                                                                                                                                   cGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGATTTATTTAGAATTTTACCCGAATCATAAGATTCTGATATCTTAAATAAGAAAATCCCGTTCGACGTTCt                                  >  2:1857730/1‑136 (MQ=255)
                                                                                                                                   cGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGATTTATTTAGAATTTTACCCGAATCATAAGATTCTGATATCTTAAATAAGAAAATCCCGTTCGACGTTCt                                  <  1:1857730/136‑1 (MQ=255)
                                                                                                                                   cGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGATTTATTTAGAATTTTACCCGAATCATAAGATTCTGATATCTTAAATAAGAAAATCCCGTTCGACGTTCt                                  <  1:4432118/136‑1 (MQ=255)
                                                                                                                                    gAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGATTTATTTAGAATTTTACCCGAATCATAAGATTCTGATATCTTAAATAAGAAAATCCCGTTCGACGTTCTGAAAAAGAAATCATa                   >  2:2552973/1‑150 (MQ=255)
                                                                                                                                     aaGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGATTTATTTAGAATTTTACCCGAATCATAAGATTCTGATATCTTAAATAAGAAAATCCCGTTCGACGTTCTGAAAAAGAAATCATaa                  >  1:1988061/1‑150 (MQ=255)
                                                                                                                                      aGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGATTTATTTAGAATTTTACCCGAATCATAAGATTCTGATATCTTAAATAAGAAAATCCCGTTCGACGTTCTGAAAAAGAAATCATaaa                 >  1:3976489/1‑150 (MQ=255)
                                                                                                                                      aGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGATTTATTTAGAATTTTACCCGAATCATAAGATTCTGATATCTTAAATAAGAAAATCCCGTTCGACGTTCTGAAAAAGAAATCATaaa                 <  1:5104392/150‑1 (MQ=255)
                                                                                                                                      aGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAAGGACATGTTGATTTATTTAGAATTTTACCCGAATCATAAGATTCTGATATCTTAAATAAGAAAATCCCGTTCGACGTTCTGAAAAAGAAATCATaaa                 <  2:4077480/150‑1 (MQ=255)
                                                                                                                                        aCATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGATTTATTTAGAATTTTACCCGAATCATAAGATTCTGATATCTTAAATAAGAAAATCCCGTTCGACGTTCTGAAAAAGAAATCATAAAat               <  1:3592848/150‑1 (MQ=255)
                                                                                                                                         cATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGATTTATTTAGAATTTTACCCGAATCATAAGATTCTGATATCTTAAATAAGAAAATCCCGTTCGACGTTCTGAAAAAGAAATCATAAAata              <  1:4135239/150‑1 (MQ=255)
                                                                                                                                            gTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGATTTATTTAGAATTTTACCCGAATCATAAGATTCTGATATCTTAAATAAGAAAATCCCGTTCGACGTTCTGAAAAAGAAATCATAAAATATTa           >  1:3514893/1‑150 (MQ=255)
                                                                                                                                             tAGTCTCTATTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGATTTATTTAGAATTTTACCCGAATCATAAGATTCTGATATCTTAAATAAGAAAATCCCGTTCGACGTTCTGAAAAAGAAATCATAAAATATTAg          <  2:1908446/150‑1 (MQ=255)
                                                                                                                                               gTCTCTATTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGATTTATTTAGAATTTTACCCGAATCATAAGATTCTGATATCTTAAATAAGAAAATCCCGTTCGACGTTCTGAAAAAGAAATCATAAAATATTAGCa        <  1:1208674/150‑1 (MQ=255)
                                                                                                                                                 ctctGTTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGATTTATTTAGAATTTTACCCGAATCATAAGATTCTGATATCTTAAATAAGAAAATCCCGTTCGACGTTCTGAAAAAGAAATCATAAAATATTAGCAGa      >  2:3620992/1‑150 (MQ=255)
                                                                                                                                                   ctATTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGATTTATTTAGAATTTTACCCGAATCATAAGATTCTGATATCTTAAATAAGAAAATCCCGTTCGACGTTCTGAAAAAGAAATCATAAAATATTAGCAGATa    >  1:481146/1‑150 (MQ=255)
                                                                                                                                                     aTTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGATTTATTTAGAATTTTACCCGAATCATAAGATTCTGATATCTTAAATAAGAAAATCCCGTTCGACGTTCTGAAAAAGAAATCATAAAATATTAGCAGATACt  <  2:834977/150‑1 (MQ=255)
                                                                                                                                                     aTTTTAGCCGAAATAGCTAATGGATATGATGATTTTAATGACTTGTTGATTTATTTAGAATTTTACCCGAATCATAAGATTCTGATATCTTAAATAAGAAAATCCCGTTCGACGTTCTGAAAAAGAAATCATAAAATATTAGCAGATACt  <  1:3021080/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
AGGAAACAAAAATCATTGTAAACCATGCCATCTTTTTAGAGATTACCAGCAGAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGATTTATTTAGAATTTTACCCGAATCATAAGATTCTGATATCTTAAATAAGAAAATCCCGTTCGACGTTCTGAAAAAGAAATCATAAAATATTAGCAGATACT  >  CP009273/3577678‑3577976

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: