Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,396,554 C→A 6.6% A46A (GCG→GCT csrD ← targeting factor for csrBC sRNA degradation

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,396,5540CA 6.6% 974.3 / 34.7 256A46A (GCG→GCTcsrDtargeting factor for csrBC sRNA degradation
Reads supporting (aligned to +/- strand):  ref base C (115/124);  new base A (7/10);  total (122/134)
Fisher's exact test for biased strand distribution p-value = 6.24e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

Only 100 of 257 total aligned reads displayed.

CGATAACTACCATTTCTGGCCAGGGTATAAACCTGTTTATCACCATGGAGCAGGTCTACACGAACGATATCTGCCGACATCATTAATTCGGTAATTTGTGGCCTTAATACGCTGAAGTCATTCGACACAAGGTGGGTATCAATCGCCGTCGCCACCGCCTGAACGCGATGACTAAACTTATACTGAATGGCGTTGTAGAAACTTAGCGAACAGCCCAGCAAAGTCACAAAAATTGTTAACCCGGTGAGCAGCGTAACAAAGGCCGAAAATTTCGTCGTTAATCTCATCCTTGTGTTAA  >  CP009273/3396405‑3396702
                                                                                                                                                     |                                                                                                                                                    
 gATAACTACCATTTCTGGCCAGGGTATAAACCTGTTTATCACCATGGAGCAGGTCTACACGAACGATATCTGCCGACATCATTAATTCGGTAATTTGTGGCCTTAATACGCTGAAGTCATTCGACACAAGGTGGGTATCAATCGCCGTCg                                                                                                                                                     <  1:4025417/150‑1 (MQ=255)
 gATAACTACCATTTCTGGCCAGGGTATAAACCTGTTTATCACCATGGAGCAGGTCTACACGAACGATATCTGCCGACATCATTAATTCGGTAATTTGTGGCCTTAATACGCTGAAGTCATTCGACACAAGGTGGGTATCAATCGCCGTCg                                                                                                                                                     <  1:1636675/150‑1 (MQ=255)
 gATAACTACCATTTCTGGCCAGGGTATAAACCTGTTTATCACCATGGAGCAGGTCTACACGAACGATATCTGCCGACATCATTAATTCGGTAATTTGTGGCCTTAATACGCTGAAGTCATTCGACACAAGGTGGGTATCAATCGCCGTCg                                                                                                                                                     >  1:1705033/1‑150 (MQ=255)
    aaCTACCATTTCTGGCCAGGGTATAAACCTGTTTATCACCATGGAGCAGGTCTACACGAACGATATCTGCCGACATCATTAATTCGGTAATTTGTGGCCTTAATACGCTGAAGTCATTCGACACAAGGTGGGTATCAATCGCCGTCGCCa                                                                                                                                                  >  1:4728719/1‑150 (MQ=255)
    aaCTACCATTTCTGGCCAGGGTATAAACCTGTTTATCACCATGGAGCAGGTCTACACGAACGATATCTGCCGACATCATTAATTCGGTAATTTGTGGCCTTAATACGCTGAAGTCATTCGACACAAGGTGGGTATCAATCGCCGTCGCCa                                                                                                                                                  >  2:2570541/1‑150 (MQ=255)
     aCTACCATTTCTGGCCAGGGTATAAACCTGTTTATCACCATGGAGCAGGTCTACACGAACGATATCTGCCGACATCATTAATTCGGTAATTTGTGGCCTTAATACGCTGAAGTCATTCGACACAAGGTGGGTATCAATCGCCGTCGCCAc                                                                                                                                                 <  1:3479587/150‑1 (MQ=255)
     aCTACCATTTCTGGCCAGGGTATAAACCTGTTTATCACCATGGAGCAGGTCTACACGAACGATATCTGCCGACATCATTAATTCGGTAATTTGTGGCCTTAATACGCTGAAGTCATTCGACACAAGGTGGGTATCAATCGCCGTCGCCAc                                                                                                                                                 <  1:2830803/150‑1 (MQ=255)
      cTACCATTTCTGGCCAGGGTATAAACCTGTTTATCACCATGGAGCAGGTCTACACGAACGATATCTGCCGACATCATTAATTCGGTAATTTGTGGCCTTAATACGCTGAAGTCATTCGACACAAGGTGGGTATCAATCGCCGTCGCCAcc                                                                                                                                                <  2:336716/150‑1 (MQ=255)
      cTACCATTTCTGGCCAGGGTATAAACCTGTTTATCACCATGGAGCAGGTCTACACGAACGATATCTGCCGACATCATTAATTCGGTAATTTGTGGCCTTAATACGCTGAAGTCATTCGACACAAGGTGGGTATCAATCGCCGTCGCCAcc                                                                                                                                                <  1:59149/150‑1 (MQ=255)
         ccATTTCTGGCCAGGGTATAAACCTGTTTATCACCAAGGAGCAGGTCTACACGAACGATATCTGCCGACATCATTAATTCGGTAATTTGTGGCCTTAATACGCTGAAGTCATTCGACACAAGGTGGGTATCAATCGCCGTCGCCACCGcc                                                                                                                                             <  1:4297501/150‑1 (MQ=255)
          cATTTCTGGCCAGGGTATAAACCTGTTTATCACCATGGAGCAGGTCTACACGAACGATATCTGCCGACATCATTAATTCGGTAATTTGTGGCCTTAATACGCTGAAGTCATTCGACACAAGGTGGGTATCAATCGCCGTCGCCACCGCCt                                                                                                                                            <  1:4696003/150‑1 (MQ=255)
          cATTTCTGGCCAGGGTATAAACCTGTTTATCACCATGGAGCAGGTCTACACGAACGATATCTGCCGACATCATTAATTCGGTAATTTGTGGCCTTAATACGCTGAAGTCATTCGACACAAGGTGGGTATCAATCGCCGTCGCCACCGCCt                                                                                                                                            <  1:2315860/150‑1 (MQ=255)
             ttCTGGCCAGGGTATAAACCTGTTTATCACCATGGAGCAGGTCTACACGAACGATATCTGCCGACATCATTAATTCGGTAATTTGTGGCCTTAATACGCTGAAGTCATTCGACACAAGGTGGGTATCAATCGCCGTCGCCACCGCCTGaa                                                                                                                                         >  2:1398970/1‑150 (MQ=255)
             ttCTGGCCAGGGTATAAACCTGTTTATCACCATGGAGCAGGTCTACACGAACGATATCTGCCGACATCATTAATTCGGTAATTTGTGGCCTTAATACGCTGAAGTCATTCGACACAAGGTGGGTATCAATCGCCGTCGCCACCGCCTGaa                                                                                                                                         >  1:1908255/1‑150 (MQ=255)
              tCTGGCCAGGGTATAAACCTGTTTATCACCATGGAGCAGGTCTACACGAACGATATCTGCCGACATCATTAATTCGGTAATTTGTGGCCTTAATACGCTGAAGTCATTCGACACAAGGTGGGTATCAATCGCCGTCGCCACCGCCTGAAc                                                                                                                                        <  1:1204061/150‑1 (MQ=255)
                tGGCCAGGGTATAAACCTGTTTATCACCATGGAGCAGGTCTACACGAACGATATCTGCCGACATCATTAATTCGGTAATTTGTGGCCTTAATACGCTGAAGTCATTCGACACAAGGTGGGTATCAATCGCCGTCGCCACCGCCTGAAcgc                                                                                                                                      >  1:3915336/1‑150 (MQ=255)
                 ggCCAGGGTATAAACCTGTTTATCACCATGGAGCAGGTCTACACGAACGATATCTGCCGACATCATTAATTCGGTAATTTGTGGCCTTAATACGCTGAAGTCATTCGACACAAGGTGGGTATCAATCGCCGTCGCCACCGCCTGAAcgcg                                                                                                                                     <  1:3618468/150‑1 (MQ=255)
                   ccAGGGTATAAACCTGTTTATCACCATGGAGCAGGTCTACACGAACGATATCTGCCGACATCATTAATTCGGTAATTTGTGGCCTTAATACGCTGAAGTCATTCGACACAAGGTGGGTATCAATCGCCGTCGCCACCGCCTGAACGCGAt                                                                                                                                   <  2:3769878/150‑1 (MQ=255)
                   ccAGGGTATAAACCTGTTTATCACCATGGAGCAGGTCTACACGAACGATATCTGCCGACATCATTAATTCGGTAATTTGTGGCCTTAATACGCTGAAGTCATTCGACACAAGGTGGGTATCAATCGCCGTAGCCACCGCCTGAACGCGAt                                                                                                                                   <  2:1137634/150‑1 (MQ=255)
                    cAGGGTATAAACCTGTTTATCACCATGGAGCAGGTCTACACGAACGATATCTGCCGACATCATTAATTCGGTAATTTGTGGCCTTAATACGCTGAAGTCATTCGACACAAGGTGGGTATCAATCGCCGTCGCCACCGCCTGAACGCGATg                                                                                                                                  >  1:3537835/1‑150 (MQ=255)
                    cAGGGTATAAACCTGTTTATCACCATGGAGCAGGTCTACACGAACGATATCTGCCGACATCATTAATTCGGTAATTTGTGGCCTTAATACGCTGAAGTCATTCGACACAAGGTGGGTATCAATCGCCGTCGCCACCGCCTGAACGCGATg                                                                                                                                  <  2:8886/150‑1 (MQ=255)
                    cAGGGTATAAACCTGTTTATCACCATGGAGCAGGTCTACACGAACGATATCTGCCGACATCATTAATTCGGTAATTTGTGGCCTTAATACGCTGAAGTCATTCGACACAAGGTGGGTATCAATCGCCGTAGCCACCGCCTGAACGCGATg                                                                                                                                  <  1:3323476/150‑1 (MQ=255)
                     aGGGTATAAACCTGTTTATCACCATGGAGCAGGTCTACACGAACGATATCTGCCGACATCATTAATTCGGTAATTTGTGGCCTTAATACGCTGAAGTCATTCGACACAAGGTGGGTATCAATCGCCGTCGCCACCGCCTGAACGCGATGa                                                                                                                                 >  2:2130161/1‑150 (MQ=255)
                       ggTATACACCTGTTTATCACCATGGAGCAGGTCTACACGAACGATATCTGCCGACATCATTAATTCGGTAATTTGTGGCCTTAATACGCTGAAGTCATTCGACACAAGGTGGGTATCAATCGCCGTCGCCACCGCCTGAACGCGATGACt                                                                                                                               <  2:1090488/150‑1 (MQ=255)
                        gTATAAACCTGTTTATCACCATGGAGCAGGTCTACACGAACGATATCTGCCGACATCATTAATTCGGTAATTTGTGGCCTTAATACGCTGAAGTCATTCGACACAAGGTGGGTATCAATCGCCGTCGCCACCGCCTGAACGCGATGACTa                                                                                                                              >  2:3485895/1‑150 (MQ=255)
                          ataAACCTGTTTATCACCATGGAGCAGGTCTACACGAACGATATCTGCCGACATCATTAATTCGGTAATTTGTGGCCTTAATACGCTGAAGTCATTCGACACAAGGTGGGTATCAATCGCCGTCGCCACCGCCTGAACGCGATGACTaaa                                                                                                                            <  1:2168918/150‑1 (MQ=255)
                          ataAACCTGTTTATCACCATGGAGCAGGTCTACACGAACGATATCTGCCGACATCATTAATTCGGTAATTTGTGGCCTTAATACGCTGAAGTCATTCGACACAAGGTGGGTATCAATCGCCGTCGCCACCGCCTGAACGCGATGACTaaa                                                                                                                            >  2:3522796/1‑150 (MQ=255)
                          ataAACCTGTTTATCACCATGGAGCAGGTCTACACGAACGATATCTGCCGACATCATTAATTCGGTAATTTGTGGCCTTAATACGCTGAAGTCATTCGACACAAGGTGGGTATCAATCGCCGTCGCCACCGCCTGAACGCGATGACTaaa                                                                                                                            <  2:2979042/150‑1 (MQ=255)
                          ataAACCTGTTTATCACCATGGAGCAGGTCTACACGAACGATATCTGCCGACATCATAAATTCGGTAATTTGTGGCCTTAATACGCTGAAGTCATTCGACACAAGGTGGGTATCAATCGCCGTCGCCACCGCCTGAACGCGATGACTaaa                                                                                                                            <  2:36330/150‑1 (MQ=255)
                            aaaCCTGTTTATCACCATGGAGCAGGTCTACACGAACGATATCTGCCGACATCATTAATTCGGTAATTTGTGGCCTTAATACGCTGAAGTCATTCGACACAAGGTGGGTATCAATCGCCGTAGCCACCGCCTGAACGCGATGACTAAAca                                                                                                                          <  1:2023931/150‑2 (MQ=255)
                             aaCCTGTTTATCACCATGGAGCAGGTCTACACGAACGATATCTGCCGACATCATTAATTCGGTAATTTGTGGCCTTAATACGCTGAAGTCATTCGACACAAGGTGGGTATCAATCGCCGTCGCCACCGCCTGAACGCGATGACTAAACtt                                                                                                                         <  1:3767358/150‑1 (MQ=255)
                              aCCTGTTTATCACCATGGAGCAGGTCTACACGAACGATATCTGCCGACATCATTAATTCGGTAATTTGTGGCCTTAATACGCTGAAGTCATTCGACACAAGGTGGGTATCAATCGCCGTCGCCACCGCCTGAACGCGATGACTAAACTta                                                                                                                        <  2:3049358/150‑1 (MQ=255)
                                cTGTTTATCACCATGGAGCAGGTCTACACGAACGATATCTGCCGACATCATTAATTCGGTAATTTGTGGCCTTAATACGCTGAAGTCATTCGACACAAGGTGGGTATCAATCGCCGTCGCCACCGCCTGAACGCGATGACTAAACTtata                                                                                                                      <  2:292578/150‑1 (MQ=255)
                                  gTTTATCACCATGGAGCAGGTCTACACGAACGATATCTGCCGACATCATTAATTCGGTAATTTGTGGCCTTAATACGCTGAAGTCATTCGACACAAGGTGGGTATCAATCGCCGTCGCCACCGCCTGAACGCGATGACTAAACTTATACt                                                                                                                    <  1:731894/150‑1 (MQ=255)
                                    ttATCACCATGGAGCAGGTCTACACGAACGATATCTGCCGACATCATTAATTCGGTAATTTGTGGCCTTAATACGCTGAAGTCATTCGACACAAGGTGGGTATCAATCGCCGTCGCCACCGCCTGAACGCGATGACTAAACTTATACTGa                                                                                                                  >  1:4415940/1‑150 (MQ=255)
                                    ttATCACCATGGAGCAGGTCTACACGAACGATATCTGCCGACATCATTAATTCGGTAATTTGTGGCCTTAATACGCTGAAGTCATTCGACACAAGGTGGGTATCAATCGCCGTCGCCACCGCCTGAACGCGATGACTAAACTTATACTGa                                                                                                                  <  1:1118600/150‑1 (MQ=255)
                                      aTCACCATGGAGCAGGTCTACACGAACGATATCTGCCGACATCATTAATTCGGTAATTTGTGGCCTTAATACGCTGAAGTCATTCGACACAAGGTGGGTATCAATCGCCGTCGCCACCGCCTGAACGCGATGACTAAACTTATACTGAAt                                                                                                                <  2:2731527/150‑1 (MQ=255)
                                      aTCACCATGGAGCAGGTCTACACGAACGATATCTGCCGACATCATTAATTCGGTAATTTGTGGCCTTAATACGCTGAAGTCATTCGACACAAGGTGGGTATCAATCGCCGTCGCCACCGCCTGAACGCGATGACTAAACTTATACTGAAt                                                                                                                <  2:2731907/150‑1 (MQ=255)
                                         aCCATGGAGCAGGTCTACACGAACGATATCTGCCGACATCATTAATTCGGTAATTTGTGGCCTTAATACGCTGAAGTCATTCGACACAAGGTGGGTATCAATCGCCGTCGCCACCGCCTGAACGCGATGACTAAACTTATACTGAATGGc                                                                                                             <  2:829519/150‑1 (MQ=255)
                                           cATGGAGCAGGTCTACACGAACGATATCTGCCGACATCATTAATTCGGTAATTTGTGGCCTTAATACGCTGAAGTCATTCGACACAAGGTGGGTATCAATCGCCGTCGCCACCGCCTGAACGCGATGACTAAACTTATACTGAATGGCGt                                                                                                           >  2:241619/1‑150 (MQ=255)
                                           cATGGAGCAGGTCTACACGAACGATATCTGCCGACATCATTAATTCGGTAATTTGTGGCCTTAATACGCTGAAGTCATTCGACACAAGGTGGGTATCAATCGCCGTCGCCACCGCCTGAACGCGATGACTAAACTTATACTGAATGGCGt                                                                                                           >  1:974859/1‑150 (MQ=255)
                                           cATGGAGCAGGTCTACACGAACGATATCTGCCGACATCATTAATTCGGTAATTTGTGGCCTTAATACGCTGAAGTCATTCGACACAAGGTGGGTATCAATCGCCGTCGCCACCGCCTGAACGCGATGACTAAACTTATACTGAATGGCGt                                                                                                           <  2:1747158/150‑1 (MQ=255)
                                             tGGAGCAGGTCTACACGAACGATATCTGCCGACATCATTAATTCGGTAATTTGTGGCCTTAATACGCTGAAGTCATTCGACACAAGGTGGGTATCAATCGCCGTCGCCACCGCCTGAACGCGATGACTAAACTTATACTGAATGGCGTTg                                                                                                         <  2:4145401/150‑1 (MQ=255)
                                                aGCAGGTCTACACGAACGATATCTGCCGACATCATTAATTCGGTAATTTGTGGCCTTAATACGCTGAAGTCATTCGACACAAGGTGGGTATCAATCGCCGTCGCCACCGCCTGAACGCGATGACTAAACTTATACTGAATGGCGTTGTAg                                                                                                      <  1:2527373/150‑1 (MQ=255)
                                                  cAGGTCTACACGAACGATATCTGCCGACATCATTAATTCGGTAATTTGTGGCCTTAATACGCTGAAGTCATTCGACACAAGGTGGGTATCAATCGCCGTCGCCACCGCCTGAACGCGATGACTAAACTTATACTGAATGGCGTTGTAGaa                                                                                                    >  2:1301442/1‑150 (MQ=255)
                                                  cAGGTCTACACGAACGATATCTGCCGACATCATTAATTCGGTAATTTGTGGCCTTAATACGCTGAAGTCATTCGACACAAGGTGGGTATCAATCGCCGTCGCCACCGCCTGAACGCGATGACTAAACTTATACTGAATGGCGTTGTAGaa                                                                                                    >  2:3696334/1‑150 (MQ=255)
                                                     gTCTACACGAACGATATCTGCCGACATCATTAATTCGGTAATTTGTGGCCTTAATACGCTGAAGTCATTCGACACAAGGTGGGTATCAATCGCCGTCGCCACCGCCTGAACGCGATGACTAAACTTATACTGAATGGCGTTGTAGAAACt                                                                                                 <  1:1781462/150‑1 (MQ=255)
                                                       cTACACGAACGATATCTGCCGACATCATTAATTCGGTAATTTGTGGCCTTAATACGCTGAAGTCATTCGACACAAGGTGGGTATCAATCGCCGTCGCCACCGCCTGAACGCGATGACTAAACTTATACTGAATGGCGTTGTAGAAACTTa                                                                                               <  1:2915336/150‑1 (MQ=255)
                                                         acacGAACGATATCTGCCGACATCATTAATTCGGTAATTTGTGGCCTTAATACGCTGAAGTCATTCGACACAAGGTGGGTATCAATCGCCGTCGCCACCGCCTGAACGCGATGACTAAACTTATACTGAATGGCGTTGTAGAAACTTAGc                                                                                             >  1:2964535/1‑150 (MQ=255)
                                                         acacGAACGATATCTGCCGACATCATTAATTCGGTAATTTGTGGCCTTAATACGCTGAAGTCATTCGACACAAGGTGGGTATCAATCGCCGTCGCCACCGCCTGAACGCGATGACTAAACTTATACTGAATGGCGTTGTAGAAACTTAGc                                                                                             >  1:2424768/1‑150 (MQ=255)
                                                           acgaacgaTATCTGCCGACATCATTAATTCGGTAATTTGTGGCCTTAATACGCTGAAGTCATTCGACACAAGGTGGGTATCAATCGCCGTCGCCACCGCCTGAACGCGATGACTAAACTTATACTGAATGGCGTTGTAGAAACTTAGCGa                                                                                           <  1:3593782/150‑1 (MQ=255)
                                                              aacgaTATCTGCCGACATCATTAATTCGGTAATTTGTGGCCTTAATACGCTGAAGTCATTCGACACAAGGTGGGTATCAATCGCCGTCGCCACCGCCTGAACGCGATGACTAAACTTATACTGAATGGCGTTGTAGAAACTTAGCGAACa                                                                                        >  1:264048/1‑150 (MQ=255)
                                                              aacgaTATCTGCCGACATCATTAATTCGGTAATTTGTGGCCTTAATACGCTGAAGTCATTCGACACAAGGTGGGTATCAATCGCCATCGCCACCGCCTGAACGCGATGACTAAACTTATACTGAATGGCGTTGTAGAAACTTAGCGAACa                                                                                        <  2:4506586/150‑1 (MQ=255)
                                                                   tatCTGCCGACATCATTAATTCGGTAATTTGTGGCCTTAATACGCTGAAGTCATTCGACACAAGGTGGGTATCAATCGCCGTCGCCACCGCCTGAACGCGATGACTAAACTTATACTGAATGGCGTTGTAGAAACTTAGCGAACAGCCCa                                                                                   <  1:4619818/150‑1 (MQ=255)
                                                                   tatCTGCCGACATCATTAATTCGGTAATTTGTGGCCTTAATACGCTGAAGTCATTCGACACAAGGTGGGTATCAATCGCCGTCGCCACCGCCTGAACGCGATGACTAAACTTATACTGAATGGCGTTGTAGAAACTTAGCGAACAGCCCa                                                                                   <  1:937205/150‑1 (MQ=255)
                                                                   tatCTGCCGACATCATTAATTCGGTAATTTGTGGCCTTAATACGCTGAAGTCATTCGACACAAGGTGGGTATCAATCGCCGTCGCCACCGCCTGAACGCGATGACTAAACTTATACTGAATGGCGTTGTAGAAACTTAGCGAACAGCCCa                                                                                   <  1:3314274/150‑1 (MQ=255)
                                                                    atCTGCCGACATCATTAATTCGGTAATTTGTGGCCTTAATACGCTGAAGTCATTCGACACAAGGTGGGTATCAATCGCCGTCGCCACCGCCTGAACGCGATGACTAAACTTATACTGAATGGCGTTGTAGAAACTTAGCGAACAGCCCAg                                                                                  <  1:1986028/150‑1 (MQ=255)
                                                                      cTGCCGACATCATTAATTCGGTAATTTGTGGCCTTAATACGCTGAAGTCATTCGACACAAGGTGGGTATCAATCGCCGTCGCCACCGCCTGAACGCGATGACTAAACTTATACTGAATGGCGTTGTAGAAACTTAGCGAACAGCCCAGCa                                                                                <  1:1009434/150‑1 (MQ=255)
                                                                      cTGCCGACATCATTAATTCGGTAATTTGTGGCCTTAATACGCTGAAGTCATTCGACACAAGGTGGGTATCAATCGCCGTCGCCACCGCCTGAACGCGATGACTAAACTTATACTGAATGGCGTTGTAGAAACTTAGCGAACAGCCCAGCa                                                                                <  1:4290604/150‑1 (MQ=255)
                                                                          cGACATCATTAATTCGGTAATTTGTGGCCTTAATACGCTGAAGTCATTCGACACAAGGTGGGTATCAATCGCCGTAGCCACCGCCTGAACGCGATGACTAAACTTATACTGAATGGCGTTGTAGAAACTTAGCGAACAGCCCAGCAAAGt                                                                            >  1:2750206/1‑150 (MQ=255)
                                                                                 atTAATTCGGTAATTTGTGGCCTTAATACGCTGAAGTCATTCGACACAAGGTGGGTATCAATCGCCGTCGCCACCGCCTGAACGCGATGACTAAACTTATACTGAATGGCGTTGTAGAAACTTAGCGAACAGCCCAGCAAAGTCACaaaa                                                                     >  1:172426/1‑150 (MQ=255)
                                                                                 atTAATTCGGTAATTTGTGGCCTTAATACGCTGAAGTCATTCGACACAAGGTGGGTATCAATCGCCGTAGCCACCGCCTGAACGCGATGACTAAACTTATACTGAATGGCGTTGTAGAAACTTAGCGAACAGCCCAGCAAAGTCACaaaa                                                                     <  2:2120644/150‑1 (MQ=255)
                                                                                   tAATTCGGTAATTTGTGGCCTTAATACGCTGAAGTCATTCGACACAAGGTGGGTATCAATCGCCGTCGCCACCGCCTGAACGCGATGACTAAACTTATACTGAATGGCGTTGTAGAAACTTAGCGAACAGCCCAGCAAAGTCACAAAAAt                                                                   >  1:1035119/1‑150 (MQ=255)
                                                                                     aTTCGGTAATTTGTGGCCTTAATACGCTGAAGTCATTCGACACAAGGTGGGTATCAATCGCCGTCGCCACCGCCTGAACGCGATGACTAAACTTATACTGAATGGCGTTGTAGAAACTTAGCGAACAGCCCAGCAAAGTCACAAAAATTg                                                                 <  2:982983/150‑1 (MQ=255)
                                                                                       tCGGTAATTTGTGGCCTTAATACGCTGAAGTCATTCGACACAAGGTGGGTATCAATCGCCGTCGCCACCGCCTGAACGCGATGACTAAACTTATACTGAATGGCGTTGTAGAAACTTAGCGAACAGCCCAGCAAAGTCACAAAAATTGtt                                                               >  1:3401316/1‑150 (MQ=255)
                                                                                         ggTAATTTGTGGCCTTAATACGCTGAAGTCATTCGACACAAGGTGGGTATCAATCGCCGTCGCCACCGCCTGAACGCGATGACTAAACTTATACTGAATGGCGTTGTAGAAACTTAGCGAACAGCCCAGCAAAGTCACAAAAATTGTTaa                                                             >  1:519080/1‑150 (MQ=255)
                                                                                             aTTTGTGGCCTTAATACGCTGAAGTCATTCGACACAAGGTGGGTATCAATCGCCGTCGCCACCGCCTGAACGCGATGACTAAACTTATACTGAATGGCGTTGTAGAAACTTAGCGAACAGCCCAGCAAAGTCACAAAAATTGTTAACCCg                                                         <  2:342000/150‑1 (MQ=255)
                                                                                                tgtgGCCTTAATACGCTGAAGTCATTCGACACAAGGTGGGTATCAATCGCCGTCGCCACCGCCTGAACGCGATGACTAAACTTATACTGAATGGCGTTGTAGAAACTTAGCGAACAGCCCAGCAAAGTCACAAAAATTGTTAACCCGGTg                                                      >  2:2491954/1‑150 (MQ=255)
                                                                                                tgtgGCCTTAATACGCTGAAGTCATTCGACACAAGGTGGGTATCAATCGCCGTCGCCACCGCCTGAACGCGATGACTAAACTTATACTGAATGGCGTTGTAGAAACTTAGCGAACAGCCCAGCAAAGTCACAAAAATTGTTAAACCGGTg                                                      >  2:2490675/1‑150 (MQ=255)
                                                                                                 gtgGCCTTAATACGCTGAAGTCATTCGACACAAGGTGGGTATCAATCGCCGTCGCCACCGCCTGAACGCGATGACTAAACTTATACTGAATGGCGTTGTAGAAACTTAGCGAACAGCCCAGCAAAGTCACAAAAATTGTTAACCCGGTGa                                                     >  2:647883/1‑150 (MQ=255)
                                                                                                 gtgGCCTTAATACGCTGAAGTCATTCGACACAAGGTGGGTATCAATCGCCGTCGCCACCGCCTGAACGCGATGACTAAACTTATACTGAATGGCGTTGTAGAAACTTAGCGAACAGCCCAGCAAAGTCACAAAAATTGTTAACCCGGTGa                                                     <  1:1962717/150‑1 (MQ=255)
                                                                                                      cTTAATACGCTGAAGTCATTCGACACAAGGTGGGTATCAATCGCCGTCGCCACCGCCTGAACGCGATGACTAAACTTATACTGAATGGCGTTGTAGAAACTTAGCGAACAGCCCAGCAAAGTCACAAAAATTGTTAACCCGGTGagcagc                                                <  2:1603199/150‑1 (MQ=255)
                                                                                                      cTTAATACGCTGAAGTCATTCGACACAAGGTGGGTATCAATCGCCGTCGCCACCGCCTGAACGCGATGACTAAACTTATACTGAATGGCGTTGTAGAAACTTAGCGAACAGCCCAGCAAAGTCACAAAAATTGTTAACCCGGTGagcagc                                                <  2:390002/150‑1 (MQ=255)
                                                                                                       ttAATACGCTGAAGTCATTCGACACAAGGTGGGTATCAATCGCCGTCGCCACCGCCTGAACGCGATGACTAAACTTATACTGAATGGCGTTGTAGAAACTTAGCGAACAGCCCAGCAAAGTCACAAAAATTGTTAACCCGGTGAGCAGCg                                               >  2:1994161/1‑150 (MQ=255)
                                                                                                        tAATACGCTGAAGTCATTCGACACAAGGTGGGTATCAATCGCCGTCGCCACCGCCTGAACGCGATGACTAAACTTATACTGAATGGCGTTGTAGAAACTTAGCGAACAGCCCAGCAAAGTCACAAAAATTGTTAACCCGGTGAGCAGCGt                                              >  1:4155748/1‑150 (MQ=255)
                                                                                                         aaTACGCTGAAGTCATTCGACACAAGGTGGGTATCAATCGCCGTCGCCACCGCCTGAACGCGATGACTAAACTTATACTGAATGGCGTTGTAGAAACTTAGCGAACAGCCCAGCAAAGTCACAAAAATTGTTAACCCGGTGAGCAGCGTa                                             >  1:2400064/1‑150 (MQ=255)
                                                                                                          aTACGCTGAAGACATTCGACACAAGGTGGGTATCAAACGCCGTAGCCACCGCCTGAACGCGATGACTAAACTAATACTGAATGGCGTTGTAGAAACTTAGCGAACAGCCCAGCAAAGTCACACAAATTGATAACCCGGTGAGCAGCGaaa                                            >  2:229106/1‑150 (MQ=255)
                                                                                                                 gAAGTCATTCGACACAAGGTGGGTATCAATCGCCGTCGCCACCGCCTGAACGCGATGACTAAACTTATACTGAATGGCGTTGTAGAAACTTAGCGAACAGCCCAGCAAAGTCACAAAAATTGTTAACCCGGTGAGCAGCGTAACAAAGGc                                     >  1:3061229/1‑150 (MQ=255)
                                                                                                                  aaGTCATTCGACACAAGGTGGGTATCAATCGCCGTCGCCACCGCCTGAACGCGATGACTAAACTTATACTGAATGGCGTTGTAGAAACTTAGCGAACAGCCCAGCAAAGTCACAAAAATTGTTAACCCGGTGAGCAGCGTAACAAAGGcc                                    <  1:2153526/150‑1 (MQ=255)
                                                                                                                     tCATTCGACACAAGGTGGGTATCAATCGCCGTCGCCACCGCCTGAACGCGATGACTAAACTTATACTGAATGGCGTTGTAGAAACTTAGCGAACAGCCCAGCAAAGTCACAAAAATTGTTAACCCGGTGAGCAGCGTAACAAAGGCCGaa                                 <  2:2668435/150‑1 (MQ=255)
                                                                                                                      cATTCGACACAAGGTGGGTATCAATCGCCGTCGCCACCGCCTGAACGCGATGACTAAACTTATACTGAATGGCGTTGTAGAAACTTAGCGAACAGCCCAGCAAAGTCACAAAAATTGTTAACCCGGTGAGCAGCGTAACAAAGGCCGaaa                                <  1:1610491/150‑1 (MQ=255)
                                                                                                                         tCGACACAAGGTGGGTATCAATCGCCGTCGCCACCGCCTGAACGCGATGACTAAACTTATACTGAATGGCGTTGTAGAAACTTAGCGAACAGCCCAGCAAAGTCACAAAAATTGTTAACCCGGTGAGCAGCGTAACAAAGGCCGAAAAtt                             <  1:1334537/150‑1 (MQ=255)
                                                                                                                          cGACACAAGGTGGGTATCAATCGCCGTCGCCACCGCCTGAACGCGATGACTAAACTTATACTGAATGGCGTTGTAGAAACTTAGCGAACAGCCCAGCAAAGTCACAAAAATTGTTAACCCGGTGAGCAGCGTAACAAAGGCCGAAAAttt                            <  2:4728719/150‑1 (MQ=255)
                                                                                                                                aaGGTGGGTATCAATCGCCGTCGCCACCGCCTGAACGCGATGACTAAACTTATACTGAATGGCGTTGTAGAAACTTAGCGAACAGCCCAGCAAAGTCACAAAAATTGTTAACCCGGTGAGCAGCGTAACAAAGGCCGAAAATTTcgtcgt                      >  1:1506951/1‑150 (MQ=255)
                                                                                                                                aaGGTGGGTATCAATCGCCGTCGCCACCGCCTGAACGCGATGACTAAACTTATACTGAATGGCGTTGTAGAAACTTAGCGAACAGCCCAGCAAAGTCACAAAAATTGTTAACCCGGTGAGCAGCGTAACAAAGGCCGAAAATTTcgtcgt                      >  2:2173832/1‑150 (MQ=255)
                                                                                                                                  ggTGGGTATCAATCGCCGTCGCCACCGCCTGAACGCGATGACTAAACTTATACTGAATGGCGTTGTAGAAACTTAGCGAACAGCCCAGCAAAGTCACAAAAATTGTTAACCCGGTGAGCAGCGTAACAAAGGCCGAAAATTTCGTCGTTa                    >  2:508898/1‑150 (MQ=255)
                                                                                                                                      ggTATCAATCGCCGTCGCCACCGCCTGAACGCGATGACTAAACTTATACTGAATGGCGTTGTAGAAACTTAGCGAACAGCCCAGCAAAGTCACAAAAATTGTTAACCCGGTGAGCAGCGTAACAAAGGCCGAAAATTTCGTCGTTAAtct                <  2:4083880/150‑1 (MQ=255)
                                                                                                                                         aTCAATCGCCGTCGCCACCGCCTGAACGCGATGACTAAACTTATACTGAATGGCGTTGTAGAAACTTAGCGAACAGCCCAGCAAAGTCACAAAAATTGTTAACCCGGTGAGCAGCGTAACAAAGGCCGAAAATTTCGTCGTTAATCTCAt             >  1:710912/1‑150 (MQ=255)
                                                                                                                                         aTCAATCGCCGTCGCCACCGCCTGAACGCGATGACTAAACTTATACTGAATGGCGTTGTAGAAACTTAGCGAACAGCCCAGCAAAGTCACAAAAATTGTTAACCCGGTGAGCAGCGTAACAAAGGCCGAAAATTTCGTCGTTAATCTCAt             >  1:3081847/1‑150 (MQ=255)
                                                                                                                                          tCAATCGCCGTCGCCACCGCCTGAACGCGATGACTAAACTTATACTGAATGGCGTTGTAGAAACTTAGCGAACAGCCCAGCAAAGTCACAAAAATTGTTAACCCGGTGAGCAGCGTAACAAAGGCCGAAAATTTCGTCGTTAATCTCATc            >  1:2310454/1‑150 (MQ=255)
                                                                                                                                           cAATCGCCGTCGCCACCGCCTGAACGCGATGACTAAACTTATACTGAATGGCGTTGTAGAAACTTAGCGAACAGCCCAGCAAAGTCACAAAAATTGTTAACCCGGTGAGCAGCGTAACAAAGGCCGAAAATTTCGTCGTTAATCTCATcc           >  2:4626253/1‑150 (MQ=255)
                                                                                                                                             aTCGCCGTAGCCACCGCCTGAACGCGATGACTAAACTTATACTGAATGGCGTTGTAGAAACTTAGCGAACAGCCCAGCAAAGTCACAAAAATTGTTAACCCGGTGAGCAGCGTAACAAAGGCCGAAAATTTCGTCGTTAATCTCATCCtt         <  1:1018643/150‑1 (MQ=255)
                                                                                                                                              tCGCCGTCGCCACCGCCTGAACGCGATGACTAAACTTATACTGAATGGCGTTGTAGAAACTTAGCGAACAGCCCAGCAAAGTCACAAAAATTGTTAACCCGGTGAGCAGCGTAACAAAGGCCGAAAATTTCGTCGTTAATCTCATCCTtg        <  1:1772467/150‑1 (MQ=255)
                                                                                                                                                 ccGTCGCCACCGCCTGAACGCGATGACTAAACTTATACTGAATGGCGTTGTAGAAACTTAGCGAACAGCCCAGCAAAGTCACAAAAATTGTTAACCCGGTGAGCAGCGTAACAAAGGCCGAAAATTTCGTCGTTAATCTCATCCTTgtgt     <  1:3784018/150‑1 (MQ=255)
                                                                                                                                                  cGTCGCCACCGCCTGAACGCGATGACTAAACTTATACTGAATGGCGTTGTAGAAACTTAGCGAACAGCCCAGCAAAGTCACAAAAATTGTTAACCCGGTGAGCAGCGTAACAAAGGCCGAAAATTTCGTCGTTAATCTCATCCTTGTGtt    >  2:384496/1‑150 (MQ=255)
                                                                                                                                                  cGTAGCCACCGCCTGAACGCGATGACTAAACTTATACTGAATGGCGTTGTAGAAACTTAGCGAACAGCCCAGCAAAGTCACAAAAATTGTTAACCCGGTGAGCAGCGTAACAAAGGCCGAAAATTTCGTCGTTAATCTCATCCTTGTGtt    <  1:857303/150‑1 (MQ=255)
                                                                                                                                                   gTCGCCACCGCCTGAACGCGATGACTAAACTTATACTGAATGGCGTTGTAGAAACTTAGCGAACAGCCCAGCAAAGTCACAAAAATTGTTAACCCGGTGAGCAGCGTAACAAAGGCCGAAAATTTCGTCGTTAATCTCATCCTTGTGTTa   >  1:2017083/1‑150 (MQ=255)
                                                                                                                                                    tCGCCACCGCCTGAACGCGATGACTAAACTTATACTGAATGGCGTTGTAGAAACTTAGCGAACAGCCCAGCAAAGTCACAAAAATTGTTAACCCGGTGAGCAGCGTAACAAAGGCCGAAAATTTCGTCGTTAATCTCATCCTTGTGTTaa  >  2:1501495/1‑150 (MQ=255)
                                                                                                                                                    tCGCCACCGCCTGAACGCGATGACTAAACTTATACTGAATGGCGTTGTAGAAACTTAGCGAACAGCCCAGCAAAGTCACAAAAATTGTTAACCCGGTGAGCAGCGTAACAAAGGCCGAAAATTTCGTCGTTAATCTCATCCTTGTGTTaa  >  2:1501559/1‑150 (MQ=255)
                                                                                                                                                    tCGCCACCGCCTGAACGCGATGACTAAACTTATACTGAATGGCGTTGTAGAAACTTAGCGAACAGCCCAGCAAAGTCACAAAAATTGTTAACCCGGTGAGCAGCGTAACAAAGGCCGAAAATTTCGTCGTTAATCTCATCCTTGTGTTaa  >  2:107650/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CGATAACTACCATTTCTGGCCAGGGTATAAACCTGTTTATCACCATGGAGCAGGTCTACACGAACGATATCTGCCGACATCATTAATTCGGTAATTTGTGGCCTTAATACGCTGAAGTCATTCGACACAAGGTGGGTATCAATCGCCGTCGCCACCGCCTGAACGCGATGACTAAACTTATACTGAATGGCGTTGTAGAAACTTAGCGAACAGCCCAGCAAAGTCACAAAAATTGTTAACCCGGTGAGCAGCGTAACAAAGGCCGAAAATTTCGTCGTTAATCTCATCCTTGTGTTAA  >  CP009273/3396405‑3396702

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: