Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 320,920 G→A 63.3% S52S (AGC→AGT ykgH ← putative inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273320,9200GA63.3% 211.3 / 254.4 219S52S (AGC→AGTykgHputative inner membrane protein
Reads supporting (aligned to +/- strand):  ref base G (29/50);  new base A (57/81);  total (87/132)
Fisher's exact test for biased strand distribution p-value = 5.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

Only 100 of 220 total aligned reads displayed.

TATTTCTCTGTTCTTGACATCTGGGCTGGTGTGAAAAAATGTAGTAATTCGCTGAATAACATCTTCAGTCATTGTAGTCTGACAGGCAAAATTATTATTGCTATCACAGCTAACAAGTATACCTGGCGTTGTCTGGTTTATTTTTAGGCTGTAAATCACTTTTATATTGTCTTTATTGATATACAGAAACAACAGCACCATAGCCATACATATAGCCATAATAAAAAGGATAACACGAATCTTTTTCTTCAGATAAAATAAAATCTCAATCAGATCGATATCCTGTTTGATTTGTT  >  CP009273/320773‑321068
                                                                                                                                                   |                                                                                                                                                    
tATTTCTCTGTTCTTGACATCTGGGCTGGTGTGAAAAAATGTAGTAATTCGCTGAATAACATCTTCAGTCATTGTAGTCTGACAGGCAAAATTATTATTGCTATCACAGCTAACAAGTATACCTGGCGTTGTCTGGTTTATTTTTAGGc                                                                                                                                                     <  1:2098486/149‑1 (MQ=255)
  tttCTCTGTTCTTGACATCTGGGCTGGTGTGAAAAAATGTAGTAATTCGCTGAATAACATCTTCAGTCATTGTAGTCTGACAGGCAAAATTATTATTGCTATCACAGCTAACAAGTATACCTGGCGTTGTCTGGTTTATTTTTAGACTg                                                                                                                                                   <  1:2135539/149‑1 (MQ=255)
   ttCTCTGTTCTTGACATCTGGGCTGGTGTGAAAAAATGTAGTAATTCGCTGAATAACATCTTCAGTCATTGTAGTCTGACAGGCAAAATTATTATTGCTATCACAGCTAACAAGTATACCTGGCGTTGTCTGGTTTATTTTTAGGCTGt                                                                                                                                                  <  1:2370451/149‑1 (MQ=255)
    tctcTGTTCTTGACATCTGGGCTGGTGTGAAAAAATGTAGTAATTCGCTGAATAACATCTTCAGTCATTGTAGTCTGACAGGCAAAATTATTATTGCTATCACAGCTAACAAGTATACCTGGCGTTGTCTGGTTTATTTTTAGACTGTa                                                                                                                                                 <  2:3022327/149‑1 (MQ=255)
        tGTTCTTGACATCTGGGCTGGTGTGAAAAAATGTAGTAATTCGCTGAATAACATCTTCAGTCATTGTAGTCTGACAGGCAAAATTATTATTGCTATCACAGCTAACAAGTATACCTGGCGTTGTCTGGTTTATTTTTAGACTGTAAATc                                                                                                                                             >  2:1760186/1‑149 (MQ=255)
         gTTCTTGACATCTGGGCTGGTGTGAAAAAATGTAGTAATTCGCTGAATAACATCTTCAGTCATTGTAGTCTGACAGGCAAAATTATTATTGCTATCACAGCTAACAAGTATACCTGGCGTTGTCTGGTTTATTTTTAGGCTGTAAATCa                                                                                                                                            <  2:4575806/149‑1 (MQ=255)
            cTTGACATCTGGGCTGGTGTGAAAAAATGTAGTAATTCGCTGAATAACATCTTCAGTCATTGTAGTCTGACAGGCAAAATTATTATTGCTATCACAGCTAACAAGTATACCTGGCGTTGTCTGGTTTATTTTTAGGCTGTAAATCACtt                                                                                                                                         <  1:423404/149‑1 (MQ=255)
              tGACATCTGGGCTGGTGTGAAAAAATGTAGTAATTCGCTGAATAACATCTTCAGTCATTGTAGTCTGACAGGCAAAATTATTATTGCTATCACAGCTAACAAGTATACCTGGCGTTGTCTGGTTTATTTTTAGACTGTAAATCACtttt                                                                                                                                       >  2:4122055/1‑149 (MQ=255)
                aCATCTGGGCTGGTGTGAAAAAATGTAGTAATTCGCTGAATAACATCTTCAGTCATTGTAGTCTGACAGGCAAAATTATTATTGCTATCACAGCTAACAAGTATACCTGGCGTTGTCTGGTTTATTTTTAGACTGTAAATCACTTTtat                                                                                                                                     <  1:540984/149‑1 (MQ=255)
                  aTCTGGGCTGGTGTGAAAAAATGTAGTAATTCGCTGAATAACATCTTCAGTCATTGTAGTCTGACAGGCAAAATTATTATTGCTATCACAGCTAACAAGTATACCTGGCGTTGTCTGGTTTATTTTTAGACTGTAAATCACTTTTatat                                                                                                                                   <  1:1671885/149‑1 (MQ=255)
                    cTGGGCTGGTGTGAAAAAATGTAGTAATTCGCTGAATAACATCTTCAGTCATTGTAGTCTGACAGGCAAAATTATTATTGCTATCACAGCTAACAAGTATACCTGGCGTTGTCTGGTTTATTTTTAGACTGTAAATCACTTTTATATTg                                                                                                                                 <  1:2655680/149‑1 (MQ=255)
                     tGGGCTGGTGTGAAAAAATGTAGTAATTCGCTGAATAACATCTTCAGTCATTGTAGTCTGACAGGCAAAATTATTATTGCTATCACAGCTAACAAGTATACCTGGCGTTGTCTGGTTTATTTTTAGACTGTAAATCACTTTTATAtttt                                                                                                                                >  1:437010/1‑147 (MQ=255)
                      gggCTGGTGTGAAAAAATGTAGTAATTCGCTGAATAACATCTTCAGTCATTGTAGTCTGACAGGCAAAATTATTATTGCTATCACAGCTAACAAGTATACCTGGCGTTGTCTGGTTTATTTTTAGACTGTAAATCACTTTTATATTGTc                                                                                                                               <  2:2764265/149‑1 (MQ=255)
                       ggCTGGTGTGAAAAAATGTAGTAATTCGCTGAATAACATCTTCAGTCATTGTAGTCTGACAGGCAAAATTATTATTGCTATCACAGCTAACAAGTATACCTGGCGTTGTCTGGTTTATTTTTAGACTGTAAATCACTTTTATATTGTCt                                                                                                                              <  1:4519229/149‑1 (MQ=255)
                        gCTGGTGTGAAAAAATGTAGTAATTCGCTGAATAACATCTTCAGTCATTGTAGTCTGACAGGCAAAATTATTATTGCTATCACAGCTAACAAGTATACCTGGCGTTGTCTGGTTTATTTTTAGACTGTAAATCACTTTTATATTGTCtt                                                                                                                             <  1:299881/149‑1 (MQ=255)
                         cTGGTGTGAAAAAATGTAGTAATTCGCTGAATAACATCTTCAGTCATTGTAGTCTGACAGGCAAAATTATTATTGCTATCACAGCTAACAAGTATACCTGGCGTTGTCTGGTTTATTTTTAGGCTGTAAATCACTTTTATATTGTCttt                                                                                                                            <  1:2866248/149‑1 (MQ=255)
                         cTGGTGTGAAAAAATGTAGTAATTCGCTGAATAACATCTTCAGTCATTGTAGTCTGACAGGCAAAATTATTATTGCTATCACAGCTAACAAGTATACCTGGCGTTGTCTGGTTTATTTTTAGGCTGTAAATCACTTTTATATTGTCttt                                                                                                                            <  1:3308764/149‑1 (MQ=255)
                           ggTGTGAAAAAATGTAGTAATTCGCTGAATAACATCTTCAGTCATTGTAGTCTGACAGGCAAAATTATTATTGCTATCACAGCTAACAAGTATACCTGGCGTTGTCTGGTTTATTTTTAGACTGTAAATCACTTTTATATTGTCTTTAt                                                                                                                          >  2:1819459/1‑149 (MQ=255)
                           ggTGTGAAAAAATGTAGTAATTCGCTGAATAACATCTTCAGTCATTGTAGTCTGACAGGCAAAATTATTATTGCTATCACAGCTAACAAGTATACCTGGCGTTGTCTGGTTTATTTTTAGACTGTAAATCACTTTTATATTGTCTTTAt                                                                                                                          <  1:3043306/149‑1 (MQ=255)
                              gtgAAAAAATGTAGTAATTCGCTGAATAACATCTTCAGTCATTGTAGTCTGACAGGCAAAATTATTATTGCTATCACAGCTAACAAGTATACCTGGCGTTGTCTGGTTTATTTTTAGACTGTAAATCACTTTTATATTGTCTTTATTGa                                                                                                                       <  1:1760236/149‑1 (MQ=255)
                               tgAAAAAATGTAGTAATTCGCTGAATAACATCTTCAGTCATTGTAGTCTGACAGGCAAAATTATTATTGCTATCACAGCTAACAAGTATACCTGGCGTTGTCTGGTTTATTTTTAGGCTGTAAATCACTTTTATATTGTCTTTATTGat                                                                                                                      >  2:83890/1‑149 (MQ=255)
                                gAAAAAATGTAGTAATTCGCTGAATAACATCTTCAGTCATTGTAGTCTGACAGGCAAAATTATTATTGCTATCACAGCTAACAAGTATACCTGGCGTTGTCTGGTTTATTTTTAGACTGTAAATCACTTTTATATTGTCTTTATTGata                                                                                                                     <  2:3579384/149‑1 (MQ=255)
                                 aaaaaaTGTAGTAATTCGCTGAATAACATCTTCAGTCATTGTAGTCTGACAGGCAAAATTATTATTGCTATCACAGCTAACAAGTATACCTGGCGTTGTCTGGTTTATTTTTAGACTGTAAATCACTTTTATATTGTCTTTATTGatat                                                                                                                    >  1:2388093/1‑149 (MQ=255)
                                 aaaaaaTGTAGTAATTCGCTGAATAACATCTTCAGTCATTGTAGTCTGACAGGCAAAATTATTATTGCTATCACAGCTAACAAGTATACCTGGCGTTGTCTGGTTTATTTTTAGACTGTAAATCACTTTTATATTGTCTTTATTGatat                                                                                                                    <  1:2729541/149‑1 (MQ=255)
                                 aaaaaaTGTAGTAATTCGCTGAATAACATCTTCAGTCATTGTAGTCTGACAGGCAAAATTATTATTGCTATCACAGCTAACAAGTATACCTGGCGTTGTCTGGTTTATTTTTAGACTGTAAATCACTTTTATATTGTCTTTATTGatat                                                                                                                    >  1:3986706/1‑149 (MQ=255)
                                  aaaaaTGTAGTAATTCGCTGAATAACATCTTCAGTCATTGTAGTCTGACAGGCAAAATTATTATTGCTATCACAGCTAACAAGTATACCTGGCGTTGTCTGGTTTATTTTTAGGCTGTAAATCACTTTTATATTGTCTTTATTGAtata                                                                                                                   <  1:4318240/149‑1 (MQ=255)
                                  aaaaaTGTAGTAATTCGCTGAATAACATCTTCAGTCATTGTAGTCTGACAGGCAAAATTATTATTGCTATCACAGCTAACAAGTATACCTGGCGTTGTCTGGTTTATTTTTAGGCTGTAAATCACTTTTATATTGTCTTTATTGAtata                                                                                                                   <  1:116570/149‑1 (MQ=255)
                                     aaTGTAGTAATTCGCTGAATAACATCTTCAGTCATTGTAGTCTGACAGGCAAAATTATTATTGCTATCACAGCTAACAAGTATACCTGGCGTTGTCTGGTTTATTTTTAGGCTGTAAATCACTTTTATATTGTCTTTATTGATATACAg                                                                                                                >  2:703024/1‑149 (MQ=255)
                                     aaTGTAGTAATTCGCTGAATAACATCTTCAGTCATTGTAGTCTGACAGGCAAAATTATTATTGCTATCACAGCTAACAAGTATACCTGGCGTTGTCTGGTTTATTTTTAGACAGTAAATCAATTTTATATTGTCTTTATTGATATACat                                                                                                                >  1:1992846/1‑148 (MQ=255)
                                       tGTAGTAATTCGCTGAATAACATCTTCAGTCATTGTAGTCTGACAGGCAAAATTATTATTGCTATCACAGCTAACAAGTATACCTGGCGTTGTCTGGTTTATTTTTAGACTGTAAATCACTTTTATATTGTCTTTATTGATATACAGaa                                                                                                              >  2:1132841/1‑149 (MQ=255)
                                        gtagtaATTCGCTGAATAACATCTTCAGTCATTGTAGTCTGACAGGCAAAATTATTATTGCTATCACAGCTAACAAGTATACCTGGCGTTGTCTGGTTTATTTTTAGACTGTAAATCACTTTTATATTGTCTTTATTGATATACAGaaa                                                                                                             <  1:1003129/149‑1 (MQ=255)
                                          agtaATTCGCTGAATAACATCTTCAGTCATTGTAGTCTGACAGGCAAAATTATTATTGCTATCACAGCTAACAAGTATACCTGGCGTTGTCTGGTTTATTTTTAGGCTGTAAATCACTTTTATATTGTCTTTATTGATATACAGAaaca                                                                                                           <  1:3121173/149‑1 (MQ=255)
                                            taATTCGCTGAATAACATCTTCAGTCATTGTAGTCTGACAGGCAAAATTATTATTGCTATCACAGCTAACAAGTATACCTGGCGTTGTCTGGTTTATTTTTAGACTGTAAATCACTTTTATATTGTCTTTATTGATATACAGAaacaac                                                                                                         <  2:406043/149‑1 (MQ=255)
                                             aaTTCGCTGAATAACATCTTCAGTCATTGTAGTCTGACAGGCAAAATTATTATTGCTATCACAGCTAACAAGTATACCTGGCGTTGTCTGGTTTATTTTTAGGCTGTAAATCACTTTTATATTGTCTTTATTGATATACAGAAacaaca                                                                                                        <  1:3306592/149‑1 (MQ=255)
                                             aaTTCGCTGAATAACATCTTCAGTCATTGTAGTCTGACAGGCAAAATTATTATTGCTATCACAGCTAACAAGTATACCTGGCGTTGTCTGGTTTATTTTTAGGCTGTAAATCACTTTTATATTGTCTTTATTGATATACAGAAacaaca                                                                                                        >  2:3229497/1‑149 (MQ=255)
                                                 cGCTGAATAACATCTTCAGTCATTGTAGTCTGACAGGCAAAATTATTATTGCTATCACAGCTAACAAGTATACCTGGCGTTGTCTGGTTTATTTTTAGGCTGTAAATCACTTTTATATTGTCTTTATTGATATACAGAAACAACAGCAc                                                                                                    <  2:3929635/149‑1 (MQ=255)
                                                   cTGAATAACATCTTCAGTCATTGTAGTCTGACAGGCAAAATTATTATTGCTATCACAGCTAACAAGTATACCTGGCGTTGTCTGGTTTATTTTTAGACTGTAAATCACTTTTATATTGTCTTTATTGATATACAGAAACAACAGCACCa                                                                                                  <  2:4678836/149‑1 (MQ=255)
                                                   cTGAATAACATCTTCAGTCATTGTAGTCTGACAGGCAAAATTATTATTGCTATCACAGCTAACAAGTATACCTGGCGTTGTCTGGTTTATTTTTAGACTGTAAATCACTTTTATATTGTCTTTATTGATATACAGAAACAACAGCACCa                                                                                                  <  2:2314371/149‑1 (MQ=255)
                                                    tGAATAACATCTTCAGTCATTGTAGTCTGACAGGCAAAATTATTATTGCTATCACAGCTAACAAGTATACCTGGCGTTGTCTGGTTTATTTTTAGACTGTAAATCAATTTTATATTGTCTTTATTGATATAAATTAAAAACAGCACCaa                                                                                                 >  2:4431335/1‑148 (MQ=255)
                                                          aCATCTTCAGTCATTGTAGTCTGACAGGCAAAATTATTATTGCTATCACAGCTAACAAGTATACCTGGCGTTGTCTGGTTTATTTTTAGGCTGTAAATCACTTTTATATTGTCTTTATTGATATACAGAAACAACAGCACCATAGCcat                                                                                           <  1:2075456/149‑1 (MQ=255)
                                                          aCATCTTCAGTCATTGTAGTCTGACAGGCAAAATTATTATTGCTATCACAGCTAACAAGTATACCTGGCGTTGTCTGGTTTATTTTTAGGCTGTAAATCACTTTTATATTGTCTTTATTGATATACAGAAACAACAGCACCATAGCcat                                                                                           >  1:2211212/1‑149 (MQ=255)
                                                          aCATCTTCAGTCATTGTAGTCTGACAGGCAAAATTATTATTGCTATCACAGCTAACAAGTATACCTGGCGTTGTCTGGTTTATTTTTAGACTGTAAATCACTTTTATATTGTCTTTATTGATATACAGAAACAACAGCACCATAGCcat                                                                                           >  2:821024/1‑149 (MQ=255)
                                                           cATCTTCAGTCATTGTAGTCTGACAGGCAAAATTATTATTGCTATCACAGCTAACAAGTATACCTGGCGTTGTCTGGTTTATTTTTAGGCTGTAAATCACTTTTATATTGTCTTTATTGATATACAGAAACAACAGCACCATAGCcata                                                                                          >  2:3438336/1‑149 (MQ=255)
                                                           cATCTTCAGTCATTGTAGTCTGACAGGCAAAATTATTATTGCTATCACAGCTAACAAGTATACCTGGCGTTGTCTGGTTTATTTTTAGACTGTAAATCACTTTTATATTGTCTTTATTGATATACAGAAACAACAGCACCATAGCcata                                                                                          <  1:2470062/149‑1 (MQ=255)
                                                            aTCTTCAGTCATTGTAGTCTGACAGGCAAAATTATTATTGCTATCACAGCTAACAAGTATACCTGGCGTTGTCTGGTTTATTTTTAGACTGTAAATCACTTTTATATTGTCTTTATTGATATACAGAAACAACAGCACCATAGCcatac                                                                                         >  2:3211996/1‑149 (MQ=255)
                                                            aTCTTCAGTCATTGTAGTCTGACAGGCAAAATTATTATTGCTATCACAGCTAACAAGTATACCTGGCGTTGTCTGGTTTATTTTTAGACTGTAAATCACTTTTATATTGTCTTTATTGATATACAGAAACAACAGCACCATAGCcatac                                                                                         >  1:2400318/1‑149 (MQ=255)
                                                             tCTTCAGTCATTGTAGTCTGACAGGCAAAATTATTATTGCTATCACAGCTAACAAGTATACCTGGCGTTGTCTGGTTTATTTTTAGGCTGTAAATCACTTTTATATTGTCTTTATTGATATACAGAAACAACAGCACCATAGCcataca                                                                                        >  2:1866992/1‑149 (MQ=255)
                                                             tCTTCAGTCATTGTAGTCTGACAGGCAAAATTATTATTGCTATCACAGCTAACAAGTATACCTGGCGTTGTCTGGTTTATTTTTAGACTGTAAATCACTTTTATATTGTCTTTATTGATATACAGAAACAACAGCACCATAGCcataca                                                                                        >  2:2115036/1‑149 (MQ=255)
                                                              cTTCAGTCATTGTAGTCTGACAGGCAAAATTATTATTGCTATCACAGCTAACAAGTATACCTGGCGTTGTCTGGTTTATTTTTAGACTGTAAATCACTTTTATATTGTCTTTATTGATATACAGAAACAACAGCACCATAGCcatacat                                                                                       <  1:4057119/149‑1 (MQ=255)
                                                                tCAGTCATTGTAGTCTGACAGGCAAAATTATTATTGCTATCACAGCTAACAAGTATACCTGGCGTTGTCTGGTTTATTTTTAGGCTGTAAATCACTTTTATATTGTCTTTATTGATATACAGAAACAACAGCACCATAGCCATACatat                                                                                     <  1:2618585/149‑1 (MQ=255)
                                                                tCAGTCATTGTAGTCTGACAGGCAAAATTATTATTGCTATCACAGCTAACAAGTATACCTGGCGTTGTCTGGTTTATTTTTAGACTGTAAATCACTTTTATATTGTCTTTATTGATATACAGAAACAACAGCACCATAGCCATACatat                                                                                     >  2:4013528/1‑149 (MQ=255)
                                                                  aGTCATTGTAGTCTGACAGGCAAAATTATTATTGCTATCACAGCTAAAAAGTATACCTGGCGTTGTCTGGTTTATTTTTAGACTGTAAATCACTTTTATATTGTCTTTATTGATATACAGAAACAACAGCACCATAGCCATACATATAg                                                                                   <  2:1080962/149‑1 (MQ=255)
                                                                   gTCATTGTAGTCTGACAGTCAAAATTATTATTGCTATCACAGCTAACAAGTATACCTGGCGTTGTCTGGTTTATTTTTAGGCTGTAAATCACTTTTATATTGTCTTTATTGATATACAGAAACAACAGCACCATAGCCATACATATAGc                                                                                  <  1:1976047/149‑1 (MQ=255)
                                                                   gTCATTGTAGTCTGACAGGCAAAATTATTATTGCTATCACAGCTAACAAGTATACCTGGCGTTGTCTGGTTTATTTTTAGACTGTAAATCACTTTTATATTGTCTTTATTGATATACAGAAACAACAGCACCATAGCCATACATATAGc                                                                                  <  1:4153924/149‑1 (MQ=255)
                                                                    tCATTGTAGTCTGACAGGCAAAATTATTATTGCTATCACAGCTAACAAGTATACCTGGCGTTGTCTGGTTTATTTTTAGGCTGTAAATCACTTTTATATTGTCTTTATTGATATACAGAAACAACAGCACCATAGCCATACATATAGcc                                                                                 <  1:3520766/149‑1 (MQ=255)
                                                                    tCATTGTAGTCTGACAGGCAAAATTATTATTGCTATCACAGCTAACAAGTATACCTGGCGTTGTCTGGTTTATTTTTAGACTGTAAATCACTTTTATATTGTCTTTATTGATATACAGAAACAACAGCACCATAGCCATACATATAGcc                                                                                 >  1:3388873/1‑149 (MQ=255)
                                                                       ttGTAGTCTGACAGGCAAAATTATTATTGCTATCACAGCTAACAAGTATACCTGGCGTTGTCTGGTTTATTTTTAGACTGTAAATCACTTTTATATTGTCTTTATTGATATACAGAAACAACAGCACCATAGCCATACATATAGCCata                                                                              <  1:1737978/149‑1 (MQ=255)
                                                                         gTAGTCTGACAGGCAAAATTATTATTGCTATCACAGCTAACAAGTATACCTGGCGTTGTCTGGTTTATTTTTAGGCTGTAAATCACTTTTATATTGTCTTTATTGATATACAGAAACAACAGCACCATAGCCATACATATAGCCataat                                                                            >  2:394353/1‑149 (MQ=255)
                                                                             tCTGACAGGCAAAATTATTATTGCTATCACAGCTAACAAGTATACCTGGCGTTGTCTGGTTTATTTTTAGACTGTAAATCACTTTTATATTGTCTTTATTGATATACAGAAACAACAGCACCATAGCCATACATATAGCCATAATaaaa                                                                        >  2:4324933/1‑149 (MQ=255)
                                                                              cTGACAGGCAAAATTATTATTGCTATCACAGCTAACAAGTATACCTGGCGTTGTCTGGTTTATTTTTAGACTGTAAATCACTTTTATATTGTCTTTATTGATATACAGAAACAACAGCACCATAGCCATACATATAGCCATAATaaaaa                                                                       <  1:1667413/149‑1 (MQ=255)
                                                                               tGACAGGCAAAATTATTATTGCTATCACAGCTAACAAGTATACCTGGCGTTGTCTGGTTTATTTTTAGGCTGTAAATCACTTTTATATTGTCTTTATTGATATACAGAAACAACAGCACCATAGCCATACATATAGCCATAATAAAAAg                                                                      <  2:153954/149‑1 (MQ=255)
                                                                                 aCAGGCAAAATTATTATTGCTATCACAGCTAACAAGTATACCTGGCGTTGTCTGGTTTATTTTTAGGCTGTAAATCACTTTTATATTGTCTTTATTGATATACAGAAACAACAGCACCATAGCCATACATATAGCCATAATAAAAAGGa                                                                    >  2:201405/1‑149 (MQ=255)
                                                                                   aGGCAAAATTATTATTGCTATCACAGCTAACAAGTATACCTGGCGTTGTCTGGTTTATTTTTAGACTGTAAATCACTTTTATATTGTCTTTATTGATATACAGAAACAACAGCACCATAGCCATACATATAGCCATAATAAAAAGGATa                                                                  <  1:3807189/149‑1 (MQ=255)
                                                                                    ggCAAAATTATTATTGCTATCACAGCTAACAAGTATACCTGGCGTTGTCTGGTTTATTTTTAGACTGTAAATCACTTTTATATTGTCTTTATTGATATACAGAAACAACAGCACCATAGCCATACATATAGCCATAATAAAAAGGATaa                                                                 <  1:3273423/149‑1 (MQ=255)
                                                                                     gCAAAATTATTATTGCTATCACAGCTAACAAGTATACCTGGCGTTGTCTGGTTTATTTTTAGACTGTAAATCACTTTTATATTGTCTTTATTGATATACAGAAACAACAGCACCATAGCCATACATATAGCCATAATAAAAAGGATAac                                                                >  2:2742237/1‑149 (MQ=255)
                                                                                      cAAAATTATTATTGCTATCACAGCTAACAAGTATACCTGGCGTTGTCTGGTTTATTTTTAGACTGTAAATCACTTTTATATTGTCTTTATTGATATACAGAAACAACAGCACCATAGCCATACATATAGCCATAATAAAAAGGATAaca                                                               <  2:1335272/149‑1 (MQ=255)
                                                                                      cAAAATTATTATTGCTATCACAGCTAACAAGTATACCTGGCGTTGTCTGGTTTATTTTTAGACTGTAAATCACTTTTATATTGTCTTTATTGATATACAGAAACAACAGCACCATAGCCATACATATAGCCATAATAAAAAGGATAaca                                                               <  1:1456276/149‑1 (MQ=255)
                                                                                          attattattGCTATCACAGCTAACAAGTATACCTGGCGTTGTCTGGTTTATTTTTAGACTGTAAATCACTTTTATATTGTCTTTATTGATATACAGAAACAACAGCACCATAGCCATACATATAGCCATAATAAAAAGGATAACACGaa                                                           <  1:2762590/149‑1 (MQ=255)
                                                                                            tattattGCTATCACAGCTAACAAGTATACCTGGCGTTGTCTGGTTTATTTTTAGACTGTAAATCACTTTTATATTGTCTTTATTGATATACAGAAACAACAGCACCATAGCCATACATATAGCCATAATAAAAAGGATAACACGAATc                                                         <  2:2372307/149‑1 (MQ=255)
                                                                                            tattattGCTATCACAGCTAACAAGTATACCTGGCGTTGTCTGGTTTATTTTTAGACTGTAAATCACTTTTATATTGTCTTTATTGATATACAGAAACAACAGCACCATAGCCATACATATAGCCATAATAAAAAGGATAACACGAATc                                                         >  1:602130/1‑149 (MQ=255)
                                                                                             attattTCTATCACAGCTAACAATTATACCTGGCTTTGTCTGGTTTATTTTTAGTCTGTAAATCTCTTTTATATTGTCTTTATTTATATACAGAAACAACAGCACCATATCCATACATATAGCCATAATAAAAAGGATAACACGAATCt                                                        <  1:1646558/149‑1 (MQ=255)
                                                                                              ttattGCTATCACAGCTAACAAGTATACCTGGCGTTGTCTGGTTTATTTTTAGACTGTAAATCACTTTTATATTGTCTTTATTGATATACAGAAACAACAGCACCATAGCCATACATATAGCCATAATAAAAAGGATAACACGAATCtt                                                       <  2:3493212/149‑1 (MQ=255)
                                                                                                attGCTATCACAGCTAACAAGTATACCTGGCGTTGTCTGGTTTATTTTTAGACTGTAAATCACTTTTATATTGTCTTTATTGATATACAGAAACAACAGCACCATAGCCATACATATAGCCATAATAAAAAGGATAACACGAATCtttt                                                     <  2:2432684/149‑1 (MQ=255)
                                                                                                  tGCTATCACAGCTAACAAGTATACCTGGCGTTGTCTGGTTTATTTTTAGGCTGTAAATCACTTTTATATTGTCTTTATTGATATACAGAAACAACAGCACCATAGCCATACATATAGCCATAATAAAAAGGATAACACGAATCTTTttc                                                   >  2:3594485/1‑149 (MQ=255)
                                                                                                     tATCACAGCTAACAAGTATACCTGGCGTTGTCTGGTTTATTTTTAGACTGTAAATCACTTTTATATTGTCTTTATTGATATACAGAAACAACAGCACCATAGCCATACATATAGCCATAATAAAAAGGATAACACGAATCTTTttcttc                                                <  2:2899004/149‑1 (MQ=255)
                                                                                                     tATCACAGCTAACAAGTATACCTGGCGTTGTCTGGTTTATTTTTAGACTGTAAATCACTTTTATATTGTCTTTATTGATATACAGAAACAACAGCACCATAGCCATACATATAGCCATAATAAAAAGGATAACACGAATCTTTttcttc                                                <  2:558389/149‑1 (MQ=255)
                                                                                                       tCACAGCTAACAAGTATACCTGGCGTTGTCTGGTTTATTTTTAGGCTGTAAATCACTTTTATATTGTCTTTATTGATATACAGAAACAACAGCACCATAGCCATACATATAGCCATAATAAAAAGGATAACACGAATCTTTTTCTTCAg                                              <  1:4656773/149‑1 (MQ=255)
                                                                                                          caGCTAACAAGTATACCTGGCGTTGTCTGGTTTATTTTTAGGCTGTAAATCACTTTTATATTGTCTTTATTGATATACAGAAACAACAGCACCATAGCCATACATATAGCCATAATAAAAAGGATAACACGAATCTTTTTCTTCAGata                                           <  2:487879/149‑1 (MQ=255)
                                                                                                            gCTAACAAGTATACCTGGCGTTGTCTGGTTTATTTTTAGACTGTAAATCACTTTTATATTGTCTTTATTGATATACAGAAACAACAGCACCATAGCCATACATATAGCCATAATAAAAAGGATAACACGAATCTTTTTCTTCAGataaa                                         >  1:3027285/1‑149 (MQ=255)
                                                                                                             cTAACAAGTATACCTGGCGTTGTCTGGTTTATTTTTAGACTGTAAATCACTTTTATATTGTCTTTATTGATATACAGAAACAACAGCACCATAGCCATACATATAGCCATAATAAAAAGGATAACACGAATCTTTTTCTTCAGataaaa                                        <  2:2443699/149‑1 (MQ=255)
                                                                                                             cTAACAAGTATACCTGGCGTTGTCTGGTTTATTTTTAGACTGTAAATCACTTTTATATTGTCTTTATTGATATACAGAAACAACAGCACCATAGCCATACATATAGCCATAATAAAAAGGATAACACGAATCTTTTTCTTCAGataaaa                                        <  2:364526/149‑1 (MQ=255)
                                                                                                               aaCAAGTATACCTGGCGTTGTCTGGTTTATTTTTAGACTGTAAATCACTTTTATATTGTCTTTATTGATATACAGAAACAACAGCACCATAGCCATACATATAGCCATAATAAAAAGGATAACACGAATCTTTTTCTTCAGataaaata                                      <  1:857444/149‑1 (MQ=255)
                                                                                                               aaCAAGTATACCTGGCGTTGTCTGGTTTATTTTTAGACTGTAAATCACTTTTATATTGTCTTTATTGATATACAGAAACAACAGCACCATAGCCATACATATAGCCATAATAAAAAGGATAACACGAATCTTTTTCTTCAGataaaata                                      <  1:2323919/149‑1 (MQ=255)
                                                                                                                aCAAGTATACCTGGCGTTGTCTGGTTTATTTTTAGACTGTAAATCACTTTTATATTGTCTTTATTGATATAAAGAAACAACAGCACCATAGCCATACATATAGCCATAATAAAAAGGATAACACGAATCTTTTTCTTCAGataaaataa                                     >  2:1672428/1‑149 (MQ=255)
                                                                                                                aCAAGTATACCTGGCGTTGTCTGGTTTATTTTTAGACTGTAAATAACTTTTATATTGTCTTTATTGATATACAGAAACAACAGCACCAAAGCCAAACATATAGCCAAAATAAAAAAGATAAAACGAATATTTTTATTCAGataaaataa                                     >  1:1307426/1‑149 (MQ=255)
                                                                                                                   aGTATACCTGGCGTTGTCTGGTTTATTTTTAGACTGTAAATCACTTTTATATTGTCTTTATTGATATACAGAAACAACAGCACCATAGCCATACATATAGCCATAATAAAAAGGATAACACGAATCTTTTTCTTCAGATaaaataaaat                                  >  2:2946504/1‑149 (MQ=255)
                                                                                                                    gTATACCTGGCGTTGTCTGGTTTATTTTTAGGCTGTAAATCACTTTTATATTGTCTTTATTGATATACAGAAACAACAGCACCATAGCCATACATATATCCATAATAAAAAGGATAACACGAATCTTTTTCTTCAGATaaaataaaata                                 >  2:898667/1‑148 (MQ=255)
                                                                                                                         ccTGGCGTTGTCTGGTTTATTTTTAGACTGTAAATCACTTTTATATTGTCTTTATTGATATACAGAAACAACAGCACCATAGCCATACATATAGCCATAATAAAAAGGATAACACGAATCTTTTTCTTCAGATAAAATAAAATCTCAAt                            <  1:4663637/149‑1 (MQ=255)
                                                                                                                           tGGCGTTGTCTGGTTTATTTTTAGGCTGTAAATCACTTTTATATTGTCTTTATTGATATACAGAAACAACAGCACCATAGCCATACATATAGCCATAATAAAAAGGATAACACGAATCTTTTTCTTCAGATAAAATAAAATCTCAATCa                          <  2:1909482/149‑1 (MQ=255)
                                                                                                                              cGTTGTCTGGTTTATTTTTAGGCTGTAAATCACTTTTATATTGTCTTTATTGATATACAGAAACAACAGCACCATAGCCATACATATAGCCATAATAAAAAGGATAACACGAATCTTTTTCTTCAGATAAAATAAAATCTCAATCAGAt                       <  2:2028445/149‑1 (MQ=255)
                                                                                                                                  gTCTGGTTTATTTTTAGACTGTAAATCACTTTTATATTGTCTTTATTGATATACAGAAACAACAGCACCATAGCCATACATATAGCCATAATAAAAAGGATAACACGAATCTTTTTCTTCAGATAAAATAAAATCTCAATCAGATCGat                   >  1:1114594/1‑149 (MQ=255)
                                                                                                                                    cTGGTTTATTTTTAGGCTGTAAATCACTTTTATATTGTCTTTATTGATATACAGAAACAACAGCACCATAGCCATACATATAGCCATAATAAAAAGGATAACACGAATCTTTTTCTTCAGATAAAATAAAATCTCAATCAGATCGatat                 <  1:3229497/149‑1 (MQ=255)
                                                                                                                                    cTGGTTTATTTTTAGACTGTAAATCACTTTTATATTGTCTTTATTGATATACAGAAACAACAGCACCATAGCCATACATATAGCCATAATAAAAAGGATAACACGAATCTTTTTCTTCAGATAAAATAAAATCTCAATCAGATCGatat                 >  2:2288538/1‑149 (MQ=255)
                                                                                                                                      ggTTTATTTTTAGGCTGTAAATCACTTTTATATTGTCTTTATTGATATACAGAAACAACAGCACCATAGCCATACATATAGCCATAATAAAAAGGATAACACGAATCTTTTTCTTCAGATAAAATAAAATCTCAATCAGATCGATATcc               <  1:4639336/149‑1 (MQ=255)
                                                                                                                                      ggTTTATTTTTAGACTGTAAATCACTTTTATATTGTCTTTATTGATATACAGAAACAACAGCACCATAGCCATACATATAGCCATAATAAAAAGGATAACACGAATCTTTTTCTTCAGATAAAATAAAATCTCAATCAGATCGATATcc               >  2:1430828/1‑149 (MQ=255)
                                                                                                                                      ggTTTATTTTTAGACTGTAAATCACTTTTATATTGTCTTTATTGATATACAGAAACAACAGCACCATAGCCATACATATAGCCATAATAAAAAGGATAACACGAATCTTTTTCTTCAGATAAAATAAAATCTCAATCAGATCGATATcc               >  1:745330/1‑149 (MQ=255)
                                                                                                                                        tttATTTTTAGGCTGTAAATCACTTTTATATTGTCTTTATTGATATACAGAAACAACAGCACCATAGCCATACATATAGCCATAATAAAAAGGATAACACGAATCTTTTTCTTCAGATAAAATAAAATCTCAATCAGATCGATATCCTg             >  1:4214459/1‑149 (MQ=255)
                                                                                                                                          tATTTTTAGGCTGTAAATCACTTTTATATTGTCTTTATTGATATACAGAAACAACAGCACCATAGCCATACATATAGCCATAATAAAAAGGATAACACGAATCTTTTTCTTCAGATAAAATAAAATCTCAATCAGATCGATATCCTGtt           <  2:2379040/149‑1 (MQ=255)
                                                                                                                                            tttttAGACTGTAAATCACTTTTATATTGTCTTTATTGATATACAGAAACAACAGCACCATAGCCATACATATAGCCATAATAAAAAGGATAACACGAATCTTTTTCTTCAGATAAAATAAAATCTCAATCAGATCGATATCCTGTTTg         >  2:4167007/1‑149 (MQ=255)
                                                                                                                                                tAGACTGTAAATCACTTTTATATTGTCTTTATTGATATACAGAAACAACAGCACCATAGCCATACATATAGCCATAATAAAAAGGATAACACGAATCTTTTTCTTCAGATAAAATAAAATCTCAATCAGATCGATATCCTGTTTGAttt     <  1:4100628/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
TATTTCTCTGTTCTTGACATCTGGGCTGGTGTGAAAAAATGTAGTAATTCGCTGAATAACATCTTCAGTCATTGTAGTCTGACAGGCAAAATTATTATTGCTATCACAGCTAACAAGTATACCTGGCGTTGTCTGGTTTATTTTTAGGCTGTAAATCACTTTTATATTGTCTTTATTGATATACAGAAACAACAGCACCATAGCCATACATATAGCCATAATAAAAAGGATAACACGAATCTTTTTCTTCAGATAAAATAAAATCTCAATCAGATCGATATCCTGTTTGATTTGTT  >  CP009273/320773‑321068

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: