Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,056,018 T→A 100% F138I (TTC→ATC)  yihN → MFS transporter family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,056,0180TA99.6% 731.5 / ‑4.6 220F138I (TTC→ATC) yihNMFS transporter family protein
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base A (97/122);  total (97/123)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 226 total aligned reads displayed.

TTTGACGACAATGCCGGGGTATTGGGGCATCCTCTTTGTCTGGGCGCTATTTGGCGTTACTTGCGACATGATGAACTGGCCGGTCTTGCTCAAGTCGGTAAGTCGATTGGGCAATAGCGAACAACAAGGTCGGTTGTTTGGCTTCTTCGAAACAGGGCGTGGCATTGTCGATACCGTGGTGGCATTTTCTGCGTTGGCAGTATTTACCTGGTTTGGCAGTGGCTTATTAGGTTTTAAAGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGCCGTAGGCATTATTATTT  >  CP009273/4055873‑4056162
                                                                                                                                                 |                                                                                                                                                
tttGACGACAATGCCGGGGTATTGGGGCATCCTCTTTGTCTGGGCGCTATTTGGCGTTACTTGCGACATGATGAACTGGCCGGTCTTGCTCAAGTCGGTAAGTCGATTGGGCAATAGCGAACAACAAGGTCGGTTGTTTGGCttcatc                                                                                                                                                >  2:1293260/1‑148 (MQ=255)
     cgacAATGCCGGGGTATTGGGGCATCCTCTTTGTCTGGGCGCTATTTGGCGTTACTTGCGACATGATGAACTGGCCGGTCTTGCTCAAGTCGGTAAGTCGATTGGGCAATAGCGAACAACAAGGTCGGTTGTTTGGCTTCATCGAAAc                                                                                                                                           <  1:3110643/148‑1 (MQ=255)
     cgacAATGCCGGGGTATTGGGGCATCCTCTTTGTCTGGGCGCTATTTGGCGTTACTTGCGACATGATGAACTGGCCGGTCTTGCTCAAGTCGGTAAGTCGATTGGGCAATAGCGAACAACAAGGTCGGTTGTTTGGCTTCATCGAAAc                                                                                                                                           <  1:499491/148‑1 (MQ=255)
      gacAATGCCGGGGTATTGGGGCATCCTCTTTGTCTGGGCGCTATTTGGCGTTACTTGCGACATGATGAACTGGCCGGTCTTGCTCAAGTCGGTAAGTCGATTGGGCAATAGCGAACAACAAGGTCGGTTGTTTGGCTTCATCGAAACa                                                                                                                                          <  2:2873089/148‑1 (MQ=255)
      gacAATGCCGGGGTATTGGGGCATCCTCTTTGTCTGGGCGCTATTTGGCGTTACTTGCGACATGATGAACTGGCCGGTCTTGCTCAAGTCGGTAAGTCGATTGGGCAATAGCGAACAACAAGGTCGGTTGTTTGGCTTCATCGAAACa                                                                                                                                          <  2:2750059/148‑1 (MQ=255)
        cAATGCCGGGGTATTGGGGCATCCTCTTTGTCTGGGCGCTATTTGGCGTTACTTGCGACATGATGAACTGGCCGGTCTTGCTCAAGTCGGTAAGTCGATTGGGCAATAGCGAACAACAAGGTCGGTTGTTTGGCTTCATCGAAACAgg                                                                                                                                        >  2:2436767/1‑148 (MQ=255)
        cAATGCCGGGGTATTGGGGCATCCTCTTTGTCTGGGCGCTATTTGGCGTTACTTGCGACATGATGAACTGGCCGGTCTTGCTCAAGTCGGTAAGTCGATTGGGCAATAGCGAACAACAAGGTCGGTTGTTTGGCTTCATCGAAACAgg                                                                                                                                        >  2:1901971/1‑148 (MQ=255)
         aaTGCCGGGGTATTGGGGCATCCTCTTTGTCTGGGCGCTATTTGGCGTTACTTGCGACATGATGAACGGGCCGGTCTTGCTCAAGTCGGTAAGTCGATTGGGCAATAGCGAACAACAAGGTCGGTTGTTTGGCTTCATCGAAACAggg                                                                                                                                       <  2:2138169/148‑1 (MQ=255)
          aTGCCGGGGTATTGGGGCATCCTCTTTGTCTGGGCGCTATTTGGCGTTACTTGCGACATGATGAACTGGCCGGTCTTGCTCAAGTCGGTAAGTCGATTGGGCAATAGCGAACAACAAGGTCGGTTGTTTGGCTTCATCGAAACAGGGc                                                                                                                                      >  1:1827304/1‑148 (MQ=255)
                 ggTATTGGGGCATCCTCTTTGTCTGGGCGCTATTTGGCGTTACTTGCGACATGATGAACTGGCCGGTCTTGCTCAAGTCGGTAAGTCGATTGGGCAATAGCGAACAACAAGGTCGGTTGTTTGGCTTCATCGAAACAGGGCGTGGCAt                                                                                                                               <  1:2317213/148‑1 (MQ=255)
                 ggTATTGGGGCATCCTCTTTGTCTGGGCGCTATTTGGCGTTACTTGCGACATGATGAACTGGCCGGTCTTGCTCAAGTCGGTAAGTCGATTGGGCAATAGCGAACAACAAGGTCGGTTGTTTGGCTTCATCGAAACAGGGCGTGGCAt                                                                                                                               <  2:2026273/148‑1 (MQ=255)
                   tATTGGGGCATCCTCTTTGTCTGGGCGCTATTTGGCGTTACTTGCGACATGATGAACTGGCCGGTCTTGCTCAAGTCGGTAAGTCGATTGGGCAATAGCGAACAACAAGGTCGGTTGTTTGGCTTCATCGAAACAGGGCGTGGCATTg                                                                                                                             <  1:717451/148‑1 (MQ=255)
                    aTTGGGGCATCCTCTTTGTCTGGGCGCTATTTGGCGTTACTTGCGACATGATGAACTGGCCGGTCTTGCTCAAGTCGGTAAGTCGATTGGGCAATAGCGAACAACAAGGTCGGTTGTTTGGCTTCATCGAAACAGGGCGTGGCATTGt                                                                                                                            <  1:3102101/148‑1 (MQ=255)
                    aTTGGGGCATCCTCTTTGTCTGGGCGCTATTTGGCGTTACTTGCGACATGATGAACTGGCCGGTCTTGCTCAAGTCGGTAAGTCGATTGGGCAATAGCGAACAACAAGGTCGGTTGTTTGGCTTCATCGAAACAGGGCGTGGCATTGt                                                                                                                            <  1:841921/148‑1 (MQ=255)
                    aTTGGGGCATCCTCTTTGTCTGGGCGCTATTTGGCGTTACTTGCGACATGATGAACTGGCCGGTCTTGCTCAAGTCGGTAAGTCGATTGGGCAATAGCGAACAACAAGGTCGGTTGTTTGGCTTCATCGAAACAGGGCGTGGCATTGt                                                                                                                            <  1:841943/148‑1 (MQ=255)
                      tGGGGCATCCTCTTTGTCTGGGCGCTATTTGGCGTTACTTGCGACATGATGAACTGGCCGGTCTTGCTCAAGTCGGTAAGTCGATTGGGCAATAGCGAACAACAAGGTCGGTTGTTTGGCTTCATCGAAACAGGGCGTGGCATTGTCg                                                                                                                          <  1:484052/148‑1 (MQ=255)
                        gggCATCCTCTTTGTCTGGGCGCTATTTGGCGTTACTTGCGACATGATGAACTGGCCGGTCTTGCTCAAGTCGGTAAGTCGATTGGGCAATAGCGAACAACAAGGTCGGTTGTTTGGCTTCATCGAAACAGGGCGTGGCATTGTCGAt                                                                                                                        <  1:2447604/148‑1 (MQ=255)
                               ctctTTGTCTGGGCGCTATTTGGCGTTACTTGCGACATGATGAACTGGCCGGTCTTGCTCAAGTCGGTAAGTCGATTGGGCAATAGCGAACAACAAGGTCGGTTGTTTGGCTTCATCGAAACAGGGCGTGGCATTGTCGATACCgtgg                                                                                                                 <  1:2661297/148‑1 (MQ=255)
                                tctTTGTCTGGGCGCTATTTGGCGTTACTTGCGACATGATGAACTGGCCGGTCTTGCTCAAGTCGGTAAGTCGATTGGGCAATAGCGAACAACAAGGTCGGTTGTTTGGCTTCATCGAAACAGGGCGTGGCATTGTAGATAACgtgtt                                                                                                                >  2:1030525/1‑146 (MQ=255)
                                  tttGTCTGGGCGCTATTTGGCGTTACTTGCGACATGATGAACTGGCCGGTCTTGCTCAAGTCGGTAAGTCGATTGGGCAATAGCGAACAACAAGGTCGGTTGTTTGGCTTCATCGAAACAGGGCGTGGCATTGTCGATACCGtggtgg                                                                                                              >  2:3116261/1‑148 (MQ=255)
                                   ttGTCTGGGCGCTATTTGGCGTTACTTGCGACATGATGAACTGGCCGGTCTTGCTCAAGTCGGTAAGTCGATTGGGCAATAGCGAACAACAAGGTCGGTTGTTTGGCTTCATCGAAACAGGGCGTGGCATTGTCGATACCGTGGTGGc                                                                                                             <  1:1373736/148‑1 (MQ=255)
                                   ttGTCTGGGCGCTATTTGGCGTTACTTGCGACATGATGAACTGGCCGGTCTTGCTCAAGTCGGTAAGTCGATTGGGCAATAGCGAACAACAAGGTCGGTTGTTTGGCTTCATCGAAACAGGGCGTGGCATTGTCGATACCGTGGTGGc                                                                                                             <  2:2283000/148‑1 (MQ=255)
                                          ggCGCTATTTGGCGTTACTTGCGACATGATGAACTGGCCGGTCTTGCTCAAGTCGGTAAGTCGATTGGGCAATAGCGAACAACAAGGTCGGTTGTTTGGCTTCATCGAAACAGGGCGTGGCATTGTCGATACCGTGGTGGCATTTTCt                                                                                                      >  2:572460/1‑148 (MQ=255)
                                           gcgcTATTTGGCGTTACTTGCGACATGATGAACTGGCCGGTCTTGCTCAAGTCGGTAAGTCGATAGGGCAATAGCGAACAACAAGGTCGGGTGTTTGGCTTCATCGAAACAGGGCGTGGCATTGTCGATACCGTGGTGGCATTTTCTg                                                                                                     >  1:484765/1‑148 (MQ=255)
                                           gcgcTATTTGGCGTTACTTGCGACATGATGAACGGGCCGGTCTTGCTCAAGTCGGTAAGTCGATTGGGCAATAGCGAACAACAAGGTCGGTTGTTTGGCTTCATCGAAACAGGGCGTGGCATTGTCGATACCGTGGTGGCATTTTCTg                                                                                                     <  2:2734568/148‑1 (MQ=255)
                                              cTATTTGGCGTTACTTGCGACATGATGAACTGGCCGGTCTTGCTCAAGTCGGTAAGTCGATTGGGCAATAGCGAACAACAAGGTCGGTTGTTTGGCTTCATCGAAACAGGGCGTGGCATTGTCGATACCGTGGTGGCATTTTCTGCGt                                                                                                  <  2:1175921/148‑1 (MQ=255)
                                               tATTTGGCGTTACTTGCGACATGATGAACTGGCCGGTCTTGCTCAAGTCGGTAAGTCGATTGGGCAATAGCGAACAACAAGGTCGGTTGTTTGGCTTCATCGAAACAGGGCGTGGCATTGTCGATACCGTGGTGGCATTTTCTGCGtt                                                                                                 <  1:1309896/148‑1 (MQ=255)
                                                aTTTGGCGTTACTTGCGACATGATGAACTGGCCGGTCTTGCTCAAGTCGGTAAGTCGATTGGGCAATAGCGAACAACAAGGTCGGTTGTTTGGCTTCATCGAAACAGGGCGTGGCATTGTCGATACCGTGGTGGCATTTTCTGCGTTg                                                                                                >  2:756330/1‑148 (MQ=255)
                                                aTTTGGCGTTACTTGCGACATGATGAACTGGCCGGTCTTGCTCAAGTCGGTAAGTCGATTGGGCAATAGCGAACAACAAGGTCGGTTGTTTGGCTTCATCGAAACAGGGCGTGGCATTGTCGATACCGTGGTGGCATTTTCTGAGTTg                                                                                                >  2:756563/1‑148 (MQ=255)
                                                    ggCGTTACTTGCGACATGATGAACTGGCCGGTCTTGCTCAAGTCGGTAAGTCGATTGGGCAATAGCGAACAACAAGGTCGGTTGTTTGGCTTCATCGAAACAGGGCGTGGCATTGTCGATACCGTGGTGGCATTTTCTGCGTTGGCAg                                                                                            <  1:93072/148‑1 (MQ=255)
                                                         tACTTGCGACATGATGAACTGGCCGGTCTTGCTCAAGTCGGTAAGTCGATTGGGCAATAGCGAACAACAAGGTCGGTTGTTTGGCTTCATCGAAACAGGGCGTGGCATTGTCGATACCGTGGTGGCATTTTCTGCGTTGGCAGTAttt                                                                                       >  2:1145651/1‑148 (MQ=255)
                                                         tACTTGCGACATGATGAACTGGCCGGTCTTGCTCAAGTCGGTAAGTCGATTGGGCAATAGCGAACAACAAGGTCGGTTGTTTGGCTTCATCGAAACAGGGCGTGGCATTGTCGATACCGTGGTGGCATTTTCTGCGTTGGCAGTAttt                                                                                       >  2:958101/1‑148 (MQ=255)
                                                         tACTTGCGACATGATGAACTGGCCGGTCTTGCTCAAGTCGGTAAGTCGATTGGGCAATAGCGAACAACAAGGTCGGTTGTTTGGCTTCATCGAAACAGGGCGTGGCATTGTCGATACCGTGGTGGCATTTTCTGCGTTGGCAGTAttt                                                                                       >  1:20606/1‑148 (MQ=255)
                                                          aCTTGCGACATGATGAACTGGCCGGTCTTGCTCAAGTCGGTAAGTCGATTGGGCAATAGCGAACAACAAGGTCGGTTGTTTGGCTTCATCGAAACAGGGCGTGGCATTGTCGATACCGTGGTGGCATTTTCTGCGTTGGCAGTATTTa                                                                                      >  1:1021513/1‑148 (MQ=255)
                                                          aCTTGCGACATGATGAACTGGCCGGTCTTGCTCAAGTCGGTAAGTCGATTGGGCAATAGCGAACAACAAGGTCGGTTGTTTGGCTTCATCGAAACAGGGCGTGGCATTGTCGATACCGTGGTGGCATTTTCTGCGTTGGCAGTATTTa                                                                                      >  1:2822419/1‑148 (MQ=255)
                                                           cTTGCGACATGATGAACTGGCCGGTCTTGCTCAAGTCGGTAAGTCGATTGGGCAATAGCGAACAACAAGGTCGGTTGTTTGGCTTCATCGAAACAGGGCGTGGCATTGTCGATACCGTGGTGGCATTTTCTGCGTTGGCAGTATTTAc                                                                                     <  1:2400573/148‑1 (MQ=255)
                                                           cTTGCGACATGATGAACTGGCCGGTCTTGCTCAAGTCGGTAAGTCGATTGGGCAATAGCGAACAACAAGGTCGGTTGTTTGGCTTCATCGAAACAGGGCGTGGCATTGTCGATACCGTGGTGGCATTTTCTGCGTTGGCAGTATTTAc                                                                                     <  2:1227763/148‑1 (MQ=255)
                                                            ttGCGACATGATGAACTGGCCGGTCTTGCTCAAGTCGGTAAGTCGATTGGGCAATAGCGAACAACAAGGTCGGTTGTTTGGCTTCATCGAAACAGGGCGTGGCATTGTCGATACCGTGGTGGCATTTTCTGCGTTGGCAGTATTTAcc                                                                                    >  2:1907530/1‑148 (MQ=255)
                                                                 aCATGCTGAACTGGCCGGTCTTGCTCAAGTCGGTAAGTCGATTGGGCAATAGCGAACAACAAGGTCGGTTGTTTGGCTTCATCGAAACAGGGCGTGGCATTGTCGATACCGTGGTGGCATTTTCTGCGTTGGCAGTATTTACCTGGtt                                                                               <  1:1895799/148‑1 (MQ=255)
                                                                 aCATGATGAACTGGCCGGTCTTGCTCAAGTCGGTAAGTCGATTGGGCAATAGCGAACAACAAGGTCGGTTGTTTGGCTTCATCGAAACAGGGCGTGGCATTGTCGATACCGTGGTGGCATTTTCTGCGTTGGCAGTATTTACCTGGtt                                                                               <  2:2725712/148‑1 (MQ=255)
                                                                 aCATGATGAACTGGCCGGTCTTGCTCAAGTCGGTAAGTCGATTGGGCAATAGCGAACAACAAGGTCGGTTGTTTGGCTTCATCGAAACAGGGCGTGGCATTGTCGATACCGTGGTGGCATTTTCTGCGTTGGCAGTATTTACCTGGtt                                                                               >  2:484912/1‑148 (MQ=255)
                                                                  cATGATGAACTGGCCGGTCTTGCTCAAGTCGGTAAGTCGATTGGGCAATAGCGAACAACAAGGTCGGTTGTTTGGCTTCATCGAAACAGGGCGTGGCATTGTCGATACCGTGGTGGCATTTTCTGCGTTGGCAGTATTTACCTGGttt                                                                              <  1:1297078/148‑1 (MQ=255)
                                                                  cATGATGAACTGGCCGGTCTTGCTCAAGTCGGTAAGTCGATTGGGCAATAGCGAACAACAAGGTCGGTTGTTTGGCTTCATCGAAACAGGGCGTGGCATTGTCGATACCGTGGTGGCATTTTCTGCGTTGGCAGTATTTACCTGGttt                                                                              <  1:1998503/148‑1 (MQ=255)
                                                                  cATGATGAACTGGCCGGTCTTGCTCAAGTCGGTAAGTCGATTGGGCAATAGCGAACAACAAGGTCGGTTGTTTGGCTTCATCGAAACAGGGCGTGGCATTGTCGATACCGTGGTGGCATTTTCTGCGTTGGCAGTATTTACCTGGttt                                                                              >  1:345774/1‑148 (MQ=255)
                                                                      atgaACTGGCCGGTCTTGCTCAAGTCGGTAAGTCGATTGGGCAATAGCGAACAACAAGGTCGGTTGTTTGGCTTCATCGAAACAGGGCGTGGCATTGTCGATACCGTGGTGGCATTTTCTGCGTTGGCAGTATTTACCTGGTTTGGCa                                                                          >  1:438398/1‑148 (MQ=255)
                                                                       tgaACTGGCCGGTCTTGCTCAAGTCGGTAAGTCGATTGGGCAATAGCGAACAACAAGGTCGGTTGTTTGGCTTCATCGAAACAGGGCGTGGCATTGTCGATACCGTGGTGGCATTTTCTGCGTTGGCAGTATTTACCTGGTTTGGCAg                                                                         <  2:3064628/148‑1 (MQ=255)
                                                                        gaACTGGCCGGTCTTGCTCAAGTCGGTAAGTCGATTGGGCAATAGCGAACAACAAGGTCGGTTGTTTGGCTTCATCGAAACAGGGCGTGGCATTGTCGATACCGTGGTGGCATTTTCTGCGTTGGCAGTATTTACCTGGTTTGGCAGt                                                                        >  1:1791534/1‑148 (MQ=255)
                                                                           cTGGCCGGTCTTGCTCAAGTCGGTAAGTCGATTGGGCAATAGCGAACAACAAGGTCGGTTGTTTGGCTTCATCGAAACAGGGCGTGGCATTGTCGATACCGTGGTGGCATTTTCTGCGTTGGCAGTATTTACCTGGTTTGGCAGTGGc                                                                     >  1:1329062/1‑148 (MQ=255)
                                                                           cTGGCCGGTCTTGCTCAAGTCGGTAAGTCGATTGGGCAATAGCGAACAACAAGGTCGGTTGTTTGGCTTCATCGAAACAGGGCGTGGCATTGTCGATACCGTGGTGGCATTTTCTGCGTTGGCAGTATTTACCTGGTTTGGCAGTGGc                                                                     <  1:1555752/148‑1 (MQ=255)
                                                                           cTGGCCGGTCTTGCTCAAGTCGGTAAGTCGATTGGGCAATAGCGAACAACAAGGTCGGTTGTTTGGCTTCATCGAAACAGGGCGTGGCATTGTCGATACCGTGGTGGCATTTTCTGCGTTGGCAGTATTTACCTGGTTTGGCAGTGGc                                                                     <  1:1307332/148‑1 (MQ=255)
                                                                              gCCGGTCTTGCTCAAGTCGGTAAGTCGATTGGGCAATAGCGAACAACAAGGTCGGTTGTTTGGCTTCATCGAAACAGGGCGTGGCATTGTCGATACCGTGGTGGCATTTTCTGCGTTGGCAGTATTTACCTGGTTTGGCAGTGGCtta                                                                  >  1:463284/1‑148 (MQ=255)
                                                                                 ggTCTTGCTCAAGTCGGTAAGTCGATTGGGCAATAGCGAACAACAAGGTCGGTTGTTTGGCTTCATCGAAACAGGGCGTGGCATTGTCGATACCGTGGTGGCATTTTCTGCGTTGGCAGTATTTACCTGGTTTGGCAGTGGCttatta                                                               <  1:2886713/148‑1 (MQ=255)
                                                                                 ggTCTTGCTCAAGTCGGTAAGTCGATTGGGCAATAGCGAACAACAAGGTCGGTTGTTTGGCTTCATCGAAACAGGGCGTGGCATTGTCGATACCGTGGTGGCATTTTCTGCGTTGGCAGTATTTACCTGGTTTGGCAGTGGCttatta                                                               <  1:2340128/148‑1 (MQ=255)
                                                                                 ggTCTTGCTCAAGTCGGTAAGTCGATTGGGCAATAGCGAACAACAAGGTCGGTTGTTTGGCTTCATCGAAACAGGGCGTGGCATTGTCGATACCGTGGTGGCATTTTCTGCGTTGGCAGTATTTACCTGGTTTGGCAGTGGCttatta                                                               <  1:678564/148‑1 (MQ=255)
                                                                                   tCTTGCTCAAGTCGGTAAGTCGATTGGGCAATAGCGAACAACAAGGTCGGTTGTTTGGCTTCATCGAAACAGGGCGTGGCATTGTCGATACCGTGGTGGCATTTTCTGCGTTGGCAGTATTTACCTGGTTTGGCAGTGGCTTATTAgg                                                             <  1:2517255/148‑1 (MQ=255)
                                                                                   tCTTGCTCAAGTCGGTAAGTCGATTGGGCAATAGCGAACAACAAGGTCGGTTGTTTGGCTTCATCGAAACAGGGCGTGGCATTGTCGATACCGTGGTGGCATTTTCTGCGTTGGCAGTATTTACCTGGTTTGGCAGTGGCTTATTAgg                                                             >  1:621713/1‑148 (MQ=255)
                                                                                    cTTGCTCAAGTCGGTAAGTCGATTGGGCAATAGCGAACAACAAGGTCGGTTGTTTGGCTTCATCGAAACAGGGCGTGGCATTGTCGATACCGTGGTGGCATTTTCTGCGTTGGCAGTATTTACCTGGTTTGGCAGTGGCTTATTAGGt                                                            <  1:2145155/148‑1 (MQ=255)
                                                                                    cTTGCTCAAGTCGGTAAGTCGATTGGGCAATAGCGAACAACAAGGTCGGTTGTTTGGCTTCATCGAAACAGGGCGTGGCATTGTCGATACCGTGGTGGCATTTTCTGCGTTGGCAGTATTTACCTGGTTTGGCAGTGGCTTATTAGGt                                                            <  2:2365948/148‑1 (MQ=255)
                                                                                    cTTGCTCAAGTCGGTAAGTCGATTGGGCAATAGCGAACAACAAGGTCGGTTGTTTGGCTTCATCGAAACAGGGCGTGGCATTGTCGATACCGTGGTGGCATTTTCTGCGTTGGCAGTATTTACCTGGTTTGGCAGTGGCTTATTAGGt                                                            >  1:3081902/1‑148 (MQ=255)
                                                                                      tGCTCAAGTCGGTAAGTCGATTGGGCAATAGCGAACAACAAGGTCGGTTGTTTGGCTTCATCGAAACAGGGCGTGGCATTGTCGATACCGTGGTGGCATTTTCTGCGTTGGCAGTATTTACCTGGTTTGGCAGTGGCTTATTAGGttt                                                          <  1:845278/148‑1 (MQ=255)
                                                                                      tGCTCAAGTCGGTAAGTCGATTGGGCAATAGCGAACAACAAGGTCGGTTGTTTGGCTTCATCGAAACAGGGCGTGGCATTGTCGATACCGTGGTGGCATTTTCTGCGTTGGCAGTATTTACCTGGTTTGGCAGTGGCTTATTAGGttt                                                          >  1:2362542/1‑148 (MQ=255)
                                                                                      tGCTCAAGTCGGTAAGTCGATTGGGCAATAGCGAACAACAAGGTCGGTTGTTTGGCTTCATCGAAACAGGGCGTGGCATTGTCGATACCGTGGTGGCATTTTCTGCGTTGGCAGTATTTACCTGGTTTGGCAGTGGCTTATTAGGttt                                                          <  2:1036280/148‑1 (MQ=255)
                                                                                         tCAAGTCGGTAAGTCGATTGGGCAATAGCGAACAACAAGGTCGGTTGTTTGGCTTCATCGAAACAGGGCGTGGCATTGTCGATACCGTGGTGGCATTTTCTGCGTTGGCAGTATTTACCTGGTTTGGCAGTGGCTTATTAGGTTTTaa                                                       <  2:1763018/148‑1 (MQ=255)
                                                                                          cAAGTCGGTAAGTCGATTGGGCAATAGCGAACAACAAGGTCGGTTGTTTGGCTTCATCGAAACAGGGCGTGGCATTGTCGATACCGTGGTGGCATTTTCTGCGTTGGCAGTATTTACCTGGTTTGGCAGTGGCTTATTAGGTTTTaaa                                                      <  2:51784/148‑1 (MQ=255)
                                                                                           aaGTCGGTAAGTCGATTGGGCAATAGCGAACAACAAGGTCGGTTGTTTGGCTTCATCGAAACAGGGCGTGGCATTGTCGATACCGTGGTGGCATTTTCTGCGTTGGCAGTATTTACCTGGTTTGGCAGTGGCTTATTAGGTTTTAAAg                                                     >  1:1609667/1‑148 (MQ=255)
                                                                                            aGTCGGTAAGTCGATTGGGCAATAGCGAACAACAAGGTCGGTTGTTTGGCTTCATCGAAACAGGGCGTGGCATTGTCGATACCGTGGTGGCATTTTCTGCGTTGGCAGTATTTACCTGGTTTGGCAGTGGCTTATTAGGTTTTAAAGc                                                    <  2:1467383/148‑1 (MQ=255)
                                                                                             gTCGGTAAGTCGATTGGGCAATAGCGAACAACAAGGTCGGTTGTTTGGCTTCATCGAAACAGGGCGTGGCATTGTCGATACCGTGGTGGCATTTTCTGCGTTGGCAGTATTTACCTGGTTTGGCAGTGGCTTATTAGGTTTTAAAGCa                                                   >  2:1368648/1‑148 (MQ=255)
                                                                                                ggTAAGTCGATTGGGCAATAGCGAACAACAAGGTCGGTTGTTTGGCTTCATCGAAACAGGGCGTGGCATTGTCGATACCGTGGTGGCATTTTCTGCGTTGGCAGTATTTACCTGGTTTGGCAGTGGCTTATTAGGTTTTAAAGCAGGc                                                >  2:221051/1‑148 (MQ=255)
                                                                                                 gTAAGTCGATTGGGCAATAGCGAACAACAAGGTCGGTTGTTTGGCTTCATCGAAACAGGGCGTGGCATTGTCGATACCGTGGTGGCATTTTCTGCGTTGGCAGTATTTACCTGGTTTGGCAGTGGCTTATTAGGTTTTAAAGCAGGCa                                               <  2:2025169/148‑1 (MQ=255)
                                                                                                    aGTCGATTGGGCAATAGCGAACAACAAGGTCGGTTGTTTGGCTTCATCGAAACAGGGCGTGGCATTGTCGATACCGTGGTGGCATTTTCTGCGTTGGCAGTATTTACCTGGTTTGGCAGTGGCTTATTAGGTTTTAAAGCAGGCATCt                                            >  1:207626/1‑148 (MQ=255)
                                                                                                     gTCGATTGGGCAATAGCGAACAACAAGGTCGGTTGTTTGGCTTCATCGAAACAGGGCGTGGCATTGTCGATACCGTGGTGGCATTTTCTGCGTTGGCAGTATTTACCTGGTTTGGCAGTGGCTTATTAGGTTTTAAAGCAGGCATCTg                                           <  2:2420679/148‑1 (MQ=255)
                                                                                                      tCGATTGGGCAATAGCGAACAACAAGGTCGGTTGTTTGGCTTCATCGAAACAGGGCGTGGCATTGTCGATACCGTGGTGGCATTTTCTGCGTTGGCAGTATTTACCTGGTTTGGCAGTGGCTTATTAGGTTTTAAAGCAGGCATCTgg                                          >  2:2256958/1‑148 (MQ=255)
                                                                                                         aTTGGGCAATAGCGAACAACAAGGTCGGTTGTTTGGCTTCATCGAAACAGGGCGTGGCATTGTCGATACCGTGGTGGCATTTTCTGCGTTGGCAGTATTTACCTGGTTTGGCAGTGGCTTATTAGGTTTTAAAGCAGGCATCTGGTTc                                       >  2:3147040/1‑148 (MQ=255)
                                                                                                         aTTGGGCAATAGCGAACAACAAGGTCGGTTGTTTGGCTTCATCGAAACAGGGCGTGGCATTGTCGATACCGTGGTGGCATTTTCTGCGTTGGCAGTATTTACCTGGTTTGGCAGTGGCTTATTAGGTTTTAAAGCAGGCATCTGGTTc                                       <  2:1335177/148‑1 (MQ=255)
                                                                                                           tGGGCAATAGCGAACAACAAGGTCGGTTGTTTGGCTTCATCGAAACAGGGCGTGGCAGTGTCGATACCGTGGTGGCATTTTCTGCGTTGGCAGTATTTACCTGGTTTGGCAGTGGCTTATTAGGTTTTAAAGCAGGCATCTGGTTCTa                                     <  1:3193278/148‑1 (MQ=255)
                                                                                                              gCAATAGCGAACAACAAGGTCGGTTGTTTGGCTTCATCGAAACAGGGCGTGGCATTGTCGATACCGTGGTGGCATTTTCTGCGTTGGCAGTATTTACCTGGTTTGGCAGTGGCTTATTAGGTTTTAAAGCAGGCATCTGGTTCTATTc                                  >  2:3040042/1‑148 (MQ=255)
                                                                                                                 aTAGCGAACAACAAGGTCGGTTGTTTGGCTTCATCGAAACAGGGCGTGGCATTGTCGATACCGTGGTGGCATTTTCTGCGTTGGCAGTATTTACCTGGTTTGGCAGTGGCTTATTAGGTTTTAAAGCAGGCATCTGGTTCTATTCCCt                               >  1:1086843/1‑148 (MQ=255)
                                                                                                                     cGAACAACAAGGTCGGTTGTTTGGCTTCATCGAAACAGGGCGTGGCATTGTCGATACCGTGGTGGCATTTTCTGCGTTGGCAGTATTTACCTGGTTTGGCAGTGGCTTATTAGGTTTTAAAGCAGGCATCTGGTTCTATTCCCTTAtt                           >  1:2043363/1‑148 (MQ=255)
                                                                                                                        acaacaAGGTCGGTTGTTTGGCTTCATCGAAACAGGGCGTGGCATTGTCGATACCGTGGTGGCATTTTCTGCGTTGGCAGTATTTACCTGGTTTGGCAGTGGCTTATTAGGTTTTAAAGCAGGCATCTGGTTCTATTCCCTTATtgtg                        <  2:1996063/148‑1 (MQ=255)
                                                                                                                        acaacaAGGTCGGTTGTTTGGCTTCATCGAAACAGGGCGTGGCATTGTCGATACCGTGGTGGCATTTTCTGCGTTGGCAGTATTTACCTGGTTTGGCAGTGGCTTATTAGGTTTTAAAGCAGGCATCTGGTTCTATTCCCTTATtgtg                        >  2:2687414/1‑148 (MQ=255)
                                                                                                                         caacaaGGTCGGTTGTTTGGCTTCATCGAAACAGGGCGTGGCATTGTCGATACCGTGGTGGCATTTTCTGCGTTGGCAGTATTTACCTGGTTTGGCAGTGGCTTATTAGGTTTTAAAGCAGGCATCTGGTTCTATTCCCTTATTGTGa                       <  1:13448/148‑1 (MQ=255)
                                                                                                                           acaaGGTCGGTTGTTTGGCTTCATCGAAACAGGGCGTGGCATTGTCGATACCGTGGTGGCATTTTCTGCGTTGGCAGTATTTACCTGGTTTGGCAGTGGCTTATTAGGTTTTAAAGCAGGCATCTGGTTCTATTCCCTTATTGTGAtt                     >  2:2564824/1‑148 (MQ=255)
                                                                                                                            caaGGTCGGTTGTTTGGCTTCATCGAAACAGGGCGTGGCATTGTCGATACCGTGGTGGCATTTTCTGCGTTGGCAGTATTTACCTGGTTTGGCAGTGGCTTATTAGGTTTTAAAGCAGGCATCTGGTTCTATTCCCTTATTGTGATTg                    >  2:46828/1‑148 (MQ=255)
                                                                                                                            caaGGTCGGTTGTTTGGCTTCATCGAAACAGGGCGTGGCATTGTCGATACCGTGGTGGCATTTTCTGCGTTGGCAGTATTTACCTGGTTTGGCAGTGGCTTATTAGGTTTTAAAGCAGGCATCTGGTTCTATTCCCTTATTGTGATTg                    <  1:468745/148‑1 (MQ=255)
                                                                                                                            caaGGTCGGTTGTTTGGCTTCATCGAAACAGGGCGTGGCAGTGTCGATACCGTGGTGGCAATTTCTGCGTTGGCAGTATTTAACTGGTTTTGCAGTGGCTTATTAGGTTTTAAAACAGGCATCTGGTTGTATTCCCTTATTGTGATTg                    >  1:1470441/1‑148 (MQ=255)
                                                                                                                               ggTCGGTTGTTTGGCTTCATCGAAACAGGGCGTGGCATTGTCGATACCGTGGTGGCATTTTCTGCGTTGGCAGTATTTACCTGGTTTGGCAGTGGCTTATTAGGTTTTAAAGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGCCg                 >  2:2063684/1‑148 (MQ=255)
                                                                                                                               ggTCGGTTGTTTGGCTTCATCGAAACAGGGCGTGGCATTGTCGATACCGTGGTGGCATTTTCTGCGTTGGCAGTATTTACCTGGTTTGGCAGTGGCTTATTAGGTTTTAAAGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGCCg                 >  2:395905/1‑148 (MQ=255)
                                                                                                                                gTCGGTTGTTTGGCTTCATCGAAACAGGGCGTGGCATTGTCGATACCGTGGTGGCATTTTCTGCGTTGGCAGTATTTACCTGGTTTGGCAGTGGCTTATTAGGTTTTAAAGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGCCGt                >  2:1840973/1‑148 (MQ=255)
                                                                                                                                gTCGGTTGTTTGGCTTCATCGAAACAGGGCGTGGCATTGTCGATACCGTGGTGGCATTTTCTGCGTTGGCAGTATTTACCTGGTTTGGCAGTGGCTTATTAGGTTTTAAAGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGCCGt                <  1:2920347/148‑1 (MQ=255)
                                                                                                                                    gttgttTGGCTTCATCGAAACAGGGCGTGGCATTGTCGATACCGTGGTGGCATTTTCTGCGTTGGCAGTATTTACCTGGTTTGGCAGTGGCTTATTAGGTTTTAAAGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGCCGTAGGc            >  2:3231217/1‑148 (MQ=255)
                                                                                                                                     ttgttTGGCTTCATCGAAACAGGGCGTGGCATTGTCGATACCGTGGTGGCATTTTCTGCGTTGGCAGTATTTACCTGGTTTGGCAGTGGCTTATTAGGTTTTAAAGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGCCGTAGGCa           >  2:1631672/1‑148 (MQ=255)
                                                                                                                                     ttgttTGGCTTCATCGAAACAGGGCGTGGCATTGTCGATACCGTGGTGGCATTTTCTGCGTTGGCAGTATTTACCTGGTTTGGCAGTGGCTTATTAGGTTTTAAAGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGCCGTAGGCa           >  2:1543127/1‑148 (MQ=255)
                                                                                                                                      tgttTGGCTTCATCGAATCAGGGCGTGGCATTGTCGATTCCGTGGTGGCATTTTCTGCGTTGGCAGTATTTACCTGGTTTGGCAGTGGCTTATTAGGTTTTAAAGCAGGCATCTGGTTCTATTCCGTTATTGTGATTGCCGTAGGCat          <  2:1387197/148‑1 (MQ=255)
                                                                                                                                      tgttTGGCTTCATCGAAACAGGGCGTGGCATTGTCGATACCGTGGTGGCATTTTCTGCGTTGGCAGTATTTACCTGGTTTGGCAGTGGCTTATTAGGTTTTAAAGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGCCGTAGGCat          <  1:1228929/148‑1 (MQ=255)
                                                                                                                                      tgttTGGCTTCATCGAAACAGGGCGTGGCATTGTCGATACCGTGGTGGCATTTTCTGCGTTGGCAGTATTTACCTGGTTTGGCAGTGGCTTATTAGGTTTTAAAGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGCCGTAGGCat          <  1:395929/148‑1 (MQ=255)
                                                                                                                                        tttGGCTTCATCGAAACAGGGCGTGGCATTGTCGATACCGTGGTGGCATTTTCTGCGTTGGCAGTATTTACCTGGTTTGGCAGTGGCTTATTAGGTTTTAAAGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGCCGTAGGCatta        <  2:3026501/148‑1 (MQ=255)
                                                                                                                                         ttGGCTTCATCGAAACAGGGCGTGGCATTGTCGATACCGTGGTGGCATTTTCTGCGTTGGCAGTATTTACCTGGTTTGGCAGTGGCTTATTAGGTTTTAAAGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGCCGTAGGCattat       <  1:775797/148‑1 (MQ=255)
                                                                                                                                             cttcatCGAAACAGGGCGTGGCATTGTCGATACCGTGGTGGCATTTTCTGCGTTGGCAGTATTTACCTGGTTTGGCAGTGGCTTATTAGGTTTTAAAGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGCCGTAGGCattattatt   <  1:1212078/148‑1 (MQ=255)
                                                                                                                                              ttcatcGAAACAGGGCGTGGCATTGTCGATACCGTGGTGGCATTTTCTGCGTTGGCAGTATTTACCTGGTTTGGCAGTGGCTTATTAGGTTTTAAAGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGCCGTAGGCATTATTAttt  <  2:1078857/148‑1 (MQ=255)
                                                                                                                                              ttcatcGAAACAGGGCGTGGCATTGTCGATACCGTGGTGGCATTTTCTGCGTTGGCAGTATTTACCTGGTTTGGCAGTGGCTTATTAGGATTTAAAGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGCCGTAGGCATTATTAttt  <  2:1512098/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
TTTGACGACAATGCCGGGGTATTGGGGCATCCTCTTTGTCTGGGCGCTATTTGGCGTTACTTGCGACATGATGAACTGGCCGGTCTTGCTCAAGTCGGTAAGTCGATTGGGCAATAGCGAACAACAAGGTCGGTTGTTTGGCTTCTTCGAAACAGGGCGTGGCATTGTCGATACCGTGGTGGCATTTTCTGCGTTGGCAGTATTTACCTGGTTTGGCAGTGGCTTATTAGGTTTTAAAGCAGGCATCTGGTTCTATTCCCTTATTGTGATTGCCGTAGGCATTATTATTT  >  CP009273/4055873‑4056162

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: