Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,816,882 C→A 10.8% A56A (GCC→GCA nadE → NAD synthetase, NH3/glutamine‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,816,8820CA10.8% 943.4 / 56.7 263A56A (GCC→GCAnadENAD synthetase, NH3/glutamine‑dependent
Reads supporting (aligned to +/- strand):  ref base C (135/99);  new base A (9/20);  total (144/119)
Fisher's exact test for biased strand distribution p-value = 9.31e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 266 total aligned reads displayed.

ATAATAAAGGCGCTGGGCGCAAAACCGCAGATTAATGCTGAAGAGGAAATTCGTCGTAGTGTCGATTTTCTGAAAAGCTACCTGCAAACTTATCCGTTCATTAAATCACTGGTGCTCGGGATCAGCGGCGGTCAGGACTCCACGCTTGCCGGAAAGCTGTGCCAGATGGCGATTAATGAGCTGCGCCTGGAAACCGGCAACGAATCACTGCAATTTATTGCCGTACGCCTGCCCTATGGTGTTCAGGCCGACGAACAAGATTGCCAGGATGCCATTGCCTTTATTCAACCGGATCGCG  >  CP009273/1816733‑1817030
                                                                                                                                                     |                                                                                                                                                    
ataataAAGGCGCTGGGCGCAAAACCGCAGATTAATGCTGAAGAGGAAATTCGTCGTAGTGTCGATTTTCTGAAAAGCTACCTGCAAACTTATCCGTTCATTAAATCACTGGTGCTCGGGATCAGCGGCGGTCAGGACTCCACGCTTGcc                                                                                                                                                      >  1:100220/1‑150 (MQ=255)
 taataaAGGCGCTGGGCGCAAAACCGCAGATTAATGCTGAAGAGGAAATTCGTCGTAGTGTCGATTTTCTGAAAAGCTACCTGCAAACTTATCCGTTCATTAAATCACTGGTGCTCGGGATCAGCGGCGGTCAGGACTCCACGCTTGCCg                                                                                                                                                     >  2:3444033/1‑150 (MQ=255)
 taataaAGGCGCTGGGCGCAAAACCGCAGATTAATGCTGAAGAGGAAATTCGTCGTAGTGTCGATTTTCTGAAAAGCTACCTGCAAACTTATCCGTTCATTAAATCACTGGTGCTCGGGATCAGCGGCGGTCAGGACTCCACGCTTGCCg                                                                                                                                                     <  2:2222857/150‑1 (MQ=255)
  aataaAGGCGCTGGGCGCAAAACCGCAGATTAATGCTGAAGAGGAAATTCGTCGTAGTGTCGATTTTCTGAAAAGCTACCTGCAAACTTATCCGTTCATTAAATCACTGGTGCTCGGGATCAGCGGCGGTCAGGACTCCACGCTTGCCgg                                                                                                                                                    >  2:2794304/1‑150 (MQ=255)
   ataaAGGCGCTGGGCGCAAAACCGCAGATTAATGCTGAAGAGGAAATTCGTCGTAGTGTCGATTTTCTGAAAAGCTACCTGCAAACTTATCCGTTCATTAAATCACTGGTGCTCGGGATCAGCGGCGGTCAGGACTCCACGCTTGCCGGa                                                                                                                                                   >  2:530089/1‑150 (MQ=255)
     aaaGGCGCTGGGCGCAAAACCGCAGATTAATGCTGAAGAGGAAATTCGTCGTAGTGTCGATTTTCTGAAAAGCTACCTGCAAACTTATCCGTTCATTAAATCACTGGTGCTCGGGATCAGCGGCGGTCAGGACTCCACGCTTGCCGGaaa                                                                                                                                                 >  2:3160226/1‑150 (MQ=255)
         gcgcTGGGCGCAAAACCGCAGATTAATGCTGAAGAGGAAATTCGTCGTAGTGTCGATTTTCTGAAAAGCTACCTGCAAACTTATCCGTTCATTAAATCACTGGTGCTCGGGATCAGCGGCGGTCAGGACTCCACGCTTGCCGGAAAGCtg                                                                                                                                             <  2:3338936/150‑1 (MQ=255)
           gcTGGGCGCAAAACCGCAGATTAANGCTGAAGANGAAATTCGTNGTAGTGTCGATTTTCTGAAAAGCTACCTGCAAACTTATCCGTTCATTTAATCACTGTTGCTCGGGATCAGCGGCGGTCAGGACTCCACGCTTGCAGGAAAGCtgtg                                                                                                                                           <  1:1347638/150‑1 (MQ=255)
              gggCGCAAAACCGCAGATTAATGCTGAAGAGGAAATTCGTCGTAGTGTCGATTTTCTGAAAAGCTACCTGCAAACTTATCCGTTCATTAAATCACTGGTGCTCGGGATCAGCGGCGGTCAGGACTCCACGCTTGCCGGAAAGCTGTGCCa                                                                                                                                        <  1:1134468/150‑1 (MQ=255)
               ggCGCAAAACCGCAGATTAATGCTGAAGAGGAAATTCGTCGTAGTGTCGATTTTCTGAAAAGCTACCTGCAAACTTATCCGTTCATTAAATCACTGGTGCTCGGGATCAGCGGCGGTCAGGACTCCACGCTTGCCGGAAAGCTGTGCCAg                                                                                                                                       >  2:2979002/1‑150 (MQ=255)
                      aaCCGCAGATTAATGCTGAAGAGGAAATTCGTCGTAGTGTCGATTTTCTGAAAAGCTACCTGCAAACTTATCCGTTCATTAAATCACTGGTGCTCGGGATCAGCGGCGGTCAGGACTCCACGCTTGCCGGAAAGCTGTGCCAGATGGCGa                                                                                                                                <  2:1522875/150‑1 (MQ=255)
                         cGCAGATTAATGCTGAAGAGGAAATTCGTCGTAGTGTCGATTTTCTGAAAATCTACCTGCAAACGTATCCGTTCATTAAATCACTGGTGCTCGGGATCAGCGGCGGTCAGGACTCCACGCTTGCCGGCAAGCTGTGCCAGATGGCGAATa                                                                                                                             >  1:3557158/1‑150 (MQ=255)
                          gCAGATTAATGCTGAAGAGGAAATTCGTCGTAGTGTCGATTTTCTGAAAAGCTACCTGCAAACTTATCCGTTCATTAAATCACTGGTGCTCGGGATCAGCGGCGGTCAGGACTCCACGCTTGCCGGAAAGCTGTGCCAGATGGCGATTaa                                                                                                                            <  1:1604156/150‑1 (MQ=255)
                           cAGATTAATGCTGAAGAGGAAATTCGTCGTAGTGTCGATTTTCTGAAAAGCTACCTGCAAACTTATCCGTTCATTAAATCACTGGTGCTCGGGATCAGCGGCGGTCAGGACTCCACGCTTGCCGGAAAGCTGTGCCAGATGGCGATTAAt                                                                                                                           <  1:2052946/150‑1 (MQ=255)
                             gATTAATGCTGAAGAGGAAATTCGTCGTAGTGTCGATTTTCTGAAAAGCTACCTGCAAACTTATCCGTTCATTAAATCACTGGTGCTCGGGATCAGCGGCGGTCAGGACTCCACGCTTGCCGGAAAGCTGTGCCAGATGGCGATTAATGa                                                                                                                         <  2:139359/150‑1 (MQ=255)
                                 aaTGCTGAAGAGGAAATTCGTCGTAGTGTCGATTTTCTGAAAAGCTACCTGCAAACTTATCCGTTCATTAAATCACTGGTGCTCGGGATCAGCGGCGGTCAGGACTCCACGCTTGCCGGAAAGCTGTGCCAGATGGCGATTAATGAGCTg                                                                                                                     <  2:4165012/150‑1 (MQ=255)
                                   tGCTGAAGAGGAAATTCGTCGTAGTGTCGATTTTCTGAAAAGCTACCTGCAAACTTATCCGTTCATTAAATCACTGGTGCTCGGGATCAGCGGCGGTCAGGACTCCACGCTTGCCGGAAAGCTGTGCCAGATGGCGATTAATGAGCTgcg                                                                                                                   <  1:1775881/150‑1 (MQ=255)
                                     cTGAAGAGGAAATTCGTCGTAGTTTCGATTTTCTGAAAAGCTACCTGCAAAATTATCCGTTCATTAAATCACTGGTGCTCGGGATCAGCGGCGGTCAGGAAGACACGCTTGACGGAAAGCTGTGCCACATGGCGATTAATGAGCTGCGcc                                                                                                                 >  2:1275461/1‑150 (MQ=255)
                                      tGAAGAGGAAATTCGTCGTAGTGTCGATTTTCTGAAAAGCTACCTGCAAACTTATCCGTTCATTAAATCACTGGTGCTCGGGATCAGCGGCGGTCAGGACTCCACGCTTGCAGGAAAGCTGTGCCAGATGGCGATTAATGAGCTGCGCCt                                                                                                                <  2:2275468/150‑1 (MQ=255)
                                       gAAGAGGAAATTCGTCGTAGTGTCGATTTTCTGAAAAGCTACCTGCAAACTTATCCGTTCATTAAATCACTGGTGCTCGGGATCAGCGGCGGTCAGGACTCCACGCTTGCCGGAAAGCTGTGCCAGATGGCGATTAATGAGCTGCGCCTg                                                                                                               <  1:3229771/150‑1 (MQ=255)
                                             gAAATTCGTCGTAGTGTCGATTTTCTGAAAAGCTACCTGCAAACTTATCCGTTCATTAAATCACTGGTGCTCGGGATCAGCGGCGGTCAGGACTCCACGCTTGCCGGAAAGCTGTGCCAGATGGCGATTAATGAGCTGCGCCTGGAAAcc                                                                                                         >  2:3579166/1‑150 (MQ=255)
                                             gAAATTCGTCGTAGTGTCGATTTTCTGAAAAGCTACCTGCAAACTTATCCGTTCATTAAATCACTGGTGCTCGGGATCAGCGGCGGTCAGGACTCCACGCTTGCAGGAAAGCTGTGCCAGATGGCGATTAATGAGCTGCGCCTGGAAAcc                                                                                                         <  1:1230943/150‑1 (MQ=255)
                                              aaaTTCGTCGTAGTGTCGATTTTCTGAAAAGCTACCTGCAAACTTATCCGTTCATTAAATCACTGGTGCTCGGGATCAGCGGCGGTCAGGACTCCACGCTTGCAGGAAAGCTGTGCCAGATGGCGATTAATGAGCTGCGCCTGGAAACCg                                                                                                        <  2:2912576/150‑1 (MQ=255)
                                               aaTTCGTCGTAGTGTCGATTTTCTGAAAAGCTACCTGCAAACTTATCCGTTCATTAAATCACTGGTGCTCGGGATCAGCGGCGGTCAGGACTCCACGCTTGCCGGAAAGCTGTGCCAGATGGCGATTAATGAGCTGCGCCTGGAAACCgg                                                                                                       >  2:119579/1‑150 (MQ=255)
                                                aTTCGTCGTAGTGTCGATTTTCTGAAAAGCTACCTGCAAACTTATCCGTTCATTAAATCACTGGTGCTCGGGCTAAGCGGCGGTCAGGACTCCACGCTTGCCGGAAAGCTGTGCCAGATGGCGATTAATGAGCTGCGCCTGGAAACCGGc                                                                                                      >  1:1997836/1‑150 (MQ=255)
                                                aTTCGTCGTAGTGTCGATTTTCTGAAAAGCTACCTGCAAACTTATCCGTTCATTAAATAACTGGTGCTCGGGATCAGCGGCGGTCAGGACTCCACGCTTGCCGGAAAGCTGTGCCAGATGGCGATTAATGAGCTGCGCCTGGAAACCGGc                                                                                                      >  1:4098656/1‑150 (MQ=255)
                                                 ttCGTCGTAGTGTCGATTTTCTGAAAAGCTACCTGCAAACTTATCCGTTCATTAAATCACTGGTGCTCGGGATCAGCGGCGGTCAGGACTCCACGCTTGCCGGAAAGCTGTGCCAGATGGCGATTAATGAGCTGCGCCTGGAAACCGGCa                                                                                                     >  2:2829425/1‑150 (MQ=255)
                                                  tcgtcgTAGTGTCGATTTTCTGAAAAGCTACCTGCAAACTTATCCGTTCATTAAATCACTGGTGCTCGGGATCAGCGGCGGTCAGGACTCCACGCTTGCCGGAAAGCTGTGCCAGATGGCGATTAATGAGCTGCGCCTGGAAACCGGCaa                                                                                                    <  2:3887110/150‑1 (MQ=255)
                                                  tcgtcgTAGTGTCGATTTTCTGAAAAGCTACCTGCAAACTTATCCGTTCATTAAATCACTGGTGCTCGGGATCAGCGGCGGTCAGGACTCCACGCTTGCCGGAAAGCTGTGCCAGATGGCGATTAATGAGCTGCGCCTGGAAACCGGCaa                                                                                                    <  2:1702025/150‑1 (MQ=255)
                                                  tcgtcgTAGTGTCGATTTTCTGAAAAGCTACCTGCAAACTTATCCGTTCAATAAATCACTGGTGCTCGGGATCAGCGGCGGTCAGGACTCCACGCTTGCCGGAAAGCTGTGCCAGATGGCGATTAATGAGCTGCGCCTGGAAACCGGCaa                                                                                                    >  2:1614534/1‑150 (MQ=255)
                                                   cgtcgtAGTGTCGATTTTCTGAAAAGCAACCTGCAAACTTATCCGTTCATTAAATCACAGGTGCTCGGGATCAGCGGCGGTCAGGACTCCACGCTTGCATGAAAGCTGTGCCAGATGGCGATTAATTAGCTGCGCCAGGAAACCGGCAAc                                                                                                   >  1:320291/1‑150 (MQ=255)
                                                    gtcgtAGTGTCGATTTTCTGAAAAGCTACCTGCAAACTTATCCGTTCATTAAATCACTGGTGCTCGGGATCAGCGGCGGTCAGGACTCCACGCTTGCCGGAAAGCTGTGCCAGATGGCGATTAATGAGCTGCGCCTGGAAACCGGCAACg                                                                                                  >  1:2309785/1‑150 (MQ=255)
                                                     tcgtAGTGTCGATTTTCTGAAAAGCTACCAGCAAACTTATCCGTTCATTAAATCACTGGTGCTCGGGATCAGCGGCGGTCAGGACTCCACGCTTGCCGGAAAGCTGTGCCAGATGGCGATTAATGAGCTGCGCCTGGAAACCGGCAACGa                                                                                                 >  2:3452541/1‑150 (MQ=255)
                                                       gtAGTGTCGATTTTCTGAAAAGCTACCTGCAAACTTATCCGTTCATTAAATCACTGGTGCTCGGGATCAGCGGCGGTCAGGACTCCACGCTTGCCGGAAAGCTGTGCCAGATGGCGATTAATGAGCTGCGCCTGGAAACCGGCAACGAAt                                                                                               <  2:1915593/150‑1 (MQ=255)
                                                       gtAGTGTCGATTTTCTGAAAAGCTACCTGCAAACTTATCCGTTCATTAAATCACTGGTGCTCGGGATCAGCGGCGGTCAGGACTCCACGCTTGCCGGAAAGCTGTGCCAGATGGCGATTAATGAGCTGCGCCTGGAAACCGGCAACGAAt                                                                                               <  2:2303606/150‑1 (MQ=255)
                                                        tAGTGTCGATTTTCTGAAAAGCTACCTGCAAACTTATCCGTTCATTAAATCACTGGTGCTCGGGATCAGCGGCGGTCAGGACTCCACGCTTGCCGGAAAGCTGTGCCAGATGGCGATTAATGAGCTGCGCCTGGAAACCGGCAACGAATc                                                                                              <  2:672783/150‑1 (MQ=255)
                                                        tAGTGTCGATTTTCTGAAAAGCTACCTGCAAACTTATCCGTTCATTAAATCACTGGTGCTCGGGATCAGCGGCGGTCAGGACTCCACGCTTGCCGGAAAGCTGTGCCAGATGGCGATTAATGAGCTGCGCCTGGAAACCGGCAACGAATc                                                                                              <  2:212976/150‑1 (MQ=255)
                                                         aGTGTCGATTTTCTGAAAAGCTACCTGCAAACTTATCCGTTCATTAAATCACTGGTGCTCGGGATCAGCGGCGGTCAGGACTCCACGCTTGCCGGAAAGCTGTGCCAGATGGCGATTAATGAGCTGCGCCTGGAAACCGGCAACGAATCa                                                                                             <  1:512946/150‑1 (MQ=255)
                                                         aGTGTCGATTTTCTGAAAAGCTACCTGCAAACTTATCCGTTCATTAAATCACTGGTGCTCGGGATCAGCGGCGGTCAGGACTCCACGCTTGCAGGAAAGCTGTGCCAGATGGCGATTAATGAGCTGCGCCTGGAAACCGGCAACGAATCa                                                                                             <  2:1010991/150‑1 (MQ=255)
                                                               gATTTTCTGAAAAGCTACCTGCAAACTTATCCGTTCATTAAATCACTGGTGCTCGGGATCAGCGGCGGTAAGGACTCCACGCTTGCAGGAAAGCTGTGCCAGATGGCGATTAATGAGCTGCGCCTGGAAACCTGCAACGAATCAATCCaa                                                                                       >  1:2945900/1‑150 (MQ=255)
                                                                aTTTTCTGAAAAGCTACCTGCAAACTTATGCGTTCATTAAATCACTGGTGCTCGGGATCAGCGGCGGTCAGGACTCCACGCTTGCCGGAAAGCTGTGCCAGATGGCGATTAATGAGCTGCGCCTGGAAACCGGCAACGAATCACTGCAAt                                                                                      <  1:2446666/150‑1 (MQ=255)
                                                                aTTTTCTGAAAAGCTACCTGCAAACTTATCCGTTCATTAAATCACTGGTGCTCGGGATCATCGGCGGTCAGGACTCCACGCTTGCCGGAAAGCTGTGCCAGATGGCGATTAATGAGCTGCGCCTGGAAACCGGCAACGAATCACTGCAAt                                                                                      >  1:3571402/1‑150 (MQ=255)
                                                                  tttCTGAAAAGCTACCTGCAAACTTATCCGTTCATTAAATCACTGGTGCTCGGGATCAGCGGCGGTCAGGACTCCACGCTTGCCGGAAAGCTGTGCCAGATGGCGATTAATGAGCTGCGCCTGGAAACCGGCAACGAATCACTGCAAttt                                                                                    >  2:2009824/1‑150 (MQ=255)
                                                                  tttCTGAAAAGCTACCTGCAAACTTATCCGTTCATTAAATCACTGGTGCTCGGGATCAGCGGCGGTCAGGACTCCACGCTTGCCGGAAAGCTGTGCCAGATGGCGATTAATGAGCTGCGCCTGGAAACCGGCAACGAATCACTGCAAttt                                                                                    >  2:2449325/1‑150 (MQ=255)
                                                                   ttCTGAAAAGCTACCTGCAAACTTATCCGTTCATTAAATCACTGGTGCTCGGGATCAGCGGCGGTCAGGACTCCACGCTTGCCGGAAAGCTGTGCCAGATGGCGATTAATGAGCTGCGCCTGGAAACCGGCAACGAATCACTGCAATTTa                                                                                   >  2:2324035/1‑150 (MQ=255)
                                                                      tGAAAAGCTACCTGCAAACTTATCCGTTCATTAAATCACTGGTGCTCGGGATCAGCGGCGGTCAGGACTCCACGCTTGCCGGAAAGCTGTGCCAGATGGCGATTAATGAGCTGCGCCTGGAAACCGGCAACGAATCACTGCAATTTATTg                                                                                >  1:3909302/1‑150 (MQ=255)
                                                                      tGAAAAGCTACCTGCAAACTTATCCGTTCATTAAATCACTGGTGCTCGGGATCAGCGGCGGTCAGGACTCCACGCTTGCCGGAAAGCTGTGCCAGATGGCGATTAATGAGCTGCGCCTGGAAACCGGCAACGAATCACTGCAATTTATTg                                                                                >  2:3712780/1‑150 (MQ=255)
                                                                      tGAAAAGCTACCTGCAAACTTATCCGTTCATTAAATCACTGGTGCTCGGGATCAGCGGCGGTCAGGACTCCACGCTTGCCGGAAAGCTGTGCCAGATGGCGATTAATGAGCTGCGCCTGGAAACCGGCAACGAATCACTGCAATTTATTg                                                                                >  2:121268/1‑150 (MQ=255)
                                                                       gAAAAGCTACCTGCAAACTTATCCGTTCATTAAATCACTGGTGCTCGGGATCAGCGGCGGTCAGGACTCCACGCTTGCCGGAAAGCTGTGCCAGATGGCGATTAATGAGCTGCGCCTGGAAACCGGCAACGAATCACTGCAATTTATTGc                                                                               <  1:821648/150‑1 (MQ=255)
                                                                        aaaaGCTACCTGCAAACTTATCCGTTCATTAAATCACTGGTGCTCGGGATCAGCGGCGGTCAGGACTCCACGCTTGCCGGAAAGCTGTGCCAGATGGCGATTAATGAGCTGCGCCTGGAAACCGGCAACGAATCACTGCAATTTATTGcc                                                                              >  2:1688725/1‑150 (MQ=255)
                                                                        aaaaGCTACCTGCAAACTTATCCGTTCATTAAATCACTGGTGCTCGGGATCAGCGGCGGTCAGGACTCCACGCTTGCCGGAAAGCTGTGCCAGATGGCGATTAATGAGCTGCGCCTGGAAACCGGCAACGAATCACTGCAATTTATTGcc                                                                              <  2:3765240/150‑1 (MQ=255)
                                                                         aaaGCTACCTGCAAACTTATCCGTTCATTAAATCACTGGTGCTCGGGATCAGCGGCGGTCAGGACTCCACGCTTGCCGGAAAGCTGTGCCAGATGGCGATTAATGAGCTGCGCCTGGAAACCGGCAACGAATCACTGCAATTTATTGCCg                                                                             >  1:3789825/1‑150 (MQ=255)
                                                                         aaaGCTACCTGCAAACTTATCCGTTCATTAAATCACTGGTGCTCGGGATCAGCGGCGGTCAGGACTCCACGCTTGCCGGAAAGCTGTGCCAGATGGCGATTAATGAGCTGCGCCTGGAAACCGGCAACGAATCACTGCAATTTATTGCCg                                                                             >  2:2523254/1‑150 (MQ=255)
                                                                          aaGCTACCTGCAAACTTATCCGTTCATTAAATCACTGGTGCTCGGGATCAGCGGCGGTCAGGACTCCACGCTTGCCGGAAAGCTGTGCCAGATGGCGATTAATGAGCTGCGCCTGGAAACCGGCAACGAATCACTGCAATTTATTGCCGt                                                                            >  2:3260700/1‑150 (MQ=255)
                                                                           aGCTACCTGCAAACTTATCCGTTCATTAAATCACTGGTGCTCGGGATCAGCGGCGGTCAGGACTCCACGCTTGCCGGAAAGCTGTGCCAGATGGCGATTAATGAGCTGCGCCTGGAAACCGGCAACGAATCACTGCAATTTATTGCCGTa                                                                           <  1:1625545/150‑1 (MQ=255)
                                                                             cTACCTGCAAACTTATCCGTTCATTAAATCACTGGTGCTCGGGATCAGCGGCGGTCAGGACTCCACGCTTGCCGGAAAGCTGTGCCAGATGGCGATTAATGAGCTGCGCCTGGAAACCGGCAACGAATCACTGCAATTTATTGCCGTACg                                                                         <  1:3993752/150‑1 (MQ=255)
                                                                              tACCTGCAAACTTATCCGTTCATTAAATCACTGGTGCTCGGGATCAGCGGCGGTCAGGACTCCACGCTTGCCGGAAAGCTGTGCCAGATGGCGATTAATGAGCTGCGCCTGGAAACCGGCAACGAATCACTGCAATTTATTGCCGTACGc                                                                        >  1:2569368/1‑150 (MQ=255)
                                                                                 cTGCAAACTTATCCGTTCATTAAATCACTGGTGCTCGGGATCAGCGGCGGTCAGGACTCCACGCTTGCCGGAAAGCTGTGCCAGATGGCGATTAATGAGCTGCGCCTGGAAACCGGCAACGAATCACTGCAATTTATTGCCGTACGCCTg                                                                     <  1:4021075/150‑1 (MQ=255)
                                                                                   gCAAACTTATCCGTTCATTAAATCACTGGTGCTCGGGATCAGCGGCGGTCAGGACTCCACGCTTGCCGGAAAGCTGTGCCAGATGGCGATTAATGAGCTGCGCCTGGAAACCGGCAACGAATCACTGCAATTTATTGCCGTACGCCTGcc                                                                   >  2:3911007/1‑150 (MQ=255)
                                                                                       aCTTATCCGTTCATTAAATCACTGGTGCTCGGGATCAGCGGCGGTCAGGACTCCACGCTTGCCGGAAAGCTGTGCCAGATGGCGATTAATGAGCTGCGCCTGGAAACCGGCAACGAATCACTGCAATTTATTGCCGTACGCCTGCCCTAt                                                               >  1:2773102/1‑150 (MQ=255)
                                                                                         ttATCCGTTCATTAAATCACTGGTGCTCGGGATCAGCGGCGGTCAGGACTCCACGCTTGCCGGAAAGCTGTGCCAGATGGCGATTAATGAGCTGCGCCTGGAAACCGGCAACGAATCACTGCAATTTATTGCCGTACGCCTGCCCTATgg                                                             >  2:2879454/1‑150 (MQ=255)
                                                                                         ttATCCGTTCATTAAATCACTGGTGCTCGGGATCAGCGGCGGTCAGGACTCCACGCTTGCCGGAAAGCTGTGCCAGATGGCGATTAATGAGCTGCGCCTGGAAACCGGCAACGAATCACTGCAATTTATTGCCGTACGCCTGCCCTATgg                                                             >  2:2258935/1‑150 (MQ=255)
                                                                                          tATCCGTTCATTAAATCACTGGTGCTCGGGATCAGCGGCGGTCAGGACTCCACGCTTGCCGGAAAGCTGTGCCAGATGGCGATTAATGAGCTGCGCCTGGAAACCGGCAACGAATCACTGCAATTTATTGCCGTACGCCTGCCCTATGgt                                                            >  2:435703/1‑150 (MQ=255)
                                                                                               gTTCATTAAATCACTGGTGCTCGGGATCAGCGGCGGTCAGGACTCCACGCTTGCCGGAAAGCTGTGCCAGATGGCGATTAATGAGCTGCGCCTGGAAACCGGCAACGAATCACTGCAATTTATTGCCGTACGCCTGCCCTATGGTGTTCa                                                       >  1:2142998/1‑150 (MQ=255)
                                                                                               gTTCATTAAATCACTGGTGCTCGGGATCAGCGGCGGTCAGGACTCCACGCTTGCCGGAAAGCTGTGCCAGATGGCGATTAATGAGCTGCGCCTGGAAACCGGCAACGAATCACTGCAATTTATTGCCGTACGCCTGCCCTATGGTGTTCa                                                       <  2:3962183/150‑1 (MQ=255)
                                                                                                 tCATTAAATCACTGGTGCTCGGGATCAGCGGCGGTCAGTACTCCACGCTTGCCGGAAAACTGTGCCAGATGGCGATTAATGATCTGCGCCTGGAAACCGGCAACGAATCACTGCAATTTATTGCCGTACGCATGCCCTATGGTGTTCAgg                                                     >  2:2536318/1‑150 (MQ=255)
                                                                                                      aaaTCACTGGTGCTCGGGATCAGCGGCGGTCAGGACTCCACGCTTGCCGGAAAGCTGTGCCAGATGGCGATTAATGAGCTGCGCCTGGAAACCGGCAACGAATCACTGCAATTTATTGCCGTACGCCTGCCCTATGGTGTTCAGGCcgac                                                >  1:3697727/1‑150 (MQ=255)
                                                                                                       aaTCACTGGTGCTCGGGATCAGCGGCGGTCAGGACTCCACGCTTGCCGGAAAGCTGTGCCAGATGGCGATTAATGAGCTGCGCCTGGAAACCGGCAACGAATCACTGCAATTTATTGCCGTACGCCTGCCCTATGGTGTTCAGGCcgacg                                               <  1:2858218/150‑1 (MQ=255)
                                                                                                       aaTCACTGGTGCTCGGGATCAGCGGCGGTCAGGACTCCACGCTTGCCGGAAAGCTGTGCCAGATGGCGATTAATGAGCTGCGCCTGGAAACCGGCAACGAATCACTGCAATTTATTGCCGTACGCCTGCCCTATGGTGTTCAGGCcgacg                                               >  1:461062/1‑150 (MQ=255)
                                                                                                       aaTCACTGGTGCTCGGGATCAGCGGCGGTCAGGACTCCACGCTTGCCGGAAAGCTGTGCCAGATGGCGATTAATGAGCTGCGCCTGGAAACCGGCAACGAATCACTGCAATTTATTGCCGTACGCCTGCCCTATGGTGTTCAGGCcgacg                                               <  1:769777/150‑1 (MQ=255)
                                                                                                        aTCACTGGTGCTCGGGATCAGCGGCGGTCAGGACTCCACGCTTGCCGGAAAGCTGTGCCAGATGGCGATTAATGAGCTGCGCCTGGAAACCGGCAACGAATCACTGCAATTTATTGCCGTACGCCTGCCCTATGGTGTTCAGGCcgacga                                              >  2:3334528/1‑150 (MQ=255)
                                                                                                         tCACTGGTGCTCGGGATCAGCGGCGGTCAGGACTCCACGCTTGCCGGAAAGCTGTGCCAGATGGCGATTAATGAGCTGCGCCTGGAAACCGGCAACGAATCACTGCAATTTATTGCCGTACGCCTGCCCTATGGTGTTCAGGCCGACGaa                                             >  2:2787490/1‑150 (MQ=255)
                                                                                                            cTGGTGCTCGGGATCAGCGGCGGTCAGGACTCCACGCTTGCCGGAAAGCTGTGCCAGATGGCGATTAATGAGCTGCGCCTGGAAACCGGCAACGAATCACTGCAATTTATTGCCGTACGCCTGCCCTATGGTGTTCAGGCCGACGAACaa                                          >  1:2319127/1‑150 (MQ=255)
                                                                                                             tGGTGCTCGGGATCAGCGGCGGTCAGGACTCCACGCTTGCAGGAAAGCTGTGCCAGATGGCGATTAATGAGCTGCGCCTGGAAACCGGCAACGAATCACTGCAATTTATTGCCGTACGCCTGCCCTATGGTGTTCAGGCCGACGAACAAg                                         >  1:2131484/1‑150 (MQ=255)
                                                                                                                tGCTCGGGATCAGCGGCGGTCAGGACTCCACGCTTGCCGGAAAGCTGTGCCAGATGGCGATTAATGAGCTGCGCCTGGAAACCGGCAACGAATCACTGCAATTTATTGCCGTACGCCTGCCCTATGGTGTTCAGGCCGACGAACAAGAtt                                      >  1:2194035/1‑150 (MQ=255)
                                                                                                                tGCTCGGGATCAGCGGCGGTCAGGACTCCACGCTTGCCGGAAAGCTGTGCCAGATGGCGATTAATGAGCTGCGCCTGGAAACCGGCAACGAATCACTGCAATTTATTGCCGTACGCCTGCCCTATGGTGTTCAGGCCGACGAACAAGAtt                                      >  1:1323906/1‑150 (MQ=255)
                                                                                                                   tCGGGATCAGCGGCGGTCAGGACTCCACGCTTGCCGGAAAGCTGTGCCAGATGGCGATTAATGAGCTGCGCCTGGAAACCGGCAACGAATCACTGCAATTTATTGCCGTACGCCTGCCCTATGGTGTTCAGGCCGACGAACAAGATTGcc                                   >  1:2404668/1‑150 (MQ=255)
                                                                                                                    cGGGATCAGCGGCGGTCAGGACTCCACGCTTGCCGGAAAGCTGTGCCAGATGGCGATTAATGAGCTGCGCCTGGAAACCGGCAACGAATCACTGCAATTTATTGCCGTACGCCTGCCCTATGGTGTTCAGGCCGACGAACAAGATTGCCa                                  <  2:936810/150‑1 (MQ=255)
                                                                                                                     gggATCAGCGGCGGTCAGGACTCCACGCTTGCCGGAAAGCTGTGCCAGATGGCGATTAATGAGCTGCGCCTGGAAACCGGCAACGAATCACTGCAATTTATTGCCGTACGCCTGCCCTATGGTGTTCAGGCCGACGAACAAGATTGCCAg                                 <  1:2811075/150‑1 (MQ=255)
                                                                                                                     gggATCAGCGGCGGTCAGGACTCCACGCTTGCCGGAAAGCTGTGCCAGATGGCGATTAATGAGCTGCGCCTGGAAACCGGCAACGAATCACTGCAATTTATTGCCGTACGCCTGCCCTATGGTGTTCAGGCCGACGAACAAGATTGCCAg                                 >  1:344353/1‑150 (MQ=255)
                                                                                                                      ggATCAGCGGCGGTCAGGACTCCACGCTTGCCGGAAAGCTGTGCCAGATGGCGATTAATGAGCTGCGCCTGGAAACCGGCAACGAATCACTGCAATTTATTGCCGTACGCCTGCCCTATGGTGTTCAGGCCGACGAACAAGATTGCCAgg                                <  1:1113502/150‑1 (MQ=255)
                                                                                                                             cggcggTCAGGACTCCACGCTTGCCGGAAAGCTGTGCCAGATGGCGATTAATGAGCTGCGCCTGGAAACCGGCAACGAATCACTGCAATTTATTGCCGTACGCCTGCCCTATGGTGTTCAGGCCGACGAACAAGATTGCCAGGATGCCAt                         <  2:2084644/150‑1 (MQ=255)
                                                                                                                              ggcggTCAGGACTCCACGCTTGCCGGAAAGCTGTGCCAGATGGCGATTAATGAGCTGCGCCTGGAAACCGGCAACGAATCACTGCAATTTATTGCCGTACGCCTGCCCTATGGTGTTCAGGCCGACGAACAAGATTGCCAGGATGCCAtt                        <  1:1561197/150‑1 (MQ=255)
                                                                                                                              ggcggTCAGGACTCCACGCTTGCCGGAAAGCTGTGCCAGATGGCGATTAATGAGCTGCGCCTGGAAACCGGCAACGAATCACTGCAATTTATTGCCGTACGCCTGCCCTATGGTGTTCAGGCCGACGAACAAGATTGCCAGGATGCCAtt                        <  1:2682417/150‑1 (MQ=255)
                                                                                                                              ggcggTCAGGACTCCACGCTTGCAGGAAAGCTGTGCCAGATGGCGATTAATGAGCTGCGCCTGGAAACCGGCAACGAATCACTGCAATTTATTGCCGTACGCCTGCCCTATGGTGTTCAGGCCGACGAACAAGATTGCCAGGATGCCAtt                        <  1:4138713/150‑1 (MQ=255)
                                                                                                                               gcggTCAGGACTCCACGCTTGCCGGAAAGCTGTGCCAGATGGCGATTAATGAGCTGCGCCTGGAAACCGGCAACGAATCACTGCAATTTATTGCCGTACGCCTGCCCTATGGTGTTCAGGCCGACGAACAAGATTGCCAGGATGCCATTg                       <  1:2267845/150‑1 (MQ=255)
                                                                                                                                   tCATGACTCCACGCTTGCCGGAAAGCTGTGCCAGATGGCGATTAATGAGCTGCGCCTGGAAAGCGGCAACGAATCACTGCAATTTATTGCCGTACGCCTGCCCTATGGTGTTCAGGCCGACGAACAAGATTGCCAGGATGCCATTGCCtt                   <  1:1421500/150‑1 (MQ=255)
                                                                                                                                   tCAGGACTCCACGCTTGCCGGAAAGCTGTGCCAGATGGCGATTAATGAGCTGCGCCTGGAAACCGGCAACGAATCACTGCAATTTATTGCCGTACGCCTGCCCTATGGTGTTCAGGCCGACGAACAAGATTGCCAGGATGCCATTGCCtt                   >  1:1133601/1‑150 (MQ=255)
                                                                                                                                   tCAGGACTCCACGCTTGCCGGAAAGCTGTGCCAGATGGCGATTAATGAGCTGCGCCTGGAAACCGGCAACGAATCACTGCAATTTATTGCCGTACGCCTGCCCTATGGTGTTCAGGCCGACGAACAAGATTGCCAGGATGCCATTGCCtt                   <  1:2965780/150‑1 (MQ=255)
                                                                                                                                    cAGGACTCCACGCTTGCCGGAAAGCTGTGCCAGATGGCGATTAATGAGCTGCGCCTGGAAACCGGCAACGAATCACTGCAATTTATTGCCGTACGCCTGCCCTATGGTGTTCAGGCCGACGAACAAGATTGCCAGGATGCCATTGCCttt                  <  2:2938468/150‑1 (MQ=255)
                                                                                                                                    cAGGACTCCACGCTTGCAGGAAAGCTGTGCCAGATGGCGATTAATGAGCTGCGCCTGGAAACCGGCAACGAATCACTGCAATTTATTGCCGTACGCCTGCCCTATGGTGTTCAGGCCGACGAACAAGATTGCCAGGATGCCATTGCCttt                  <  1:1524184/150‑1 (MQ=255)
                                                                                                                                      ggACTCCACGCTTGCCGGAAAGCTGTGCCAGATGGCGATTAATGAGCTGCGCCTGGAAACCGGCAACGAATCACTGCAATTTATTGCCGTACGCCTGCCCTATGGTGTTCAGGCCGACGAACAAGATTGCCAGGATGCCATTGCCTTTAt                <  1:41106/150‑1 (MQ=255)
                                                                                                                                        aCTCCACGCTTGCCGGAAAGCTGTGCCAGATGGCGATTAATGAGCTGCGCCTGGAAACCGGCAACGAATCACTGCAATTTATTGCCGTACGCCTGCCCTATGGTGTTCAGGCCGACGAACAAGATTGCCAGGATGCCATTGCCTTTATTc              <  1:1472669/150‑1 (MQ=255)
                                                                                                                                             aCGCTTGCCGGAAAGGTGTGCCAGATGGCGATTAATGAGCTGCGTCTGGAAACCGGCTACGAATCTCTGCTTTTTTTTGCCGTACGCCTGACCTATGGTGTTCAGGGCGTCGAACAAGATTGCCAGGATGACATTGCCTTTATTCAACCg         <  1:3001913/150‑1 (MQ=255)
                                                                                                                                               gCTTGCCGGAAAGCTGTGCCAGATGGCGATTAATGAGCTGCGCCTGGAAACCGGCAACGAATCACTGCAATTTATTGCCGTACGCCTGCCCTATGGTGTTCAGGCCGACGAACAAGATTGCCAGGATGCCATTGCCTTTATTCAACCGGa       <  2:3874821/150‑1 (MQ=255)
                                                                                                                                               gCTTGCCGGAAAGCTGTGCCAGATGGCGATTAATGAGCTGCGCCTGGAAACCGGCAACGAATCACTGCAATTTATTGCCGTACGCCTGCCCTATGGTGTTCAGGCCGACGAACAAGATTGCCAGGATGCCATTGCCTTTATTCAACCGGa       <  1:3396694/150‑1 (MQ=255)
                                                                                                                                                 ttGCCGGAAAGCTGTGCCAGATGGCGATTAATGAGCTGCGCCTGGAAACCGGCAACGAATCACTGCAATTTATTGCCGTACGCCTGCCCTATGGTGTTCAGGCCGACGAACAAGATTGCCAGGATGCCATTGCCTTTATTCAACCGGATc     >  2:2652088/1‑150 (MQ=255)
                                                                                                                                                 ttGCCGGAAAGCTGTGCCAGATGGCGATTAATGAGCTGCGCCTGGAAACCGGCAACGAATCACTGCAATTTATTGCCGTACGCCTGCCCTATGGTGTTCAGGCCGACGAACAAGATTGCCAGGATGCCATTGCCTTTATTCAACCGGATc     <  2:3408939/150‑1 (MQ=255)
                                                                                                                                                 ttGCCGGAAAGCTGTGCCAGATGGCGATTAATGAGCTGCGCCTGGAAACCGGCAACGAATCACTGCAATTTATTGCCGTACGCCTGCCCTATGGTGTTCAGGCCGACGAACAAGATTGCCAGGATGCCATTGCCTTTATTCAACCGGATc     <  1:638878/150‑1 (MQ=255)
                                                                                                                                                  tGCCGGAAAGCTGTGCCAGATGGCGATTAATGAGCTGCGCCTGGAAACCGGCAACGAATCACTGCAATTTATTGCCGTACGCCTGCCCTATGGTGTTCAGGCCGACGAACAAGATTGCCAGGATGCCATTGCCTTTATTCAACCGGATcg    <  2:808148/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
ATAATAAAGGCGCTGGGCGCAAAACCGCAGATTAATGCTGAAGAGGAAATTCGTCGTAGTGTCGATTTTCTGAAAAGCTACCTGCAAACTTATCCGTTCATTAAATCACTGGTGCTCGGGATCAGCGGCGGTCAGGACTCCACGCTTGCCGGAAAGCTGTGCCAGATGGCGATTAATGAGCTGCGCCTGGAAACCGGCAACGAATCACTGCAATTTATTGCCGTACGCCTGCCCTATGGTGTTCAGGCCGACGAACAAGATTGCCAGGATGCCATTGCCTTTATTCAACCGGATCGCG  >  CP009273/1816733‑1817030

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: