Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,825,702 C→A 100% K732N (AAG→AAT yicI ← putative alpha‑glucosidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,825,7020CA99.7% 945.3 / ‑4.5 291K732N (AAG→AATyicIputative alpha‑glucosidase
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base A (143/147);  total (144/147)
Fisher's exact test for biased strand distribution p-value = 4.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.84e-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 299 total aligned reads displayed.

CGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTCACAACATTGCGCAGGCACAGTGTCCAGTTCTTCGCCTCGCCCGCACCAGTCACAGTAATCGTGTTGCCAGTACGTGCTGCTTTTAAAGTAAAGATCACCGATCCGTCAGCAGCGGGCACTTCACAGACGGCTTCATGCCCGTCTTGCAGATTGAAGAGGTGGAATGCCGTGCCTTC  >  CP009273/3825555‑3825848
                                                                                                                                                   |                                                                                                                                                  
cGTACCCGTCATTCCATCAATAGATTACACCGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTCACAACATTGCGCAGGCACAGTGTCCAGttatt                                                                                                                                                  <  2:4972921/150‑1 (MQ=255)
cGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCAGTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTCACAACATTGCGCAGGCACAGTGTCCAGttatt                                                                                                                                                  <  1:5214938/150‑1 (MQ=255)
 gTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTCACAACATTGCGCAGGCACAGTGTCCAGttattc                                                                                                                                                 <  2:4419082/150‑1 (MQ=255)
  tACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCTCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCTTCTTGCAGACCATTTACTTTCACAACATTGCGCAGGCACAGTGTCCAGTTATtcg                                                                                                                                                <  1:4390422/150‑1 (MQ=255)
  tACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTCACAACATTGCGCAGGCACAGTGTCCAGTTATtcg                                                                                                                                                >  2:4533796/1‑150 (MQ=255)
  tACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTCACAACATTGCGCAGGCACAGTGTCCAGTTATtcg                                                                                                                                                >  1:3187547/1‑150 (MQ=255)
   aCCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTCACAACATTGCGCAGGCACAGTGTCCAGATATtcgc                                                                                                                                               >  2:1338143/1‑150 (MQ=255)
     ccGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTCACAACATTGCGCAGGCACAGTGTCCAGTTATtcgcct                                                                                                                                             >  2:2676575/1‑150 (MQ=255)
       ggCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTCACAACATTGCGCAGGCACAGTGTCCAGTTATtcgcctcg                                                                                                                                           <  1:238540/150‑1 (MQ=255)
       ggCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTCACAACATTGCGCAGGCACAGTGTCCAGTTATtcgcctcg                                                                                                                                           >  1:5126347/1‑150 (MQ=255)
         cATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTCACAACATTGCGCAGGCACAGTGTCCAGTTATtcgcctcgcc                                                                                                                                         >  2:700891/1‑150 (MQ=255)
             ccATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTCACAACATTGCGCAGGCACAGTGTCCAGTTATTCGCCTCGCCCGCa                                                                                                                                     >  2:1918743/1‑150 (MQ=255)
             ccATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTCACAACATTGCGCAGGCACAGTGTCCAGTTATTCGCCTCGCCCGCa                                                                                                                                     >  2:3134863/1‑150 (MQ=255)
              cATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTCACAACATTGCGCAGGCACAGTGTCCAGTTATTCGCCTCGCCCGCAc                                                                                                                                    >  2:3829808/1‑150 (MQ=255)
              cATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTCACAACATTGCGCAGGCACAGTGTCCAGTTATTCGCCTCGCCCGCAc                                                                                                                                    >  2:1405297/1‑150 (MQ=255)
                tCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTCACAACATTGCGCAGGCACAGTGTCCAGTTATTCGCCTCGCCCGCACCa                                                                                                                                  >  1:329741/1‑150 (MQ=255)
                tCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTCACAACATTGCGCAGGCACAGTGTCCAGTTATTCGCCTCGCCCGCACCa                                                                                                                                  >  2:2614931/1‑150 (MQ=255)
                 cAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTCACAACATTGCGCAGGCACAGTGTCCAGTTATTCGCCTCGCCCGCACCAg                                                                                                                                 >  1:1593097/1‑150 (MQ=255)
                  aaTAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTCACAACATTGCGCAGGCACAGTGTCCAGTTATTCGCCTCGCCCGCACCAGt                                                                                                                                <  2:5313002/150‑1 (MQ=255)
                  aaTAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTCACAACATTGCGCAGGCACAGTGTCCAGTTATTCGCCTCGCCCGCACCAGt                                                                                                                                >  1:1638782/1‑150 (MQ=255)
                      gATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTACTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTCACAACATTGCGCAGGCACAGTGTCCAGTTATTCGCCTCGCCCGCACCAGTcaca                                                                                                                            <  2:166523/150‑1 (MQ=255)
                       aTTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTCACAACATTGCTCAGGCACAGTGTCCAGTTATTCGCCTCGCCCGCACCAGTCACAg                                                                                                                           >  2:569128/1‑150 (MQ=255)
                       aTTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTCACAACATTGCGCAGGCACAGTGTCCAGTTATTCGCCTCGCCCGCACCAGTCACAg                                                                                                                           >  2:5215126/1‑150 (MQ=255)
                       aTTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTCACAACATTGCGCAGGCACAGTGTCCAGTTATTCGCCTCGCCCGCACCAGTCACAg                                                                                                                           >  1:3218111/1‑150 (MQ=255)
                        ttACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTCACAACATTGCGCAGGCACAGTGTCCAGTTATTCGCCTCGCCCGCACCAGTCACAGt                                                                                                                          <  2:1528229/150‑1 (MQ=255)
                         tACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTCACAACATTGCGCAGGCACAGTGTCCAGTTATTCGCCTCGCCCGCACCAGTCACAGTa                                                                                                                         >  2:2283467/1‑150 (MQ=255)
                         tACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTCACAACATTGCGCAGGCACAGTGTCCAGTTATTCGCCTCGCCCGCACCAGTCACAGTa                                                                                                                         >  1:2218482/1‑150 (MQ=255)
                          aCAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTCACAACATTGCGCAGGCACAGTGTCCAGTTATTCGCCTCGCCCGCACCAGTCACAGTaa                                                                                                                        <  2:3739987/150‑1 (MQ=255)
                           cAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTCACAACATTGCGCAGGCACAGTGTCCAGTTATTCGCCTCGCCCGCACCAGTCACAGTAAt                                                                                                                       <  1:352640/150‑1 (MQ=255)
                            aaCGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTCACAACATTGCGCAGGCACAGTGTCCAGTTATTCGCCTCGCCCGCACCAGTCACAGTAATc                                                                                                                      <  2:32894/150‑1 (MQ=255)
                                 aaTTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTCACAACATTGCGCAGGCACAGTGTCCAGTTATTCGCCTCGCCCGCACCAGTCACAGTAATCGTGtt                                                                                                                 <  1:2490965/150‑1 (MQ=255)
                                  aTTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTCACAACATTGCGCAGGCACAGTGTCCAGTTATTCGCCTCGCCCGCACCAGTCACAGTAATCGTGTTg                                                                                                                >  2:4448519/1‑150 (MQ=255)
                                     gTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTCACAACATTGCGCAGGCACAGTGTCCAGTTATTCGCCTCGCCCGCACCAGTCACAGTAATCGTGTTGCCa                                                                                                             >  2:2259231/1‑150 (MQ=255)
                                        aGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTCACAACATTGCGCAGGCACAGTGTCCAGTTATTCGCCTCGCCCGCACCAGTCACAGTAATCGTGTTGCCAGTa                                                                                                          >  1:481087/1‑150 (MQ=255)
                                           gcATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTCACAACATTGCGCAGGCACAGTGTCCAGTTATTCGCCTCGCCCGCACCAGTCACAGTAATCGTGTTGCCAGTACGt                                                                                                       >  1:4137386/1‑150 (MQ=255)
                                            cATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTCACAACATTGCGCAGGCACAGTGTCCAGTTATTCGCCTCGCCCGCACCAGTCACAGTAATCGTGTTGCCAGTACGtg                                                                                                      >  2:1822612/1‑150 (MQ=255)
                                             aTTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTCACAACATTGCGCAGGCACAGTGTCCAGTTATTCGCCTCGCCCGCACCAGTCACAGTAATCGTGTTGCCAGTACGtgc                                                                                                     >  1:533837/1‑150 (MQ=255)
                                               tCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTCACAACATTGCGCAGGCACAGTGTCCAGTTATTCGCCTCGCCCGCACCAGTCACAGTAATCGTGTTGCCAGTACGtgctg                                                                                                   <  1:106241/150‑1 (MQ=255)
                                                 ccTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTCACAACATTGCGCAGGCACAGTGTCCAGTTATTCGCCTCGCCCGCACCAGTCACAGTAATCGTGTTGCCAGTACGTgctgct                                                                                                 >  1:5114262/1‑150 (MQ=255)
                                                   ttGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTCACAACATTGCGCAGGCACAGTGTCCAGTTATTCGCCTCGCCCGCACCAGTCACAGTAATCGTGTTGCCAGTACGTGCTGCttt                                                                                               <  2:1157098/150‑1 (MQ=255)
                                                     gAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTCACAACATTGCGCAGGCACAGTGTCCAGTTATTCGCCTCGCCCGCACCAGTCACAGTAATCGTGTTGCCAGTACGTGCTGCTTTta                                                                                             <  1:4206590/150‑1 (MQ=255)
                                                     gAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTCACAACATTGCGCAGGCACAGGGTCCAGTTATTCGCCTCGCCCGCACCAGTCACAGTCATCGTGTGGCCAGTACGTGCGGCtgttc                                                                                             >  1:4165929/1‑149 (MQ=255)
                                                        gCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTCACAACATTGCGCAGGCACAGTGTCCAGTTATTCGCCTCGCCCGCACCAGTCACAGTAATCGTGTTGCCAGTACGTGCTGCTTTtaaag                                                                                          <  2:2500668/150‑1 (MQ=255)
                                                        gCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTCACAACATTGCGCAGGCACAGTGTCCAGTTATTCGCCTCGCCCGCACCAGTCACAGTAATCGTGTTGCCAGTACGTGCTGCTTTtaaag                                                                                          <  2:4076925/150‑1 (MQ=255)
                                                            caccacCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTCACAACATTGCGCAGGCACAGTGTCCAGTTATTCGCCTCGCCCGCACCAGTCACAGTAATCGTGTTGCCAGTACGTGCTGCTTTtaaagtaaa                                                                                      >  2:4772542/1‑150 (MQ=255)
                                                             accaccAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTCACAACATTGCGCAGGCACAGTGTCCAGTTATTCGCCTCGCCCGCACCAGTCACAGTAATCGTGTTGCCAGTACGTGCTGCTTTtaaagtaaag                                                                                     <  2:4468460/150‑1 (MQ=255)
                                                                accaGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTCACAACATTGCGCAGGCACAGTGTCCAGTTATTCGCCTCGCCCGCACCAGTCACAGTAATCGTGTTGCCAGTACGTGCTGCTTTTAAAGTAAAGATc                                                                                  <  2:880278/150‑1 (MQ=255)
                                                                  caGCCCCTGCTCTCTTGCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTCACATCATTGCGCAGGCACTGTGTCCAGTTATTCGCCTCGCCCGCACCAGTCACAGTAATCGTGTTGCCAGTACGTGCTGCTTTTAAAGTAAAGATCAc                                                                                <  2:2836999/150‑1 (MQ=255)
                                                                  caGCCCCTGCTCACTATCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTCACAACATTGCGCAGGCACAGTGTCCAGTTATTCGCCTCGCCCGCACCAGTCACAGTAATCGTGTTGCCAGTACGTGCTGCTTTTAAAGTAAAGATCAc                                                                                <  2:803898/150‑1 (MQ=255)
                                                                   aGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTCACAACATTGCGCAGGCACAGTGTCCAGTTATTCGCCTCGCCCGCACCAGTCACAGTAATCGTGTTGCCAGTACGTGCTGCTTTTAAAGTAAAGATCAcc                                                                               <  2:215433/150‑1 (MQ=255)
                                                                   aGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTCACAACATTGCGCAGGCACAGTGTCCAGTTATTCGCCTCGCCCGCACCAGTCACAGTAATCGTGTTGCCAGTACGTGCTGCTTTTAAAGTAAAGATCAcc                                                                               >  2:4030571/1‑150 (MQ=255)
                                                                      cccTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTCACAACATTGCGCAGGCACAGTGTCCAGTTATTCGCCTCGCCCGCACCAGTCACAGTAATCGTGTTGCCAGTACGTGCTGCTTTTAAAGTAAAGATCACCGAt                                                                            <  2:2670742/150‑1 (MQ=255)
                                                                      cccTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTCACAACATTGCGCAGGCACAGTGTCCAGTTATTCGCCTCGCCCGCACCAGTCACAGTAATCGTGTTGCCAGTACGTGCTGCTTTTAAAGTAAAGATCACCGAt                                                                            <  2:4124755/150‑1 (MQ=255)
                                                                         tGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTCACAACATTGCGCAGGCACAGTGTCCAGTTATTCGCCTCGCCCGCACCAGTCACAGTAATCGTGTTGCCAGTACGTGCTGCTTTTAAAGTAAAGATCACCGATCCg                                                                         >  1:1546760/1‑150 (MQ=255)
                                                                          gCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTCACAACATTGCGCAGGCACAGTGTCCAGTTATTCGCCTCGCCCGCACCAGTCACAGTAATCGTGTTGCCAGTACGTGCTGCTTTTAAAGTAAAGATCACCGATCCGt                                                                        <  2:4289484/150‑1 (MQ=255)
                                                                            tCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTCACAACATTGCGCAGGCACAGTGTCCAGTTATTCGCCTCGCCCGCACCAGTCACAGTAATCGTGTTGCCAGTACGTGCTGCTTTTAAAGTAAAGATCACCGATCCGTca                                                                      <  2:3536704/150‑1 (MQ=255)
                                                                            tCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTCACAACATTGCGCAGGCACAGTGTCCAGTTATTCGCCTCGCCCGCACCAGTCACAGTAATCGTGTTGCCAGTACGTGCTGCTTTTAAAGTAAAGATCACCGATCCGTca                                                                      <  1:1936998/150‑1 (MQ=255)
                                                                                          cGAACCGTCTTGCAGACCATTTACTTTCACCACATTGCGCAGGCACAGTGTCCAGTTATTCGCCTCGCCCGCACCAGTCACAGTAATCGTGTTGCCAGTACGTGCTGCTTTTAAAGTAAAGATCACCGATCCGTCAGCAGCGGGCACTTc                                                        <  2:274155/150‑1 (MQ=255)
                                                                                           gAACCGTCTTGCAGACCATTTACTTTCACAACATTGCGCAGGCACAGTGTCCAGTTATTCGCCTCGCCCGTACCAGTCACAGTAATCGTGTTGCCAGTACGTGCTGCTTTTAAAGTAAAGATCACCGATCCGTCAGCAGCGGGCACTAca                                                       >  2:593920/1‑150 (MQ=255)
                                                                                           gAACCGTCTTGCAGACCATTTACTTTCACAACATTGCGCAGGCACAGTGTCCAGTTATTCGCCTCGCCCGCACCAGTCACAGTAATCGTGTTGCCAGTACGTGCTGCTTTTAAAGTAAAGATCACCGATCCGTCAGCAGCGGGCACTTca                                                       >  1:4109846/1‑150 (MQ=255)
                                                                                           gAACCGTCTTGCAGACCATTTACTTTCACAACATTGCGCAGGCACAGCGTCCAGTTATTCGCCTCGCCCGCACCAGTCACAGTAATCGTGTTGCCAGTACGTGCTGCTTTTAAAGTAAAGATCACCGATCCGTCAGCAGCGGGCACTTca                                                       >  1:1017324/1‑150 (MQ=255)
                                                                                              ccGTCTTGCAGACCATTTACTTTCACAACATTGCGCAGGCACAGTGTCCAGTTATTCGCCTCGCCCGCACCAGTCACAGTAATCGTGTTGCCAGTACGTGCTGCTTTTAAAGTAAAGATCACCGATCCGTCAGCAGCGGGCACTTCACAg                                                    >  2:4648630/1‑150 (MQ=255)
                                                                                                  cTTGCAGACCATTTACTTTCACAACATTGCGCAGGCACAGTGTCCAGTTATTCGCCTCGCCCGCACCAGTCACAGTAATCGTGTTGCCAGTACGTGCTGCTTTTAAAGTAAAGATCACCGATCCGTCAGCAGCGGGCACTTCACAGACgg                                                <  1:3862770/150‑1 (MQ=255)
                                                                                                     gCAGACCATTTACTTTCACAACATTGCGCAGGCACAGTGTCCAGTTATTCGCCTCGCCCGCACCAGTCACAGTAATCGTGTTGCCAGTACGTGCTGCTTTTAAAGTAAAGATCACCGATCCGTCAGCAGCGGGCACTTCACAGACGGCtt                                             >  2:1614760/1‑150 (MQ=255)
                                                                                                       aGACCATTTACTTTCACAACATTGCGCAGGCACAGTGTCCAGTTATTCGCCTCGCCCGCACCAGTCACAGTAATCGTGTTGCCAGTACGTGCTGCTTTTAAAGTAAAGATCACCGATCCGTCAGCAGCGGGCACTTCACAGACGGCTTCa                                           >  1:1301215/1‑150 (MQ=255)
                                                                                                       aGACCATTTACTTTCACAACATTGCGCAGGCACAGTGTCCAGTTATTCGCCTCGCCCGCACCAGTCACAGTAATCGTGTTGCCAGTACGTGCTGCTTTTAAAGTAAAGATCACCGATCCGTCAGCAGCGGGCACTTCACAGACGGCTTCa                                           >  2:561619/1‑150 (MQ=255)
                                                                                                        gACCATTTACTTTCACAACATTGCGCAGGCACAGTGTCCAGTTATTCGCCTCGCCCGCACCAGTCACAGTAATCGTGTTGCCAGTACGTGCTGCTTTTAAAGTAAAGATCACCGATCCGTCAGCAGCGGGCACTTCACAGACGGCTTCAt                                          <  1:4643958/150‑1 (MQ=255)
                                                                                                         aCCATTTACTTTCACAACATTGCGCAGGCACAGTGTCCAGTTATTCGCCTCGCCCGCACCAGTCACAGTAATCGTGTTGCCAGTACGTGCTGCTTTTAAAGTAAAGATCACCGATCCGTCAGCAGCGGGCACTTCACAGACGGCTTCATg                                         >  1:4035137/1‑150 (MQ=255)
                                                                                                           cATTTACTTTCACAACATTGCGCAGGCACAGTGTCCAGTTATTCGCCTCGCCCGCACCAGTCACAGTAATCGTGTTGCCAGTACGTGCTGCTTTTAAAGTAAAGATCACCGATCCGTCAGCAGCGGGCACTTCACAGACGGCTTCATGcc                                       >  2:3950878/1‑150 (MQ=255)
                                                                                                            aTTTACTTTCACAACATTGCGCAGGCACAGTGTCCAGTTATTCGCCTCGCCCGCACCAGTCACAGTAATCGTGTTGCCAGTACGTGCTGCTTTTAAAGTAAAGATCACCGATCCGTCAGCAGCGGGCACTTCACAGACGGCTTCATGccc                                      >  1:1485096/1‑150 (MQ=255)
                                                                                                            aTTTACTTTCACAACATTGCGCAGGCACAGTGTCCAGTTATTCGCCTCGCCCGCACCAGTCACAGTAATCGTGTTGCCAGTACGTGCTGCTTTTAAAGTAAAGATCACCGATCCGTCAGCAGCGGGCACTTCACAGACGGCTTCATGccc                                      <  1:951405/150‑1 (MQ=255)
                                                                                                             tttACTTTCACAACATTGCGCAGGCACAGTGTCCAGTTATTCGCCTCGCCCGCACCAGTCACAGTAATCGTGTTGCCAGTACGTGCTGCTTTTAAAGTAAAGATCACCGATCCGTCAGCAGCGGGCACTTCACAGACGGCTTCATGCCCg                                     >  2:5089037/1‑150 (MQ=255)
                                                                                                              ttACTTTCACAACATTGCGCAGGCACAGTGTCCAGTTATTCGCCTCGCCCGCACCAGTCACAGTAATCGTGTTGCCAGTACGTGCTGCTTTTAAAGTAAAGATCACCGATCCGTCAGCAGCGGGCACTTCACAGACGGCTTCATGCCCGt                                    >  2:4355998/1‑150 (MQ=255)
                                                                                                              ttACTTTCACAAAATTGCGCAGGCACAGTGTCAAGTTATTCGCATCGACCGCAACAGTCACAGTAATCGTGTTGAAAGTACGTGCTGCTTTTAAAGTAAAGATCACCGATCCGTCAGCAGCGGGCACTACACAGACTGCTTCATGCCCga                                    >  2:1828513/1‑149 (MQ=255)
                                                                                                                aCTTTCACAACATTGCGCAGGCACAGTGTCCAGTTATTCGCCTCGCCCGCACCAGTCACAGTAATCGTGTTGCCAGTACGTGCTGCTTTTAAAGTAAAGATCACCGATCCGTCAGCAGCGGGCACTTCACAGACGGCTTCATGCCCGTCt                                  >  1:4718439/1‑150 (MQ=255)
                                                                                                                 cTTTCACAACATTGCGCAGGCACAGTGTCCAGTTATTCGCCTCGCCCGCACCAGTCACAGTAATCGTGTTGCCAGTACGTGCTGCTTTTAAAGTAAAGATCACCGATCCGTCAGCAGCGGGCACTTCACAGACGGCTTCATGCCCGTCtt                                 >  1:1705235/1‑150 (MQ=255)
                                                                                                                  tttCACAACATTGCGCAGGCACAGTGTCCAGTTATTCGCCTCGCCCGCACCAGTCACAGTAATCGTGTTGCCAGTACGTGCTGCTTTTAAAGTAAAGATCACCGATCCGTCAGCAGCGGGCACTTCACAGACGGCTTCATGCCCGTCTTg                                <  1:2045143/150‑1 (MQ=255)
                                                                                                                      acaacaTTGCGCAGGCACAGTGTCCAGTTATTCGCCTCGCCCGCACCAGTCACAGTAATCGTGTTGCCAGTACGTGCTGCTTTTAAAGTAAAGATCACCGATCCGTCAGCAGCGGGCACTTCACAGACGGCTTCATGCCCGTCTTGCAGa                            <  1:3671252/150‑1 (MQ=255)
                                                                                                                       caacaTTGCGCAGGCACAGTGTCCAGTTATTCGCCTCGCCCGCACCAGTCACAGTAATCGTGTTGCCAGTACGTGCTGCTTTTAAAGTAAAGATCACCGATCCGTCAGCAGCGGGCACTTCACAGACGGCTTCATGCCCGTCTTGCAGAt                           <  1:2361612/150‑1 (MQ=255)
                                                                                                                       caacaTTGCGCAGGCACAGTGTCCAGTTATTCGCCTCGCCCGCACCAGTCACAGTAATCGTGTTGCCAGTACGTGCTGCTTTTAAAGTAAAGATCACCGATCCGTCAGCAGCGGGCACTTCACAGACGGCTTCATGCCCGTCTTGCAGAt                           >  1:427510/1‑150 (MQ=255)
                                                                                                                       caacaTTGCGCAGGCACAGTGTCCAGTTATTCGCCTCGCCCGCACCAGTCACAGTAATCGTGTTGCCAGTACGTGCTGCTTTTAAAGTAAAGATCACCGATCCGTCAGCAGCGGGCACTTCACAGACGGCTTCATGCCCGTCTTGCAGAt                           >  1:3384098/1‑150 (MQ=255)
                                                                                                                       caacaTTGCGCAGGCACAGTGTCCAGTTATTCGCCTCGCCCGCACCAGTCACAGTAATCGTGTTGCCAGTACGTGCTGCTTTTAAAGTAAAGATCACCGATCCGTCAGCAGCGGGCACTTCACAGACGGCTTCATGCCCGTCTTGCAGAt                           <  2:2030836/150‑1 (MQ=255)
                                                                                                                         acaTTGCGCAGGCACAGTGTCCAGTTATTCGCCTCGCCCGCACCAGTCACAGTAATCGTGTTGCCAGTACGTGCTGCTTTTAAAGTAAAGATCACCGATCCGTCAGCAGCGGGCACTTCACAGACGGCTTCATGCCCGTCTTGCAGATTg                         >  1:1762059/1‑150 (MQ=255)
                                                                                                                            ttGCGCAGGCACAGTGTCCAGTTATTCGCCTCGCCCGCACCAGTCACAGTAATCGTGTTGCCAGTACGTGCTGCTTTTAAAGTAAAGATCACCGATCCGTCAGCAGCGGGCACTTCACAGACGGCTTCATGCCCGTCTTGCAGATTGAag                      <  1:4532240/150‑1 (MQ=255)
                                                                                                                            ttGCGCAGGCACAGTGTCCAGTTATTCGCCTCGCCCGCACCAGTCACAGTAATCGTGTTGCCAGTACGTGCTGCTTTTAAAGTAAAGATCACCGATCCGTCAGCAGCGGGCACTTCACAGACGGCTTCATGCCCGTCTTGCAGATTGAag                      <  2:3468923/150‑1 (MQ=255)
                                                                                                                               cgcAGGCACAGTGTCCAGTTATTCGCCTCGCCCGCACCAGTCACAGTAATCGTGTTGCCAGTACGTGCTGCTTTTAAAGTAAAGATCACCGATCCGTCAGCAGCGGGCACTTCACAGACGGCTTCATGCCCGTCTTGCAGATTGAAGAgg                   >  2:2538963/1‑150 (MQ=255)
                                                                                                                                gcAGGCACAGTGTCCAGTTATTCGCCTCGCCCGCACCAGTCACAGTAATCGTGTTGCCAGTACGTGCTGCTTTTAAAGTAAAGATCACCGATCCGTCAGCAGCGGGCACTTCACAGACGGCTTCATGCCCGTCTTGCAGATTGAAGAGGt                  >  1:3020413/1‑150 (MQ=255)
                                                                                                                                gcAGGCACAGTGTCCAGTTATTCGCCTCGCCCGCACCAGTCACAGTAATCGTGTTGCCAGTACGTGCTGCTTTTAAAGTAAAGATCACCGATCCGTCAGCAGCGGGCACTTCACAGACGGCTTCATGCCCGTCTTGCAGATTGAAGAGGt                  <  2:3239219/150‑1 (MQ=255)
                                                                                                                                 cAGGCACAGTGTCCAGTTATTCGCCTCGCCCGCACCAGTCACAGTAATCGTGTTGCCAGTACGTGCTGCTTTTAAAGTAAAGATCACCGATCCGTCAGCAGCGGGCACTTCACAGACGGCTTCATGCCCGTCTTGCAGATTGAAGAGGTg                 <  1:2973046/150‑1 (MQ=255)
                                                                                                                                    gCACAGTGTCCAGTTATTCGCCTCGCCCGCACCAGTCACAGTAATCGTGTTGCCAGTACGTGCTGCTTTTAAAGTAAAGATCACCGATCCGTCAGCAGCGGGCACTTCACAGACGGCTTCATGCCCGTCTTGCAGATTGAAGAGGTGGaa              >  1:3830518/1‑150 (MQ=255)
                                                                                                                                     cacaGTGTCCAGTTATTCGCCTCGCCCGCACCAGTCACAGTAATCGTGTTGCCAGTACGTGCTGCTTTTAAAGTAAAGATCACCGATCCGTCAGCAGCGGGCACTTCACAGACGGCTTCATGCCCGTCTTGCAGATTGAAGAGGTGGAAt             >  1:1183705/1‑150 (MQ=255)
                                                                                                                                      acaGTGTCCAGTTATTCGCCTCGCCCGCACCAGTCACAGTAATCGTGTTGCCAGTACGTGCTGCTTTTAAAGTAAAGATCACCGATCCGTCAGCAGCGGGCACTTCACAGACGGCTTCATGCCCGTCTTGCAGATTGAAGAGGTGGAATg            >  2:2932092/1‑150 (MQ=255)
                                                                                                                                      acaGTGTCCAGTTATTCGCCTCGCCCGCACCAGTCACAGTAATCGTGTTGCCAGTACGTGCTGCTTTTAAAGTAAAGATCACCGATCCGTCAGCAGCGGGCACTTCACAGACGGCTTCATGCCCGTCTTGCAGATTGAAGAGGTGGAATg            <  1:774675/150‑1 (MQ=255)
                                                                                                                                        aGTGTCCAGTTATTCGCCTCGCCCGCACCAGTCACAGTAATCGTGTTGCCAGTACGTGCTGCTTTTAAAGTAAAGATCACCGATCCGTCAGCAGCGGGCACTTCACAGACGGCTTCATGCCCGTCTTGCAGATTGAAGAGGTGGAATGcc          >  2:3600538/1‑150 (MQ=255)
                                                                                                                                           gtCCAGTTATTCGCCTCGCCCGCACCAGTCACAGTAATCGTGTTGCCAGTACGTGCTGCTTTTAAAGTAAAGATCACCGATCCGTCAGCAGCGGGCACTTCACAGACGGCTTCATGCCCGTCTTGCAGATTGAAGAGGTGGAATGCCGTg       >  1:272839/1‑150 (MQ=255)
                                                                                                                                            tCCAGTTATTCGCCTCGCCCGCACCAGTCACAGTAATCGTGTTGCCAGTACGTGCTGCTTTTAAAGTAAAGATCACCGATCCGTCAGCAGCGGGCACTTCACAGACGGCTTCATGCCCGTCTTGCAGATTGAAGAGGTGGAATGCCGTGc      >  1:2521522/1‑150 (MQ=255)
                                                                                                                                             ccAGTTATTCGCCTCGCCCGCACCAGTCACAGTAATCGTGTTGCCAGTACGTGCTGCTTTTAAAGTAAAGATCACCGATCCGTCAGCAGCGGGCACTTCACAGACGGCTTCATGCCCGTCTTGCAGATTGAAGAGGTGGAATGCCGTGcc     <  1:2836848/150‑1 (MQ=255)
                                                                                                                                              cAGTTATTCGCCTCGCCCGCACCAGTCACAGTAATCGTGTTGCCAGTACGTGCTGCTTTTAAAGTAAAGATCACCGATCCGTCAGCAGCGGGCACTTCACAGACGGCTTCATGCCCGTCTTGCAGATTGAAGAGGTGGAATGCCGTGCCt    >  1:5182550/1‑150 (MQ=255)
                                                                                                                                                gTTATTCGCCTCGCCCGCACCAGTCACAGTTATAGTGTTGCCAGTACGTGCTGCTTTTAAAGTAAAGTTCACCGATCCGTCAGCAGCGGGCACTTCACAGACGGCTTCATGCCCGTCTTGCAGATTGAAGAGGTGGAATGCCGTGCCTTc  <  1:1728859/150‑1 (MQ=255)
                                                                                                                                                gTTATTCGCCTCGCCCGCACCAGTCACAGTAATCGTGTTGCCAGTACGTGCTGCTTTTAAAGTAAAGATCACCGATCCGTCAGCAGCGGGCACTTCACAGACGGCTTCATGCCCGTCTTGCAGATTGAAGAGGTGGAATGCCGTGCCTTc  <  2:1062226/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
CGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTCACAACATTGCGCAGGCACAGTGTCCAGTTCTTCGCCTCGCCCGCACCAGTCACAGTAATCGTGTTGCCAGTACGTGCTGCTTTTAAAGTAAAGATCACCGATCCGTCAGCAGCGGGCACTTCACAGACGGCTTCATGCCCGTCTTGCAGATTGAAGAGGTGGAATGCCGTGCCTTC  >  CP009273/3825555‑3825848

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: