Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,084,305 G→T 93.3% K200N (AAG→AAT endA → DNA‑specific endonuclease I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,084,3050GT93.3% 490.8 / 27.8 178K200N (AAG→AATendADNA‑specific endonuclease I
Reads supporting (aligned to +/- strand):  ref base G (7/5);  new base T (80/86);  total (87/91)
Fisher's exact test for biased strand distribution p-value = 5.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.12e-01

Only 100 of 180 total aligned reads displayed.

CCATGAAGGTCGATTTCAAAGAAAAAGCTGCCGAACCACCAGCGCGTGCACGCGGTGCCATTGCGCGCACCTACTTCTATATGCGCGACCAATACAACCTGACACTCTCTCGCCAGCAAACGCAGCTGTTCAACGCATGGAACAAGATGTATCCGGTTACCGACTGGGAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATAACCCGTATGTGCAACGCGCTTGCCAGGCGCGAAAGAGCTAACCTACACTAGCGGGATTCTTTTTGTTAACCCCTACC  >  CP009273/3084160‑3084449
                                                                                                                                                 |                                                                                                                                                  
ccATGAAGGTCGATTTCAAAGAAAAAGCTGCCGAACCACCAGCGCGTGCACGCGGTGCCATTGCGCGCACCTACTTCTATATGCGCGACCAATACAACCTGACACTCTCTCGCCAGCAAACGCAGCTGTTCAACGCATGGAACAatat                                                                                                                                                  <  1:3324941/148‑1 (MQ=255)
ccATGAAGGTCGATTTCAAAGAAAAAGCTGCCGAACCACCAGCGCGTGCACGCGGTGCCATTGCGCGCACCTACTTCTATATGCGCGACCAATACAACCTGACACTCTCTCGCCAGCAAACGCAGCTGTTCAACGCATGGAACAatat                                                                                                                                                  >  1:708322/1‑148 (MQ=255)
   tGAAGGTCGATTTCAAAGAAAAAGCTGCCGAACCACCAGCGCGTGCACGCGGTGCCATTGCGCGCACCTACTTCTATATGCGCGACCAATACAACCTGACACTCTCTCGCCAGCAAACGCAGCTGTTCAACGCATGGAACAATATGTa                                                                                                                                               >  2:2006432/1‑148 (MQ=255)
   tGAAGGTCGATTTCAAAGAAAAAGCTGCCGAACCACCAGCGCGTGCACGCGGTGCCATTGCGCGCACCTACTTCTATATGCGCGACCAATACAACCTGACACTCTCTCGCCAGCAAACGCAGCTGTTCAACGCATGGAACAATATGTa                                                                                                                                               >  2:1366296/1‑148 (MQ=255)
   tGAAGGTCGATTTCAAAGAAAAAGCTGCCGAACCACCAGCGCGTGCACGCGGTGCCATTGCGCGCACCTACTTCTATATGCGCGACCAATACAACCTGACACTCTCTCGCCAGCAAACGCAGCTGTTCAACGCATGGAACAATATGTa                                                                                                                                               <  1:2848414/148‑1 (MQ=255)
       ggTCGATTTCAAAGAAAAAGCTGCCGAACCACCAGCGCGTGCACGCGGTGCCATTGCGCGCACCTACTTCTATATGCGCGACCAATACAACCTGACACTCTCTCGCCAGCAAACGCAGCTGTTCAACGCATGGAACAATATGTATCCg                                                                                                                                           >  1:2234550/1‑148 (MQ=255)
           gATTTCAAAGAAAAAGCTGCCGAACCACCAGCGCGTGCACGCGGTGCCATTGCGCGCACCTACTTCTATATGCGCGACCAATACAACCTGACACTCTCTCGCCAGCAAACGCAGCTGTTCAACGCATGGAACAATATGTATCCGGTTa                                                                                                                                       <  1:556004/148‑1 (MQ=255)
           gATTTCAAAGAAAAAGCTGCCGAACCACCAGCGCGTGCACGCGGTGCCATTGCGCGCACCTACTTCTATATGCGCGACCAATACAACCTGACACTCTCTCGCCAGCAAACGCAGCTGTTCAACGCATGGAACAATATGTATCCGGTTa                                                                                                                                       <  2:2700250/148‑1 (MQ=255)
             tttCAAAGAAAAAGCTGCCGAACCACCAGCGCGTGCACGCGGTGCCATTGCGCGCACCTACTTCTATATGCGCGACCAATACAACCTGACACTCTCTCGCCAGCAAACGCAGCTGTTCAACGCATGGAACAATATGTATCCGGTTAcc                                                                                                                                     >  1:976867/1‑148 (MQ=255)
              ttCAAAGAAAAAGCTGCCGAACCACCAGCGCGTGCACGCGGTGCCATTGCGCGCACCTACTTCTATATGCGCGACCAATACAACCTGACACTCTCTCGCCAGCAAACGCAGCTGTTCAACGCATGGAACAATATGTATCCGGTTACCg                                                                                                                                    >  2:585005/1‑148 (MQ=255)
               tCAAAGAAAAAGCTGCCGAACCACCAGCGCGTGCACGCGGTGCCATATCGCGCACCTCCTTCTATATGCGCCACCAATACAACCTGACACTCTATCGCCCGCAAACCCAGCTGTTCAACGCATGGAACAATATGAAGCCGCTTACCGa                                                                                                                                   >  1:2203726/1‑148 (MQ=255)
                 aaaGAAAAAGCTGCCGAACCACCAGCGCGTGCACGCGGTGCCATTGCGCGCACCTACTTCTATATGCGCGACCAATACAACCTGACACTCTCTCGCCAGCAAACGCAGCTGTTCAACGCATGGAACAATATGTATCCGGTTACCGACt                                                                                                                                 <  2:2347037/148‑1 (MQ=255)
                 aaaGAAAAAGCTGCCGAACCACCAGCGCGTGCACGCGGTGCCATTGCGCGCACCTACTTCTATATGCGCGACCAATACAACCTGACACTCTCTCGCCAGCAAACGCAGCTGTTCAACGCATGGAACAATATGTATCCGGTTACCGACt                                                                                                                                 <  2:624696/148‑1 (MQ=255)
                  aaGAAAAAGCTGCCGAACCACCAGCGCGTGCACGCGGTGCCATTGCGCGCACCTACTTCTATATGCGCGACCAATACAACCTGACACTCTCTCGCCAGCAAACGCAGCTGTTCAACGCATGGAACAATATGTATCCGGTAACCGACTg                                                                                                                                >  1:343475/1‑148 (MQ=255)
                   aGAAAAAGCTGCCGAACCACCAGCGCGTGCACGCGGTGCCATTGCGCGCACCTACTTCTATATGCGCGACCAATACAACCTGACACTCTCTCGCCAGCAAACGCAGCTGTTCAACGCATGGAACAATATGTATCCGGTTACCGACTgg                                                                                                                               <  2:2113257/148‑1 (MQ=255)
                    gAAAAAGCTGCCGAACCACCAGCGCGTGCACGCGGTGCCATTGCGCGCACCTACTTCTATATGCGCGACCAATACAACCTGACACTCTCTCGCCAGCAAACGCAGCTGTTCAACGCATGGAACAATATGTATCCGGTTACCGACTggg                                                                                                                              >  1:3305729/1‑148 (MQ=255)
                    gAAAAAGCTGCCGAACCACCAGCGCGTGCACGCGGTGCCATTGCGCGCACCTACTTCTATATGCGCGACCAATACAACCTGACACTCTCTCGCCAGCAAACGCAGCTGTTCAACGCATGGAACAATATGTATCCGGTTACCGACTggg                                                                                                                              >  2:104310/1‑148 (MQ=255)
                      aaaaGCTGCCGAACCACCAGCGCGTGCACGCGGTGCCATTGCGCGCACCTACTTCTATATGCGCGACCAATACAACCTGACACTCTCTCGCCAGCAAACGCAGCTGTTCAACGCATGGAACAATATGTATCCGGTTACCGACTGGGAg                                                                                                                            <  2:3219580/148‑1 (MQ=255)
                        aaGCTGCCGAACCACCAGCGCGTGCACGCGGTGCCATTGCGCGCACCTACTTCTATATGCGCGACCAATACAACCTGACACTCTCTCGCCAGCAAACGCAGCTGTTCAACGCATGGAACAATATGTATCCGGTTACCGACTGGGAGTg                                                                                                                          <  2:2927053/148‑1 (MQ=255)
                         aGCTGCCGAACCACCAGCGCGTGCACGCGGTGCCATTGCGCGCACCTACTTCTATATGCGCGACCAATACAACCTGACACTCTCTCGCCAGCAAACGCAGCTGTTCAACGCATGGAACAATATGTATCCGGTTACCGACTGGGAGTGc                                                                                                                         >  2:1210524/1‑148 (MQ=255)
                           cTGCCGAACCACCAGCGCGTGCACGCGGTGCCATTGCGCGCACCTACTTCTATATGCGCGACCAATACAACCTGACACTCTCTCGCCAGCAAACGCAGCTGTTCAACGCATGGAACAATATGTATCCGGTTACCGACTGGGAGTGCGa                                                                                                                       <  2:2132551/148‑1 (MQ=255)
                               cGAACCACCAGCGCGTGCACGCGGTGCCATTGCGCGCACCTACTTCTATATGCGCGACCAATACAACCTGACACTCTCTCGCCAGCAAACGCAGCTGTTCAACGCATGGAACAATATGTATCCGGTTACCGACTGGGAGTGCGAgcgc                                                                                                                   <  1:3235363/148‑1 (MQ=255)
                               cGAACCACCAGCGCGTGCACGCGGTGCCATTGCGCGCACCTACTTCTATATGCGCGACCAATACAACCTGACACTCTCTCGCCAGCAAACGCAGCTGTTCAACGCATGGAACAATATGTATCCGGTTACCGACTGGGAGTGCGAgcgc                                                                                                                   <  1:3235348/148‑1 (MQ=255)
                                 aaCCACCAGCGCGTGCACGCGGTGCCATTGCGCGCACCTACTTCTATATGCGCGACCAATACAACCTGACACTCTCTCGCCAGCAAACGCAGCTGTTCAACGCATGGAACAATATGTATCCGGTTACCGACTGGGAGTGCGAGCGCGa                                                                                                                 >  1:1119235/1‑148 (MQ=255)
                                  accaccAGCGCGTGCACGCGGTGCCATTGCGCGCACCTACTTCTATATGCGCGACCAATACAACCTGACACTCTCTCGCCAGCAAACGCAGCTGTTCAACGCATGGAACAATATGTATCCGGTTACCGACTGGGAGTGCGAGCGCGAt                                                                                                                >  2:1453850/1‑148 (MQ=255)
                                  accaccAGCGCGTGCACGCGGTGCCATTGCGCGCACCTACTTCTATATGCGCGACCAATACAACCTGACACTCTCTCGCCAGCAAACGCAGCTGTTCAACGCATGGAACAATATGTATCCGGTTACCGACTGGGAGTGCGAGCGCGAt                                                                                                                >  2:2305842/1‑148 (MQ=255)
                                     accaGCGCGTGCACGCGGTGCCATTGCGCGCACCTACTTCTATATGCGCGACCAATACAACCTGACACTCTCTCGCCAGCAAACGCAGCTGTTCAACGCATGGAACAATATGTATCCGGTTACCGACTGGGAGTGCGAGCGCGATGaa                                                                                                             >  2:3337498/1‑148 (MQ=255)
                                       caGCGCGTGCACGCGGTGCCATTGCGCGTACCTACTTCTATCTGCGCGACCAATACAACCTGACACTCTCTCGCCAGCAAACGCAGCTGTTCAACGCATGTAACAATATGTATCCGGTTACCGACTGGGAGTGCGAGCGCGATGAACg                                                                                                           <  1:2586915/148‑1 (MQ=255)
                                       caGCGCGTGCACGCGGTGCCATTGCGCGCACCTACTTCTATATGCGCGACCAATACAACCTGACACTCTCTCGCCAGCAAACGCAGCTGTTCAACGCATGGAACAATCTGTATCCGGTTACCGACTGGGAGTGCGAGCGCGATGAACg                                                                                                           >  1:77041/1‑148 (MQ=255)
                                       caGCGCGTGCACGCGGTGCCATTGCGCGCACCTACTTCTATATGCGCGACCAATACAACCTGACACTCTCTCGCCAGCAAACGCAGCTGTTCAACGCATGGAACAATATGTATCCGGTTACCGACTGGGAGTGCGAGCGCGATGAACg                                                                                                           >  1:1747873/1‑148 (MQ=255)
                                          cgcgTGCACGCGGTGCCATTGCGCGCACCTACTTCTATATGCGCGACCAATACAACCTGACACTCTCTCGCCAGCAAACGCAGCTGTTCAACGCATGGAACAATATGTATCCGGTTACCGACTGGGAGTGCGAGCGCGATGAACGCAt                                                                                                        >  1:1991263/1‑148 (MQ=255)
                                               gCACGCGGTGCCATTGCGCGCACCTACTTCTATATGCGCGACCAATACAACCTGACACTCTCTCGCCAGCAAACGCAGCTGTTCAACGCATGGAACAATATGTATCCGGTTACCGACTGGGAGTGCGAGCGCGATGAACGCATCGCGa                                                                                                   <  2:724931/148‑1 (MQ=255)
                                                cACGCGGTGCCAGTGCGCGCACCTACTTCTATAGGCGCGACCAATACAACCTGACACTCTCTCGCCAGCAAACGCAGCTGTTCAACGCATGGAACAATATGTATCCGGTTACCGACTGGGAGTGCGAGCGCGATGAACGCATCGCGaa                                                                                                  <  1:716097/148‑1 (MQ=255)
                                                      gTGCCATTGCGCGCACCTACTTCTATATGCGCGACCAATACAACCTGACACTCTCTCGCCAGCAAACGCAGCTGTTCAACGCATGGAACAATATGTATCCGGTTACCGACTGGGAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCa                                                                                            <  1:2870545/148‑1 (MQ=255)
                                                      gTGCCATTGCGCGCACCTACTTCTATATGCGCGACCAATACAACCTGACACTCTCTCGCCAGCAAACGCAGCTGTTCAACGCATGGAACAATATGTATCCGGTTACCGACTGGGAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCa                                                                                            <  2:2275872/148‑1 (MQ=255)
                                                      gTGCCATTGCGCGCACCTACTTCTATATGCGCGACCAATACAACCTGACACTCTCTCGCCAGCAAACGCAGCTGTTCAACGCATGGAACAATATGTATCCGGTTACCGACTGGGAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCa                                                                                            <  2:663863/148‑1 (MQ=255)
                                                       tGCCATTGCGCGCACCTACTTCTATATGCGCGACCAATACAACCTGACACTCTCTCGCCAGCAAACGCAGCTGTTCAACGCATGGAACAATATGTATCCGGTTACCGACTGGGAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAg                                                                                           >  1:2941538/1‑148 (MQ=255)
                                                         ccATTGCGCGCACCTACTTCTATATGCGCGACCAATACAACCTGACACTCTCTCGCCAGCAAACGCAGCTGTTCAACGCATGGAACAATATGTATCCGGTTACCGACTGGGAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAggg                                                                                         >  1:2356956/1‑148 (MQ=255)
                                                         ccATTGCGCGCACCTACTTCTATATGCGCGACCAATACAACCTGACACTCTCTCGCCAGCAAACGCAGCTGTTCAACGCATGGAACAATATGTATCCGGTTACCGACTGGGAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAggg                                                                                         >  1:2356953/1‑148 (MQ=255)
                                                         ccATTGCGCGCACCTACTTCTATATGCGCGACCAATACAACCTGACACTCTCTCGCCAGCAAACGCAGCTGTTCAACGCATGGAACAATATGTATCCGGTTACCGACTGGGAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAggg                                                                                         >  1:2356925/1‑148 (MQ=255)
                                                           aTTGCGCGCACCTACTTCTATATGCGCGACCAATACAACCTGACACTCTCTCGCCAGCAAACGCAGCTGTTCAACGCATGGAACAATATGTATCCGGTTACCGACTGGGAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCa                                                                                       <  2:3428258/148‑1 (MQ=255)
                                                              gcgcgcACCTACTTCTATATGCGCGACCAATACAACCTGACACTCTCTCGCCAGCAAACGCAGCTGTTCAACGCATGGAACAATATGTATCCGGTTACCGACTGGGAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATc                                                                                    <  1:1429080/148‑1 (MQ=255)
                                                               cgcgcACCTACTTCTATATGCGCGACCAATACAACCTGACACTCTCTCGCCAGCAAACGCAGCTGTTCAACGCATGGAACAATATGTATCCGGTTACCGACTGGGAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCa                                                                                   <  2:2039523/148‑1 (MQ=255)
                                                                 cgcACCTACTTCTATATGCGCGACCAATACAACCTGACACTCTCTCGCCAGCAAACGCAGCTGTTCAACGCATGGAACAATATGTATCCGGTTACCGACTGGGAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATa                                                                                 >  2:2609070/1‑148 (MQ=255)
                                                                    aCCTACTTCTATATGCGCGACCAATACAACCTGACACTCTCTCGCCAGCAAACGCAGCTGTTCAACGCATGGAACAATATGTATCCGGTTACCGACTGGGAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATAAcc                                                                              <  2:1802651/148‑1 (MQ=255)
                                                                     ccTACTTCTATATGCGCGACCAATACAACCTGACACTCTCTCGCCAGCAAACGCAGCTGTTCAACGCATGGAACAATATGTATCCGGTTACCGACTGGGAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATAAccc                                                                             >  1:105589/1‑148 (MQ=255)
                                                                      cTACGTCTATATGCGCGACCAATACAACCTGACACTCTCTCGCCAGCAAACGCAGCTGTTCAACGCATGGAACAATATGTATCCGGTTACCGACTGGGAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATAACCCg                                                                            <  2:1106468/148‑1 (MQ=255)
                                                                       tACTTCTATATGCGCGACCAATACAACCTGACACTCTCTCGCCAGCAAACGCAGCTGTTCAACGCATGGAACAATATGTATCCGGTTACCGACTGGGAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATAACCCGt                                                                           <  2:167095/148‑1 (MQ=255)
                                                                        aCTTCTATATGCGCGACCAATACAACCTGACACTCTCTCGCCAGCAAACGCAGCTGTTCAACGCATGGAACAATATGTATCCGGTTACCGACTGGGAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATAACCCGTa                                                                          <  1:2192151/148‑1 (MQ=255)
                                                                             tataTGCGCGACCAATACAACCTGACACTCTCTCGCCAGCAAACGCAGCTGTTCAACGCATGGAACAATATGTATCCGGTTACCGACTGGGAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATAACCCGTATGTGc                                                                     <  2:3507801/148‑1 (MQ=255)
                                                                              atatGCGCGACCAATACAACCTGACACTCTCTCGCCAGCAAACGCAGCTGTTCAACGCATGGAACAATATGTATCCGGTTACCGACTGGGAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATAACCCGTATGTGCa                                                                    >  2:1552658/1‑148 (MQ=255)
                                                                               tatGCGCGACCAATACAACCTGACACTCTCTCGCCAGCAAACGCAGCTGTTCAACGCATGGAACAATATGTATCCGGTTACCGACTGGGAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATAACCCGTATGTGCaa                                                                   >  1:204388/1‑148 (MQ=255)
                                                                               tatGCGCGACCAATACAACCTGACACTCTCTCGCCAGCAAACGCAGCTGTTCAACGCATGGAACAATATGTATCCGGTTACCGACTGGGAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATAACCCGTATGTGCaa                                                                   <  1:1265968/148‑1 (MQ=255)
                                                                                 tGCGCGACCAATACAACCTGACACTCTCTCGCCAGCAAACGCAGCTGTTCAACGCATGGAACAATATGTATCCGGTTACCGACTGGGAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATAACCCGTATGTGCAAcg                                                                 >  2:902837/1‑148 (MQ=255)
                                                                                 tGCGCGACCAATACAACCTGACACTCTCTCGCCAGCAAACGCAGCTGTTCAACGCATGGAACAATATGTATCCGGTTACCGACTGGGAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATAACCCGTATGTGCAAcg                                                                 <  2:113610/148‑1 (MQ=255)
                                                                                   cgcgACCAATACAACCTGACACTCTCTCGCCAGCAAACGCAGCTGTTCAACGCATGGAACAATATGTATCCGGTTACCGACTGGGAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATAACCCGTATGTGCAAcgcg                                                               <  2:240325/148‑1 (MQ=255)
                                                                                    gcgACCAATACAACCTGACACTCTCTCGCCAGCAAACGCAGCTGTTCAACGCATGGAACAATATGTATCCGGTTACCGACTGGGAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATAACCCGTATGTGCAACgcgc                                                              <  1:78806/148‑1 (MQ=255)
                                                                                    gcgACCAATACAACCTGACACTCTCTCGCCAGCAAACGCAGCTGTTCAACGCATGGAACAATATGTATCCGGTTACCGACTGGGAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATAACCCGTATGTGCAACgcgc                                                              <  2:688507/148‑1 (MQ=255)
                                                                                    gcgACCAATACAACCTGACACTCTCTCGCCAGCAAACGCAGCTGTTCAACGCATGGAACAATATGTATCCGGTTACCGACTGGGAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATAACCCGTATGTGCAACgcgc                                                              <  2:688474/148‑1 (MQ=255)
                                                                                     cgACCAATACAACCTGACACTCTCTCGCCAGCAAACGCAGCTGTTCAACGCATGGAACAATATGTATCCGGTTACCGACTGGGAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATAACCCGTATGTGCAACGCGCt                                                            <  2:2973794/148‑1 (MQ=255)
                                                                                     cgACCAATACAACCTGACACTCTCTCGCCAGCAAACGCAGCTGTTCAACGCATGGAACAAGATGTATCCGGTTACCGACTGGGAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATAACCCGTATGTGCAACGCGCt                                                            >  2:347550/1‑148 (MQ=255)
                                                                                         cAATACAACCTGACACTCTCTCGCCAGCAAACGCAGCTGTTCAACGCATGGAACAAGATGTATCCGGTTACCGACTGGGAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATAACCCGTATGTGCAACGCGCTTGcc                                                        >  1:339270/1‑148 (MQ=255)
                                                                                            tACAACCTGACACTCTCTCGCCAGCAAACGCAGCTGTTCAACGCATGGAACAAGATGTATCCGGTTACCGACTGGGAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATAACCCGTATGTGCAACGCGCTTGCCAgg                                                     >  1:541803/1‑148 (MQ=255)
                                                                                               aaCCTGACACTCTCTCGCCAGCAAACGCAGCTGTTCAACGCATGGAACAATATGTATCCGGTTACCGACTGGGAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATAACCCGTATGTGCAACGCGCTTGCCAGgcgc                                                  >  1:471187/1‑148 (MQ=255)
                                                                                                aCCTGACACTCTCTCGCCAGCAAACGCAGCTGTTCAACGCATGGAACAATATGTATCCGGTTACCGACTGGGAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATAACCCGTATGTGCAACGCGCTTGCCAGGcgcg                                                 >  2:2267743/1‑148 (MQ=255)
                                                                                                    gACACTCTCTCGCCAGCAAACGCAGCTGTTCAACGCATGGAACAATATGTATCCGGTTACCGACTGGGAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATAACCCGTATGTGCAACGCGCTTGCCAGGCGCGAAag                                             >  2:1104135/1‑148 (MQ=255)
                                                                                                    gACACTCTCTCGCCAGCAAACGCAGCTGTTCAACGCATGGAACAATATGTATCCGGTTACCGACTGGGAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATAACCCGTATGTGCAACGCGCTTGCCAGGCGCGAAag                                             <  1:61583/148‑1 (MQ=255)
                                                                                                    gACACTCTCTCGCCAGCAAACGCAGCTGTTCAACGCATGGAACAATATGTATCCGGTTACCGACTGGGAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATAACCCGTATGTGCAACGCGCTTGCCAGGCGCGAAag                                             >  2:465760/1‑148 (MQ=255)
                                                                                                     acacTCTCTCGCCAGCAAACGCAGCTGTTCAACGCATGGAACAATATGTATCCGGTTACCGACTGGGAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATAACCCGTATGTGCAACGCGCTTGCCAGGCGCGAAaga                                            <  2:1964530/148‑1 (MQ=255)
                                                                                                      cacTCTCTCGCCAGCAAACGCAGCTGTTCAACGCATGGAACAATATGTATCCGGTTACCGACTGGGAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATAACCCGTATGTGCAACGCGCTTGCCAGGCGCGAAagag                                           <  2:1709486/148‑1 (MQ=255)
                                                                                                       acTCTCTCGCCAGCAAACGCAGCTGTTCAACGCATGGAACAATATGTATCCGGTTACCGACTGGGAGTGCGAGCGCGATGAACGCATCGCGACGGTGCAGGGCAATCATAACCCGTATGTGCAACGCGCTTGCCAGGCGCGAAAGAGc                                          >  1:3180767/1‑148 (MQ=255)
                                                                                                       acTCTCTCGCCAGCAAACGCAGCTGTTCAACGCATGGAACAAGATGTATCCGGTTACCGACTGGGAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATAACCCGTATGTGCAACGCGCTTGCCAGGCGCGAAAGAGc                                          <  1:1068057/148‑1 (MQ=255)
                                                                                                          ctctcGCCAGCAAACGCAGCTGTTCAACGCATGGAACAATATGTATCCGGTTACCGACTGGGAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATAACCCGTATGTGCAACGCGCTTGCCAGGCGCGAAAGAGCTaa                                       >  1:2195702/1‑148 (MQ=255)
                                                                                                                ccAGCAAACGCAGCTGTTCAACGCATGGAACAATATGTATCCGGTTACCGACTGGGAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATAACCCGTATGTGCAACGCGCTTGCCAGGCGCGAAAGAGCTCACCTaca                                 >  2:1387227/1‑148 (MQ=255)
                                                                                                                ccAGCAAACGCAGCTGTTCAACGCATGGAACAATATGTATCCGGTTACCGACTGGGAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATAACCCGTATGTGCAACGCGCTTGCCAGGCGCGAAAGAGCTAACCTaca                                 <  1:2641910/148‑1 (MQ=255)
                                                                                                                  aGCAAACGCAGCTGTTCAACGCATGGAACAATATGTATCCGGTTACCGACTGGGAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATAACCCGTATGTGCAACGCGCTTGCCAGGCGCGAAAGAGCTAACCTACACt                               <  1:2346799/148‑1 (MQ=255)
                                                                                                                   gCAAACGCAGCTGTTCAACGCATGGAACAATATGTATCCGGTTACCGACTGGGAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATAACCCGTATGTGCAACGCGCTTGCCAGGCGCGAAAGAGCTAACCTACACTa                              >  1:2780791/1‑148 (MQ=255)
                                                                                                                   gCAAACGCAGCTGTTCAACGCATGGAACAATATGTATCCGGTTACCGACTGGGAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATAACCCGTATGTGCAACGCGCTTGCCAGGCGCGAAAGAGCTAACCTACACTa                              >  1:2780825/1‑148 (MQ=255)
                                                                                                                       aCGCAGCTGTTCAACGCATGGAACAATATGTATCCGGTTACCGACTGGGAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATAACCCGTATGTGCAACGCGCTTGCCAGGCGCGAAAGAGCTAACCTACACTAGCgg                          >  1:1081780/1‑148 (MQ=255)
                                                                                                                       aCGCAGCTGTTCAACGCATGGAACAATATGTATCCGGTTACCGACTGGGAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATAACCCGTATGTGCAACGCGCTTGCCAGGCGCGAAAGAGCTAACCTACACTAGCgg                          >  1:3400996/1‑148 (MQ=255)
                                                                                                                       aCGCAGCTGTTCAACGCATGGAACAATATGTATCCGGTTACCGACTGGGAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATAACCCGTATGTGCAACGCGCTTGCCAGGCGCGAAAGAGCTAACCTACACTAGCgg                          >  1:3401026/1‑148 (MQ=255)
                                                                                                                       aCGCAGCTGTTCAACGCATGGAACAATATGTATCCGGTTACCGACTGGGAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATAACCCGTATGTGCAACGCGCTTGCCAGGCGCGAAAGAGCTAACCTACACTAGCgg                          >  1:912364/1‑148 (MQ=255)
                                                                                                                         gCAGCTGTTCAACGCATGGAACAATATGTATCCGGTTACCGACTGGGAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATAACCCGTATGTGCAACGCGCTTGCCAGGCGCGAAAGAGCTAACCTACACTAGCGGGa                        <  1:1996805/148‑1 (MQ=255)
                                                                                                                         gCAGCTGTTCAACGCATGGAACAATATGTATCCGGTTACCGACTGGGAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATAACCCGTATGTGCAACGCGCTTGCCAGGCGCGAAAGAGCTAACCTACACTAGCGGGa                        <  1:2305842/148‑1 (MQ=255)
                                                                                                                          cAGCTGTTCAACGCATGGAACAATATGTATCCGGTTACCGACTGGGAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATAACCCGTATGTGCAACGCGCTTGCCAGGCGCGAAAGAGCTAACCTACACTAGCGGGAt                       <  1:284167/148‑1 (MQ=255)
                                                                                                                           aGCTGTTCAACGCATGGAACAATATGTATCCGGTTACCGACTGGGAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATAACCCGTATGTGCAACGCGCTTGCCAGGCGCGAAAGAGCTAACCTACACTAGCGGGAtt                      <  2:1468916/148‑1 (MQ=255)
                                                                                                                            gCTGTTCAACGCATGGAACAATATGTATCCGGTTACCGACTGGGAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATAACCCGTATGTGCAACGCGCTTGCCAGGCGCGAAAGAGCTAACCTACACTAGCGGGATTc                    >  1:3049445/1‑148 (MQ=255)
                                                                                                                              tGTTCAACGCATGGAACAATATGTATCCGGTTACCGACTGGGAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATAACCCGTATGTGCAACGCGCTTGCCAGGCGCGAAAGAGCTAACCTACACTAGCGGGATTCtt                  <  1:3398430/148‑1 (MQ=255)
                                                                                                                               gTTCAACGCATGGAACAATATGTATCCGGTTACCGACTGGGAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATAACCCGTATGTGCAACGCGCTTGCCAGGCGCGAAAGAGCTAACCTACACTAGCGGGATTCttt                 <  2:2234550/148‑1 (MQ=255)
                                                                                                                                ttCAACGCATGGAACAATATGTATCCGGTTACCGACTGGGAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATAACCCGTATGTGCAACGCGCTTGCCAGGCGCGAAAGAGCTAACCTACACTAGCGGGATTCtttt                >  1:2965939/1‑148 (MQ=255)
                                                                                                                                     cGCATGGAACAATATGTATCCGGTTACCGACTGGGAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATAACCCGTATGTGCAACGCGCTTGCCAGGCGCGAAAGAGCTAACCTACACTAGCGGGATTCTTTTTGTTa           >  1:115762/1‑148 (MQ=255)
                                                                                                                                     cGCATGGAACAATATGTATCCGGTTACCGACTGGGAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATAACCCGTATGTGCAACGCGCTTGCCAGGCGCGAAAGAGCTAACCTACACTAGCGGGATTCTTTTTGTTa           <  1:585005/148‑1 (MQ=255)
                                                                                                                                        aTGGAACAATATGTATCCGGTTACCGACTGGGAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATAACCCGTATGTGCAACGCGCTTGCCAGGCGCGAAAGAGCTAACCTACACTAGCGGGATTCTTTTTGTTAAcc        <  2:976867/148‑1 (MQ=255)
                                                                                                                                        aTGGAACAATATGTATCCGGTTACCGACTGGGAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATAACCCGTATGTGCAACGCGCTTGCCAGGCGCGAAAGAGCTAACCTACACTAGCGGGATTCTTTTTGTTAAcc        <  2:2625121/148‑1 (MQ=255)
                                                                                                                                        aTGGAACAATATGTATCCGGTTACCGACTGGGAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATAACCCGTATGTGCAACGCGCTTGCCAGGCGCGAAAGAGCTAACCTACACTAGCGGGATTCTTTTTGTTAAcc        <  2:1119199/148‑1 (MQ=255)
                                                                                                                                          ggAACAATATGTATCCGGTTACCGACTGGGAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATAACCCGTATGTGCAACGCGCTTGCCAGGCGCGAAAGAGCTAACCTACACTAGCGGGATTCTTTTTGTTAAcccc      <  1:1210524/148‑1 (MQ=255)
                                                                                                                                          ggAACAAGATGTATCCGGTTACCGACTGGGAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATAACCCGTATGTGCAACGCGCTTGCCAGGCGCGAAAGAGCTAACCTACACTAGCGGGATTCTTTTTGTTAAcccc      <  1:1159325/148‑1 (MQ=255)
                                                                                                                                            aaCAATATGTATCCGGTTACCGACTGGGAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATAACCCGTATGTGCAACGCGCTTGCCAGGCGCGAAAGAGCTAACCTACACTAGCGGGATTCTTTTTGTTAACCCCTa    >  2:2891845/1‑148 (MQ=255)
                                                                                                                                             aCAATATGTATCCGGTTACCGACTGGGAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGCCAAACAAACCCCTATGTGCAACGCCCATTCCCGGGCCGGAAGTAGCAACCTACCCGAGGGGGCTTAGTTTTGGTTAAcccgaac   >  2:3236400/1‑144 (MQ=255)
                                                                                                                                              cAATATGTATCCGGTTACCGACTGGGAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATAACCCGTATGTGCAACGCGCTTGCCAGGCGCGAAAGAGCTAACCTACACTAGCGGGATTCTTTTTGTTAACCCCTAcc  >  1:2612020/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
CCATGAAGGTCGATTTCAAAGAAAAAGCTGCCGAACCACCAGCGCGTGCACGCGGTGCCATTGCGCGCACCTACTTCTATATGCGCGACCAATACAACCTGACACTCTCTCGCCAGCAAACGCAGCTGTTCAACGCATGGAACAAGATGTATCCGGTTACCGACTGGGAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATAACCCGTATGTGCAACGCGCTTGCCAGGCGCGAAAGAGCTAACCTACACTAGCGGGATTCTTTTTGTTAACCCCTACC  >  CP009273/3084160‑3084449

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: