Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,991,843 G→A 13.1% L40F (CTC→TTC)  yecS ← ABC family putative inner membrane permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,991,8430GA13.1% 761.9 / 80.4 251L40F (CTC→TTC) yecSABC family putative inner membrane permease
Reads supporting (aligned to +/- strand):  ref base G (109/109);  new base A (12/21);  total (121/130)
Fisher's exact test for biased strand distribution p-value = 1.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.58e-01

Only 100 of 252 total aligned reads displayed.

GAATGGGATCTAATTCGATACCAAACTGCGGCAGGCCGTAATAGATCATAAACAGTTGTGCGATGAGCGGCGTACCACGAAAAATAGAGATATAAAAACGCGCCAGCCAGCGCACGGGCCAGATTGGCGACAGACGCATCAGCGCGAGAATAAACCCGAGCAGTAAGCCAAAAAACATGCCGCCGATACTGAGTTGCAGCGTATACCCAGCCCCTTTCAGCAGGAACGGCAAGGAATCAATAACCAGTTGTATACTTTCTTGCATGAGCGTTTCCGGAGCTAAACGTGGGGATGAT  >  CP009273/1991696‑1991991
                                                                                                                                                   |                                                                                                                                                    
gAATGGGATCTAATTCGATACCAAACTGCGGCAGGCCGTAATAGATCATAAACAGTTGTGCGATGAGCGGCGTACCACGAAAAATAGAGATATAAAAACGCGCCAGCCAGCGCACGGGCCAGATTGGCGACAGACGCATCAGCGCGAGaa                                                                                                                                                    <  1:732345/150‑1 (MQ=255)
   tGGGATCTAATTCGATACCAAACTGCGGCAGGCCGTAATAGATCATAAACAGTTGTGCGATGAGCGGCGTACCACGAAAAATAGAGATATAAAAACGCGCCAGCCAGCGCACGGGCCAGATTGGCGACAGACGCATCAGCGCGAGAATaa                                                                                                                                                 >  1:1682951/1‑150 (MQ=255)
    gggATCTAATTCGATACCAAACTGCGGCAGGCCGTAATAGATCATAAACAGTTGTGCGATGAGCGGCGTACCACGAAAAATAGAGATATAAAAACGCGCCAGCCAGCGCACGGGCCAGATTGGCGACAGACGCATCAGCGCGAGAATaaa                                                                                                                                                <  1:786779/150‑1 (MQ=255)
    gggATCTAATTCGATACCAAACTGCGGCAGGCCGTAATAGATCATAAACAGTTGTGCGATGAGCGGCGTACCACGAAAAATAGAGATATAAAAACGCGCCAGCCAGCGCACGGGCCAGATTGGCGACAGACGCATCAGCGCGAGAATaaa                                                                                                                                                >  1:3390672/1‑150 (MQ=255)
     ggATCTAATTCGATACCAAACTGCGGCAGGCCGTAATAGATCATAAACAGTTGTGCGATGAGCGGCGTACCACGAAAAATAGAGATATAAAAACGCGCCAGCCAGCGCACGGGCCAGATTGGCGACAGACGCATCAGCGCGAGAATAAAc                                                                                                                                               <  1:1244312/150‑1 (MQ=255)
     ggATCTAATTCGATACCAAACTGCGGCAGGCCGTAATAGATCATAAACAGTTGTGCGATGAGCGGCGTACCACGAAAAATAGAGATATAAAAACGCGCCAGCCAGAGCACGGGCCAGATTGGCGACAGACGCATCAGCGCGAAAATAAAc                                                                                                                                               >  2:641951/1‑150 (MQ=255)
      gATCTAATTCGATACCAAACTGCGGCAGGCCGTAATAGATCATAAACAGTTGTGCGATGAGCGGCGTACCACGAAAAATAGAGATATAAAAACGCGCCAGCCAGCGCACGGGCCAGATTGGCGACAGACGCATCAGCGCGAGAATAAAcc                                                                                                                                              <  1:2332428/150‑1 (MQ=255)
          tAATTCGATACCAAACTGCGGCAGGCCGTAATAGATCATAAACAGTTGTGCGATGAGCGGCGTACCACGAAAAATAGAGATATAAAAACGCGCCAGCCAGCGCACGGGCCAGATTGGCGACAGACGCATCAGCGCGAGAATAAACCCGAg                                                                                                                                          <  1:525614/150‑1 (MQ=255)
          tAATTCGATACCAAACTGCGGCAGGCCGTAATAGATCATAAACAGTTGTGCGATGAGCGGCGTACCACGAAAAATAGAGATATAAAAACGCGCCAGCCAGCGCACGGGCCAGATTGGCGACAGACGCATCAGCGCGAGAATAAACCCGAg                                                                                                                                          <  1:161326/150‑1 (MQ=255)
          tAATTCGATACCAAACTGCGGCAGGCCGTAATAGATCATAAACAGTTGTGCGATGAGCGGCGTACCACGAAAAATAGAGATATAAAAACGCGCCAGCCAGCGCACGGGCCAGATTGGCGACAGACGCATCAGCGCGAGAATAAACCCGAg                                                                                                                                          <  1:2100544/150‑1 (MQ=255)
           aaTTCGATACCAAACTGCGGCAGGCCGTAATAGATCATAAACAGTTGTGCGATGAGCGGCGTACCACGAAAAATAGAGATATAAAAACGCGCCAGCCAGCGCACGGGCCAGATTGGCGACAGACGCATCAGCGCGAGAATAAACCCGAGc                                                                                                                                         >  2:2201812/1‑150 (MQ=255)
            aTTCGATACCAAACTGCGGCAGGCCGTAATAGATCATAAACAGTTGTGCGATGAGCGGCGTACCACGAAAAATAGAGATATAAAAACGCGCCAGCCAGCGCACGGGCCAGATTGGCGACAGACGCATCAGCGCGAGACTAAACCCGAGCa                                                                                                                                        >  2:2363609/1‑150 (MQ=255)
                gATACCAAACTGCGGCAGGCCGTAATAGATCATAAACAGTTGTGCGATGAGCGGCGTACCACGAAAAATAGAGATATAAAAACGCGCCAGCCAGCGCACGGGCCAGATTGGCGACAGACGCATCAGCGCGAGAATAAACCCGAGCAGTaa                                                                                                                                    <  2:3422927/150‑1 (MQ=255)
                gATACCAAACTGCGGCAGGCCGTAATAGATCATAAACAGTTGTGCGATGAGCGGCGTACCACGAAAAATAGAGATATAAAAACGCGCCAGCCAGCGCACGGGCCAGATTGGCGACAGACGCATCAGCGCGAGAATAAACCCGAGCAGTaa                                                                                                                                    >  1:2447852/1‑150 (MQ=255)
                gATACCAAACTGCGGCAGGCCGTAATAGATCATAAACAGTTGTGCGATGAGCGGCGTACCACGAAAAATAGAGATATAAAAACGCGCCAGCCAGCGCACGGGCCAGATTGGCGACAGACGCATCAGCGCGAAAATAAACCCGAGCAGTaa                                                                                                                                    <  1:995762/150‑1 (MQ=255)
                 aTACCAAACTGCGGCAGGCCGTAATAGATCATAAACAGTTGTGCGATGAGCGGCGTACCACGAAAAATAGAGATATAAAAACGCGCCAGCCAGCGCACGGGCCAGATTGGCGACAGACGCATCAGCGCGAGAATAAACCCGAGCAGTAAg                                                                                                                                   <  2:1237112/150‑1 (MQ=255)
                  tACCAAACTGCGGCAGGCCGTAATAGATCATAAACAGTTGTGCGATGAGCGGCGTACCACGAAAAATAGAGATATAAAAACGCGCCAGCCAGCGCACGGGCCAGATTGGCGACAGACGCATCAGCGCGAGAATAAACCCGAGCAGTAAGc                                                                                                                                  <  1:587116/150‑1 (MQ=255)
                  tACCAAACTGCGGCAGGCCGTAATAGATCATAAACAGTTGTGCGATGAGCGGCGTACCACGAAAAATAGAGATATAAAAACGCGCCAGCCAGCGCACGGGCCAGATTGGCGACAGACGCATCAGCGCGAAAATAAACCCGAGCAGTAAGc                                                                                                                                  <  2:2803102/150‑1 (MQ=255)
                    ccAAACTGCGGCAGGCCGTAATAGATCATAAACAGTTGTGCGATGAGCGGCGTACCACGAAAAATAGAGATATAAAAACGCGCCAGCCAGCGCACGGGCCAGATTGGCGACAGACGCATCAGCGCGAGAATAAACCCGAGCAGTAAGCCa                                                                                                                                <  2:1384047/150‑1 (MQ=255)
                    ccAAACTGCGGCAGGCCGTAATAGATCATAAACAGTTGTGCGATGAGCGGCGTACCACGAAAAATAGAGATATAAAAACGCGCCAGCCAGCGCACGGGCCAGATTGGCGACAGACGCATCAGCGCGAGAATAAACCCGAGCAGTAAGCCa                                                                                                                                <  1:1934456/150‑1 (MQ=255)
                       aaCTGCGGCAGGCCGTAATAGATCATAAACAGTTGTGCGATGAGCGGCGTACCACGAAAAATAGAGATATAAAAACGCGCCAGCCAGCGCACGGGCCAGATTGGCGACAGACGCATCAGCGCGAGAATAAACCCGAGCAGTAAGCCaaaa                                                                                                                             >  1:4035723/1‑150 (MQ=255)
                       aaCTGCGGCAGGCCGTAATAGATCATAAACAGTTGTGCGATGAGCGGCGTACCACGAAAAATAGAGATATAAAAACGCGCCAGCCAGCGCACGGGCCAGATTGGCGACAGACGCATCAGCGCGAGAATAAACCCGAGCAGTAAGCCaaaa                                                                                                                             <  1:2332820/150‑1 (MQ=255)
                       aaCTGCGGCAGGCCGTAATAGATCATAAACAGTTGTGCGATGAGCGGCGTACCACGAAAAATAGAGATATAAAAACGCGCCAGCCAGCGCACGGGCCAGATTGGCGACAGACGCATCAGCGCGAGAATAAACCCGAGCAGTAAGCCaaaa                                                                                                                             >  2:1812904/1‑150 (MQ=255)
                        aCTGCGGCAGGCCGTAATAGATCATAAACAGTTGTGCGATGAGCGGCGTACCACGAAAAATAGAGATATAAAAACGCGCCAGCCAGCGCACGGGCCAGATTGGCGACAGACGCATCAGCGCGAGAATAAACCCGAGCAGTAAGCCaaaaa                                                                                                                            <  1:2714887/150‑1 (MQ=255)
                          tGCGGCAGGCCGTAATAGATCATAAACAGTTGTGCGATGAGCGGCGTACCACGAAAAATAGAGATATAAAAACGCGCCAGCCAGCGCACGGGCCAGATTGGCGACAGACGCATCAGCGCGAGAATAAACCCGAGCAGTAAGCCAAAAAAc                                                                                                                          <  1:1491463/150‑1 (MQ=255)
                          tGCGGCAGGCCGTAATAGATCATAAACAGTTGTGCGATGAGCGGCGTACCACGAAAAATAGAGATATAAAAACGCGCCAGCCAGCGCACGGGCCAGATTGGCGACAGACGCATCAGCGCGAAAATAAACCCGAGCAGTAAGCCAAAAAAc                                                                                                                          >  2:1655407/1‑150 (MQ=255)
                                aGGCCGTAATAGATCATAAACAGTTGTGCGATGAGCGGCGTACCACGAAAAATAGAGATATAAAAACGCGCCAGCCAGCGCACGGGCCAGATTGGCGACAGACGCATCAGCGCGAGAATAAACCCGAGCAGTAAGCCAAAAAACATgccg                                                                                                                    >  2:401875/1‑150 (MQ=255)
                                  gCCGTAATAGATCATAAACAGTTGTGCGATGAGCGGCGTACCACGAAAAATAGAGATATAAAAACGCGCCAGCCAGCGCACGGGCCAGATTGGCGACAGACGCATCAGCGCGAGAATAAACCCGAGCAGTAAGCCAAAAAACATgccgcc                                                                                                                  <  1:2142506/150‑1 (MQ=255)
                                  gCCGTAATAGATCATAAACAGTTGTGCGATGAGCGGCGTACCACGAAAAATAGAGATATAAAAACGCGCCAGCCAGCGCACGGGCCAGATTGGCGACAGACGCATCAGCGCGAGAATAAACCCGAGCAGTAAGCCAAAAAACATgccgcc                                                                                                                  >  1:4367756/1‑150 (MQ=255)
                                  gCCGTAATAGATCATAAACAGTTGTGCGATGAGCGGCGTACCACGAAAAATAGAGATATAAAAACGCGCCAGCCAGCGCACGGGCCAGATTGGCGACAGACGCATCAGCGCGAGAATAAACCCGAGCAGTAAGCCAAAAAACATgccgcc                                                                                                                  <  1:3133713/150‑1 (MQ=255)
                                     gTAATAGATCATAAACAGTTGTGCGATGAGCGGCGTACCACGAAAAATAGAGATATAAAAACGCGCCAGCCAGCGCACGGGCCAGATTGGCGACAGACGCATCAGCGCGAGAATAAACCCGAGCAGTAAGCCAAAAAACATGCCGCCGAt                                                                                                               >  1:1601412/1‑150 (MQ=255)
                                     gTAATAGATCATAAACAGTTGTGCGATGAGCGGCGTACCACGAAAAATAGAGATATAAAAACGCGCCAGCCAGCGCACGGGCCAGATTGGCGACAGACGCATCAGCGCGAGAATAAACCCGAGCAGTAAGCCAAAAAACATGCCGCCGAt                                                                                                               >  1:3436518/1‑150 (MQ=255)
                                     gTAATAGATCATAAACAGTTGTGCGATGAGCGGCGTACCACGAAAAATAGAGATATAAAAACGCGCCAGCCAGCGCACGGGCCAGATTGGCGACAGACGCATCAGCGCGAGAATAAACCCGAGCAGTAAGCCAAAAAACATGCCGCCGAt                                                                                                               >  1:194513/1‑150 (MQ=255)
                                     gTAATAGATCATAAACAGTTGTGCGATGAGCGGCGTACCACGAAAAATAGAGATATAAAAACGCGCCAGCCAGCGCACGGGCCAGATTGGCGACAGACGCATCAGCGCGAAAATAAACCCGAGCAGTAAGCCAAAAAACATGCCGCCGAt                                                                                                               <  1:3397188/150‑1 (MQ=255)
                                     gTAATAGATCATAAACAGTTGTGCGATGAGCGGCGTACCACGAAAAATAGAGATATAAAAACGCGCCAGCCAGCGCACGGGCCAGATTGGCGACAGACGCATCAGCGCGAAAATAAACCCGAGCAGTAAGCCAAAAAACATGCCGCCGAt                                                                                                               >  1:4599321/1‑150 (MQ=255)
                                        aTAGATCATAAACAGTTGTGCGATGAGCGGCGTACCACGAAAAATAGAGATATAAAAACGCGCCAGCCAGCGCACGGGCCAGATTGGCGACAGACGCATCAGCGCGAGAATAAACCCGAGCAGTAAGCCAAAAAACATGCCGCCGATACt                                                                                                            >  1:737025/1‑150 (MQ=255)
                                           gATCATAAACAGTTGTGCGATGAGCGGCGTACCACGAAAAATAGAGATATAAAAACGCGCCAGCCAGCGCACGGGCCAGATTGGCGACAGACGCATCAGCGCGAGAATAAACCCGAGCAGTAAGCCAAAAAACATGCCGCCGATACTGAg                                                                                                         >  1:37105/1‑150 (MQ=255)
                                           gATCATAAACAGTTGTGCGATGAGCGGCGTACCACGAAAAATAGAGATATAAAAACGCGCCAGCCAGCGCACGGGCCAGATTGGCGACAGACGCATCAGCGCGAAAATAAACCCGAGCAGTAAGCCAAAAAACATGCCGCCGATACTGAg                                                                                                         <  1:3900278/150‑1 (MQ=255)
                                            aTCATAAACAGTTGTGCGATGAGCGGCGTACCACGAAAAATAGAGATATAAAAACGCGCCAGCCAGCGCACGGGCCAGATTGGCGACAGACGCATCAGCGCGAGAATAAACCCGAGCAGTAAGCCAAAAAACATGCCGCCGATACTGAGt                                                                                                        >  2:1856056/1‑150 (MQ=255)
                                              cATAAACAGTTGTGCGATGAGCGGCGTACCACGAAAAATAGAGATATAAAAACGCGCCAGCCAGCGCACGGGCCAGATTGGCGACAGACGCATCAGCGCGAGAATAAACCCGAGCAGTAAGCCAAAAAACATGCCGCCGATACTGAGTTg                                                                                                      <  1:2683770/150‑1 (MQ=255)
                                              cATAAACAGTTGTGCGATGAGCGGCGTACCACGAAAAATAGAGATATAAAAACGCGCCAGCCAGCGCACGGGCCAGATTGGCGACAGACGCATCAGCGCGAGAATAAACCCGAGCAGTAAGCCAAAAAACATGCCGCCGATACTGAGTTg                                                                                                      >  2:2675156/1‑150 (MQ=255)
                                              cATAAACAGTTGTGCGATGAGCGGCGTACCACGAAAAATAGAGATATAAAAACGCGCCAGCCAGCGCACGGGCCAGATTGGCGACAGACGCATCAGCGCGAGAATAAACCCGAGCAGTAAGCCAAAAAACATGCCGCCGATACTGAGTTg                                                                                                      <  2:4483165/150‑1 (MQ=255)
                                                tAAACAGTTGTGCGATGAGCGGCGTACCACGAAAAATAGAGATATAAAAACGCGCCAGCCAGCGCACGGGCCAGATTGGCGACAGACGCATCAGCGCGAGAATAAACCCGAGCAGTAAGCCAAAAAACATGCCGCCGATACTGAGTTGCa                                                                                                    >  2:745955/1‑150 (MQ=255)
                                                 aaaCAGTTGTGCGATGAGCGGCGTACCACGAAAAATAGAGATATAAAAACGCGCCAGCCAGCGCACGGGCCAGATTGGCGACAGACGCATCAGCGCGAGAATAAACCCGAGCAGTAAGCCAAACAACATGCCGCCGATACTGAGTTGCAg                                                                                                   >  2:2365952/1‑150 (MQ=255)
                                                   aCAGTTGTGCGATGAGCGGCGTACCACGAAAAATAGAGATATAAAAACGCGCCAGCCAGCGCACGGGCCAGATTGGCGACAGACGCATCAGCGCGAGAATAAACCCGAGCAGTAAGCCAAAAAACATGCCGCCGATACTGAGTTGCAGCg                                                                                                 <  1:4421941/150‑1 (MQ=255)
                                                     aGTTGTGCGATGAGCGGCGTACCACGAAAAATAGAGATATAAAAACGCGCCAGCCAGCGCACGGGCCAGATTGGCGACAGACGCATCAGCGCGAGAATAAACCCGAGCAGTAAGCCAAAAAACATGCCGCCGATACTGAGTTGCAGCGta                                                                                               >  2:1671556/1‑150 (MQ=255)
                                                            cGATGAGCGGCGTACCACGAAAAATAGAGATATAAAAACGCGCCAGCCAGCGCACGGGCCAGATTGGCGACAGACGCATCAGCGCGAGAATAAACCCGAGCAGTAAGCCAAAAAACATGCCGCCGATACTGAGTTGCAGCGTATACCCAg                                                                                        <  1:3795752/150‑1 (MQ=255)
                                                            cGATGAGCGGCGTACCACGAAAAATAGAGATATAAAAACGCGCCAGCCAGCGCACGGGCCAGATTGGCGACAGACGCATCAGCGCGAGAATAAACCCGAGCAGTAAGCCAAAAAACATGCCGCCGATACTGAGTTGCAGCGTATACCCAg                                                                                        >  2:4233406/1‑150 (MQ=255)
                                                              aTGAGCGGCGTACCACGAAAAATAGAGATATAAAAACGCGCCAGCCAGCGCACGGGCCAGATTGGCGACAGACGCATCAGCGCGAGAATAAACCCGAGCAGTAAGCCAAAAAACATGCCGCCGATACTGAGTTGCAGCGTATACCCAGcc                                                                                      <  2:4391209/150‑1 (MQ=255)
                                                               tGAGCGGCGTACCACGAAAAATAGAGATATAAAAACGCGCCAGCCAGCGCACGGGCCAGATTGGCGACAGACGCATCAGCGCGAGAATAAACCCGAGCAGTAAGCCAAAAAACATGCCGCCGATACTGAGTTGCAGCGTATACCCAGccc                                                                                     <  1:1046256/150‑1 (MQ=255)
                                                                 aGCGGCGTACCACGAAAAATAGAGATATAAAACCGCGCCAGCCAGCGCACGGGCCAGATTGGCGACAGACGCATCAGCGCGAGAATAAACCCGAGCAGTAAGCCAAAAAACATGCCGCCGATACTGAGTTGCAGCGTATACCCAGCCCCt                                                                                   <  1:4159749/150‑1 (MQ=255)
                                                                    ggcgTACCACGAAAAATAGAGATATAAAAACGCGCCAGCCAGCGCACGGGCCAGATTGGCGACAGACGCATCAGCGCGAAAATAAACCCGAGCAGTAAGCCAAAAAACATGCCGCCGATACTGAGTTGCAGCGTATACCCAGCCCCTTTc                                                                                >  1:4219662/1‑150 (MQ=255)
                                                                     gcgTACCACGAAAAATAGAGATATAAAAACGCGCCAGCCAGCGCACGGGCCAGATTGGCGACAGACGCATCAGCGCGAGAATAAACCCGAGCAGTAAGCCAAAAAACATGCCGCCGATACTGAGTTGCAGCGTATACCCAGCCCCTTTca                                                                               <  2:3878822/150‑1 (MQ=255)
                                                                         aCCACGAAAAATAGAGATATAAAAACGCGCCAGCCAGCGCACGGGCCAGATTGGCGACAGACGCATCAGCGCGAGAATAAACCCGAGCAGTAAGCCAAAAAACATGCCGCCGATACTGAGTTGCAGCGTATACCCAGCCCCTTTcagcag                                                                           <  2:497703/150‑1 (MQ=255)
                                                                              gAAAAATAGAGATATAAAAACGCGCCAGCCAGCGCACGGGCCAGATTGGCGACAGACGCATCAGCGCGAGAATAAACCCGAGCAGTAAGCCAAAAAACATGCCGCCGATACTGAGTTGCAGCGTATACCCAGCCCCTTTCAGCAGGAACg                                                                      <  2:3804924/150‑1 (MQ=255)
                                                                              gAAAAATAGAGATATAAAAACGCGCCAGCCAGCGCACGGGCCAGATTGGCGACAGACGCATCAGCGCGAGAATAAACCCGAGCAGTAAGCCAAAAAACATGCCGCCGATACTGAGTTGCAGCGTATACCCAGCCCCTTTCAGCAGGAACg                                                                      <  2:1620281/150‑1 (MQ=255)
                                                                               aaaaaTAGAGATATAAAAACGCGCCAGCCAGCGCACGGGCCAGATTGGCGACAGACGCATCAGCGCGAGAATAAACCCGAGCAGTAAGCCAAAAAACATGCCGCCGATACTGAGTTGCAGCGTATACCCAGCCCCTTTCAGCAGGAACgg                                                                     <  2:3648124/150‑1 (MQ=255)
                                                                               aaaaaTAGAGATATAAAAACGCGCCAGCCAGCGCACGGGCCAGATTGGCGACAGACGCATCAGCGCGAAAATAAACCCGAGCAGTAAGCCAAAAAACATGCCGCCGATACTGAGTTGCAGCGTATACCCAGCCCCTTTCAGCAGGAACgg                                                                     <  2:797121/150‑1 (MQ=255)
                                                                                aaaaTAGAGATATAAAAACGCGCCAGCCAGCGCACGGGCCAGATTGGCGACAGACGCATCAGCGCGAGAATAAACCCGAGCAGTAAGCCAAAAAACATGCCGCCGATACTGAGTTGCAGCGTATACCCAGCCCCTTTCAGCAGGAACGGc                                                                    >  1:3349054/1‑150 (MQ=255)
                                                                                aaaaTAGAGATATAAAAACGCGCCAGCCAGCGCACGGGCCAGATTGGCGACAGACGCATCAGCGCGAGAATAAACCCGAGCAGTAAGCCAAAAAACATGCCGCCGATACTGAGTTGCAGCGTATACCCAGCCCCTTTCAGCAGGAACGGc                                                                    >  2:98363/1‑150 (MQ=255)
                                                                                  aaTAGAGATATAAAAACGCGCCAGCCAGCGCACGGGCCAGATTGGCGACAGACGCATCAGCGCGAGAATAAACCCGAGCAGTAAGCCAAAAAACATGCCGCCGATACTGAGTTGCAGCGTATACCCAGCCCCTTTCAGCAGGAACGGCaa                                                                  >  2:4551419/1‑150 (MQ=255)
                                                                                  aaTAGAGATATAAAAACGCGCCAGCCAGCGCACGGGCCAGATTGGCGACAGACGCATCAGCGCGAGAATAAACCCGAGCAGTAAGCCAAAAAACATGCCGCCGATACTGAGTTGCAGCGTAGACCCAGCCCCTTTCAGCAGGAACGGCaa                                                                  >  2:2891279/1‑150 (MQ=255)
                                                                                    tAGAGATATAAAAACGCGCCAGCCAGCGCACGGGCCAGATTGGCGACAGACGCATCAGCGCGAGAATAAACCCGAGCAGTAAGCCAAAAAACATGCCGCCGATACTGAGTTGCAGCGTATACCCAGCCCCTTTCAGCAGGAACGGCAAgg                                                                <  2:1447354/150‑1 (MQ=255)
                                                                                    tAGAGATATAAAAACGCGCCAGCCAGCGCACGGGCCAGATTGGCGACAGACGCATCAGCGCGAGAATAAACCCGAGCAGTAAGCCAAAAAACATGCCGCCGATACTGAGTTGCAGCGTATACCCAGCCCCTTTCAGCAGGAACGGCAAgg                                                                <  1:153373/150‑1 (MQ=255)
                                                                                       agaTATAAAAACGCGCCAGCCAGCGCACGGGCCAGATTGGCGACAGACGCATCAGCGCGAGAATAAACCCGAGCAGTAAGCCAAAAAACATGCCGCCGATACTGAGTTGCAGCGTATACCCAGCCCCTTTCAGCAGGAACGGCAAGGAAt                                                             <  2:3904619/150‑1 (MQ=255)
                                                                                        gaTATAAAAACGCGCCAGCCAGCGCACGGGCCAGATTGGCGACAGACGCATCAGCGCGAAAATAAACCCGAGCAGTAAGCCAAAAAACATGCCGCCGATACTGAGTTGCAGCGTATACCCAGCCCCTTTCAGCAGGAACGGCAAGGAATc                                                            >  2:2046730/1‑150 (MQ=255)
                                                                                           ataAAACCGCGCCAGCCAGCGCACGGGCCAGATTGGCGACAGACGCATCAGCGCGAAAATAAACCCGAGCAGTAAGCCAAAAAACATGCCGCCGATACTGAGTTGCAGCGTATACCCAGCCCCTTTCAGCAGGAACGGCAAGGAATCTAt                                                         <  1:3578836/150‑1 (MQ=255)
                                                                                           ataAAAACGCGCCAGCCAGCGCACGGGCCAGATTGGCGACAGACGCATCAGCGCGAAAATAAACCCGAGCAGTAAGCCAAAAAACATGCCGCCGATACTGAGTTGCAGCGTATACCCAGCCCCTTTCAGCAGGAACGGCAAGGAATCAAt                                                         <  1:2578133/150‑1 (MQ=255)
                                                                                             aaaaaCGCGCCAGCCAGCGCACGGGCCAGATTGGCGACAGACGCATCAGCGCGAGAATAAACCCGAGCAGTAAGCCAAAAAACATGCCGCCGATACTGAGTTGCAGCGTATACCCAGCCCCTTTCAGCAGGAACGGCAAGGAATCAATaa                                                       >  2:1738039/1‑150 (MQ=255)
                                                                                                aaCGCGCCAGCCAGCGCACGGGCCAGATTGGCGACAGACGCATCAGCGCGAGAATAAACCCGAGCAGTAAGCCAAAAAACATGCCGCCGATACTGAGTTGCAGCGTATACCCAGCCCCTTTCAGCAGGAACGGCAAGGAATCAATAACCa                                                    >  2:2750406/1‑150 (MQ=255)
                                                                                                  cgcgCCAGCCAGCGCACGGGCCAGATTGGCGACAGACGCATCAGCGCGAGAATAAACCCGAGCAGTAAGCCAAAAAACATGCCGCCGATACTGAGTTGCAGCGTATACCCAGCCCCTTTCAGCAGGAACGGCAAGGAATCAATAACCAGt                                                  <  2:1837093/150‑1 (MQ=255)
                                                                                                        agccagcGCACGGGCCAGATTGGCGACAGACGCATCAGCGCGAGAATAAACCCGAGCAGTAAGCCAAAAAACATGCCGCCGATACTGAGTTGCAGCGTATACCCAGCCCCTTTCAGCAGGAACGGCAAGGAATCAATAACCAGTTGtata                                            >  1:1356255/1‑150 (MQ=255)
                                                                                                             gcgcACGGGCCAGATTGGCGACAGACGCATCAGCGCGAGAATAAACCCGAGCAGTAAGCCAAAAAACATGCCGCCGATACTGAGTTGCAGCGTATACCCAGCCCCTTTCAGCAGGAACGGCAAGGAATCAATAACCAGTTGTATACTTTc                                       <  2:3005661/150‑1 (MQ=255)
                                                                                                              cgcACGGGCCAGATTGGCGACAGACGCATCAGCGCGAGAATAAACCCGAGCAGTAAGCCAAAAAACATGCCGCCGATACTGAGTTGCAGCGTATACCCAGCCCCTTTCAGCAGGAACGGCAAGGAATCAATAACCAGTTGTATACTTTCt                                      <  2:1086822/150‑1 (MQ=255)
                                                                                                               gcACGGGCCAGATTGGCGACAGACGCATCAGCGCGAGAATAAACCCGAGCAGTAAGCCAAAAAACATGCCGCCGATACTGAGTTGCAGCGTATACCCAGCCCCTTTCAGCAGGAACGGCAAGGAATCAATAACCAGTTGTATACTTTCtt                                     <  1:1896799/150‑1 (MQ=255)
                                                                                                                cACGGGCCAGATTGGCGACAGACGCATCAGCGCGAGAATAAACCCGAGCAGTAAGCCAAAAAACATGCCGCCGATACTGAGTTGCAGCGTATACCCAGCCCCTTTCAGCAGGAACGGCAAGGAATCAATAACCAGTTGTATACTTTCTTg                                    >  2:2503622/1‑150 (MQ=255)
                                                                                                                cACGGGCCAGATTGGCGACAGACGCATCAGCGCGAGAATAAACCCGAGCAGTAAGCCAAAAAACATGCCGCCGATACTGAGTTGCAGCGTATACCCAGCCCCTTTCAGCAGGAACGGCAAGGAATCAATAACCAGTTGTATACTTTCTTg                                    <  1:720363/150‑1 (MQ=255)
                                                                                                                 aCGGGCCAGATTGGCGACAGACGCATCAGCGCGAGAATAAACCCGAGCAGTAAGCCAAAAAACATGCCGCCGATACTGAGTTGCAGCGTATACCCAGCCCCTTTCAGCAGGAACGGCAAGGAATCAATAACCAGTTGTATACTTTCTTGc                                   >  2:2704889/1‑150 (MQ=255)
                                                                                                                 aCGGGCCAGATTGGCGACAGACGCATCAGCGCGAGAATAAACCCGAGCAGTAAGCCAAAAAACATGCCGCCGATACTGAGTTGCAGCGTATACCCAGCCCCTTTCAGCAGGAACGGCAAGGAATCAATAACCAGTTGTATACTTTCTTGc                                   >  2:3295721/1‑150 (MQ=255)
                                                                                                                  cGGGCCAGATTGGCGACAGACGCATCAGCGCGAGAATAAACCCGAGCAGTAAGCCAAAAAACATGCCGCCGATACTGAGTTGCAGCGTATACCCAGCCCCTTTCAGCAGGAACGGCAAGGAATCAATAACCAGTTGTATACTTTCTTGCa                                  <  1:928596/150‑1 (MQ=255)
                                                                                                                   gggCCAGATTGGCGACAGACGCATCAGCGCGAGAATAAACCCGAGCAGTAAGCCAAAAAACATGCCGCCGATACTGAGTTGCAGCGTATACCCAGCCCCTTTCAGCAGGAACGGCAAGGAATCAATAACCAGTTGTATACTTTCTTGCAt                                 >  2:3611236/1‑150 (MQ=255)
                                                                                                                     gCCAGATTGGCGACAGACGCATCAGCGCGAGAATAAACCCGAGCAGTAAGCCAAAAAACATGCCGCCGATACTGAGTTGCAGCGTATACCCAGCCCCTTTCAGCAGGAACGGCAAGGAATCAATAACCAGTTGTATACTTTCTTGCATGa                               <  1:1429596/150‑1 (MQ=255)
                                                                                                                      ccAGATTGGCGACAGACGCATCAGCGCGAGAATAAACCCGAGCAGTAAGCCAAAAAACATGCCGCCGATACTGAGTTGCAGCGTATACCCAGCCCCTTTCAGCAGGAACGGCAAGGAATCAATAACCAGTTGTATACTTTCTTGCATGAg                              <  2:3390672/150‑1 (MQ=255)
                                                                                                                       cAGATTGGCGACAGACGCATCAGCGCGAGAATAAACCCGAGCAGTAAGCCAAAAAACATGCCGCCGATACTGAGTTGCAGCGTATACCCAGCCCCTTTCAGCAGGAACGGCAAGGAATCAATAACCAGTTGTATACTTTCTTGCATGAGc                             >  1:325800/1‑150 (MQ=255)
                                                                                                                             ggCGACAGACGCATCAGCGCGAGAATAAACCCGAGCAGTAAGCCAAAAAACATGCCGCCGATACTGAGTTGCAGCGTATACCCAGCCCCTTTCAGCAGGAACGGCAAGGAATCAATAACCAGTTGTATACTTTCTTGCATGAGCGTTTcc                       <  1:1983710/150‑1 (MQ=255)
                                                                                                                               cGACAGACGCATCAGCGCGAGAATAAACCCGAGCAGTAAGCCAAAAAACATGCCGCCGATACTGAGTTGCAGCGTATACCCAGCCCCTTTCAGCAGGAACGGCAAGGAATCAATAACCAGTTGTATACTTTCTTGCATGAGCGTTTCCgg                     <  1:1658755/150‑1 (MQ=255)
                                                                                                                                 aCAGACGCATCAGCGCGAGAATAAACCCGAGCAGTAAGCCAAAAAACATGCCGCCGATACTGAGTTGCAGCGTATACCCAGCCCCTTTCAGCAGGAACGGCAAGGAATCAATAACCAGTTGTATACTTTCTTGCATGAGCGTTTCCGGAg                   <  2:2626722/150‑1 (MQ=255)
                                                                                                                                 aCAGACGCATCAGCGCGAGAATAAACCCGAGCAGTAAGCCAAAAAACATGCCGCCGATACTGAGTTGCAGCGTATACCCAGCCCCTTTCAGCAGGAACGGCAAGGAATCAATAACCAGTTGTATACTTTCTTGCATGAGCGTTTCCGGAg                   >  1:3928435/1‑150 (MQ=255)
                                                                                                                                   aGACGCATCAGCGCGAGAATAAACCCGAGCAGTAAGCCAAAAAACATGCCGCCGATACTGAGTTGCAGCGTATACCCAGCCCCTTTCAGCAGGAACGGCAAGGAATCAATAACCAGTTGTATACTTTCTTGCATGAGCGTTTCCGGAGCt                 <  1:1019547/150‑1 (MQ=255)
                                                                                                                                    gACGCATCAGCGCGAGAATAAACCCGAGCAGTAAGCCAAAAAACATGCCGCCGATACTGAGTTGCAGCGTATACCCAGCCCCTTTCAGCAGGAACGGCAAGGAATCAATAACCAGTTGTATACTTTCTTGCATGAGCGTTTCCGGAGCTa                <  1:3669717/150‑1 (MQ=255)
                                                                                                                                      cGCATCAGCGCGAGAATAAACCCGAGCAGTAAGCCAAAAAACATGCCGCCGATACTGAGTTGCAGCGTATACCCAGCCCCTTTCAGCAGGAACGGCAAGGAATCAATAACCAGTTGTATACTTTCTTGCATGAGCGTTTCCGGAGCTaaa              >  1:1111435/1‑150 (MQ=255)
                                                                                                                                        cATCAGCGCGAAAATAAACCCGAGCAGTAAGCCAAAAAACATGCCGCCGATACTGAGTTGCAGCGTATACCCAGCCCCTTTCAGCAGGAACGGCAAGGAATCAATAACCAGTTGTATACTTTCTTGCATGAGCGTTTCCGGAGCTAAACg            <  1:80168/150‑1 (MQ=255)
                                                                                                                                          tCAGCGCGAGAATAAACCCGAGCAGTAAGCCAAAAAACATGCCGCCGATACTGAGTTGCAGCGTATACCCAGCCCCTTTCAGCAGGAACGGCAAGGAATCAATAACCAGTTGTATACTTTCTTGCATGAGCGTTTCCGGAGCTAAACGTg          >  1:3483175/1‑150 (MQ=255)
                                                                                                                                           cAGCGCGAGAATAAACCCGAGCAGTAAGCCAAAAAACATGCCGCCGATACTGAGTTGCAGCGTATACCCAGCCCCTTTCAGCAGGAACGGCAAGGAATCAATAACCAGTTGTATACTTTCTTGCATGAGCGTTTCCGGAGCTAAACGTgg         >  2:879942/1‑150 (MQ=255)
                                                                                                                                            aGCGCGAGAATAAACCCGAGCAGTAAGCCAAAAAACATGCCGCCGATACTGAGTTGCAGCGTATACCCAGCCCCTTTCAGCAGGAACGGCAAGGAATCAATAACCAGTTGTATACTTTCTTGCATGAGCGTTTCCGGAGCTAAACGTggg        >  2:352988/1‑150 (MQ=255)
                                                                                                                                            aGCGCGAGAATAAACCCGAGCAGTAAGCCAAAAAACATGCCGCCGATACTGAGTTGCAGCGTATACCCAGCCCCTTTCAGCAGGAACGGCAAGGAATCAATAACCAGTTGTATACTTTCTTGCATGAGCGTTTCCGGAGCTAAACGTggg        >  2:3364443/1‑150 (MQ=255)
                                                                                                                                            aGCGCGAGAATAAACCCGAGCAGTAAGCCAAAAAACATGCCGCCGATACTGAGTTGCAGCGTATACCCAGCCCCTTTCAGCAGGAACGGCAAGGAATCAATAACCAGTTGTATACTTTCTTGCATGAGCGTTTCCGGAGCTAAACGTggg        >  1:4567542/1‑150 (MQ=255)
                                                                                                                                             gcgcGAGAATAAACCCGAGCAGTAAGCCAAAAAACATGCCGCCGATACTGAGTTGCAGCGTATACCCAGCCCCTTTCAGCAGGAACGGCAAGGAATCAATAACCAGTTGTATACTTTCTTGCATGAGCGTTTCCGGAGCTAAACGTgggg       <  1:3008209/150‑1 (MQ=255)
                                                                                                                                             gcgcGAGAATAAACCCGAGCAGTAAGCCAAAAAACATGCCGCCGATACTGAGTTGCAGCGTATACCCAGCCCCTTTCAGCAGGAACGGCAAGGAATCAATAACCAGTTGTATACTTTCTTGCATGAGCGTTTCCGGAGCTAAACGTgggg       <  1:890590/150‑1 (MQ=255)
                                                                                                                                                cgAGAATAAACCCGAGCAGTAAGCCAAAAAACATGCCGCCGATACTGAGTTGCAGCGTATACCCAGCCCCTTTCAGCAGGAACGGCAAGGAATCAATAACCAGTTGTATACTTTCTTGCATGAGCGTTTCCGGAGCTAAACGTGGGgatg    <  2:3931465/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
GAATGGGATCTAATTCGATACCAAACTGCGGCAGGCCGTAATAGATCATAAACAGTTGTGCGATGAGCGGCGTACCACGAAAAATAGAGATATAAAAACGCGCCAGCCAGCGCACGGGCCAGATTGGCGACAGACGCATCAGCGCGAGAATAAACCCGAGCAGTAAGCCAAAAAACATGCCGCCGATACTGAGTTGCAGCGTATACCCAGCCCCTTTCAGCAGGAACGGCAAGGAATCAATAACCAGTTGTATACTTTCTTGCATGAGCGTTTCCGGAGCTAAACGTGGGGATGAT  >  CP009273/1991696‑1991991

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: ‑