Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,655,084 A→C 72.4% L136R (CTG→CGG)  iscR ← isc operon transcriptional repressor; suf operon transcriptional activator; oxidative stress‑ and iron starvation‑inducible; autorepressor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,655,0840AC72.4% 419.8 / 129.2 199L136R (CTG→CGG) iscRisc operon transcriptional repressor; suf operon transcriptional activator; oxidative stress‑ and iron starvation‑inducible; autorepressor
Reads supporting (aligned to +/- strand):  ref base A (31/24);  new base C (82/62);  total (113/86)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

Only 100 of 200 total aligned reads displayed.

GCTGGATGTACGACCGTGTTTACGGAGTATTTAGCACTCCGGCCTGATTCTGAATTCTTTTTATTAAGCGCGTAACTTAACGTCGATCGCGTCTTGTGTGCGGGTGCGTGGCGCGTCGTGAGTATGCTGACGACCAGACACATCCAGCACTTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTA  >  CP009273/2654939‑2655231
                                                                                                                                                 |                                                                                                                                                   
  tGGATGTACGACCGTGTTTACGGAGTATTTAGCACTCCGGCCTGATTCTGAATTCTTTTTATTAAGCGCGTAACTTAACGTCGATCGCGTCTTGTGTGCGGGTGCGTGGCGCGTCGTGAGTATGCTGACGACCAGACACATCCAGCAc                                                                                                                                                 <  2:1678268/148‑1 (MQ=255)
   ggATGTACGACCGTGTTTACGGAGTATTTAGCACTCCGGCCTGATTCTGAATTCTTTTTATTAAGCGCGTAACTTAACGTCGATCGCGTCTTGTGTGCGGGTGCGTGGCGCGTCGTGAGTATGCTGACGACCAGACACATCCAGCACt                                                                                                                                                <  2:663404/148‑1 (MQ=255)
    gATGTACGACCGTGTTTACGGAGTATTTAGCACTCCGGCCTGATTCTGAATTCTTTTTATTAAGCGCGTAACTTAACGTCGATCGCGTCTTGTGTGCGGGTGCGTGGCGCGTCGTGAGTATGCTGACGACCAGACACATCCCGCACtt                                                                                                                                               <  1:1788724/148‑1 (MQ=255)
     aTGTACGACCGTGTTTACGGAGTATTTAGCACTCCGGCCTGATTCTGAATTCTTTTTATTAAGCGCGTAACTTAACGTCGATCGCGTCTTGTGTGCGGGTGCGTGGCGCGTCGTGAGTATGCTGACGACCAGACACATCCCGCACTTc                                                                                                                                              >  2:993203/1‑148 (MQ=255)
     aTGTACGACCGTGTTTACGGAGTATTTAGCACTCCGGCCTGATTCTGAATTCTTTTTATTAAGCGCGTAACTTAACGTCGATCGCGTCTTGTGTGCGGGTGCGTGGCGCGTCGTGAGTATGCTGACGACCAGACACATCCCGCACTTc                                                                                                                                              >  1:2340772/1‑148 (MQ=255)
       gTACGACCGTGTTTACGGAGTATTTAGCACTCCGGCCTGATTCTGAATTCTTTTTATTAAGCGCGTAACTTAACGTCGATCGCGTCTTGTGTGCGGGTGCGTGGCGCGTCGTGAGTATGCTGACGACCAGACACATCCCGCACTTCCt                                                                                                                                            <  2:304109/148‑1 (MQ=255)
       gTACGACCGTGTTTACGGAGTATTTAGCACTCCGGCCTGATTCTGAATTCTTTTTATTAAGCGCGTAACTTAACGTCGATCGCGTCTTGTGTGCGGGTGCGTGGCGCGTCGTGAGTATGCTGACGACCAGACACATCCCGCACTTCCt                                                                                                                                            <  2:160750/148‑1 (MQ=255)
        tACGACCGTGTTTACGGAGTATTTAGCACTCCGGCCTGATTCTGAATTCTTTTTATTAAGCGCGTAACTTAACGTCGATCGCGTCTTGTGTGCGGGTGCGTGGCGCGTCGTGAGTATGCTGACGACCAGACACATCCCGCACTTCCTg                                                                                                                                           <  2:2800010/148‑1 (MQ=255)
          cGACCGTGTTTACGGAGTATTTAGCACTCCGGCCTGATTCTGAATTCTTTTTATTAAGCGCGTAACTTAACGTCGATCGCGTCTTGTGTGCGGGTGCGTGGCGCGTCGTGAGTATGCTGACGACCAGACACATCCCGCACTTCCTGGt                                                                                                                                         >  1:24825/1‑148 (MQ=255)
             ccGTGTTTACGGAGTATTTAGCACTCCGGCCTGATTCTGAATTCTTTTTATTAAGCGCGTAACTTAACGTCGATCGCGTCTTGTGTGCGGGTGCGTGGCGCGTCGTGAGTATGCTGACGACCAGACACATCCCGCACTTCCTGGttat                                                                                                                                      >  1:12052/1‑148 (MQ=255)
             ccGTGTTTACGGAGTATTTAGCACTCCGGCCTGATTCTGAATTCTTTTTATTAAGCGCGTAACTTAACGTCGATCGCGTCTTGTGTGCGGGTGCGTGGCGCGTCGTGAGTATGCTGACGACCAGACACATCCCGCACTTCCTGGttat                                                                                                                                      >  1:3387052/1‑148 (MQ=255)
              cGTGTTTACGGAGTATTTAGCACTCCGGCCTGATTCTGAATTCTTTTTATTAAGCGCGTAACTTAACGTCGATCGCGTCTTGTGTGCGGGTGCGTGGCGCGTCGTGAGTATGCTGACGACCAGACACATCCCGCACTTCCTGGttatt                                                                                                                                     >  2:1833663/1‑148 (MQ=255)
               gtgtTTACGGAGTATTTAGCACTCCGGCCTGATTCTGAATTCTTTTTATTAAGCGCGTAACTTAACGTCGATCGCGTCTTGTGTGCGGGTGCGTGGCGCGTCGTGAGTATGCTGACGACCAGACACATCCAGCACTTCCTGGttatta                                                                                                                                    >  1:2543641/1‑148 (MQ=255)
               gtgtTTACGGAGTATTTAGCACTCCGGCCTGATTCTGAATTCTTTTTATTAAGCGCGTAACTTAACGTCGATCGCGTCTTGTGTGCGGGTGCGTGGCGCGTCGTGAGTATGCTGACGACCAGACACATCCAGCACTTCCTGGttatta                                                                                                                                    >  1:2543684/1‑148 (MQ=255)
                  tttACGGAGTATTTAGCACTCCGGCCTGATTCTGAATTCTTTTTATTAAGCGCGTAACTTAACGTCGATCGCGTCTTGTGTGCGGGTGCGTGGCGCGTCGTGAGTATGCTGACGACCAGACACATCCAGCACTTCCTGGTTATTAAcc                                                                                                                                 <  1:3347676/148‑1 (MQ=255)
                   ttACGGAGTATTTAGCACTCCGGCCTGATTCTGAATTCTTTTTATTAAGCGCGTAACTTAACGTCGATCGCGTCTTGTGTGCGGGTGCGTGGCGCGTCGTGAGTATGCTGACGACCAGACACATCCCGCACTTCCTGGTTATTAACCa                                                                                                                                <  2:1022910/148‑1 (MQ=255)
                      cGGAGTATTTAGCACTCCGGCCTGATTCTGAATTCTTTTTATTAAGCGCGTAACTTAACGTCGATCGCGTCTTGTGTGCGGGTGCGTGGCGCGTCGTGAGTATGCTGACGACCAGACACATCCCGCACTTCCTGGTTATTAACCAGtt                                                                                                                             <  2:3247805/148‑1 (MQ=255)
                      cGGAGTATTTAGCACTCCGGCCTGATTCTGAATTCTTTTTATTAAGCGCGTAACTTAACGTCGATCGCGTCTTGTGTGCGGGTGCGTGGCGCGTCGTGAGTATGCTGACGACCAGACACATCCCGCACTTCCTGGTTATTAACCAGtt                                                                                                                             <  2:1985476/148‑1 (MQ=255)
                      cGGAGTATTTAGCACTCCGGCCTGATTCTGAATTCTTTTTATTAAGCGCGTAACTTAACGTCGATCGCGTCTTGTGTGCGGGTGCGTGGCGCGTCGTGAGTATGCTGACGACCAGACACATCCCGCACTTCCTGGTTATTAACCAGtt                                                                                                                             <  1:2719302/148‑1 (MQ=255)
                         aGTATTTAGCACTCCGGCCTGATTCTGAATTCTTTTTATTAAGCGCGTAACTTAACGTCGATCGCGTCTTGTGTGCGGGTGCGTGGCGCGTCGTGAGTATGCTGACGACCAGACACATCCCGCACTTCCTGGTTATTAACCAGTTCGc                                                                                                                          <  1:2705523/148‑1 (MQ=255)
                           tATTTAGCACTCCGGCCTGATTCTGAATTCTTTTTATTAAGCGCGTAACTTAACGTCGATCGCGTCTTGTGTGCGGGTGCGTGGCGCGTCGTGAGTATGCTGACGACCAGACACATCCCGCACTTCCTGGTTATTAACCAGTTCGCCt                                                                                                                        >  2:1851771/1‑148 (MQ=255)
                             tttAGCACTCCGGCCTGATTCTGAATTCTTTTTATTAAGCGCGTAACTTAACGTCGATCGCGTCTTGTGTGCGGGTGCGTGGCGCGTCGTGAGTATGCTGACGACCAGACACATCCCGCACTTCCTGGTTATTAACCAGTTCGCCTaa                                                                                                                      >  1:44081/1‑148 (MQ=255)
                              ttAGCACTCCGGCCTGATTCTGAATTCTTTTTATTAAGCGCGTAACTTAACGTCGATCGCGTCTTGTGTGCGGGTGCGTGGCGCGTCGTGAGTATGCTGACGACCAGACACATCCCGCACTTCCTGGTTATTAACCAGTTCGCCTaaa                                                                                                                     >  1:815085/1‑148 (MQ=255)
                              ttAGCACTCCGGCCTGATTCTGAATTCTTTTTATTAAGCGCGTAACTGAACGTCGATCGCGTCTTGTGTGCGGGTGCGTGGCGCGTCGTGAGTATGCTGACGACCAGACACATCCCGCACTTCCTGGTTATTAACCAGTTCGCCTaaa                                                                                                                     <  2:2723800/148‑1 (MQ=255)
                               tAGCACTCCGGCCTGATTCTGAATTCTTTTTATTAAGCGCGTAACTTAACGTCGATCGCGTCTTGTGTGCGGGTGCGTGGCGCGTCGTGAGTATGCTGACGACCAGACACATCCCGCACTTCCTGGTTATTAACCAGTTCGCCTAAAg                                                                                                                    <  2:1733877/148‑1 (MQ=255)
                                 gCACTCCGGCCTGATTCTGAATTCTTTTTATTAAGCGCGTAACTTAACGTCGATCGCGTCTTGTGTGCGGGTGCGTGGCGCGTCGTGAGTATGCTGACGACCAGACACATCCCGCACTTCCTGGTTATTAACCAGTTCGCCTAAAGTa                                                                                                                  >  1:1545873/1‑148 (MQ=255)
                                 gCACTCCGGCCTGATTCTGAATTCTTTTTATTAAGCGCGTAACTTAACGTCGATCGCGTCTTGTGTGCGGGTGCGTGGCGCGTCGTGAGTATGCTGACGACCAGACACATCCAGCACTTCCTGGTTATTAACCAGTTCGCCTAAAGTa                                                                                                                  >  1:1227160/1‑148 (MQ=255)
                                  cACTCCGGCCTGATTCTGAATTCTTTTTATTAAGCGCGTAACTTAACGTCGATCGCGTCTTGTGTGCGGGTGCGTGGCGCGTCGTGAGTATGCTGACGACCAGACACATCCCGCACTTCCTGGTTATTAACCAGTTCGCCTAAAGTaa                                                                                                                 >  2:292067/1‑148 (MQ=255)
                                  cACTCCGGCCTGATTCTGAATTCTTTTTATTAAGCGCGTAACTTAACGTCGATCGCGTCTTGTGTGCGGGTGCGTGGCGCGTCGTGAGTATGCTGACGACCAGACACATCCAGCACTTCCTGGTTATTAACCAGTTCGCCTAAAGTaa                                                                                                                 <  2:3460726/148‑1 (MQ=255)
                                    cTCCGGCCTGATTCTGAATTCTTTTTATTAAGCGCGTAACTTAACGTCGATCGCGTCTTGTGTGCGGGTGCGTGGCGCGTCGTGAGTATGCTGACGACCAGACACATCCAGCACTTCCTGGTTATTAACCAGTTCGCCTAAAGTAAtg                                                                                                               >  1:276805/1‑148 (MQ=255)
                                      ccGGCCTGATTCTGAATTCTTTTTATTAAGCGCGTAACTTAACGTCGATCGCGTCTTGTGTGCGGGTGAGTGGCGCGTCGTGAGTAAGCTGACGACCAGACACATCCCGCACTTCCTGGTTATTAACCAGTTCGCCTAAAGTAAtgtt                                                                                                             >  2:2388552/1‑148 (MQ=255)
                                        ggCCTGATTCTGAATTCTTTTTATTAAGCGCGTAACTTAACGTCGATCGCGTCTTGTGTGCGGGTGCGTGGCGCGTCGTGAGTATGCTGACGACCAGACACATCCCGCACTTCCTGGTTATTAACCAGTTCGCCTAAAGTAAtgttgt                                                                                                           <  1:2964185/148‑1 (MQ=255)
                                          ccTGATTCTGAATTCTTTTTATTAAGCGCGTAACTTAACGTCGATCGCGTCTTGTGTGCGGGTGCGTGGCGCGTCGTGAGTATGCTGACGACCAGACACATCCCGCACTTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGttgttg                                                                                                         >  2:1933015/1‑148 (MQ=255)
                                          ccTGATTCTGAATTCTTTTTATTAAGCGCGTAACTTAACGTCGATCGCGTCTTGTGTGCGGGTGCGTGGCGCGTCGTGAGTATGCTGACGACCAGACACATCCAGCACTTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGttgttg                                                                                                         >  2:121151/1‑148 (MQ=255)
                                            tGATTCTGAATTCTTTTTATTAAGCGCGTAACTTAACGTCGATCGCGTCTTGTGTGCGGGTGCGTGGCGCGTCGTGAGTATGCTGACGACCAGACACATCCCGCACTTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTgag                                                                                                       >  1:3209026/1‑148 (MQ=255)
                                                tCTGAATTCTTTTTATTAAGCGCGTAACTTAACGTCGATCGCGTCTTGTGTGCGGGTGCGTGGCGCGTCGTGAGTATGCTGACGACCAGACACATCCCGCACTTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGaaaa                                                                                                   >  1:1347353/1‑148 (MQ=255)
                                                    aaTTCTTTTTATTAAGCGCGTAACTTAACGTCGATCGCGTCTTGTGTGCGGGTGCGTGGCGCGTCGTGAGTATGCTGACGACCAGACACATCCAGCACTTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCgg                                                                                               >  2:2287795/1‑148 (MQ=255)
                                                      ttCTTTTTATTAAGCGCGTAACTTAACGTCGATCGCGTCTTGTGTGCGGGTGCGTGGCGCGTCGTGAGTATGCTGACGACCAGACACATCCAGCACTTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTg                                                                                             >  2:3314460/1‑148 (MQ=255)
                                                                 aaGCGCGTAACTTACCGTCGATCGCGTCTTGTGTGCGGGTGCGTGGCGCGTCGTGAGTATGCTGACGACCAGACACATCCAGCACTTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCt                                                                                  <  1:928010/148‑1 (MQ=255)
                                                                 aaGCGCGTAACTTAACGTCGATCGCGTCTTGTGTGCGGGTGCGTGGCGCGTCGTGAGTATGCTGACGACCAGACACATCCCGCACTTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCt                                                                                  >  1:1835458/1‑148 (MQ=255)
                                                                 aaGCGCGTAACTTAACGTCGATCGCGTCTTGTGTGCGGGTGCGTGGCGCGTCGTGAGTATGCTGACGACCAGACACATCCCGCACTTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCt                                                                                  >  1:54616/1‑148 (MQ=255)
                                                                 aaGCGCGTAACTTAACGTCGATCGCGTCTTGTGTGCGGGTGCGTGGCGCGTCGTGAGTATGCTGACGACCAGACACATCCAGCACTTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCt                                                                                  <  1:2362122/148‑1 (MQ=255)
                                                                  aGCGCGTAACTTAACGTCGATCGCGTCTTGTGTGCGGGTGCGTGGCGCGTCGTGAGTATGCTGACGACCAGACACATCCCGCACTTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTc                                                                                 >  1:2112682/1‑148 (MQ=255)
                                                                    cgcgTAACTTAACGTCGATCGCGTCTTGTGTGCGGGTGCGTGGCGCGTCGTGAGTATGCTGACGACCAGACACATCCCGCACTTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCaa                                                                               >  1:872570/1‑148 (MQ=255)
                                                                       gTAACTTAACGTCGATCGCGTCTTGTGTGCGGGTGCGTGGCGCGTCGTGAGTATGCTGACGACCAGACACATCCCGCACTTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATc                                                                            <  1:1933015/148‑1 (MQ=255)
                                                                          aCTTAACGTCGATCGCGTCTTGTGTGCGGGTGCGTGGCGCGTCGTGAGTATGCTGACGACCAGACACATCCCGCACTTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACg                                                                         <  2:1003075/148‑1 (MQ=255)
                                                                          aCTTAACGTCGATCGCGTCTTGTGTGCGGGTGCGTGGCGCGTCGTGAGTATGCTGACGACCAGACACATCCCGCACTTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACg                                                                         <  2:1005514/148‑1 (MQ=255)
                                                                           cTTAACGTCGATCGCGTCTTGTGTGCGGGTGCGTGGCGCGTCGTGAGTATGCTGACGACCAGACACATCCCGCACTTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGc                                                                        <  2:802361/148‑1 (MQ=255)
                                                                           cTTAACGTCGATCGCGTCTTGTGTGCGGGTGCGTGGCGCGTCGTGAGTATGCTGACGACCAGACACATCCCGCACTTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGc                                                                        >  2:2895775/1‑148 (MQ=255)
                                                                             tAACGTCGATCGCGTCTTGTGTGCGGGTGCGTGGCGCGTCGTGAGTATGCTGACGACCAGACACATCCCGCACTTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCca                                                                      >  2:736984/1‑148 (MQ=255)
                                                                                cGTCGATCGCGTCTTGTGTGCGGGTGCGTGGCGCGTCGTGAGTATGCTGACGACCAGACACATCCAGCACTTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAg                                                                   >  1:3166726/1‑148 (MQ=255)
                                                                                cGTCGATCGCGTCTTGTGTGCCGGTGCGTGGCGCGTCGTGAGTATGCTGACGACCAGACACATACCGCACATACTGGTTATTAACCAGTTCGCCTAAACTAATGTTGTTTAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAg                                                                   >  2:407499/1‑148 (MQ=255)
                                                                                  tCGATCGCGTCTTGTGTGCGGGTGCGTGGCGCGTCGTGAGTATGCTGACGACCAGACACATCCCGCACTTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAgcg                                                                 >  1:1758459/1‑148 (MQ=255)
                                                                                  tCGATCGCGTCTTGTGTGCGGGTGCGTGGCGCGTCGTGAGTATGCTGACGACCAGACACATCCCGCACTTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAgcg                                                                 >  1:3495830/1‑148 (MQ=255)
                                                                                  tCGATCGCGTCTTGTGTGCGGGTGCGTGGCGCGTCGTGAGTATGCTGACGACCAGACACATCCCGCACTTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAgcg                                                                 >  1:1430287/1‑148 (MQ=255)
                                                                                  tCGATCGCGTCTTGTGTGCGGGTGCGTGGCGCGTCGTGAGTATGCTGACGACCAGACACATCCAGCACTTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAgcg                                                                 >  1:2545067/1‑148 (MQ=255)
                                                                                  tCGATCGCGTCTTGTGTGCGGGTGCGTGGCGCGTCGTGAGTATGCTGACGACCAGACACATCCAGCACTTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAgcg                                                                 >  2:561601/1‑148 (MQ=255)
                                                                                   cGATCGCGTCTTGTGTGCGGGTGCGTGGCGCGTCGTGAGTATGCTGACGACCAGACACATCCCGCACTTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAgcgc                                                                >  2:3445856/1‑148 (MQ=255)
                                                                                      tCGCGTCTTGTGTGCGGGTGCGTGGCGCGTCGTGAGTATGCTGACGACCAGACACATCCCGCACTTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTg                                                             >  2:1632856/1‑148 (MQ=255)
                                                                                      tCGCGTCTTGTGTGCGGGTGCGTGGCGCGTCGTGAGTATGCTGACGACCAGACACATCCCGCACTTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACAGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTg                                                             >  2:3152016/1‑148 (MQ=255)
                                                                                      tCGCGTCTTGTGTGCGGGTGCGTGGCGCGTCGTGAGTATGCTGACGACCAGACACATCCAGCACTTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTg                                                             >  2:2421084/1‑148 (MQ=255)
                                                                                         cgTCTTGTGTGCGGGTGCGTGGCGCGTCGTGAGTATGCTGACGACCAGACACATCCAGCACTTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGt                                                          <  1:1645219/148‑1 (MQ=255)
                                                                                         cgTCTTGTGTGCGGGTGCGTGGCGCGTCGTGAGTATGCTGACGACCAGACACATCCAGCACTTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGt                                                          <  1:2481821/148‑1 (MQ=255)
                                                                                            cTTGTGTGCGGGTGCGTGGCGCGTCGTGAGTATGCTGACGACCAGACACATCCCGCACTTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAg                                                       >  2:2685337/1‑148 (MQ=255)
                                                                                            cTTGTGTGCGGGTGCGTGGCGCGTCGTGAGTATGCTGACGACCAGACACATCCCGCACTTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAg                                                       >  1:3327086/1‑148 (MQ=255)
                                                                                            cTTGTGTGCGGGTGCGTGGCGCGTCGTGAGTATGCTGACGACCAGACACATCCCGCACTTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAg                                                       >  2:2084774/1‑148 (MQ=255)
                                                                                            cTTGTGTGCGGGTGCGTGGCGCGTCGTGAGTATGCTGACGACCAGACACATCCCGCACTTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAg                                                       >  1:643862/1‑148 (MQ=255)
                                                                                            cTTGTGTGCGGGTGCGTGGCGCGTCGTGAGTATGCTGACGACCAGACACATCCAGCACTTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAg                                                       >  2:2143478/1‑148 (MQ=255)
                                                                                              tgtgtgCGGGTGCGTGGCGCGTCGTGAGTATGCTGACGACCAGACACATCCCGCACTTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGc                                                     >  2:1466584/1‑148 (MQ=255)
                                                                                               gtgtgCGGGTGCGTGGCGCGTCGTGAGTATGCTGACGACCAGACACATCCCGCACTTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGCa                                                    <  1:3147244/148‑1 (MQ=255)
                                                                                                tgtgCGGGTGCGTGGCGCGTCGTGAGTATGCTGACGACCAGACACATCCCGCACTTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGCAt                                                   >  1:614791/1‑148 (MQ=255)
                                                                                                tgtgCGGGTGCGTGGCGCGTCGTGAGTATGCTGACGACCAGACACATCCCGCACTTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGCAt                                                   <  2:44081/148‑1 (MQ=255)
                                                                                                  tgCGGGTGCGTGGCGCGTCGTGAGTATGCTGACGACCAGACACATCCCGCACTTCCTGGTTATTCACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGCAttt                                                 <  2:2739285/148‑1 (MQ=255)
                                                                                                  tgCGGGTGCGTGGCGCGTCGTGAGTATGCTGACGACCAGACACATCCCGCACTTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGCAttt                                                 <  1:2412167/148‑1 (MQ=255)
                                                                                                            tGGCGCGTCGTGAGTATGCTGACGACCAGACACATCCCGCACTTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGccc                                       <  1:1058829/148‑1 (MQ=255)
                                                                                                                 cgtcgtGAGTATGCTGACGACCAGACACATCCCGCACTTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCa                                  >  1:1520410/1‑148 (MQ=255)
                                                                                                                 cgtcgtGAGTATGCTGACGACCAGACACATCCCGCACTTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCa                                  >  1:66721/1‑148 (MQ=255)
                                                                                                                   tcgtGAGTATGCTGACGACCAGACACATCCAGCACTTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAgc                                <  1:123882/148‑1 (MQ=255)
                                                                                                                      tGAGTATGCTGACGACCAGACACATCCCGCACTTCCTGGTTATTAACCAGTTCGCCTAATGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCGCAAATCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAgccgc                             >  2:1547928/1‑148 (MQ=255)
                                                                                                                       gAGTATGCTGACGACCAGACACATCCAGCACTTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAgccgcc                            <  2:2021871/148‑1 (MQ=255)
                                                                                                                          tATGCTGACGACCAGACACATCCAGCACTTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCttt                         <  2:941377/148‑1 (MQ=255)
                                                                                                                          tATGCTGACGACCAGACACATCCAGCACTTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCttg                         <  2:941174/148‑2 (MQ=255)
                                                                                                                             gCTGACGACCAGACACATCCCGCACTTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTAcc                      <  2:1347353/148‑1 (MQ=255)
                                                                                                                             gCTGACGACCAGACACATCCAGCACTTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTAcc                      <  2:276805/148‑1 (MQ=255)
                                                                                                                              cTGACGACCAGACACATCCCGCACTTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTAccc                     <  1:2937399/148‑1 (MQ=255)
                                                                                                                               tGACGACCAGACACATCCCGCACTTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCt                    >  2:3102229/1‑148 (MQ=255)
                                                                                                                                 acgacCAGACACATCCCGCACTTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGCCTTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGa                  >  1:2591581/1‑148 (MQ=255)
                                                                                                                                 acgacCAGACACATCCAGCACTTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGa                  >  2:114895/1‑148 (MQ=255)
                                                                                                                                   gacCAGACACATCCCGCACTTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACa                <  1:2872005/148‑1 (MQ=255)
                                                                                                                                   gacCAGACACATCCAGCACTTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACa                <  1:3268804/148‑1 (MQ=255)
                                                                                                                                     ccAGACACATCCCGCACTTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAAc              >  2:631040/1‑148 (MQ=255)
                                                                                                                                     ccAGACACATCCAGCACTTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAAc              <  2:1597244/148‑1 (MQ=255)
                                                                                                                                      cAGACACATCCCGCACTTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACg             <  2:2502327/148‑1 (MQ=255)
                                                                                                                                      cAGACACATCCAGCACTTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACg             <  1:561601/148‑1 (MQ=255)
                                                                                                                                       aGACACATCCCGCACTTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAATTCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACgg            >  2:1783030/1‑148 (MQ=255)
                                                                                                                                          cacaTCCAGCACTTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCAGACAACGGGTg         >  1:2288682/1‑148 (MQ=255)
                                                                                                                                           acaTCCCGCACTTCCTGGGTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTgg        <  1:3201427/148‑1 (MQ=255)
                                                                                                                                              tCCCGCACTTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCAt     >  2:3323575/1‑148 (MQ=255)
                                                                                                                                              tCCCGCACTTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCAt     <  2:3270274/148‑1 (MQ=255)
                                                                                                                                              tCCCGCACTTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCAt     <  2:270926/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
GCTGGATGTACGACCGTGTTTACGGAGTATTTAGCACTCCGGCCTGATTCTGAATTCTTTTTATTAAGCGCGTAACTTAACGTCGATCGCGTCTTGTGTGCGGGTGCGTGGCGCGTCGTGAGTATGCTGACGACCAGACACATCCAGCACTTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTA  >  CP009273/2654939‑2655231

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: