Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,577,827 A→G 26.1% 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,8270AG26.1% 755.8 / inf 429intergenic (+158/‑292)yrhA/yrhBpseudogene, interrupted by IS1E/stable heat shock chaperone
Reads supporting (aligned to +/- strand):  ref base A (162/154);  new base G (45/67);  total (208/221)
Fisher's exact test for biased strand distribution p-value = 4.80e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 432 total aligned reads displayed.

GAAACAAAAATCATTGTAAACCATGCCATCTTTTTAGAGATTACCAGCAGAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGATTTATTTAGAATTTTACCCGAATCATAAGATTCTGATATCTTAAATAAGAAAATCCCGTTCGACGTTCTGAAAAAGAAATCATAAAATATTAGCAGATA  >  CP009273/3577680‑3577974
                                                                                                                                                   |                                                                                                                                                   
gAAACAAAAATCATTGTAAACCATGCCATCTTTTTAGAGATTACCAGCAGAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTa                                                                                                                                                     >  2:1588916/1‑148 (MQ=255)
    cAAAAATCATTGTAAACCATGCCATCTTTTTAGAGATTACCAGCAGAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACCACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTAtttt                                                                                                                                                 <  1:404186/148‑1 (MQ=255)
      aaaaTCATTGTAAACCATGCCATCTTTTTAGAGATTACCAGCAGAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAg                                                                                                                                               <  1:604062/148‑1 (MQ=255)
       aaaTCATTGTAAACCATGCCATCTTTTTAGAGATTACCAGCAGAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTGTTTTAGc                                                                                                                                              >  2:1790266/1‑148 (MQ=255)
         aTCATTGTAAACCATGCCATCTTTTTAGAGATTACCAGCAGAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCg                                                                                                                                            >  2:496911/1‑148 (MQ=255)
          tCATTGTAAACCATGCCATCTTTTTAGAGATTACCAGCAGAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTGTTTTAGCCGa                                                                                                                                           >  2:1682804/1‑148 (MQ=255)
          tCATTGTAAACCATGCCATCTTTTTAGAGATTACCAGCAGAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGa                                                                                                                                           <  2:954032/148‑1 (MQ=255)
          tCATTGAAAACCATGCCATCTTTTTAGAGATTACCAGCAGAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTGTTTTAGCCGa                                                                                                                                           <  2:1284409/148‑1 (MQ=255)
           cATTGTAAACCATGCCATCTTTTTAGAGATTACCAGCAGAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTGTTTTAGCCGaa                                                                                                                                          >  2:421230/1‑148 (MQ=255)
           cATTGTAAACCATGCCATCTTTTTAGAGATTACCAGCAGAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTGTTTTAGCCGaa                                                                                                                                          <  1:2478266/148‑1 (MQ=255)
             ttGTAAACCATGCCATCTTTTTAGAGATTACCAGCAGAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTGTTTTAGCCGAAAt                                                                                                                                        <  2:1543722/148‑1 (MQ=255)
             ttGTAAACCATGCCATCTTTTTAGAGATTACCAGCAGAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAAt                                                                                                                                        <  1:2029454/148‑1 (MQ=255)
             ttGTAAACCATGCCATCTTTTTAGAGATTACCAGCAGAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAAt                                                                                                                                        >  2:233338/1‑148 (MQ=255)
             ttGTAAACCATGCCATCTTTTTAGAGATTACCAGCAGAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGACGAAAt                                                                                                                                        >  2:1194561/1‑148 (MQ=255)
               gTAAACCATGCCATCTTTTTAGAGATTACCAGCAGAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAg                                                                                                                                      <  2:2425713/148‑1 (MQ=255)
                 aaaCCATGCCATCTTTTTAGAGATTACCAGCAGAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCt                                                                                                                                    >  1:1922786/1‑148 (MQ=255)
                 aaaCCATGCCATCTTTTTAGAGATTACCAGCAGAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCt                                                                                                                                    >  1:1018670/1‑148 (MQ=255)
                  aaCCATGCCATCTTTTTAGAGATTACCAGCAGAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTGTTTTAGCCGAAATAGCTa                                                                                                                                   <  2:2689167/148‑1 (MQ=255)
                  aaCCATGCCATCTTTTTAGAGATTACCAGCAGAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTa                                                                                                                                   >  1:1390130/1‑148 (MQ=255)
                  aaCCATGCCATCTTTTTAGAGATTACCAGCAGAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTa                                                                                                                                   <  2:657948/148‑1 (MQ=255)
                     cATGCCATCTTTTTAGAGATTACCAGCAGAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATg                                                                                                                                <  2:2728365/148‑1 (MQ=255)
                      aTGCCATCTTTTTAGAGATTACCAGCAGAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTGTTTTAGCCGAAATAGCTAATgg                                                                                                                               >  1:3393470/1‑148 (MQ=255)
                       tGCCATCTTTTTAGAGATTACCAGCAGAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGc                                                                                                                              >  1:1350700/1‑148 (MQ=255)
                           aTCTTTTTAGAGATTACCAGCAGAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTGTTTTAGCCGAAATAGCTAATGGCTatg                                                                                                                          >  2:3133285/1‑148 (MQ=255)
                            tCTTTTTAGAGATTACCAGCAGAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTatga                                                                                                                         <  1:3709117/148‑1 (MQ=255)
                             cTTTTTAGAGATTACCAGCAGAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTGTTTTAGCCGAAATAGCTAATGGCTatgat                                                                                                                        >  2:1971570/1‑148 (MQ=255)
                             cTTTTTAGAGATTACCAGCAGAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTatgat                                                                                                                        <  1:1038561/148‑1 (MQ=255)
                             cTTTTTAGAGATTACCAGCAGAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTatgat                                                                                                                        >  1:2213247/1‑148 (MQ=255)
                                tttAGAGATTACCAGCAGAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTATgatgat                                                                                                                     <  2:544136/148‑1 (MQ=255)
                                    gagaTTACCAGCAGAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTGTTTTAGCCGAAATAGCTAATGGCTATGATGATTTTa                                                                                                                 <  2:3038550/148‑1 (MQ=255)
                                    gagaTTACCAGCAGAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTGTTTTAGCCGAAATAGCTAATGGCTATGATGATTTTa                                                                                                                 <  1:1252136/148‑1 (MQ=255)
                                    gagaTTACCAGCAGAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTATGATGATTTTa                                                                                                                 >  2:2793789/1‑148 (MQ=255)
                                     agaTTACCAGCAGAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTGTTTTAGCCGAAATAGCTAATGGCTATGATGATTTTaa                                                                                                                <  2:3393470/148‑1 (MQ=255)
                                     agaTTACCAGCAGAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTATGATGATTTTaa                                                                                                                >  1:2393273/1‑148 (MQ=255)
                                       aTTACCAGCAGAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATg                                                                                                              >  1:2097306/1‑148 (MQ=255)
                                           ccAGCAGAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACAt                                                                                                          <  1:3322418/148‑1 (MQ=255)
                                              gcagAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGtt                                                                                                       >  2:1391747/1‑148 (MQ=255)
                                                agAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGa                                                                                                     >  1:731265/1‑148 (MQ=255)
                                                agAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGa                                                                                                     <  2:1158354/148‑1 (MQ=255)
                                                  aaaaTTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTGTTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGatt                                                                                                   >  1:2555255/1‑148 (MQ=255)
                                                  aaaaTTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTGTTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGatt                                                                                                   >  1:637696/1‑148 (MQ=255)
                                                  aaaaTTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTGTTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGatt                                                                                                   <  1:2412509/148‑1 (MQ=255)
                                                  aaaaTTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGatt                                                                                                   >  2:1741966/1‑148 (MQ=255)
                                                   aaaTTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGattt                                                                                                  <  2:756834/148‑1 (MQ=255)
                                                   aaaTTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGattt                                                                                                  >  2:3156009/1‑148 (MQ=255)
                                                    aaTTCCTTTTTCTAAATTTAGCAGTGAAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGattta                                                                                                 <  1:3456168/148‑1 (MQ=255)
                                                    aaTTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGattta                                                                                                 <  2:1933649/148‑1 (MQ=255)
                                                     aTTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTGTTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGatttat                                                                                                <  1:113222/148‑1 (MQ=255)
                                                     aTTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGatttat                                                                                                <  2:3290564/148‑1 (MQ=255)
                                                      ttCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTCTGATGATTTTAATGACATGTTGatttatt                                                                                               >  1:3677024/1‑148 (MQ=255)
                                                      ttCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGatttatt                                                                                               <  1:1865599/148‑1 (MQ=255)
                                                       tCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGatttattt                                                                                              <  2:1054225/148‑1 (MQ=255)
                                                          ttGTTCGAAATTTAGCAGTGCAAGCTATGTAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGATTTATTTAGa                                                                                           <  1:717743/148‑1 (MQ=255)
                                                             ttCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGATTTATTTAGAAtt                                                                                        <  1:2967143/148‑1 (MQ=255)
                                                             ttCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGATTTATTTAGAAtt                                                                                        <  1:3775119/148‑1 (MQ=255)
                                                                 aaaTTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTGTTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGATTTATTTAGAATTTTAc                                                                                    >  2:2961910/1‑148 (MQ=255)
                                                                    tttAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGATTTATTTAGAATTTTACCCg                                                                                 >  2:2573875/1‑148 (MQ=255)
                                                                        gCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGATTTATTTAGAATTTTACCCGAATc                                                                             >  1:281127/1‑148 (MQ=255)
                                                                         cAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGATTTATTTAGAATTTTACCCGAATCa                                                                            >  2:374617/1‑148 (MQ=255)
                                                                                   tATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGATTTATTTAGAATTTTACCCGAATCATAAGATTCTg                                                                  <  1:946448/148‑1 (MQ=255)
                                                                                   tATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGATTTATTTAGAATTTTACCCGAATCATAAGATTCTg                                                                  <  1:3753898/148‑1 (MQ=255)
                                                                                    aTGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGATTTATTTAGAATTTTACCCGAATCATAAGATTCTGa                                                                 <  1:3039675/148‑1 (MQ=255)
                                                                                    aTGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGATTTATTTAGAATTTTACCCGAATCATAAGATTCTGa                                                                 >  2:1830275/1‑148 (MQ=255)
                                                                                            aCTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTGTTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGATTTATTTAGAATTTTACCCGAATCATAAGATTCTGATATCTTaa                                                         >  2:3719214/1‑148 (MQ=255)
                                                                                            aCTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGATTTATTTAGAATTTTACCCGAATCATAAGATTCTGATATCTTaa                                                         >  2:3700469/1‑148 (MQ=255)
                                                                                               tACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGATTTATTTAGAATTTTACCCGAATCATAAGATTCTGATATCTTAAATa                                                      >  2:3505937/1‑148 (MQ=255)
                                                                                                   aCGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGATTTATTTAGAATTTTACCCGAATCATAAGATTCTGATATCTTAAATAAGaa                                                  >  2:2941202/1‑148 (MQ=255)
                                                                                                   aCGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGATTTATTTAGAATTTTACCCGAATCATAAGATTCTGATATCTTAAATAAGaa                                                  <  2:731103/148‑1 (MQ=255)
                                                                                                     ggAACTACACTCCTTAAGTGAGATGGAACGAAGACATGTAGTCTCTGTTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGATTTATTTAGAATTTTACCCGAATCATAAGATTCTGATATCTTAAATAAGaaaa                                                <  2:2199240/148‑1 (MQ=255)
                                                                                                      gAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTGTTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGATTTATTTAGAATTTTACCAGAATCATAAGATTCTGATATCTTAAATAAGAAAAt                                               >  1:2332831/1‑148 (MQ=255)
                                                                                                       aaCTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTGTTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGATTTATTTAGAATTTTACCCGAATCATAAGATTCTGATATCTTAAATAAGAAAATc                                              <  2:2160940/148‑1 (MQ=255)
                                                                                                          tACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGATTTATTTAGAATTTTACCCGAATCATAAGATTCTGATATCTTAAATAAGAAAATCCCg                                           >  1:855067/1‑148 (MQ=255)
                                                                                                           acacTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGATTTATTTAGAATTTTACCCGAATCATAAGATTCTGATATCTTAAATAAGAAAATCCCGt                                          >  2:261991/1‑148 (MQ=255)
                                                                                                           acacTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGATTTATTTAGAATTTTACCCGAATCATAAGATTCTGATATCTTAAATAAGAAAATCCAGt                                          >  2:3210759/1‑148 (MQ=255)
                                                                                                            cacTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGATTTATTTAGAATTTTACCCGAATCATAAGATTCTGATATCTTAAATAAGAAAATCCCGtt                                         <  1:17473/148‑1 (MQ=255)
                                                                                                               tCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGATTTATTTAGAATTTTACCCGAATCATAAGATTCTGATATCTTAAATAAGAAAATCCCGTTCGa                                      >  1:2632576/1‑148 (MQ=255)
                                                                                                                ccTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCTAAATAGCTAATGGCTATGATTATTTTAATGAAATGATTATTTATTTACAATATTAACCGAATAATAACAATATTAAATAATAAATAATAAAATACATTTCGaa                                     >  1:1658456/1‑147 (MQ=255)
                                                                                                                 cTTAAGTGAGATGGACCGAAGACATGTAGTCTCTGTTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATTACATGTTGATTTATTTAGAATTTTACCCGAATCATAAGATTCTGATATATTAAATAAGAAAATCCCGTTCGACg                                    >  2:3330837/1‑148 (MQ=255)
                                                                                                                 cTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGATTTATTTAGAATTTTACCCGAATCATAAGATTCTGATATCTTAAATAAGAAAATCCCGTTCGACg                                    <  2:2190014/148‑1 (MQ=255)
                                                                                                                    aaGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGATTTATTTAGAATTTTACCCGAATCATAAGATTCTGATATCTTAAATAAGAAAATCCCGTTCGACGTTc                                 <  1:527118/148‑1 (MQ=255)
                                                                                                                        gagaTGGACCGAAGACATGTAGTCTCTGTTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGATTTATTTAGAATTTTACCCGAATCATAAGATTCTGATATCTTAAATAAGAAAATCCCGTTCGACGTTCTGaa                             >  1:3107746/1‑148 (MQ=255)
                                                                                                                         agaTGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGATTTATTTAGAATTTTACCCGAATCATAAGATTCTGATATCTTAAATAAGAAAATCCCGTTCGACGTTCTGaaa                            >  2:1362487/1‑148 (MQ=255)
                                                                                                                          gaTGGACCGAAGACATGTAGTCTCTGTTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGATTTATTTAGAATTTTACCCGAATCATAAGATTCTGATATCTTAAATAAGAAAATCCCGTTCGACGTTCTGaaaa                           >  2:1452641/1‑148 (MQ=255)
                                                                                                                          gaTGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGATTTATTTAGAATTTTACCCGAATCATAAGATTCTGATATCTTCAATAAGAAAATCCCGTTCGACGTTCTGaaaa                           >  2:821468/1‑148 (MQ=255)
                                                                                                                           aTGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGATTTATTTAGAATTTTACCCGAGTCATAAGATTCTGATATCTTAAATAAGAAAATCCCGTTCGACGTTCTGaaaaa                          >  1:3582348/1‑148 (MQ=255)
                                                                                                                            tGGACCGAAGACATGTAGTCTCTGTTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGATTTATTTAGAATTTTACCCGAATCATAAGATTCTGATATCTTAAATAAGAAAATCCCGTTCGACGTTCTGAAAAAg                         <  1:2858264/148‑1 (MQ=255)
                                                                                                                               aCCGAAGACATGTAGTCTCTGTTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGATTTATTTAGAATTTTACCCGAATCATAAGATTCTGATATCTTAAATAAGAAAATCCCGTTCGACGTTCTGAAAAAGaaa                      <  2:3410731/148‑1 (MQ=255)
                                                                                                                                ccGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGATTTATTTAGAATTTTACCCGAATCATAAGATTCTGATATCTTAAATAAGAAAATCCCGTTCGACGTTCTGAAAAAGAAAt                     >  1:2566965/1‑148 (MQ=255)
                                                                                                                                   aaGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGATTTATTTAGAATTTTACCCGAATCATAAGATTCTGATATCTTAAATAAGAAAATCCCGTTCGACGTTCTGAAAAAGAAATCAt                  >  2:2486369/1‑148 (MQ=255)
                                                                                                                                    aGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGATTTATTTAGAATTTTACCCGAATCATAAGATTCTGATATCTTAAATAAGAAAATCCCGTTCGACGTTCTGAAAAAGAAATCATa                 >  1:3209793/1‑148 (MQ=255)
                                                                                                                                       cATGTAGTCTCTGTTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGATTTATTTAGAATTTTACCCGaa                                                                               <  1:788552/83‑1 (MQ=255)
                                                                                                                                       cATGTAGTCTCTGTTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGATTTATTTAGAATTTTACCCGaa                                                                               >  2:1322627/1‑83 (MQ=255)
                                                                                                                                       cATGTAGTCTCTGTTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGATTTATTTAGAATTTTACCCGaa                                                                               <  1:1322627/83‑1 (MQ=255)
                                                                                                                                        aTGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGATTTATTTAGAATTTTACCCGAATCATAAGATTCTGATATCTTAAATAAGAAAATCCCGTTCGACGTTCTGAAAAAGAAATCATAAAat             >  1:1595244/1‑148 (MQ=255)
                                                                                                                                         tGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGATTTATTTAGAATTTTACCCGAATCATAAGATTCTGATATCTTAAATAAGAAAATCCCGTTCGACGTTCTGAAAAAGAAATCATAAAata            <  2:3556445/148‑1 (MQ=255)
                                                                                                                                            aGTCTCTGTTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGATTTATTTAGAATTTTACCCGAATCATAAGATTCTGATATCTTAAATAAGAAAATCCCGTTCGACGTTCTGAAAAAGAAATCATAAAATATTa         <  1:1790266/148‑1 (MQ=255)
                                                                                                                                             gTCTCTATTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGATTTATTTAGAATTTTACCCGAATCATAAGATTCTGATATCTTAAATAAGAAAATCCCGTTCGACGTTCTGACAAAGAAATCATAAAATATTAg        <  1:679641/148‑1 (MQ=255)
                                                                                                                                             gTCTCTATTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGATTTATTTAGAATTTTACCCGAATCATAAGATTCTGATATCTTAAATAAGAAAATCCCGTTCGACGTTCTGAAAAAGAAATCATAAAATATTAg        <  2:788789/148‑1 (MQ=255)
                                                                                                                                              tctcTATTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGATTTATTTAGAATTTTACCCGAATCATAAGATTCTGATATCTTAAATAAGAAAATCCCGTTCGACGTTCTGAAAAAGAAATCATAATATATTAGc       >  2:2066348/1‑148 (MQ=255)
                                                                                                                                                tctATTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGATTTATTTAGAATTTTACCCGAATCATAAGATTCTGATATCTTAAATAAGAAAATCCCGTTCGACGTTCTGAAAAAGAAATCATAAAATATTAGCAg     <  2:1846550/148‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
GAAACAAAAATCATTGTAAACCATGCCATCTTTTTAGAGATTACCAGCAGAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGATTTATTTAGAATTTTACCCGAATCATAAGATTCTGATATCTTAAATAAGAAAATCCCGTTCGACGTTCTGAAAAAGAAATCATAAAATATTAGCAGATA  >  CP009273/3577680‑3577974

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: