Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,404,951 A→T 5.5% intergenic (+25/‑61) fis → / → yhdJ global DNA‑binding transcriptional dual regulator/DNA adenine methyltransferase, SAM‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,404,9510AT 5.5% 739.9 / 24.1 234intergenic (+25/‑61)fis/yhdJglobal DNA‑binding transcriptional dual regulator/DNA adenine methyltransferase, SAM‑dependent
Reads supporting (aligned to +/- strand):  ref base A (106/115);  new base T (4/9);  total (110/124)
Fisher's exact test for biased strand distribution p-value = 2.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

Only 100 of 246 total aligned reads displayed.

AGAACAGCCCCTGTTGGACATGGTGATGCAATACACCCGTGGTAACCAGACCCGTGCTGCGCTGATGATGGGCATCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATGAAGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATCC  >  CP009273/3404803‑3405099
                                                                                                                                                    |                                                                                                                                                    
 gAACAGCCCCTGTTGGACATGGTGATGCAATACACCCGTGGTAACCAGACCCGTGCTGCGCTGATGATGGGCATCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTaa                                                                                                                                                     >  2:8280/1‑149 (MQ=255)
    cAGCCCCTGTTGGACATGGTGATGCAATACACCCGTGGTAACCAGACCCGTGCTGCGCTGATGATGGGCATCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTaaaaa                                                                                                                                                  >  2:3360238/1‑149 (MQ=255)
    cAGCCCCTGTTGGACATGGTGATGCAATACACCCGTGGTAACCAGACCCGTGCTGCGCTGATGATGGGCATCAACCGTGGTACGCGGCGTAAAAAATTGAAAAAATACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTaaaaa                                                                                                                                                  >  1:669703/1‑149 (MQ=255)
     aGCCCCTGTTGGACATGGTGATGCAATACACCCGTGGTAACCAGACCCGTGCTGCGCTGATGATGGGCATCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAg                                                                                                                                                 <  2:451089/149‑1 (MQ=255)
        cccTGTTGGACATGGTGATGCAATACACCCGTGGTAACCAGACCCGTGCTGCGCTGATGATGGGCATCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGgcg                                                                                                                                              <  1:1672111/149‑1 (MQ=255)
            gTTGGACATGGTGATGCAATACACCCGTGGTAACCAGACCCGTGCTGCGCTGATGATGGGCATCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTTAAAAGGCGCTAc                                                                                                                                          <  1:840527/149‑1 (MQ=255)
             ttGGACATGGTGATGCAATACACCCGTGGTAACCAGACCCGTGCTGCGCTGATGATGGGCATCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACt                                                                                                                                         >  2:1560918/1‑149 (MQ=255)
                  cATGGTGATGCAATACACCCGTGGTAACCAGACCCGTGCTGCGCTGATGATGGGCATCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCa                                                                                                                                    <  1:3022216/149‑1 (MQ=255)
                   aTGGTGATGCAATACACCCGTGGTAACCAGACCCGTGCTGCGCTGATGATGGGCATCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCAt                                                                                                                                   >  1:3188724/1‑149 (MQ=255)
                    tGGTGATGCAATACACCCGTGGTAACCAGACCCGTGCTGCGCTGATGATGGGCATCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACGCGGCATg                                                                                                                                  <  2:1893055/149‑1 (MQ=255)
                      gTGATGCAATACACCCGTGGTAACCAGACCCGTGCTGCGCTGATGATGGGCATCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATggg                                                                                                                                <  2:1699963/149‑1 (MQ=255)
                          tGCAATACACCCGTGGTAACCAGACCCGTGCTGCGCTGATGATGGGCATCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGacg                                                                                                                            >  1:2235988/1‑147 (MQ=255)
                          tGCAATACACCCGTGGTAACCAGACCCGTGCTGCGCTGATGATGGGCATCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAg                                                                                                                            >  1:3406019/1‑149 (MQ=255)
                           gCAATACACCCGTGGTAACCAGACCCGTGCTGCGCTGATGATGGGCATCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAgc                                                                                                                           >  2:109557/1‑149 (MQ=255)
                            cAATACACCCGTGGTAACCAGACCCGTGCTGCGCTGATGATGGGCATCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATAAAAAAGGCGCTACTCGGCATGGGGAAgcg                                                                                                                          <  1:1548431/149‑1 (MQ=255)
                                  acCCGTGGTAACCAGACCCGTGCTGCGCTGATGATGGGCATCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCtttt                                                                                                                    <  2:4048510/149‑1 (MQ=255)
                                   cccGTGGTAACCAGACCCGTGCTGCGCTGATGATGGGCATCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCttttt                                                                                                                   <  2:1011148/149‑1 (MQ=255)
                                      gTGGTAACCCGACCCGTGCTGCGCTGATGATGGGCATCAACCGTGGTACGATGCGTAAAAAATTGAAAAAATACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTtat                                                                                                                <  1:3784110/149‑1 (MQ=255)
                                      gTGGTAACCAGACCCGTGCTGCGCTGATGATGGGCATCACCCGTGGTACGCTGCGTAAAAAATTGAAAAAATACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTtat                                                                                                                <  2:1883584/149‑1 (MQ=255)
                                           aaCCAGACCCGTGCTGCGCTGATGATGGGCATCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGgtg                                                                                                           >  2:1478597/1‑149 (MQ=255)
                                             ccAGACCCGTGCTGCGCTGATGATGGGCATCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTc                                                                                                         >  2:1671414/1‑149 (MQ=255)
                                                 aCCCGTGCTGCGCTGATGATGGGCATCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGACGCGCCTTTTTTATAGGTGTCACaa                                                                                                     >  1:3101166/1‑149 (MQ=255)
                                                   ccGTGCTGCGCTGATGATGGGCATCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAg                                                                                                   >  1:2671433/1‑149 (MQ=255)
                                                   ccGTGCTGCGCTGATGATGGGCATCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAg                                                                                                   >  2:1292565/1‑149 (MQ=255)
                                                    cGTGCTGCGCTGATGATGGGCATCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAgg                                                                                                  <  1:2199671/149‑1 (MQ=255)
                                                     gTGCTGCGCTGATGATGGGCATCAACCGTGGTCCGCTGCGTAAAAAATTGAAAAAATACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGACGCGCCTTTTTTATAGGTGTCACCAAggg                                                                                                 >  2:2127507/1‑149 (MQ=255)
                                                     gTGCTGCGCTGATGATGGGCATCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAggg                                                                                                 <  1:2900547/149‑1 (MQ=255)
                                                     gTGCTGCGCTGATGATGGGCATCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAggg                                                                                                 <  1:1071932/149‑1 (MQ=255)
                                                     gTGCTGCGCTGATGATGGGCATCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCGACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAggg                                                                                                 <  1:2461375/149‑1 (MQ=255)
                                                      tgctgcGCTGATGATGGGCATCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGa                                                                                                <  1:781261/149‑1 (MQ=255)
                                                       gctgcGCTGATGATGGGCATCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATACGGCATGATCTAAGTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAg                                                                                               <  1:1008937/149‑1 (MQ=255)
                                                       gctgcGCTGATGATGGGCATCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAg                                                                                               <  2:2335894/149‑1 (MQ=255)
                                                       gctgcGCTGATGATGGGCATCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAg                                                                                               <  1:1008976/149‑1 (MQ=255)
                                                        ctgcGCTGATGATGGGCATCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGt                                                                                              <  2:2268699/149‑1 (MQ=255)
                                                          gcgcTGATGATGGGCATCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGa                                                                                            >  2:2192341/1‑149 (MQ=255)
                                                          gcgcTGATGATGGGCATCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGa                                                                                            <  1:1936774/149‑1 (MQ=255)
                                                           cgcTGATGATGGGCATCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGAc                                                                                           <  1:3798797/149‑1 (MQ=255)
                                                           cgcTGATGATGGGCATCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGAc                                                                                           <  2:3327578/149‑1 (MQ=255)
                                                                    tgGGCATCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGaa                                                                                   <  1:1322608/148‑1 (MQ=255)
                                                                    tgGGCATCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAAc                                                                                  <  1:3628140/149‑1 (MQ=255)
                                                                         aTCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGAt                                                                             >  2:893516/1‑149 (MQ=255)
                                                                          tCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGAtg                                                                            >  2:240623/1‑149 (MQ=255)
                                                                          tCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGAtg                                                                            >  2:4087500/1‑149 (MQ=255)
                                                                          tCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCGACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGAtg                                                                            <  1:1496847/149‑1 (MQ=255)
                                                                           cAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGAtgt                                                                           <  2:2614449/149‑1 (MQ=255)
                                                                           cAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGAtgt                                                                           >  1:3273120/1‑149 (MQ=255)
                                                                           cAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATACGGCATGAACTAAGTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAACGGGAGTGACCATGAGAACAGGAtgt                                                                           <  1:1464164/149‑1 (MQ=255)
                                                                               cGTGGTACGCTGCGTAAAAAATTGAAAAAATACGGCATGAACTAATTGAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAAc                                                                       <  2:475961/149‑1 (MQ=255)
                                                                               cGTGGTACGCTGCGTAAAAAATTGAAAAAATACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAAc                                                                       <  1:2539304/149‑1 (MQ=255)
                                                                               cGTGGTACGCTGCGTAAAAAATTGAAAAAATACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAAc                                                                       <  2:1320977/149‑1 (MQ=255)
                                                                                  ggTACGCTGCGTAAAAAATTGAAAAAATACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGACGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGa                                                                    >  1:3900504/1‑149 (MQ=255)
                                                                                    tACGCTGCGTAAAAAATTGAAAAAATACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGAcc                                                                  >  2:3741418/1‑149 (MQ=255)
                                                                                    tACGCTGCGTAAAAAATTGAAAAAATACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGAcc                                                                  >  2:3008206/1‑149 (MQ=255)
                                                                                    tACGCTGCGTAAAAAATTGAAAAAATACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGAcc                                                                  >  2:3741609/1‑149 (MQ=255)
                                                                                      cGCTGCGTAAAAAATTGAAAAAATACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCg                                                                >  2:1766900/1‑149 (MQ=255)
                                                                                          gCGTAAAAAATTGAAAAAATACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGACGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGttt                                                            >  1:3430739/1‑149 (MQ=255)
                                                                                            gaaaaaaaTTGAAAAAATACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTgg                                                          <  1:2738943/147‑1 (MQ=255)
                                                                                            gTAAAAAATTGAAAAAATACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGACGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTgg                                                          >  1:2959471/1‑149 (MQ=255)
                                                                                            gTAAAAAATTGAAAAAATACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTgg                                                          >  1:900276/1‑149 (MQ=255)
                                                                                             tAAAAAATTGAAAAAATACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGt                                                         >  1:652934/1‑149 (MQ=255)
                                                                                               aaaaaTTGAAAAAATACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTaa                                                       >  1:3986913/1‑149 (MQ=255)
                                                                                               aaaaaTTGAAAAAATACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTaa                                                       >  2:2013320/1‑149 (MQ=255)
                                                                                                  aaTTTAAAAAATACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCATTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATGa                                                    <  1:369877/149‑1 (MQ=255)
                                                                                                  aaTTGAAAAAATACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATGa                                                    >  2:2821964/1‑149 (MQ=255)
                                                                                                   aTTGAAAAAATACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATGaa                                                   >  2:2098601/1‑149 (MQ=255)
                                                                                                    ttGAAAAAATACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTGTTTATTGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATGAAg                                                  >  2:1411703/1‑149 (MQ=255)
                                                                                                      gAAAAAATACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATGAAGCt                                                >  1:1859747/1‑149 (MQ=255)
                                                                                                      gAAAAAATACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATGAAGCt                                                >  1:2730114/1‑149 (MQ=255)
                                                                                                       aaaaaaTACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTTAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATGAAGCTa                                               <  2:77296/149‑1 (MQ=255)
                                                                                                           aaTACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTTAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATGAAGCTAAGAc                                           <  2:3136632/149‑1 (MQ=255)
                                                                                                           aaTACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATGAAGCTAAGAc                                           >  1:2442000/1‑149 (MQ=255)
                                                                                                            aTACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATGAAGCTAAGAcc                                          >  2:3681986/1‑149 (MQ=255)
                                                                                                             tACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGACGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATGAAGCTAAGACCa                                         >  2:3787852/1‑149 (MQ=255)
                                                                                                              aCGGCATTAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATGAAGCTAAGACCat                                        <  2:21514/149‑1 (MQ=255)
                                                                                                              aCGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATGAAGCTAAGACCat                                        <  2:3971749/149‑1 (MQ=255)
                                                                                                                  cATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCGACTCGGCATGGGTTAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATGAAGCTAAGACCattatt                                    <  2:342292/149‑1 (MQ=255)
                                                                                                                   aTGAACTAATTCATGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATGAAGCTAAGACCATTATTc                                   <  2:669703/149‑1 (MQ=255)
                                                                                                                    tGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATGAAGCTAAGACCATTATTCa                                  >  2:1822880/1‑149 (MQ=255)
                                                                                                                    tGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATGAAGCTAAGACCATTATTCa                                  >  1:943388/1‑149 (MQ=255)
                                                                                                                     gAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATGAAGCTAAGACCATTATTCAc                                 <  1:533360/149‑1 (MQ=255)
                                                                                                                         tAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGACGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGCTGTGAACCGACCCGGTTTGGTAATGAAGCTAAGACCATTATTCACGGTg                             >  2:35432/1‑149 (MQ=255)
                                                                                                                          aaTTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATGAAGCTAAGACCATTATTCACGGTGa                            >  1:348862/1‑149 (MQ=255)
                                                                                                                               aGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATGAAGCTAAGACCATTATTCACGGTGATgcgc                       <  2:2317608/149‑1 (MQ=255)
                                                                                                                                ggTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATGAAGCTAAGACCATTATTCACGGTGATGCGCt                      <  2:4083314/149‑1 (MQ=255)
                                                                                                                                   tAGCTACCTGCTTGATTACAAAGGCGCTCCTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGCACAGGATGTGAACCGACCCGGTTTGGTAATGAAGCTAAGACCATTATTCACGGTGATGCGCTTGc                   <  1:3752688/149‑1 (MQ=255)
                                                                                                                                   tAGCTAAATGCTTGATTTAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATGAAGCTAAGACCATTATTCACGGTGATGCGCTTGc                   <  1:3960019/149‑1 (MQ=255)
                                                                                                                                    aGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATGAAGCTAAGACCATTATTCACGGTGATGCGCTTGcc                  <  1:438563/149‑1 (MQ=255)
                                                                                                                                       tAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATGAAGCTACGACCATTATTCCCGGTGATGCGCTTGCCGaa               >  2:1992790/1‑149 (MQ=255)
                                                                                                                                         aaTGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATGAAGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACt             <  2:3188724/149‑1 (MQ=255)
                                                                                                                                         aaTGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATGAAGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACt             >  1:2129081/1‑149 (MQ=255)
                                                                                                                                          aTGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATGAAGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACtt            <  2:1733733/149‑1 (MQ=255)
                                                                                                                                            gCTTGATTTAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATGAAGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTaa          >  2:265933/1‑149 (MQ=255)
                                                                                                                                            gCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAGAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATGAAGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTaa          >  1:120344/1‑149 (MQ=255)
                                                                                                                                              ttGATTTAAAAGGCGCGCCTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATGAAGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTaaaa        <  1:448925/149‑1 (MQ=255)
                                                                                                                                              ttGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATGAAGCTCAGACCATTATTCACGGTGATGCGCTTGCCGAACTTaaaa        >  2:1974383/1‑149 (MQ=255)
                                                                                                                                              ttGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGGGAACCGACCCGGTTTGGTAATGAAGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTaaaa        >  1:974610/1‑149 (MQ=255)
                                                                                                                                                gATTAAACAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATGAAGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAg      <  2:3624742/149‑1 (MQ=255)
                                                                                                                                                gATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATGAAGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAg      >  2:2631078/1‑149 (MQ=255)
                                                                                                                                                 aTTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATGAAGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGa     <  2:118686/149‑1 (MQ=255)
                                                                                                                                                 aTTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATGAAGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGa     <  2:3030270/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
AGAACAGCCCCTGTTGGACATGGTGATGCAATACACCCGTGGTAACCAGACCCGTGCTGCGCTGATGATGGGCATCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATGAAGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATCC  >  CP009273/3404803‑3405099

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: