Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,009,187 G→T 100% intergenic (+90/‑232) mocA → / → ygfK CTP:molybdopterin cytidylyltransferase/putative Fe‑S subunit oxidoreductase subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,009,1870GT97.1% 796.3 / 11.6 255intergenic (+90/‑232)mocA/ygfKCTP:molybdopterin cytidylyltransferase/putative Fe‑S subunit oxidoreductase subunit
Reads supporting (aligned to +/- strand):  ref base G (3/5);  new base T (117/130);  total (120/135)
Fisher's exact test for biased strand distribution p-value = 7.26e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 272 total aligned reads displayed.

ATATTGATACTCCGGACGACTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAGAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGAGAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCGAGCCGCCTCACATTTTTTTATATTTCCCGCCAAACCTGGCAAGAGT  >  CP009273/3009042‑3009331
                                                                                                                                                 |                                                                                                                                                
  atTGATACTCCGGACGACTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAGAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTgagaaaagataaaa                                                                                                                                              >  1:3455089/1‑148 (MQ=255)
  atTGATACTCCGGACGACTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAGAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTgagaaaagataaaa                                                                                                                                              >  2:2035422/1‑148 (MQ=255)
  atTGATACTCCGGACGACTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAGAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTgagaaaagataaaa                                                                                                                                              <  1:1537687/148‑1 (MQ=255)
    tGATACTCCGGACGACTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAGAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAAt                                                                                                                                            <  1:2785449/148‑1 (MQ=255)
      aTACTCCTGACGACTTTATTACATCGAAGGAAAGGTATACTGAAATTTAGAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCTGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTg                                                                                                                                          <  2:2547457/148‑1 (MQ=255)
      aTACTCCGGACGACTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAGAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTg                                                                                                                                          <  2:1080576/148‑1 (MQ=255)
       tACTCCGGACGACTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAGAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGa                                                                                                                                         >  2:370217/1‑148 (MQ=255)
          tCCGGACGACTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAGAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAATGATAAAAATTGATaa                                                                                                                                      >  1:778229/1‑148 (MQ=255)
          tCCGGACGACTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAGAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATaa                                                                                                                                      <  2:276856/148‑1 (MQ=255)
           ccGGACGACTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAGAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATaaa                                                                                                                                     <  1:256824/148‑1 (MQ=255)
           ccGGACGACTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAGAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATaaa                                                                                                                                     <  2:354209/148‑1 (MQ=255)
            cGGACGACTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAGAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAAAt                                                                                                                                    <  2:440624/148‑1 (MQ=255)
            cGGACGACTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAGAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAAAt                                                                                                                                    <  2:440836/148‑1 (MQ=255)
               acgacTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAGAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGAGAAAAATTGATAAATggg                                                                                                                                 >  1:1802108/1‑148 (MQ=255)
                  acTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAGAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTg                                                                                                                              <  2:2869511/148‑1 (MQ=255)
                  acTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAGAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTg                                                                                                                              >  1:437287/1‑148 (MQ=255)
                   cTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAGAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTgg                                                                                                                             >  1:3519927/1‑148 (MQ=255)
                         taCAGCGAAGGAAAGGTATACTGAAATTTAGAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCaataa                                                                                                                       >  2:1880280/1‑148 (MQ=255)
                            aGCGAAGGAAAGGTATACTGAAATTTAGAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCaataataa                                                                                                                    >  1:3677842/1‑148 (MQ=255)
                              cGAAGGAAAGGTATACTGAAATTTAGAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAAtt                                                                                                                  >  1:863511/1‑148 (MQ=255)
                                 agaaaaGGTATACTGAAATTTAGAAAATGTAGTTAAACGATTGCGTTCAAATATTTAACCCCTCCAGCACTCTTTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGAGAAAAAGTGATAAATGGGCTGGCAATAATAATTTTc                                                                                                               <  1:888937/145‑1 (MQ=255)
                                 aGGAAAGGTATACTGAAATTTAGAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTc                                                                                                               <  2:407558/148‑1 (MQ=255)
                                      aGGTATACTGAAATTTAGAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAg                                                                                                          >  1:3584310/1‑148 (MQ=255)
                                        gTATACTGAAATTTAGAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTg                                                                                                        <  1:3547233/148‑1 (MQ=255)
                                        gTATACTGAAATTTAGAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCGTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTg                                                                                                        <  2:2990602/148‑1 (MQ=255)
                                            aCTGAAATTTAGAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTTCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTTATAAATGGGCTGGCAATAATAATTTTCAGGAGTGAttt                                                                                                    >  1:496678/1‑148 (MQ=255)
                                            aCTGAAATTTAGAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTTCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGAttt                                                                                                    >  1:496078/1‑148 (MQ=255)
                                              tGAAATTTAGAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTa                                                                                                  <  2:589819/148‑1 (MQ=255)
                                               gAAATTTAGAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTaa                                                                                                 >  1:3803495/1‑148 (MQ=255)
                                               gAAATTTAGAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTaa                                                                                                 <  1:3670184/148‑1 (MQ=255)
                                                aaaTTTAGAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAAt                                                                                                <  1:190738/148‑1 (MQ=255)
                                                  aTTTAGAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATtc                                                                                              <  1:3696948/148‑1 (MQ=255)
                                                    ttAGAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATtctc                                                                                            <  1:1186270/148‑1 (MQ=255)
                                                    ttAGAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATtctc                                                                                            <  2:2882781/148‑1 (MQ=255)
                                                    ttAGAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATtctc                                                                                            >  2:3058134/1‑148 (MQ=255)
                                                    ttAGAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATtctc                                                                                            >  1:1807048/1‑148 (MQ=255)
                                                        aaaaCGTAGTTAAACGATTGCGTTCAAATATTTAATCCTACCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCAttt                                                                                        <  2:761/148‑1 (MQ=255)
                                                         aaaCGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTa                                                                                       <  2:3113971/148‑1 (MQ=255)
                                                         aaaCGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTa                                                                                       >  2:1237167/1‑148 (MQ=255)
                                                          aaCGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAg                                                                                      >  2:3771792/1‑148 (MQ=255)
                                                           aCGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGc                                                                                     >  1:1987509/1‑148 (MQ=255)
                                                             gTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCAc                                                                                   <  2:1144691/148‑1 (MQ=255)
                                                                   aaaCGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGa                                                                             >  1:3050917/1‑148 (MQ=255)
                                                                    aaCGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAg                                                                            >  2:2243515/1‑148 (MQ=255)
                                                                    aaCGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAg                                                                            <  2:1814545/148‑1 (MQ=255)
                                                                     aCGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGt                                                                           >  1:1185240/1‑148 (MQ=255)
                                                                       gATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGttc                                                                         <  2:2590705/148‑1 (MQ=255)
                                                                        aTTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGttca                                                                        <  1:1448876/148‑1 (MQ=255)
                                                                          tGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGttcatt                                                                      >  1:2147098/1‑148 (MQ=255)
                                                                                  aaTATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACaaaaa                                                              >  1:1595223/1‑148 (MQ=255)
                                                                                    tatTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGc                                                            <  2:3455089/148‑1 (MQ=255)
                                                                                     atTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGcc                                                           >  1:3162428/1‑148 (MQ=255)
                                                                                       ttAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCt                                                         <  1:269029/148‑1 (MQ=255)
                                                                                       ttAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCt                                                         >  1:3190486/1‑148 (MQ=255)
                                                                                          aTCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTaa                                                      >  1:2947734/1‑148 (MQ=255)
                                                                                          aTCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTaa                                                      <  2:3664679/148‑1 (MQ=255)
                                                                                            ccTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGt                                                    >  2:1446554/1‑148 (MQ=255)
                                                                                              ttCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGttt                                                  >  2:3602958/1‑148 (MQ=255)
                                                                                              ttCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGttt                                                  >  1:593412/1‑148 (MQ=255)
                                                                                                  ggCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCGAg                                              <  2:2021811/148‑1 (MQ=255)
                                                                                                   gCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCGAgc                                             >  2:2372627/1‑148 (MQ=255)
                                                                                                   gCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCGAgc                                             <  2:1204327/148‑1 (MQ=255)
                                                                                                       tCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCGAgccgcc                                         <  2:474667/148‑1 (MQ=255)
                                                                                                       tCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCGAgccgcc                                         >  1:17092/1‑148 (MQ=255)
                                                                                                       tCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCGAgccgcc                                         >  2:3564140/1‑148 (MQ=255)
                                                                                                         aTTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCGAGCCGCCTc                                       <  2:2683088/148‑1 (MQ=255)
                                                                                                          ttATCATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCGAGCCGCCTca                                      >  2:324067/1‑148 (MQ=255)
                                                                                                          ttATCATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCGAGCCGCCTca                                      >  1:1101511/1‑148 (MQ=255)
                                                                                                          ttATCATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCGAGCCGCCTca                                      >  2:1717431/1‑148 (MQ=255)
                                                                                                          ttATCATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCGAGCCGCCTca                                      <  2:1595223/148‑1 (MQ=255)
                                                                                                           tATCATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCGAGCCGCCTcac                                     >  1:2759090/1‑148 (MQ=255)
                                                                                                             tCATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCGAGCCGCCTCACAt                                   >  2:3730847/1‑148 (MQ=255)
                                                                                                              cATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCGAGCCGCCTCACAtt                                  <  1:709252/148‑1 (MQ=255)
                                                                                                              cATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCGAGCCGCCTCACAtt                                  <  2:1089896/148‑1 (MQ=255)
                                                                                                                  aTGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCGAGCCGCCTCACAtttttt                              <  2:432019/148‑1 (MQ=255)
                                                                                                                   tGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCGAGCCGCCTCACAttttttt                             >  1:1021564/1‑148 (MQ=255)
                                                                                                                   tGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCGAGCCGCCTCACAttttttt                             >  2:1342224/1‑148 (MQ=255)
                                                                                                                    gATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCGAGCCGCCTCACATTTTTTta                            <  2:3273665/148‑1 (MQ=255)
                                                                                                                    gATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCGAGCCGCCTCACATTTTTTta                            <  1:2659960/148‑1 (MQ=255)
                                                                                                                     aTAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCGAGCCGCCTCACATTTTTTtat                           <  1:2061441/148‑1 (MQ=255)
                                                                                                                      tAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCGAGCCGCCTCACATTTTTTtata                          >  1:730435/1‑148 (MQ=255)
                                                                                                                         ggAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCGAGCCGCCTCACATTTTTTTATAttt                       <  2:1497524/148‑1 (MQ=255)
                                                                                                                         ggAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCGAGCCGCCTCACATTTTTTTATAttt                       >  2:2916989/1‑148 (MQ=255)
                                                                                                                          gAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCGAGCCGCCTCACATTTTTTTATATTTc                      <  2:27502/148‑1 (MQ=255)
                                                                                                                          gAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCGAGCCGCCTCACATTTTTTTATATTTc                      <  2:3181775/148‑1 (MQ=255)
                                                                                                                             tCACTGATTTTGAGAAAAGAGAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCGAGCCGCCTCACATTTTTTTATATTTCCCg                   >  1:3200529/1‑148 (MQ=255)
                                                                                                                              cACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCGAGCCGCCTCACATTTTTTTATATTTCCCGc                  >  2:3539111/1‑148 (MQ=255)
                                                                                                                              cACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCGAGCCGCCTCACATTTTTTTATATTTCCCGc                  >  2:2697774/1‑148 (MQ=255)
                                                                                                                              cACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCAACTGTGAGTTCATTAACAAAAAGCCATTAAGTTTCGAGCCGACTCACATTTTTTTATATATCACga                  >  2:3208974/1‑147 (MQ=255)
                                                                                                                               aCTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCGAGCCGCCTCACATTTTTTTATATTTCCCGcc                 >  1:3400943/1‑148 (MQ=255)
                                                                                                                               aCTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCGAGCCGCCTCACATTTTTTTATATTTCCCGcc                 >  1:3259768/1‑148 (MQ=255)
                                                                                                                                  gATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCGAGCCGCCTCACATTTTTTTATATTTCCCGCCaaa              >  1:1146105/1‑148 (MQ=255)
                                                                                                                                  gATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCGAGCCGCCTCACATTTTTTTATATTTCCCGCCaaa              <  2:863511/148‑1 (MQ=255)
                                                                                                                                     tttGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCGAGCCGCCTCACATTTTTTTATATTTCCCGCCAAACCt           >  2:1477126/1‑148 (MQ=255)
                                                                                                                                     tttGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCGAGCCGCCTCACATTTTTTTATATTTCCCGCCAAACCt           <  2:939876/148‑1 (MQ=255)
                                                                                                                                      ttGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCGAGCCGCCTCACATTTTTTTATATTTCCCGCCAAACCTg          <  1:2774601/148‑1 (MQ=255)
                                                                                                                                      ttGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCGAGCCGCCTCACATTTTTTTATATTTCCCGCCAAACCTg          <  2:730435/148‑1 (MQ=255)
                                                                                                                                       tGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCGAGCCGCCTCACATTTTTTTATATTTCCCGCCAAACCTgg         <  1:215889/148‑1 (MQ=255)
                                                                                                                                           aaaagataaaaATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCGAGCCGCCTCACATTTTTTTATATTTCCCGCCAAACCTGGCAag     >  2:3349325/1‑148 (MQ=255)
                                                                                                                                              agataaaaATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCGAGCCGCCTCACATTTTTTTATATTTCCCGCCAAACCTGGCAAGAgt  <  1:1346319/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
ATATTGATACTCCGGACGACTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAGAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGAGAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCGAGCCGCCTCACATTTTTTTATATTTCCCGCCAAACCTGGCAAGAGT  >  CP009273/3009042‑3009331

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: