Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,347,644 G→A 8.5% intergenic (‑235/‑440) yhcC ← / → gltB putative Fe‑S oxidoreductase, Radical SAM superfamily protein/glutamate synthase, large subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,347,6440GA 8.5% 831.0 / 42.9 248intergenic (‑235/‑440)yhcC/gltBputative Fe‑S oxidoreductase, Radical SAM superfamily protein/glutamate synthase, large subunit
Reads supporting (aligned to +/- strand):  ref base G (112/115);  new base A (11/10);  total (123/125)
Fisher's exact test for biased strand distribution p-value = 8.23e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.28e-01

Only 100 of 252 total aligned reads displayed.

GCTTTTATCAACTGCATAATCAATCAAAATTACCGAAATTTCATGCATAATCACATAAATCACTTTTGCTTATCTTGTGTCAGATTTTTTTATCTCCTGATGGATTTTAGGCAAAAACAGTAGCATGAAACGTCATTACCAATTAAGGCAGTATAAAATGCTGGTTTTGTCGTCAGTTCAAGGCAGGATAAGGGTTAACACACCTTTATGACAGTCAGGAATTGACTGTTTCTCTAACGACTTCCCTTTTAGCCTTAAAGATAAAATCCATTTTAATTTCAGTCATTTAATAAAG  >  CP009273/3347497‑3347791
                                                                                                                                                   |                                                                                                                                                   
 cTTTTATCAACTGCATAATCAATCAAAATTACCGAAATTTCATGCATAATCACATAAATCACTTTTGCTTATCTTGTGTCAGATTTTTTTATCTCCTGATGGATTTTAGGCAAAAACAGTAGCATGAAACGTCATTACCAATTAAGGc                                                                                                                                                    <  2:856607/148‑1 (MQ=255)
 cTTTTATCAACTGCATAATCAATCAAAATTACCGAAATTTCATGCATAATCACATAAATCACTTTTGCTTATCTTGTGTCAGATTTTTTTATCTCCTGATGGATTTTAGGCAAAAACAGTAGCATGAAACGTCATTACCAATTAAGGc                                                                                                                                                    <  1:3510548/148‑1 (MQ=255)
    ttATCAACTGCATAATCAATCAAAATTACCGAAATTTCATGCATAATCACATAAATCACTTTTGCTTATCTTGTGTCAGATTTTTTTATCTCCTGATGGATTTTAGGCAAAAACAGTAGCATGAAACGTCATTACCAATTAAGGCAGt                                                                                                                                                 <  2:3098448/148‑1 (MQ=255)
        cAACTGCATAATCAATCAAAATTACCGAAATTTCATGCATAATCACATAAATCACTTTTGCTTATCTTGTGTCAGATTTTTTTATCTCCTGATGGATTTTAGGCAAAAACAGTAGCATGAAACGTCATTACCAATTAAGGCAGTATaa                                                                                                                                             >  1:1434775/1‑148 (MQ=255)
            tGCATAATCAATCAAAATTACCGAAATTTCATGCATAATCACATAAATCACTTTTGCTTATCTTGTGTCAGATTTTTTTATCTCCTGATGGATTTTAGGCAAAAACAGTAGCATGAAACGTCATTACCAATTAAGGCAGTATAAAATg                                                                                                                                         <  2:103800/148‑1 (MQ=255)
             gCATAATCAATCAAAATTACCGAAATTTCATGCATAATCACATAAATCACTTTTGCTTATCTTGTGTCAGATTTTTTTATCTCCTGATGGATTTTAGGCAAAAACAGTAGCATGAAACGTCATTACCAATTAAGGCAGTATAAAATGc                                                                                                                                        >  2:2040453/1‑148 (MQ=255)
             gCATAATCAATCAAAATTACCGAAATTTCATGCATAATCACATAAATCACTTTTGCTTATCTTGTGTCAGATTTTTTTATCTCCTGATGGATTTTAGGCAAAAACAGTAGCATGAAACGTCATTACCAATTAAGGCAGTATAAAATGc                                                                                                                                        >  1:229502/1‑148 (MQ=255)
              cATAATCAATCAAAATTACCGAAATTTCATGCATAATCACATAAATCACTTTTGCTTATCTTGTGTCAGATTTTTTTATCTCCTGATGGATTTTAGGCAAAAACAGTAGCATGAAACGTCATTACCAATTAAGGCAGTATAAAATGCt                                                                                                                                       <  1:2522765/148‑1 (MQ=255)
               aTAATCAATCAAAATTACCGAAATTTCATGCATAATCACATAAATCACTTTTGCTTATCTTGTGTCAGATTTTTTTATCTCCTGATGGATTTTAGGCAAAAACAGTAGCATGAAACGTCATTACCAATTAAGGCAGTATAAAATGCTg                                                                                                                                      >  2:2959837/1‑148 (MQ=255)
               aTAATCAATCAAAATTACCGAAATTTCATGCATAATCACATAAATCACTTTTGCTTATCTTGTGTCAGATTTTTTTATCTCCTGATGGATTTTAGGCAAAAACAGTAGCATGAAACGTCATTACCAATTAAGGCAGTATAAAATGCTg                                                                                                                                      <  1:545173/148‑1 (MQ=255)
                tattcaatataAATTACCGATTTTTAATGGATATTCACATAAATCACTTTTGCTTATGTTTTGTCAGATTTTTTTTTCTCCTGTTTGATTTTAGGCAAATAGAGTATCATGAATCGTCATTACCTATTAAGGCAGTATAAAATGCTgg                                                                                                                                     <  1:15177/145‑1 (MQ=255)
                tAATCAATCAAAATTACCGAAATTTCATGCATAATCACATAAATCACTTTTGCTTATCTTGTGTCAGATTTTTTTATCTCCTGATGGATTTTAGGCAAAAACAGTAGCATGAAACGTCATTACCAATTAAGGCAGTATAAAATGCTgg                                                                                                                                     <  2:1382051/148‑1 (MQ=255)
                tAATCAATCAAAATTACCGAAATTTCATGCATAATCACATAAATCACTTTTGCTTATCTTGTGTCAGATTTTTTTATCTCCTGATGGATTTTAGGCAAAAACAGTAGCATGAAACGTCATTACCAATTAAGGCAGTATAAAATGCTgg                                                                                                                                     <  2:1093345/148‑1 (MQ=255)
                 aatcaatcAAAATTACCGAAATTTCATGCATAATCACATAAATCACTTTTGCTTATCTTGTGTCAGATTTTTTTATCTCCTGATGGATTTTAGGCAAAAACAGTAGCATGAAACGTCATTACCAATTAAGGCAGTATAAAATGCTGGt                                                                                                                                    >  1:2786280/1‑148 (MQ=255)
                 aatcaatcAAAATTACCGAAATTTCATGCATAATCACATAAATCACTTTTGCTTATCTTGTGTCAGATTTTTTTATCTCCTGATGGATTTTAGGCAAAAACAGTAGCATGAAACGTCATTACCAATTAAGGCAGTATAAAATGCTGGt                                                                                                                                    >  2:10657/1‑148 (MQ=255)
                  atcaatcaAAATTACCGAAATTTCATGCATAATCACATAAATCACTTTTGCTTATCTTGTGTCAGATTTTTTTATCTCCTGATGGATTTTAGGCAAAAACAGTAGCATGAAACGTCATTACCAATTAATGCAGTATAAAATGCTGGtt                                                                                                                                   >  1:1549600/1‑148 (MQ=255)
                  atcaatcaAAATTACCGAAATTTCATGCATAATCACATAAATCACTTTTGCTTATCTTGTGTCAGATTTTTTTATCTCCTGATGGATTTTAGGCAAAAACAGTAGCATGAAACGTCATTACCAATTAAGGCAGTATAAAATGCTGGtt                                                                                                                                   <  1:1478365/148‑1 (MQ=255)
                   tcaatcaaAATTACCGAAATTTCATGCATAATCACATAAATCACTTTTGCTTATCTTGTGTCAGATTTTTTTATCTCCTGATGGATTTTAGGCAAAAACAGTAGCATGAAACGTCATTACCAATTAAGGCAGTATAAAATGCTGGttt                                                                                                                                  >  2:158973/1‑148 (MQ=255)
                   tcaatcaaAATTACCGAAATTTCATGCATAATCACATAAATCACTTTTGCTTATCTTGTGTCAGATTTTTTTATCTCCTGATGGATTTTAGGCAAAAACAGTAGCATGAAACGTCATTACCAATTAAGGCAGTATAAAATGCTGGttt                                                                                                                                  >  2:3560811/1‑148 (MQ=255)
                   tcaatcaaAATTACCGAAATTTCATGCATAATCACATAAATCACTTTTGCTTATCTTGTGTCAGATTTTTTTATCTCCTGATGGATTTTAGGCAAAAACAGTAGCATGAAACGTCATTACCAATTAAGACAGTATAAAATGCTGGttt                                                                                                                                  >  1:1307077/1‑148 (MQ=255)
                    caatcaaAATTACCGAAATTTCATGCATAATCACATAAATCACTTTTGCTTATCTTGTGTCAGATTTTTTTATCTCCTGATGGATTTTAGGCAAAAACAGTAGCATGAAACGTCATTACCAATTAAGGCAGTATAAAATGCTGGtttt                                                                                                                                 <  1:245315/148‑1 (MQ=255)
                      atcaaAATTACCGAAATTTCATGCATAATCACATAAATCACTTTTGCTTATCTTGTGTCAGATTTTTTTATCTCCTGATGGATTTTAGGCAAAAACAGTAGCATGAAACGTCATTACCAATTAAGGCAGTATAAAATGCTGGTTTTgt                                                                                                                               <  2:2798709/148‑1 (MQ=255)
                         aaaaTTACCGAAATTTCATGCATAATCACATAAATCACTTTTGCTTATCTTGTGTCAGATTTTTTTATCTCCTGATGGATTTTAGGCAAAAACAGTAGCATGAAACGTCATTACCAATTAAGGCAGTATAAAATGCTGGTTTTgtcgt                                                                                                                            <  1:2224845/148‑1 (MQ=255)
                             ttACCGAAATTTCATGCATAATCACATAAATCACTTTTGCTTATCTTGTGTCAGATTTTTTTATCTCCTGATGGATTTTAGGCAAAAACAGTAGCATGAAACGTCATTACCAATTAAGGCAGTATAAAATGCTGGTTTTGTCGTCAGt                                                                                                                        <  2:3197580/148‑1 (MQ=255)
                             ttACCGAAATTTCATGCATAATCACATAAATCACTTTTGCTTATCTTGTGTCAGATTTTTTTATCTCCTGATGGATTTTAGGCAAAAACAGTAGCATGAAACGTCATTACCAATTAAGGCAGTATAAAATGCTGGTTTTGTCGTCAGt                                                                                                                        <  2:2661733/148‑1 (MQ=255)
                              tACCGAAATTTCATGCATAATCACATAAATCACTTTTGCTTATCTTGTGTCAGATTTTTTTATCTCCTGATGGATTTTAGGCAAAAACAGTAGCATGAAACGTCATTACCAATTAAGGCAGTATAAAATGCTGGTTTTGTCGTCAGtt                                                                                                                       >  1:3404439/1‑148 (MQ=255)
                                   aaaTTTCATGCATAATCACATAAATCACTTTTGCTTATCTTGTGTCAGATTTTTTTATCTCCTGATGGATTTTAGGCAAAAACAGTAGCATGAAACGTCATTACCAATTAAGGCAGTATAAAATGCTGGTTTTGTCGTCAGTTCAAgg                                                                                                                  >  2:290098/1‑148 (MQ=255)
                                     aTTTCATGCATAATCACATAAATCACTTTTGCTTATCTTGTGTCAGATTTTTTTATCTCCTGATGGATTTTAGGCAAAAACAGTAGCATGAAACGTCATTACCAATTAAGGCAGTATAAAATGCTGGTTTTGTCGTCAGTTCAAGGCa                                                                                                                <  1:2343853/148‑1 (MQ=255)
                                      tttCATGCATAATCACATAAATCACTTTTGCTTATCTTGTGTCAGATTTTTTTATCTCCTGATGGATTTTAGGCAAAAACAGTAGCATGAAACGTCATTACCAATTAAGGCAGTATAAAATGCTGGTTTTGTCGTCAGTTCAAGGCAg                                                                                                               >  1:747283/1‑148 (MQ=255)
                                       ttCATGCATAATCACATAAATCACTTTTGCTTATCTTGTGTCAGATTTTTTTATCTCCTGATGGATTTTAGGCAAAAACAGTAGCATGAAACGTCATTACCAATTAAGGCAGTATAAAATGCTGGTTTTTTCGTCAGTTCAAGGCAgg                                                                                                              >  2:1105591/1‑148 (MQ=255)
                                       ttCATGCATAATCACATAAATCACTTTTGCTTATCTTGTGTCAGATTTTTTTATCTCCTGATGGATTTTAGGCAAAAACAGTAGCATGAAACGTCATTACCAATTAAGGCAGTATAAAATGCTGGTTTTGTCGTCAGTTCAAGGCAgg                                                                                                              <  1:158973/148‑1 (MQ=255)
                                        tCATGCATAATCACATAAATCACTTTTGCTTATCTTGTGTCAGATTTTTTTATCTCCTGATGGATTTTAGGCAAAAACAGTAGCATGAAACGTCATTACCAATTAAGGCAGTATAAAATGCTGGTTTTGTCGTCAGTTCAAGGCAGGa                                                                                                             <  2:2025829/148‑1 (MQ=255)
                                          aTGCATAATCACATAAATCACTTTTGCTTATCTTGTGTCAGATTTTTTTATCTCCTGATGGATTTTAGGCAAAAACAGTAGCATGAAACGTCATTACCAATTAAGGCAGTATAAAATGCTGGTTTTGTCGTCAGTTCAAGGCAGGATa                                                                                                           <  2:1962632/148‑1 (MQ=255)
                                          aTGCATAATCACATAAATCACTTTTGCTTATCTTGTGTCAGATTTTTTTATCTCCTGATGGATTTTAGGCAAAAACAGTAGCATGAAACGTCATTACCAATTAAGGCAGTATAAAATGCTGGTTTTGTCGTCAGTTCAAGGCAGGATa                                                                                                           >  2:670407/1‑148 (MQ=255)
                                           tGCATAATCACATAAATCACTTTTGCTTATCTTGTGTCAGATTTTTTTATCTCCTGATGGATTTTAGGCAAAAACAGTAGCATGAAACGTCATTACCAATTAAGGCAGAATAAAATGCTGGTTTTGTCGTCAGTTCAAGGAAGGATaa                                                                                                          >  2:2774153/1‑148 (MQ=255)
                                             cATAATCACATAAATCACTTTTGCTTATCTTGTGTCAGATTTTTTTATCTCCTGATGGATTTTAGGCAAAAACAGTAGCATGAAACGTCATTACCAATTAAGGCAGTATAAAATGCTGGTTTTGTCGTCAGTTCAAGGCAGGATAAgg                                                                                                        <  1:2995148/148‑1 (MQ=255)
                                             cATAATCACATAAATCACTTTTGCTTATCTTGTGTCAGATTTTTTTATCTCCTGATGGATTTTAGGCAAAAACAGTAGCATGAAACGTCATTACCAATTAAGACAGTATAAAATGCTGGTTTTGTCGTCAGTTCAAGGCAGGATAAgg                                                                                                        <  1:2641268/148‑1 (MQ=255)
                                                 aTCACAGAAATCACTTGTGCTTATCTGGTGTCAGATTTTTTTATCTCCTGATGGATTTTAGGCAAAAACAGTAGCATGAAACGTCATTACCAATTAAGGCAGTATAAAATGCTGGTTTTGTCGTCAGTTCAAGGCAGGATAAGGGTTa                                                                                                    <  2:3478696/148‑1 (MQ=255)
                                                     caTAAATCACTTTTGCTTATCTTGTGCCAGATTTTTTTATCTCCTGATGGATTTTAGGCAAAAACAGTAGCATGAAACGTCATTACCAATTAAGGCAGTATAAAATGCTGGTTTTGTCGTCAGTTCAAGGCAGGATAAGGGTTAacac                                                                                                >  2:1532795/1‑148 (MQ=255)
                                                        aaaTCACTTTTGCTTATCTTGTGTCAGATTTTTTTATCTCCTGATGGATTTTAGGCAAAAACAGTAGCATGAAACGTCATTACCAATTAAGGCAGTATAAAATGCTGGTTTTGTCGTCAGTTCAAGGCAGGATAAGGGTTAACACAcc                                                                                             >  1:1725266/1‑148 (MQ=255)
                                                         aaTCACTTTTGCTTATCTTGTGTCAGATTTTTTTATCTCCTGATGGATTTTAGGCAAAAACAGTAGCATGAAACGTCATTACCAATTAAGGCAGTATAAAATGCTGGTTTTGTCGTCAGTTCAAGGCAGGATAAGGGTTAACACACCt                                                                                            <  2:1732793/148‑1 (MQ=255)
                                                         aaTCACTTTTGCTTATCTTGTGTCAGATTTTTTTATCTCCTGATGGATTTTAGGCAAAAACAGTAGCATGAAACGTCATTACCAATTAAGACAGTATAAAATGCTGGTTTTGTCGTCAGTTCAAGGCAGGATAAGGGTTAACACACCt                                                                                            >  1:2115908/1‑148 (MQ=255)
                                                          aTCACTTTTGCTTATCTTGTGTCAGATTTTTTTATCTCCTGATGGATTTTAGGCAAAAACAGTAGCATGAAACGTCATTACCAATTAAGGCAGTATAAAATGCTGGTTTTGTCGTCAGTTCAAGGCAGGATAAGGGTTAACACACCtt                                                                                           <  2:579921/148‑1 (MQ=255)
                                                                tttGCTTATCTTGTGTCAGATTTTTTTATCTCCTGATGGATTTTAGGCAAAAACAGTAGCATGAAACGTCATTACCAATTAAGGCAGTATAAAATGCTGGTTTTGTCGTCAGTTCAAGGCAGGATAAGGGTTAACACACCTTTATGAc                                                                                     <  2:3177666/148‑1 (MQ=255)
                                                                  tGCTTATCTTGTGTCAGATTTTTTTATCTCCTGATGGATTTTAGGCAAAAACAGTAGCATGAAACGTCATTACCAATTAAGGCAGTATAAAATGCTGGTTTTGTCGTCAGTTCAAGGCAGGATAAGGGTTAACACACCTTTATGACAg                                                                                   <  2:2469490/148‑1 (MQ=255)
                                                                  tGCTTATCTTGTGTCAGATTTTTTTATCTCCTGATGGATTTTAGGCAAAAACAGTAGCATGAAACGTCATTACCAATTAAGACAGTATAAAATGCTGGTTTTGTCGTCAGTTCAAGGCAGGATAAGGGTTAACACACCTTTATGACAg                                                                                   >  2:441505/1‑148 (MQ=255)
                                                                    cTTATCTTGTGTCAGATTTTTTTATCTCCTGATGGATTTTAGGCAAAAACAGTAGCATGAAACGTCATTACCAATTAAGGCAGTATAAAATGCTGGTTTTGTCGTCAGTTCAAGGCAGGATAAGGGTTAACACACCTTTATGACAGTc                                                                                 >  2:1512851/1‑148 (MQ=255)
                                                                        tCTTGTGTCAGATTTTTTTATCTCCTGATGGATTTTAGGCAAAAACAGTAGCATGAAACGTCATTACCAATTAAGGCAGTATAAAATGCTGGTTTTGTCGTCAGTTCAAGGCAGGATAAGGGTTAACACACCTTTATGACAGTCAGGa                                                                             >  1:449673/1‑148 (MQ=255)
                                                                         cTTGTGTCAGATTTTTTTATCTCCTGATGGATTTTAGGCAAAAACAGTAGCATGAAACGTCATTACCAATTAAGGCAGTATAAAATGCTGGTTTTGTCGTCAGTTCAAGGCAGGATAAGGGTTAACACACCTTTATGACAGTCAGGaa                                                                            <  1:2567935/148‑1 (MQ=255)
                                                                          ttGTGTCAGATTTTTTTATCTCCTGATGGATTTTAGGCAAAAACAGTAGCATGAAACGTCATTACCAATTAAGGCAGTATAAAATGCTGGTTTTGTCGTCAGTTCAAGGTAGGATAAGGGTTAACACACCTTTATGACAGTCAGGAAt                                                                           >  1:13149/1‑148 (MQ=255)
                                                                          ttGTGTCAGATTTTTTTATCTCCTGATGGATTTTAGGCAAAAACAGTAGCATGAAACGTCATTACCAATTAAGGCAGTATAAAATGCTGGTTTTGTCGTCAGTTCAAGGCAGGATAAGGGTTAACACACCTTTATGACAGTCAGGAAt                                                                           >  2:344377/1‑148 (MQ=255)
                                                                              gtCAGATTTTTTTATCTCCTGATGGATTTTAGGCAAAAACAGTAGCATGAAACGTCATTACCAATTAAGGCAGTATAAAATGCTGGTTTTGTCGTCAGTTCAAGGCAGGATAAGGGTTAACACACCTTTATGACAGTCAGGAATTGAc                                                                       <  1:891411/148‑1 (MQ=255)
                                                                              gtCAGATTTTTTTATCTCCTGATGGATTTTAGGCAAAAACAGTAGCATGAAACGTCATTACCAATTAAGGCAGTATAAAATGCTGGTTTTGTCGTCAGTTCAAGGCAGGATAAGGGTTAACACACCTTTATGACAGTCAGGAATTGAc                                                                       >  2:1226310/1‑148 (MQ=255)
                                                                                  gATTTTTTTATCTCCTGATGGATTTTAGGCAAAAACAGTAGCATGAAACGTCATTACCAATTAAGGCAGTATAAAATGCTGGTTTTGTCGTCAGTTCAAGGCAGGATAAGGGTTAACACACCTTTATGACAGTCAGGAATTGACTGtt                                                                   >  1:2904890/1‑148 (MQ=255)
                                                                                  gATTTTTTTATCTCCTGATGGATTTTAGGCAAAAACAGTAGCATGAAACGTCATTACCAATTAAGGCAGTATAAAATGCTGGTTTTGTCGTCAGTTCAAGGCAGGATAAGGGTTAACACACCTTTATGACAGTCAGGAATTGACTGtt                                                                   >  1:361175/1‑148 (MQ=255)
                                                                                     ttttttATCTCCTGATGGATTTTAGGCAAAAACAGTAGCATGAAACGTCATTACCAATTAAGGCAGTATAAAATGCTGGTTTTGTCGTCAGTTCAAGGCAGGATAAGGGTTAACACACCTTTATGACAGTCAGGAATTGACTGTTtct                                                                <  2:3404439/148‑1 (MQ=255)
                                                                                     ttttttATCTCCTGATGGATTTTAGGCAAAAACAGTAGCATGAAACGTCATTACCAATTAAGGCAGTATAAAATGCTGGTTTTGTCGTCAGTTCAAGGCAGGATAAGGGTTAACACACCTTTATGACAGTCAGGAATTGACTGTTtct                                                                <  1:2836051/148‑1 (MQ=255)
                                                                                      tttttATCTCCTGATGGATTTTAGGCAAAAACAGTAGCATGAAACGTCATTACCAATTAAGGCAGTATAAAATGCTGGTTTTGTCTTCAGTTCAAGGCAGGATAAGGGTTAACACACCTTTATGACAGTCAGGAATTGACTGTTtctc                                                               <  1:2959837/148‑1 (MQ=255)
                                                                                        tttATCTCCTGATGGATTTTAGGCAAAAACAGTAGCATGAAACGTCATTACCAATTAAGGCAGTATAAAATGCTGGTTTTGTCGTCAGTTCAAGGCAGGATAAGGGTTAACACACCTTTATGACAGTCAGGAATTGACTGTTTCTCTa                                                             <  2:2252993/148‑1 (MQ=255)
                                                                                           aTCTCCTGATGGATTTTAGGCAAAAACAGTAGCATGAAACGTCATTACCAATTAAGGCAGTATAAAATGCTGGTTTTGTCGTCAGTTCAAGGCAGGATAAGGGTTAACACACCTTTATGACAGTCAGGAATTGACTGTTTCTCTAACg                                                          >  2:170872/1‑148 (MQ=255)
                                                                                             ctcCTGATGGATTTTAGGCAAAAACAGTAGCATGAAACGTCATTACCAATTAAGGCAGTATAAAATGCTGGTTTTGTCGTCAGTTCAAGGCAGGATAAGGGTTAACACACCTTTATGACAGTCAGGACTTGACTGTTTCTCTAACGAc                                                        >  2:2180708/1‑148 (MQ=255)
                                                                                             ctcCTGATGGATTTTAGGCAAAAACAGTAGCATGAAACGTCATTACCAATTAAGGCAGTATAAAATGCTGGTTTTGTCGTCAGTTCAAGGCAGGATAAGGGTTAACACACCTTTATGACAGTCAGGAATTGACTGTTTCTCTAACGAc                                                        >  1:424429/1‑148 (MQ=255)
                                                                                                 tGATGGATTTTAGGCAAAAACAGTAGCACGAAACGTCATTACCAATTAAGGCAGTATAAAATGCTGGTTTTGTCGTCAGTTCAAGGCAGGATAAGGGTTAACACACCTTTATGACAGTCAGGAATTGACTGTTTCTCTAACGACTTcc                                                    >  1:1181041/1‑148 (MQ=255)
                                                                                                   aTGGATTTTTGGCAAAAACAGTAGCATGAAACGTCATTACCAATTAAGGCAGTATAAAATGCTGGTTTTGTCGTCAGTTCAAGGCAGGATAAGGGTTAACACACCTTTATGACAGTCAGGACTTGACTGTTTCTCTAACGACTTCCCt                                                  >  1:2419419/1‑148 (MQ=255)
                                                                                                   aTGGATTTTAGGCAAAAACAGTAGCATGAAACGTCATTACCAATTAAGGCAGTATAAAATGCTGGTTTTGTCGTCAGTTCAAGGCAGGATAAGGGTTAACACACCTTTATGACAGTCAGGAATTGACTGTTTCTCTAACGACTTCCCt                                                  >  1:801947/1‑148 (MQ=255)
                                                                                                    tGGATTTTAGGCAAAAACAGTAGCATGAAACGTCATTACCAATTAAGGCAGTATAAAATGCTGGTTTTGTCGTCAGTTCAAGGCAGGATAAGGGTTAACACACCTTTATGACAGTCAGGAATTGACTGTTTCTCTAACGACTTCCCtt                                                 >  1:1999317/1‑148 (MQ=255)
                                                                                                      gATTTTAGGCAAAAACAGTAGCATGAAACGTCATTACCAATTAAGGCAGTATAAAATGCTGGTTTTGTCGTCAGTTCAAGGCAGGATAAGGGTTAACACACCTTTATGACAGTCAGGAATTGACTGTTTCTCTAACGACTTCCCtttt                                               >  1:177750/1‑148 (MQ=255)
                                                                                                       aTTTTAGGCAAAAACAGTAGCATGAAACGTCATTACCAATTAAGGCAGTATAAAATGCTGGTTTTGTCGTCAGTTCAAGGCAGGATAAGGGTTAACACACCTTTATGACAGTCAGGAATTGACTGTTTCTCTAACGACTTCCCTTTTa                                              >  2:660930/1‑148 (MQ=255)
                                                                                                        ttttAGGCAAAAACAGTAGCATGAAACGTCATTACCAATTAAGGCAGTATAAAATGCTGGTTTTGTCGTCAGTTCAAGGCAGGATAAGGGTTAACACACCTTTATTACAGTCAGGAATTGACTGTTTCTCTAACGACTTCCCTTTTAg                                             >  2:1237619/1‑148 (MQ=255)
                                                                                                        ttttAGGCAAAAACAGTAGCATGAAACGTCATTACCAATTAAGGCAGTATAAAATGCTGGTTTTGTCGTCAGTTCAAGGCAGGATAAGGGTTAACACACCTTTATGACAGTCAGGAATTGACTGTTTCTCTAACGACTTCCCTTTTAg                                             >  2:890881/1‑148 (MQ=255)
                                                                                                           tAGGCAAAAACAGTAGCATGAAACGTCATTACCAATTAAGGCAGTATAAAATGCTGGTTTTGTCGTCAGTTCAAGGCAGGATAAGGGTTAACACACCTTTATGACAGTCAGGAATTGACTGTTTCTCTAACGACTTCCCTTTTAGCCt                                          <  1:3448254/148‑1 (MQ=255)
                                                                                                             ggCAAAAACAGTAGCATGAAACGTCATTACCAATTAAGACAGTATAAAATGCTGGTTTTGTCGTCAGTTCAAGGCAGGATAAGGGTTAACACACCTTTATGACAGTCAGGAATTGACTGTTTCTCTAACGACTTCCCTTTTAGCCTTa                                        >  1:3140563/1‑148 (MQ=255)
                                                                                                              gCAAAAACAGTAGCATGAAACGTCATTACCAATTAAGGCAGTATAAAATGCTGGTTTTGTCGTCAGTTCAAGGCAGGATAAGGGTTAACACACCTTTATGACAGTCAGGAATTGACTGTTTCTCTAACGACTTCCCTTTTAGCCTTaa                                       >  2:2245859/1‑148 (MQ=255)
                                                                                                               cAAAAACAGTAGCATGAAACGTCATTACCAATTAAGGCAGTATAAAATGCTGGTTTTGTCGTCAGTTCAAGGCAGGATAAGGGTTAACACACCTTTATGACAGTCAGGAATTGACTGTTTCTCTAACGACTTCCCTTTTAGCCTTaaa                                      <  2:3489338/148‑1 (MQ=255)
                                                                                                                aaaaaCAGTAGCATGAAACGTCATTACCAATTAAGGCAGTATAAAATGCTGGTTTTGTCGTCAGTTCAAGGCAGGATAAGGGTTAACACACCTTTATGACAGTCAGGAATTGACTGTTTCTCTAACGACTTCCCTTTTAGCCTTAAAg                                     <  1:686750/148‑1 (MQ=255)
                                                                                                                 aaaaCAGTAGCATGAAACGTCATTACCAATTAAGGCAGTATAAAATGCTGGTTTTGTCGTCAGTTCAATGCAGGATAAGGGTTAACACACCTTTATGACAGTCAGGAATTGACTGTTTCTCTAACGACTTCCCTTTTAGCCTTAAAGa                                    <  2:1765060/148‑1 (MQ=255)
                                                                                                                 aaaaCAGTAGCATGAAACGTCATTACCAATTAAGGCAGTATAAAATGCTGGTTTTGTCGTCAGTTCAAGGCAGGATAAGGGTTAACACACCTTTATGACAGTCAGGAATTGACTGTTTCTCTAACGACTTCCCTTTTAGCCTTAAAGa                                    >  1:245186/1‑148 (MQ=255)
                                                                                                                  aaaCAGTAGCATGAAACGTCATTACCAATTAAGGCAGTATAAAATGCTGGTTTTGTCGTCAGTTCAAGGCAGGATAAGGGTTAACACACCTTTATGACAGTCAGGAATTGACTGTTTCTCTAACGACTTCCCTTTTAGCCTTAAAGAt                                   <  2:544207/148‑1 (MQ=255)
                                                                                                                  aaaCAGTAGCATGAAACGTCATTACCAATTAAGGCAGTATAAAATGCTGGTTTTGTCGTCAGTTCAAGGCAGGATAAGGGTTAACACACCTTTATGACAGTCAGGAATTGACTGTTTCTCTAACGACTTCCCTTTTAGCCTTAAAGAt                                   >  1:2612191/1‑148 (MQ=255)
                                                                                                                           cATGAAACGTCATTACCAATTAAGGCAGTATAAAATGCTGGTTTTGTCGTCAGTTCAAGGCAGGATAAGGGTTAACACACCTTTATGACAGTCAGGAATTGACTGTTTCTCTAACGACTTCCCTTTTAGCCTTAAAGATAAAATCCAt                          <  1:2072671/148‑1 (MQ=255)
                                                                                                                           cATGAAACGTCATTACCAATTAAGGCAGTATAAAATGCTGGTTTTGTCGTCAGTTCAAGGCAGGATAAGGGTTAACACACCTTTATGACAGTCAGGAATTGACTGTTTCTCTAACGACTTCCCTTTTAGCCTTAAAGATAAAATCCAt                          <  2:2507299/148‑1 (MQ=255)
                                                                                                                           cATGAAACGTCATTACCAATTAAGGCAGTATAAAATGCTGGTTTTGTCGTCAGTTCAAGGCAGGATAAGGGTTAACACACCTTTATGACAGTCAGGAATTGACTGTTTCTCTAACGACTTCCCTTTTAGCCTTAAAGATAAAATCCAt                          <  2:2731826/148‑1 (MQ=255)
                                                                                                                           cATGAAACGTCATTACCAATTAAGGCAGTATAAAATGCTGGTTTTGTCGTCAGTTCAAGGCAGGATAAGGGTTAACACACCTTTATGACAGTCAGGAATTGACTGTTTCTCTAACGACTTCCCTTTTAGCCTTAAAGATAAAATCCAt                          <  2:2559913/148‑1 (MQ=255)
                                                                                                                           cATGAAACGTCATTACCAATTAAGGCAGTATAAAATGCTGGTTTTGTCGTCAGTTCAAGGCAGGATAAGGGTTAACACACCTTTATGACAGTCAGGAATTGACTGTTTCTCTAACGACTTCCCTTTTAGCCTTAAAGATAAAATCCAt                          <  1:603911/148‑1 (MQ=255)
                                                                                                                                 aCGTCATTACCAATTAAGACAGTATAAAATGCTGGTTTTGTCGTCAGTTCAAGGCAGGATAAGGGTTAACACACCTTTATGACAGTCAGGAATTGACTGTTTCTCTAACGACTTCCCTTTTAGCCTTAAAGATAAAATCCATTTTAAt                    <  1:1633023/148‑1 (MQ=255)
                                                                                                                                     cATTACCAATTAAGGCAGTATAAAATGCTGGTTTTGTCGTCAGTTCAAGGCAGGATAAGGGTTAACACACCTTTATGACAGTCAGGAATTGACTGTTTCTCTAACGACTTCCCTTTTAGCCTTAAAGATAAAATCCATTTTAATTTCa                >  1:3520030/1‑148 (MQ=255)
                                                                                                                                       ttACCAATTAAGGCAGTATAAAATGCTGGTTTTGTCGTCAGTTCAAGGCAGGATAAGGGTTAACACACCTTTATGACAGTCAGGAATTGACTGTTTCTCTAACGACTTCCCTTTTAGCCTTAAAGATAAAATCCATTTTAATTTCaat              >  2:321881/1‑146 (MQ=255)
                                                                                                                                        tACCAATTAAGGCAGTATAAAATGCTGGTTTTGTCGTCAGTTCAAGGCAGGATAAGGGTTAACACACCTTTATGACAGTCAGGAATTGACTGTTTCTCTAACGACTTCCCTTTTAGCCTTAAAGATAAAATCCATTTTa                      <  2:2210156/139‑1 (MQ=255)
                                                                                                                                        tACCAATTAAGGCAGTATAAAATGCTGGTTTTGTCGTCAGTTCAAGGCAGGATAAGGGTTAACACACCTTTATGACAGTCAGGAATTGACTGTTTCTCTAACGACTTCCCTTTTAGCCTTAAAGATAAAATCCATTTTAATTTCAATc             >  2:1064429/1‑148 (MQ=255)
                                                                                                                                        tACCAATTAAGGCAGTATAAAATGCTGGTTTTGTCGTCAGTTCAAGGCAGGATAAGGGTTAACACACCTTTATGACAGTCAGGAATTGACTGTTTCTCTAACGACTTCCCTTTTAGCCTTAAAGATAAAATCCATTTTAATTTCAATc             <  2:2073052/148‑1 (MQ=255)
                                                                                                                                        tACCAATTAAGGCAGTATAAAATGCTGGTTTTGTCGTCAGTTCAAGGCAGGATAAGGGTTAACACACCTTTATGACAGTCAGGAATTGACTGTTTCTCTAACGACTTCCCTTTTAGCCTTAAAGATAAAATCCATTTTAATTTCAATc             >  2:1661466/1‑148 (MQ=255)
                                                                                                                                        tACCAATTAAGGCAGTATAAAATGCTGGTTTTGTCGTCAGTTCAAGGCAGGATAAGGGTTAACACACCTTTATGACAGTCAGGAATTGACTGTTTCTCTAACGACTTCCCTTTTAGCCTTAAAGATAAAATCCATTTTAATTTCAATc             >  1:2492925/1‑148 (MQ=255)
                                                                                                                                         aCCAATTAAGGCAGTATAAAATGCTGGTTTTGTCGTCAGTTCAAGGCAGGATAAGGGTTAACACACCTTTATGACAGTCAGGAATTGACTGTTTCTCTAACGACTTCCCTTTTAGCCTTAAAGATAAAATCCATTTTAATTTCAATCa            <  2:1653516/148‑1 (MQ=255)
                                                                                                                                           cAATTAAGGCAGTATAAAATGCTGGTTTTGTCGTCAGTTCAAGGCAGGATAAGGGTTAACACACCTTTATGACAGTCAGGAATTGACTGTTTCTCTAACGACTTCCCTTTTAGCCTTAAAGATAAAATCCATTTTAATTTCAATCAtt          >  1:663391/1‑148 (MQ=255)
                                                                                                                                            aaTTAAGGCAGTATAAAATGCTGGTTTTGTCGTCAGTTCAAGGCAGGATAAGGGTTAACACACCTTTATGACAGTCAGGAATTGACTGTTTCTCTAACGACTTCCCTTTTAGCCTTAAAGATAAAATCCATTTTAATTTCAATCAttt         <  1:3252605/148‑1 (MQ=255)
                                                                                                                                              ttAAGGCAGTATAAAATGCTGGTTTTGTCGTCAGTTCAAGGCAGGATAAGGGTTAACACACCTTTATGACAGTCAGGAATTGACTGTTTCTCTAACGACTTCCCTTTTAGCCTTAAAGATAAAATCCATTTTAATTTCAATCATTtaa       >  2:94938/1‑148 (MQ=255)
                                                                                                                                              ttAAGACAGTATAAAATGCTGGTTTTGTCGTCAGTTCAAGGCAGGATAAGGGTTAACACACCTTTATGACAGTCAGGAATTGACTGTTTCTCTAACGACTTCCCTTTTAGCCTTAAAGATAAAATCCATTTTAATTTCAATCATTtaa       >  2:1397434/1‑148 (MQ=255)
                                                                                                                                               tAAGGCAGTATAAAATGCTGGTTTTGTCGTCAGTTCAAGGCAGGATAAGGGTTAACACACCTTTATGACAGTCAGGAATTGACTGTTTCTCTAACGACTTCCCTTTTAGCCTTAAAGATAAAATCCATTTTAATTTCAATCATTtaat      >  1:2281759/1‑148 (MQ=255)
                                                                                                                                                 aGGCAGTATAAAATGCTGGTTTTCTCGTCAGTTCAAGGCAGGATAAGGGTTAACACACCTTTATGACAGTCAGGAATTGACTGTTTCTCTAACGACTTCCCTTTTAGCCTTAAAGATAAAATCCATTTTAATTTCAATCATTtaataa    <  1:1105591/148‑1 (MQ=255)
                                                                                                                                                   gCAGTATAAAATGCTGGTTTTGTCGTCAGTTCAAGGCAGGATAAGGGTTAACACACCTTTATGACAGTCAGGAATTGACTGTTTCTCTAACGACTTCCCTTTTAGCCTTAAAGATAAAATCCATTTTAATTTCAATCATTTAATAAAg  >  1:2697700/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
GCTTTTATCAACTGCATAATCAATCAAAATTACCGAAATTTCATGCATAATCACATAAATCACTTTTGCTTATCTTGTGTCAGATTTTTTTATCTCCTGATGGATTTTAGGCAAAAACAGTAGCATGAAACGTCATTACCAATTAAGGCAGTATAAAATGCTGGTTTTGTCGTCAGTTCAAGGCAGGATAAGGGTTAACACACCTTTATGACAGTCAGGAATTGACTGTTTCTCTAACGACTTCCCTTTTAGCCTTAAAGATAAAATCCATTTTAATTTCAGTCATTTAATAAAG  >  CP009273/3347497‑3347791

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: