Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,219,533 A→G 93.2% intergenic (+170/+202) ycgI → / ← minE pseudogene/cell division topological specificity factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,219,5330AG93.2% 857.0 / 36.2 250intergenic (+170/+202)ycgI/minEpseudogene/cell division topological specificity factor
Reads supporting (aligned to +/- strand):  ref base A (7/10);  new base G (120/113);  total (127/123)
Fisher's exact test for biased strand distribution p-value = 4.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.45e-01

Only 100 of 250 total aligned reads displayed.

GGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAACCTCATCCTATGGGAGAACAATGAATAATGAAATTGCGGGGTTATCATCTCCCAGTATATCCATACTAACAATAAGGTTATTTACTCAACCAGGCATAAACATTTTGTTTTGTGCGTGGGAACAGCCTTAAGGTGTAAAGGGGGAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAA  >  CP009273/1219388‑1219679
                                                                                                                                                 |                                                                                                                                                  
 gcgcAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAACCTCATCCTATGGGAGAACAATGAATAATGAAATTGCGGGGTTATCATCTCCCAGTATATCCATACTAACGATAAg                                                                                                                                               <  1:1659172/150‑1 (MQ=255)
    cAGTTTTCGTTTACAGGATGTTGAACGGGTAAATTCTGGGGCACAAAAAGCCCGCCAGTTACGGCGGGAAACCTCATCCTATGGGAGAACAATGAATAATGAAATTGCGGGGTTATCATCTCCCAGTATATCCATACTAACGAGACGGtt                                                                                                                                            >  2:1617962/1‑150 (MQ=255)
    cAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAACCTCATCCTATGGGAGAACAATGAATAATGAAATTGCGGGGTTATCATCTCCCAGTATATCCATACTAACGATAAGGtt                                                                                                                                            >  2:1618006/1‑150 (MQ=255)
    cAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAACCTCATCCTATGGGAGAACAATGAATAATGAAATTGCGGGGTTATCATCTCCCAGTATATCCATACTAACAATAAGGtt                                                                                                                                            <  1:775340/150‑1 (MQ=255)
       ttttCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAACCTCATCCTATGGGAGAACAATGAATAATGAAATTGCGGGGTTATCATCTCCCAGTATATCCATACTAACGATAAGGTTAtt                                                                                                                                         <  1:778301/150‑1 (MQ=255)
       ttttCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAACCTCATCCTATGGGAGAACAATGAATAATGAAATTGCGGGGTTATCATCTCCCAGTATATCCATACTAACGATAAGGTTAtt                                                                                                                                         <  1:1326112/150‑1 (MQ=255)
        tttCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAACCTCATCCTATGGGAGAACAATGAATAATGAAATTGCGGGGTTATCATCTCCCAGTATATCCATACTAACGATAAggcttatc                                                                                                                                        <  2:2489923/150‑7 (MQ=255)
        tttCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAACCTCATCCTATGGGAGAACAATGAATAATGAAATTGCGGGGTTATCATCTCCCAGTATATCCATACTAACGATAAGGTTAttt                                                                                                                                        >  1:2857892/1‑150 (MQ=255)
              ttACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAACCTCATCCTATGGGAGAACAATGAATAATGAAATTGCGGGGTTATCATCTCCCAGTATATCCATACTAACGATAAGGTTATTTTCTCaa                                                                                                                                  >  1:2370968/1‑150 (MQ=255)
               tACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAACCTCATCCTATGGGAGAACAATGAATAATGAAATTGCGGGGTTATCATCTCCCAGTATTTCCATACTAACGATAAGGTTATTTACTCAAc                                                                                                                                 <  1:3963001/150‑1 (MQ=255)
                     aTGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAACCTCATCCTATGGGAGAACAATGAATAATGAAATTGCGGGGTTATCATCTCCCAGTATATCCATACTAACGATAAGGTTATTTACTCAACCAGGCa                                                                                                                           >  1:1926377/1‑150 (MQ=255)
                         tGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAACCTCATCCTATGGGAGAACAATGAATAATGAAATTGCGGGGTTATCATCTCCCAGTATATCCATACTAACGATAAGGTTATTTACTCAACCAGGCATaaa                                                                                                                       >  1:2827562/1‑150 (MQ=255)
                          gAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAACCTCATCCTATGGGAGATCAATGAATAATGAAATTGCGGGGTTATCATCTCCGAGTATATCCATACTAACGATAAGGTTATTTACTCAACCAGGCATAAAc                                                                                                                      <  1:255908/150‑1 (MQ=255)
                          gAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAACCTCATCCTATGGGAGAACAATGAATAATGAAATTGCGGGGTTATCATCTCCCAGTATATCCATACTAACGATAAGGTTATTTACTCAACCAGGCATAAAc                                                                                                                      >  2:4254190/1‑150 (MQ=255)
                            aaGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAACCTCATCCTATGGGAGAACAATGAATAATGAAATTGCGGGGTTATCATCTCCCAGTATATCCATACTAACGATAAGGTTATTTACTCAACCAGGCATAAACAt                                                                                                                    >  2:2958663/1‑150 (MQ=255)
                              gggAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAACCTCATCCTATGGGAGAACAATGAATAATGAAATTGCGGGGTTATCATCTCCCAGTATATCCATACTAACGATAAGGTTATTTACTCAACCAGGCATAAACAttt                                                                                                                  <  1:1105726/150‑1 (MQ=255)
                               ggAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAACCTCATCCTATGGGAGAACAATGAATAATGAAATTGCGGGGTTATCATCTCCCAGTATATCCATACTAACGATAAGGTTATTTACTCAACCAGGCATAAACAtttt                                                                                                                 <  1:259370/150‑1 (MQ=255)
                                  aaaTTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAACCTCATCCTATGGGAGAACAATGAATAATGAAATTGCGGGGTTATCATCTCCCAGTATATCCATACTAACGATAAGGTTATTTACTCAACCAGGCATAAACAttttgtt                                                                                                              <  1:4098918/150‑1 (MQ=255)
                                   aaTTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAACCTCATCCTATGGGAGAACAATGAATAATGAAATTGCGGGGTTATCATCTCCCAGTATATCCATACTAACGATAAGGTTATTTACTCAACCAGGCATAAACAttttgttt                                                                                                             <  1:3392039/150‑1 (MQ=255)
                                    aTTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAACCTCATCCTATGGGAGAACAATGAATAATGAAATTGCGGGGTTATCATCTCCCAGTATATCCATACTAACGATAAGGTTATTTACTCAACCAGGCATAAACAttttgtttt                                                                                                            >  1:3657625/1‑150 (MQ=255)
                                    aTTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAACCTCATCCTATGGGAGAACAATGAATAATGAAATTGCGGGGTTATCATCTCCCAGTATATCCATACTAACGATAAGGTTATTTACTCAACCAGGCATAAACAttttgtttt                                                                                                            <  2:2476984/150‑1 (MQ=255)
                                    aTTCAGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAACCTCATCCTATGGGAGAACAATGAATAATGAAATTGCGGGGTTATCATCTCCCAGTATATCCATACTAACGATAAGGTTATTTACTCAACCAGGCATAAACAttttgtttt                                                                                                            <  1:1643565/150‑1 (MQ=255)
                                     ttCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAACCTCATCCTATGGGAGAACAATGAATCATGAAATTGCGGGGTTATCATCTCCCAGTATATCCATACTAAAGATAAGGTTATTTACTCAACCAGGCATAAACAttttgttttg                                                                                                           <  2:2793910/150‑1 (MQ=255)
                                          gggCAAAAAAAGCCCGCCAGTTACGGCGGGAAACCTCATCCTATGGGAGAACAATGAATAATGAAATTGCGGGGTTATCATCTCCCAGTATATCCATACTAACGATAAGGTTATTTACTCAACCAGGCATAAACATTTTGTTTTGTGCGt                                                                                                      >  1:904578/1‑150 (MQ=255)
                                          gggCAAAAAAAGCCCGCCAGTTACGGCGGGAAACCTCATCCTATGGGAGAACAATGAATAATGAAATTGCGGGGTTATCATCTCCCAGTATATCCATACTAACGATAAGGTTATTTACTCAACCAGGCATAAACATTTTGTTTTGTGCGt                                                                                                      >  1:3184053/1‑150 (MQ=255)
                                          gggCAAAAAAAGCCCGCCAGTTACGGCGGGAAACCTCATCCTATGGGAGAACAATGAATAATGAAATTGCGGGGTTATCATCTCCCAGTATATCCATACTAACAATAAGGTTATTTACTCAACCAGGCATAAACATTTTGTTTTGTGCGt                                                                                                      <  2:4057952/150‑1 (MQ=255)
                                               aaaaaaGCCCGCCAGTTACGGCGGGAAACCTCATCCTATGGGAGAACAATGAATAATGAAATTGCGGGGTTATCATCTCCCAGTATATCCATACTAACGATAAGGTTATTTACTCAACCAGGCATAAACATTTTGTTTTGTGCGTGGGaa                                                                                                 <  2:4368879/150‑1 (MQ=255)
                                                aaaaaGCCCGCCAGTTACGGCGGGAAACCTCATCCTATGGGAGAACAATGAATAATGAAATTGCGGGGTTATCATCTCCCAGTATATCCATACTAACGATAAGGTTATTTACTCAACCAGGCATAAACATTTTGTTTTGTGCGTGGGAAc                                                                                                <  2:3685983/150‑1 (MQ=255)
                                                aaaaaGCCCGCCAGTTACGGCGGGAAACCTCATCCTATGGGAGAACAATGAATAATGAAATTGCGGGGTTATCATCTCCCAGTATATCCATACTAACAATAAGGTTATTTACTCAACCAGGCATAAACATTTTGTTTTGTGCGTGGGAAc                                                                                                <  2:1904291/150‑1 (MQ=255)
                                                  aaaGCCCGCCAGTTACGGCGGGAAACCTCATCCTATGGGAGAACAATGAATAATGAAATTGCGGGGTTATCATCTCCCAGTATATCCATACTAACGATAAGGTTATTTACTCAACCAGGCATAAACATTTTGTTTTGTGCGTGGGAACAg                                                                                              >  1:1037679/1‑150 (MQ=255)
                                                   aaGCCCGCCAGTTACGGCGGGAAACCTCATCCTATGGGAGAACAATGAATAATGAAATTGCGGGGTTATCATCTCCCAGTATATCCATACTAACGATAAGGTTATTTACTCAACCAGGCATAAACATTTTGTTTTGTGCGTGGGAACAGc                                                                                             >  2:1414388/1‑150 (MQ=255)
                                                    aGCCCGCCAGTTACGGCGGGAAACCTCATCCTATGGGAGAACAATGAATAATGAAATTGCGGGGTTATCATCTCCCAGTATATCCATACTAACGATAAGGTTATTTACTCAACCAGGCATAAACATTTTGTTTTGTGCGTGGGAACAGcc                                                                                            <  2:1344417/150‑1 (MQ=255)
                                                      cccGCCAGTTACGTCGGGAAACCTCATACAATGGGAGAACAATGAATAAGGAAATTGCGGGGTTATCAACTCCCAGTATTTCCATACTAACGATAAGGTTATTTACTCAACCAGGCATAAACATTTGGTTTTGTGCGTGGGAACAGCCtt                                                                                          <  2:4205704/150‑1 (MQ=255)
                                                      cccGCCAGTTACGGCGGGAAACCTCATCCTATGGGAGAACAATGAATAATGAAATTGCGGGGTTATCATCTCCCAGTATATCCATACTAACAATAAGGTTATTTACTCAACCAGGCATAAACATTTTGTTTTGTGCGTGGGAACAGCCtt                                                                                          >  2:4286030/1‑150 (MQ=255)
                                                          ccAGTTACGGCGGGAAACCTCATCCTATGGGAGAACAATGAATAATGAAATTGCGGGGTTATCATCTCCCAGTATATCCATACTAACGATAAGGTTATTTACTCAACCAGGCATAAACATTTTGTTTTGTGCGTGGGAACAGCCTTAAgg                                                                                      >  2:633672/1‑150 (MQ=255)
                                                           cAGTTACGGCGGGAAACCTCATCCTATGGGAGAACAATGAATAATGAAATTGCGGGGTTATCATCTCCCAGTATATCCATACTAACGATAAGGTTATTTACTCAACCAGGCATAAACATTTTGTTTTGTGCGTGGGAACAGCCTTAAGgt                                                                                     <  1:3783937/150‑1 (MQ=255)
                                                           cAGTTACGGCGGGAAACCTCATCCTATGGGAGAACAATGAATAATGAAATTGCGGGGTTATCATCTCCCAGTATATCCATACTAACGATAAGGTTATTTACTCAACCAGGCATAAACATTTTGTTTTGTGCGTGGGAACAGCCTTAAGgt                                                                                     <  1:3783788/150‑1 (MQ=255)
                                                            aGTTACGGCGGGAAACCTCATCCTATGGGAGAACAATGAATAATGAAATTGCGGGGTTATCATCTCCCAGTATATCCATACTAACGATAAGGTTATTTACTCAACCAGGCATAAACATTTTGTTTTGTGCGTGGGAACAGCCTTAAGgtg                                                                                    >  2:3374808/1‑150 (MQ=255)
                                                               tACGGCGGGAAACCTCATCCTATGGGAGAACAATGAATAATGAAATTGCGGGGTTATCATCTCCCAGTATATCCATACTAACGATAAGGTTATTTACTCAACCAGGCATAAACATTTTGTTTTGTGCGTGGGAACAGCCTTAAGGTGTaa                                                                                 <  2:2723826/150‑1 (MQ=255)
                                                                    cggGAAACCTCATCCTATGGGAGAACAATGAATAATGAAATTGCGGGGTTATCATCTCCCAGTATATCCATACTAACGATAAGGTTATTTACTCAACCAGGCATAAACATTTTGTTTTGTGCGTGGGAACAGCCTTAAGGTGTAAAgggg                                                                            >  1:1183805/1‑150 (MQ=255)
                                                                     gggAAACCTCATCCTATGGGAGAACAATGAATAATGAAATTGCGGGGTTATCATCTCCCAGTATATCCATACTAACGATAAGGTTATTTACTCAACCAGGCCTAAACATTTTGTTTTGTGCGTGGGAACAGCCTTAAGGTGTAAAggggg                                                                           <  1:1075842/150‑1 (MQ=255)
                                                                     gggAAACCTCATCCTATGGGAGAACAATGAATAATGAAATTGCGGGGTTATCATCTCCCAGTATATCCATACTAACGATAAGGTTATTTACTCAACCAGGCATAAACATTTTGTTTTGTGCGTGGGAACAGCCTTAAGGTGTAAAggggg                                                                           <  1:305020/150‑1 (MQ=255)
                                                                      ggAAACCTCATCCTATGGGAGAACAATGAATAATGAAATTGCGGGGTTATCATCTCCCAGTATATCCATACTAACGATAAGGTTATTTACTCAACCAGGCATAAACATTTTGTTTTGTGCGTGGGAACAGCCTTAAGGTGTAAAGGGGGa                                                                          >  2:4140494/1‑150 (MQ=255)
                                                                      ggAAACCTCATCCTATGGGAGAACAATGAATAATGAAATTGCGGGGTTATCATCTCCCAGTATATCCATACTAACAATAAGGTTATTTACTCAACCAGGCATAAACATTTTGTTTTGTGCGTGGGAACAGCCTTAAGGTGTAAAGGGGGa                                                                          >  2:1185935/1‑150 (MQ=255)
                                                                        aaaCCTCATCCTATGGGAGAACAATGAATAATGAAATTGCGGGGTTATCATCTCCCAGTATATCCATACTAACGATAAGGTTATTTACTCAACCAGGCATAAACATTTTGTTTTGTGCGTGGGAACAGCCTTAAGGTGTAAAGGGGGAgg                                                                        <  2:329515/150‑1 (MQ=255)
                                                                         aaCCTCATCCTATGGGAGAACAATGAATAATGAAATTGCGGGGTTATCATCTCCCAGTATATCCATACTAACGATAAGGTTATTTACTCAACCAGGCATAAACATTTTGTTTTGTGCGTGGGAACAGCCTTAAGGTGTAAAGGGGGAGGt                                                                       >  1:859239/1‑150 (MQ=255)
                                                                            cTCATCCTATGGGAGAACAATGAATAATGAAATTGCGGGGTTATCATCTCCCAGTATATCCATACTAACGATAAGGTTATTTACTCAACCAGGCATAAACATTTTGTTTTGTGCGTGGGAACAGCCTTAAGGTGTAAAGGGGGAGGTGGa                                                                    >  1:1937085/1‑150 (MQ=255)
                                                                            cTCATCCTATGGGAGAACAATGAATAATGAAATTGCGGGGTTATCATCTCCCAGTATATCCATACTAACGATAAGGTTATTTACTCAACCAGGCATAAACATTTTGTTTTGTGCGTGGGAACAGCCTTAAGGTGTAAAGGGGGAGGTGGa                                                                    >  1:1937099/1‑150 (MQ=255)
                                                                             tCATCCTATGGGAGAACAATGAATAATGAAATTGCGGGGTTATCATCTCCCAGTATATCCATACTAACGATAAGGTTATTTACTCAACCAGGCATAAACATTTTGTTTTGTGCGTGGGAACAGCCTTAAGGTGTAAAGGGGGAGGTGGaa                                                                   <  1:1943180/150‑1 (MQ=255)
                                                                                   tATGGGAGAACAATGAATAATGAAATTGCGGGGTTATCATCTCCCAGTATATCCATACTAACGATAAGGTTATTTACTCAACCAGGCATAAACATTTTGTTTTGTGCGTGGGAACAGCCTTAAGGTGTAAAGGGGGAGGTGGAAATAGCa                                                             <  1:825139/150‑1 (MQ=255)
                                                                                   tATGGGAGAACAATGAATAATGAAATTGCGGGGTTATCATCTCCCAGTATATCCATACTAACGATAAGGTTATTTACTCAACCAGGCATAAACATTTTGTTTTGTGCGTGGGAACAGCCTTAAGGTGTAAAGGGGGAGGTGGAAATAGCa                                                             <  1:3662802/150‑1 (MQ=255)
                                                                                     tGGGAGAACAATGAATAATGAAATTGCGGGGTTATCATCTCCCAGTATATCCATACTAACGATAAGGTTATTTACTCAACCAGGCATAAACATTTTGTTTTGTGCGTGGGAACAGCCTTAAGGTGTAAAGGGGGAGGTGGAAATAGCAAt                                                           >  2:1141257/1‑150 (MQ=255)
                                                                                      gggAGAACAATGAATAATGAAATTGCGGGGTTATCATCTCCCAGTATATCCATACTAACGATAAGGTTATTTACTCAACCAGGCATAAACATTTTGTTTTGTGCGTGGGAACAGCCTTAAGGTGTAAAGGGGGAGGTGGAAATAGCAATg                                                          <  2:2332777/150‑1 (MQ=255)
                                                                                      gggAGAACAATGAATAATGAAATTGCGGGGTTATCATCTCCCAGTATATCCATACTAACGATAAGGTTATTTACTCAACCAGGCATAAACATTTTGTTTTGTGCGTGGGAACAGCCTTAAGGTGTAAAGGGGGAGGTGGAAATAGCAATg                                                          >  1:2644868/1‑150 (MQ=255)
                                                                                       ggAGAACAATGAATAATGAAATTGCGGGGTTATCATCTCCCAGTATATCCATACTAACGATAAGGTTATTTACTCAACCAGGCATAAACATTTTGTTTTGTGCGTGGGAACAGCCTTAAGGTGTAAAGGGGGAGGTGGAAATAGCAATga                                                         <  2:2330386/150‑1 (MQ=255)
                                                                                          gaACAATGAATAATGAAATTGCGGGGTTATCATCTCCCAGTATATCCATACTAACGATAAGGTTATTTACTCAACCAGGCATAAACATTTTGTTTTGTGCGTGGGAACAGCCTTAAGGTGTAAAGGGGGAGGTGGAAATAGCAATgagga                                                      >  2:2038018/1‑150 (MQ=255)
                                                                                           aaCAATGAATAATGAAATTGCGGGGTTATCATCTCCCAGTATATCCATACTAACGATAAGGTTATTTACTCAACCAGGCATAAACATTTTGTTTTGTGCGTGGGAACAGCCTTAAGGTGTAAAGGGGGAGGTGGAAATAGCAATgaggag                                                     >  2:2772271/1‑150 (MQ=255)
                                                                                           aaCAATGAATAATGAAATTGCGGGGTTATCATCTCCCAGTATATCCATACTAACGATAAGGTTATTTACTCAACCAGGCATAAACATTTTGTTTTGTGCGTGGGAACAGCCTTAAGGTGTAAAGGGGGAGGTGGAAATAGCAATgaggag                                                     <  2:2214505/150‑1 (MQ=255)
                                                                                           aaCAATGAATAATGAAATTGCGGGGTTATCATCTCCCAGTATATCCATACTAACAATAAGGTTATTTACTCAACCAGGCATAAACATTTTGTTTTGTGCGTGGGAACAGCCTTAAGGTGTAAAGGGGGAGGTGGAAATAGCAATgaggag                                                     >  1:2204328/1‑150 (MQ=255)
                                                                                            aCAATGAATAATGAAATTGCGGGGTTATCATCTCCCAGTATATCCATACTAACGATAAGGTTATTTACTCAACCAGGCATAAACATTTTGTTTTGTGCGTGGGAACAGCCTTAAGGTGTAAAGGGGGAGGTGGAAATAGCAATGAGGAGt                                                    <  2:1585369/150‑1 (MQ=255)
                                                                                            aCAATGAATAATGAAATTGCGGGGTTATCATCTCCCAGTATATCCATACTAACGATAAGGTTATTTACTCAACCAGGCATAAACATTTTGTTTTGTGCGTGGGAACAGCCTTAAGGTGTAAAGGGGGAGGTGGAAATAGCAATGAGGAGt                                                    <  2:762842/150‑1 (MQ=255)
                                                                                            aCAATGAATAATGAAATTGCGGGGTTATCATCTCCCAGTATATCCATACTAACGATAAGGTTATTTACTCAACCAGGCATAAACATTTTGTTTTGTGCGTGGGAACAGCCTTAAGGTGTAAAGGGGGAGGTGGAAATAGCAATGAGGAGt                                                    >  1:733839/1‑150 (MQ=255)
                                                                                              aaTGAATAATGAAATTGCGGGGTTATCATCTCCCAGTATATCCATACTAACGATAAGGTTATTTACTCAACCAGGCATAAACATTTTGTTTTGTGCGTGGGAACAGCCTTAAGGTGTAAAGGGGGAGGTGGAAATAGCAATGAGGAGTAt                                                  >  1:4059024/1‑150 (MQ=255)
                                                                                              aaTGAATAATGAAATTGCGGGGTTATCATCTCCCAGTATATCCATACTAACGATAAGGTTATTTACTCAACCAGGCATAAACATTTTGTTTTGTGCGTGGGAACAGCCTTAAGGTGTAAAGGGGGAGGTGGAAATAGCAATGAGGAGTAt                                                  >  1:661663/1‑150 (MQ=255)
                                                                                              aaTGAATAATGAAAGTGCGGGGTTATCATCTCCCAGTATATCCATACTAACGATAAGGTTATTTACTCAACCAGGCATAAACATTTTGTTTTGTGCGTGGGAACAGCCTTAAGGTGTAAAGGGGGAGGTGGAAATAGCAATGAGGAGTAt                                                  <  2:2838145/150‑1 (MQ=255)
                                                                                               aTGAATAATGAAATTGCGGGGTTATCATCTCCCAGTATATCCATACTAACGATAAGGTTATTTACTCAACCAGGCATAAACATTTTGTTTTGTGCGTGGGAACAGCCTTAAGGTGTAAAGGGGGAGGTGGAAATAGCAATGAGGAGTATc                                                 <  1:2082651/150‑1 (MQ=255)
                                                                                                tGAATAATGAAATTGCGGGGTTATCATCTCCCAGTATATCCATACTAACGATAAGGTTATTTACTCAACCAGGCATAAACATTTTGTTTTGTGCGTGGGAACAGCCTTAAGGTGTAAAGGGGGAGGTGGGAATAGCAATGAGGAGTATCa                                                >  1:911895/1‑150 (MQ=255)
                                                                                                tGAATAATGAAATTGCGGGGTTATCATCTCCCAGTATATCCATACTAACGATAAGGTTATTTACTCAACCAGGCATAAACATTTTGTTTTGTGCGTGGGAACAGCCTTAAGGTGTAAAGGGGGAGGTGGGAATAGCAATGAGGAGTATCa                                                >  1:911873/1‑150 (MQ=255)
                                                                                                  aataatGAAATTGCGGGGTTATCATCTCCCAGTATATCCATACTAACGATACGGTTATTTACTCAACCAGTCATAAACATTTTGTTTTGTGCGTGGGAACAGCCTTAAGGTGTAAAGGGGGAGGTGGAAATATCAATGAGGAGTATCAGc                                              >  2:3819091/1‑150 (MQ=255)
                                                                                                   ataatGAAATTGCGGGGTTATCATCTCCCAGTATATCCATACTAACGATAAGGTTATTTACTCAACCAGGCATAAACATTTTGTTTTGTGCGTGGGAACAGCCTTAAGGTGTAAAGGGGGAGGTGGAAATAGCAATGAGGAGTATCAGCa                                             >  2:3561071/1‑150 (MQ=255)
                                                                                                   ataatGAAATTGCGGGGTTATCATCTCCCAGTATATCCATACTAACGATAAGGTTATTTACTCAACCAGGCATAAACATTTTGTTTTGTGCGTGGGAACAGCCTTAAGGTGTAAAGGGGGAGGTGGAAATAGCAATGAGGAGTATCAGCa                                             >  1:1448490/1‑150 (MQ=255)
                                                                                                    taatGAAATTGCGGGGTTATCATCTCCCAGTATATCCATACTAACGATAAGGTTATTTACTCAACCAGGCATAAACATTTTGTTTTGTGCGTGGGAACAGCCTTAAGGTGTAAAGGGGGAGGTGGAAATAGCAATGAGGAGTATCAGCaa                                            <  2:342431/150‑1 (MQ=255)
                                                                                                     aatGAAATTGCGGGGTTATCATCTCCCAGTATATCCATACTAACGATAAGGTTATTTACTCAACCAGGCATAAACATTTTGTTTTGTGCGTGGGAACAGCCTTAAGGTGTAAAGGGGGAGGTGGAAATAGCAATGAGGAGTATCAGCAAg                                           <  1:1372765/150‑1 (MQ=255)
                                                                                                       tGAAATTGCGGGGTTATCATCTCCCAGTATATCCATACTAACGATAAGGTTATTTACTCAACCAGGCATAAACATTTTGTTTTGTGCGTGGGAACAGCCTTAAGGTGTAAAGGGGGAGGTGGAAATAGCAATGAGGAGTATCAGCAAGaa                                         >  1:1365733/1‑150 (MQ=255)
                                                                                                        gAAATTGCGGGGTTATCATCTCCCAGTATATCCATACTAACGATAAGGTTATTTACTCAACCAGGCATAAACATTTTGTTTTGTGCGTGGGAACAGCCTTAAGGTGTAAAGGGGGAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAAt                                        <  2:1937099/150‑1 (MQ=255)
                                                                                                        gAAATTGCGGGGTTATCATCTCCCAGTATATCCATACTAACGATAAGGTTATTTACTCAACCAGGCATAAACATTTTGTTTTGTGCGTGGGAACAGCCTTAAGGTGTAAAGGGGGAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAAt                                        <  2:1936203/150‑1 (MQ=255)
                                                                                                            ttGCGGGGTTATCATCTCCCAGTATATCCATACTAACGATAAGGTTATTTACTCAACCAGGCATAAACATTTTGTTTTGTGCGTGGGAACAGCCTTAAGGTGTAAAGGGGGAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTc                                    >  2:3843901/1‑150 (MQ=255)
                                                                                                                ggggTTATCATCTCCCAGTATATCCATACTAACGATAAGGTTATTTACTCAACCAGGCATAAACATTTTGTTTTGTGCGTGGGAACAGCCTTAAGGTGTAAAGGGGGAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTcgccg                                >  2:1627481/1‑150 (MQ=255)
                                                                                                                  ggTTATCATCTCCCAGTATATCCATACTAACGATAAGGTTATTTACTCAACCAGGCATAAACATTTTGTTTTGTGCGTGGGAACAGCCTTAAGGTGTAAAGGGGGAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCt                              >  2:4509820/1‑150 (MQ=255)
                                                                                                                        catcTCCCAGTATATCCATACTAACGATAAGGTTATTTACTCAACCAGGCATAAACATTTTGTTTTGTGCGTGGGAACAGCCTTAAGGTGTAAAGGGGGAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACca                        >  2:3582122/1‑150 (MQ=255)
                                                                                                                         atcTCCCAGTATATCCATACTAACGATAAGGTTATTTACTCAACCAGGCATAAACATTTTGTTTTGTGCGTGGGAACAGCCTTAAGGTGTAAAGGGGGAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACcac                       >  2:1828235/1‑150 (MQ=255)
                                                                                                                           ctcCCAGTATATCCATACTAACGATAAGGTTATTTACTCAACCAGGCATAAACATTTTGTTTTGTGCGTGGGAACAGCCTTAAGGTGTAAAGGGGGAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACaa                     >  2:2217878/1‑150 (MQ=255)
                                                                                                                           ctcCCAGTATATCCATACTAACGATAAGGTTATTTACTCAACCAGGCATAAACATTTTGTTTTGTGCGTGGGAACAGCCTTAAGGTGTAAAGGGGGAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACaa                     >  2:4450739/1‑150 (MQ=255)
                                                                                                                             cccAGTATATCCATACTAACGATAAGGTTATTTACTCAACCAGGCATAAACATTTTGTTTTGTGCGTGGGAACAGCCTTAAGGTGTAAAGGGGGAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACg                   >  1:2301420/1‑150 (MQ=255)
                                                                                                                              ccAGTATATCCATACTAACAATAAGGTTATTTACTCAACCAGGCATAAACATTTTGTTTTGTGCGTGGGAACAGCCTTAAGGTGTAAAGGGGGAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGt                  <  1:4286030/150‑1 (MQ=255)
                                                                                                                                   atatCCATACTAACGATAAGGTTATTTACTCAACCAGGCATAAACATTTTGTTTTGTGCGTGGGAACAGCCTTAAGGTGTAAAGGGGGAGGTGTAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATg             >  1:582391/1‑150 (MQ=255)
                                                                                                                                    tatCCATACTAACGATAAGGTTATTTACTCAACCAGGCATAAACATTTTGTTTTGTGCGTGGGAACAGCCTTAAGGTGTAAAGGGGGAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATga            <  1:2875897/150‑1 (MQ=255)
                                                                                                                                    tatCCATACTAACGATAAGGTTATTTACTCAACCAGGCATAAACATTTTGTTTTGTGCGTGGGAACAGCCTTAAGGTGTAAAGGGGGAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATga            >  2:2084880/1‑150 (MQ=255)
                                                                                                                                     atCCATACTAACGATAAGGTTATTTACTCAACCAGGCATAAACATTTTGTTTTGTGCGTGGGAACAGCCTTAAGGTGTAAAGGGGGAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATgag           <  2:842072/150‑1 (MQ=255)
                                                                                                                                       ccATACTAACAATAAGGTTATTTACTCAACCAGGCATAAACATTTTGTTTTGTGCGTGGGAACAGCCTTAAGGTGTAAAGGGGGAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGagag         <  2:1077868/150‑1 (MQ=255)
                                                                                                                                          tACTAACAATAAGGTTATTTACTCAACCAGGCATAAACATTTTGTTTTGTGCGTGGGAACAGCCTTAAGGTGTAAAGGGGGAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGa      <  1:1837305/150‑1 (MQ=255)
                                                                                                                                           aCTAACAATAAGGTTATTTACTCAACCAGGCATAAACATTTTGTTTTGTGCGTGGGAACAGCCTTAAGGTGTAAAGGGGGAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGAt     >  1:3625062/1‑150 (MQ=255)
                                                                                                                                            cTAACGATAATGTTATTTACTCAACCAGGCATAAACATTTTGTTTTGTGCGTGGGAACAGCCTTAAGGTGTAAAGGGGGAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATg    <  2:2008469/150‑1 (MQ=255)
                                                                                                                                            cTAACGATAAGGTTATTTACTCAACCAGGCATAAACATTTTGTTTTGTGCGTGGGAACAGCCTTAAGGTGTAAAGGGGGAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATg    <  2:4053532/150‑1 (MQ=255)
                                                                                                                                            cTAACGATAAGGTTATTTACTCAACCAGGCATAAACATTTTGTTTTGTGCGTGGGAACAGCCTTAAGGTGTAAAGGGGGAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATg    <  2:2717048/150‑1 (MQ=255)
                                                                                                                                             tAACGATAAGGTTATTTACTCAACCAGGCATAAACATTTTGTTTTGTGCGTGGGAACAGCCTTAAGGTGTAAAGGGGTAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGa   >  2:1818213/1‑150 (MQ=255)
                                                                                                                                             tAACGATAAGGTTATTTACTCAACCAGGCATAAACATTTTGTTTTGTGCGTGGGAACAGCCTTAAGGTGTAAAGGGGGAGGTGGAAATAGCTATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGa   <  1:3884316/150‑1 (MQ=255)
                                                                                                                                              aaCGATAAGGTTATTTACTCAACCAGGCATAAACATTTTGTTTTGTGCGTGGGAACAGCCTTAAGGTGTAAAGGGGGAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGaa  <  1:2038018/150‑1 (MQ=255)
                                                                                                                                              aaCGATAAGGTTATTTACTCAACCAGGCATAAACATTTTGTTTTGTGCGTGGGAACAGCCTTAAGGTGTAAAGGGGGAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGaa  <  1:4366764/150‑1 (MQ=255)
                                                                                                                                              aaCGATAAGGTTATTTACTCAACCAGGCATAAACATTTTGTTTTGTGCGTGGGAACAGCCTTAAGGTGTAAAGGGGGAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGaa  <  1:4372089/150‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
GGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAACCTCATCCTATGGGAGAACAATGAATAATGAAATTGCGGGGTTATCATCTCCCAGTATATCCATACTAACAATAAGGTTATTTACTCAACCAGGCATAAACATTTTGTTTTGTGCGTGGGAACAGCCTTAAGGTGTAAAGGGGGAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAA  >  CP009273/1219388‑1219679

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: