Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,081,697 T→A 16.2% C19S (TGC→AGC)  galP → D‑galactose transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,081,6970TA16.2% 782.9 / 131.6 284C19S (TGC→AGC) galPD‑galactose transporter
Reads supporting (aligned to +/- strand):  ref base T (120/118);  new base A (30/16);  total (150/134)
Fisher's exact test for biased strand distribution p-value = 7.65e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.96e-01

Only 100 of 284 total aligned reads displayed.

AAGCGATTACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAATAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCTGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCTGGATATCGGTGTAATTGCTGGCGCACTGCCGTTTATTGCAGATGAATTCCAGATTACTTCGCACACGCAAGAATGGGTCGTAAGCTCCATGATGTTCGGTGCGGCAGTCGG  >  CP009273/3081549‑3081845
                                                                                                                                                    |                                                                                                                                                    
 aGCGATTACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAATAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCTGc                                                                                                                                                    >  1:2948423/1‑150 (MQ=255)
 aGCGATTACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAATAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCTGc                                                                                                                                                    >  2:3791736/1‑150 (MQ=255)
 aGCGATTACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAATAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTAGTATGc                                                                                                                                                    >  1:3682172/1‑150 (MQ=255)
    gATTACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAATAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCAGCTTc                                                                                                                                                 >  2:2426828/1‑150 (MQ=255)
    gATTACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAATAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCAGCTTc                                                                                                                                                 >  1:4158942/1‑150 (MQ=255)
      ttACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAATAAAAAATAACCATATTGGAGGGCATCATGACTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCTGCTTCCt                                                                                                                                               >  2:1107267/1‑150 (MQ=255)
       tACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAATAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCTGCTTCCtt                                                                                                                                              <  1:2714575/150‑1 (MQ=255)
       tACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAATAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCTGCTTCCtt                                                                                                                                              >  1:4313109/1‑150 (MQ=255)
       tACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAATAAAAAATAACCAGATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCTGCTTCCtt                                                                                                                                              >  1:101311/1‑150 (MQ=255)
       cacacTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAATAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCTGCTTCCtt                                                                                                                                              >  1:2333244/2‑150 (MQ=255)
           cTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAATAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCTGCTTCCTTGCCg                                                                                                                                          >  2:3456669/1‑150 (MQ=255)
           cTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAATAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCTGCTTCCTTGCCg                                                                                                                                          >  2:3904251/1‑150 (MQ=255)
              aTGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAATAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCTGCTTCCTTGCCGctc                                                                                                                                       <  1:2018406/150‑1 (MQ=255)
              aTGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAATAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCAGCTTCCTTGCCGctc                                                                                                                                       >  2:2842526/1‑150 (MQ=255)
              aTGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAATAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCAGCTTCCTTGCCGctc                                                                                                                                       >  2:1917596/1‑150 (MQ=255)
               tgtgATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAATAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCAGCTTCCTTGCCGctct                                                                                                                                      >  1:1941162/1‑150 (MQ=255)
               tgtgATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAATAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGGTTTTCGTCTGCTTCCATGCCGctct                                                                                                                                      >  2:1872839/1‑150 (MQ=255)
                        cTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAATAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCTGCTTCCTTGCCGCTCTGGCGGGAtt                                                                                                                             >  2:33583/1‑150 (MQ=255)
                        cTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAATAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCTGCTTCCTTGCCGCTCTGGCGGGAtt                                                                                                                             >  2:4251237/1‑150 (MQ=255)
                          tCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAATAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCTGCTTCCTTGCCGCTCTGGCGGGATTAc                                                                                                                           <  2:4351386/150‑1 (MQ=255)
                             caTCTTTTTACGTCGTACTCACCTATCTTAATTCACAATAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCTGCTTCCTTGCCGCTCTGGCGGGATTActct                                                                                                                        <  2:4505924/150‑1 (MQ=255)
                              aTCTTTTTACGTCGTACTCACCTATCTTAATTCACAATAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCTGCTTCCTTGCCGCTCTGGCGGGATTACTCtt                                                                                                                       <  1:2640580/150‑1 (MQ=255)
                                    ttACGTCGTACTCACCTATCTTAATTCACAATAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCTGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCt                                                                                                                 <  2:1557785/150‑1 (MQ=255)
                                    ttACGTCGTACTCACCTATCTTAATTCACAATAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCTGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCt                                                                                                                 <  2:1975820/150‑1 (MQ=255)
                                    ttACGTCGTACTCACCTATCTTAATTCACAATAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCTGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCt                                                                                                                 <  2:2569506/150‑1 (MQ=255)
                                    ttACGTCGTACTCACCTATCTTAATTCACAATAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCAGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCt                                                                                                                 <  1:2864879/150‑1 (MQ=255)
                                        gtcgtACTCACCTATCTTAATTCACAATAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCAGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCTGGat                                                                                                             <  2:736256/150‑1 (MQ=255)
                                           gtACTCACCTATCTTAATTCACAATAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCTGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCTGGATATc                                                                                                          >  2:342171/1‑150 (MQ=255)
                                           gtACTCACCTATCTTAATTCACAATAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCTGCTTCCTTGCCGCTCTGGCGGGAGTACTCTTTGGCCTGGATATc                                                                                                          >  1:188407/1‑150 (MQ=255)
                                            tACTCACCTATCTTAATTCACAATAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCTGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCTGGATATCg                                                                                                         <  1:1632050/150‑1 (MQ=255)
                                              cTCACCTATCTTAATTCACAATAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCTGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCTGGATATCGgt                                                                                                       <  1:1347179/150‑1 (MQ=255)
                                               tCACCTATCTTAATTCACAATAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCTGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCTGGATATCGgtg                                                                                                      >  2:2266838/1‑150 (MQ=255)
                                                 aCCTATCTTAATTCACAATAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCTGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCTGGATATCGGTGTa                                                                                                    >  1:3458081/1‑150 (MQ=255)
                                                   cTATCTTAATTCACAATAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCTGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCTGGATATCGGTGTAAt                                                                                                  >  2:59809/1‑150 (MQ=255)
                                                   cTATCTTAATTCACAATAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCTGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCTGGATATCGGTGTAAt                                                                                                  <  1:2150358/150‑1 (MQ=255)
                                                      tCTTAATTCACAATAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCTGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCTGGATATCGGTGTAATTGc                                                                                               >  2:2185335/1‑150 (MQ=255)
                                                        ttAATTCACAATAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCTGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCTGGATATCGGTGTAATTGCTg                                                                                             <  1:3945003/150‑1 (MQ=255)
                                                          aaTTCACAATAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCTGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCTGGATATCGGTGTAATTGCTGgc                                                                                           <  1:3632859/150‑1 (MQ=255)
                                                          aaTTCACAATAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCTGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCTGGATATCGGTGTAATTGCTGgc                                                                                           <  1:1986860/150‑1 (MQ=255)
                                                          aaTTCACAATAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCTGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCTGGATATCGGTGTAATTGCTGgc                                                                                           >  2:3764207/1‑150 (MQ=255)
                                                           aTTCACAATAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCTGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCTGGATATCGGTGTAATTGCTGgcg                                                                                          >  2:4548293/1‑150 (MQ=255)
                                                           aTTCACAATAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCAGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCTGGATATCGGTGTAATTGCTGgcg                                                                                          >  2:2079689/1‑150 (MQ=255)
                                                              cacaATAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCTGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCTGGATATCGGTGTAATTGCTGGCGCAc                                                                                       >  1:888481/1‑150 (MQ=255)
                                                               acaATAAAAAATAACCATCTTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCTGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCTGGATATCGGTGTAATTGCTGGCGCACt                                                                                      <  2:4150189/150‑1 (MQ=255)
                                                               acaATAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCTGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCTGGATATCGGTGTAATTGCTGGCGCACt                                                                                      >  1:480014/1‑150 (MQ=255)
                                                                  aTAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCTGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCTGGATATCGGTGTAATTGCTGGCGCACTGcc                                                                                   <  1:3309503/150‑1 (MQ=255)
                                                                  aTAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCTGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCTGGATATCGGTGTAATTGCTGGCGCACTGcc                                                                                   <  2:4515425/150‑1 (MQ=255)
                                                                      aaaaTAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCTGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCTGGATATCGGTGTAATTGCTGGCGCACTGCCGttt                                                                               <  1:3536953/150‑1 (MQ=255)
                                                                       aaaTAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCTGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCTGGATATCGGTGTAATTGCTGGCGCACTGCCGTTTa                                                                              >  2:3694243/1‑150 (MQ=255)
                                                                       aaaTAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCTGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCTGGATATCGGTGTAATTGCTGGCGCACTGCCGTTTa                                                                              >  1:2923302/1‑150 (MQ=255)
                                                                         aTAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCTGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCTGGATATCGGTGTAATTGCTGGCGCACTGCCGTTTAtt                                                                            <  1:666178/150‑1 (MQ=255)
                                                                         aTAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCTGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCTGGATATCGGTGTAATTGCTGGCGCACTGCCGTTTAtt                                                                            <  2:3143738/150‑1 (MQ=255)
                                                                           aaCCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCTGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCTGGATATCGGTGTAATTGCTGGCGCACTGCCGTTTATTGc                                                                          >  1:2024164/1‑150 (MQ=255)
                                                                            aCCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCTGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCTGGATATCGGTGTAATTGCTGGCGCACTGCCGTTTATTGCa                                                                         <  1:4286711/150‑1 (MQ=255)
                                                                             ccATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCTGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCTGGATATCGGTGTAATTGCTGGCGCACTGCCGTTTATTGCAg                                                                        >  2:2996304/1‑150 (MQ=255)
                                                                             ccATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCTGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCTGGATATCGGTGTAATTGCTGGCGCACTGCCGTTTATTGCAg                                                                        >  2:470366/1‑150 (MQ=255)
                                                                              cATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCTGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCTGGATATCGGTGTAATTGCTGGCGCACTGCCGTTTATTGCAGa                                                                       <  2:4482459/150‑1 (MQ=255)
                                                                              cATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCTGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCTGGATATCGGTGTAATTGCTGGCGCACTGCCGTTTATTGCAGa                                                                       >  2:4078827/1‑150 (MQ=255)
                                                                                 atTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCTGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCTGGATATCGGTGTAATTGCTGGCGCACTGCCGTTTATTGCAGATGa                                                                    >  2:2781887/1‑150 (MQ=255)
                                                                                 atTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCTGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCTGGATATCGGTGTAATTGCTGGCGCACTGCCGTTTATTGCAGATGa                                                                    >  1:4332637/1‑150 (MQ=255)
                                                                                 atTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCAGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCTGGATATCGGTGTAATTGCTGGCGCACTGCCGTTTATTGCAGATGa                                                                    >  1:1395752/1‑150 (MQ=255)
                                                                                   tGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCTGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCTGGATATCGGTGTAATTGCTGGCGCACTGCCGTTTATTGCAGATGAAt                                                                  >  2:918598/1‑150 (MQ=255)
                                                                                   tGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCTGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCTGGATATCGGTGTAATTGCTGGCGCACTGCCGTTTATTGCAGATGAAt                                                                  >  1:3876230/1‑150 (MQ=255)
                                                                                   tGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCTGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCTGGATATCGGTGTAATTGCTGGCGCACTGCCGTTTATTGCAGATGAAt                                                                  <  1:4419737/150‑1 (MQ=255)
                                                                                       gggCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCAGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCTGGATATCGGTGTAATTGCTGGCGCACTGCCGTTTATTGCAGATGAATTCCa                                                              >  2:1739716/1‑150 (MQ=255)
                                                                                        ggCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCTGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCTGGATATCGGTGTAATTGCTGGCGCACTGCCGTTTATTGCAGATGAATTCCAg                                                             <  1:2482404/150‑1 (MQ=255)
                                                                                           atcatGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCTGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCTGGATATCGGTGTAATTGCTGGCGCACTGCCGTTTATTGCAGATGAATTCCAGAtt                                                          >  1:647613/1‑150 (MQ=255)
                                                                                             catGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCAGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCTGGATATCGGTGTAATTGCTGGCGCACTGCCGTTTATTGCAGATGAATTCCAGATTAc                                                        <  2:3571791/150‑1 (MQ=255)
                                                                                              atGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCTGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCTGGATATCGGTGTAATTGCTGGCGCACTGCCGTTTATTGCAGATGAATTCCAGATTACt                                                       <  2:4437395/150‑1 (MQ=255)
                                                                                              atGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCTGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCTGGATATCGGTGTAATTGCTGGCGCACTGCCGTTTATTGCAGATGAATTCCAGATTACt                                                       >  1:4328381/1‑150 (MQ=255)
                                                                                              atGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCTGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCTGGATATCGGTGTAATTGCTGGCGCACTGCCGTGTATTGCAGATGTATTCCAGATTACt                                                       >  2:2856665/1‑150 (MQ=255)
                                                                                              atGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCAGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCTGGATATCGGTGTAATTGCTGGCGCACTGCCGTTTATTGCAGATGAATTCCAGATTACt                                                       >  2:1358988/1‑150 (MQ=255)
                                                                                               tGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCTGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCTGGATATCGGTGTAATTGCTGGCGCACTGCCGTTTATTGCAGATGAATTCCAGATTACtt                                                      <  2:2134558/150‑1 (MQ=255)
                                                                                               tGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCTGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCTGGATATCGGTGTAATTGCTGGCGCACTGCCGTTTATTGCAGATGAATTCCAGATTACtt                                                      <  1:144953/150‑1 (MQ=255)
                                                                                               tGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCTGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCTGGATATCGGTGTAATTGCTGGCGCACTGCCGTTTATTGCAGATGAATTCCAGATTACtt                                                      <  1:3243968/150‑1 (MQ=255)
                                                                                                        cTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCTGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCTGGATATCGGTGTAATTGCTGGCGCACTGCCGTTTATTGCAGATGAATTCCAGATTACTTCGCACACGc                                             >  1:4562552/1‑150 (MQ=255)
                                                                                                         tAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCTGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCTGGATATCGGTGTAATTGCTGGCGCACTGCCGTTTATTGCAGATGAATTCCAGATTACTTCGCACACGCa                                            >  2:3517063/1‑150 (MQ=255)
                                                                                                           aaaaaCAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCTGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCTGGATATCGGTGTAATTGCTGGCGCACTGCCGTTTATTGCAGATGAATTCCAGATTACTTCGCACACGCAAg                                          <  2:3123727/150‑1 (MQ=255)
                                                                                                           aaaaaCAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCTGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCTGGATATCGGTGTAATTGCTGGCGCACTGCCGTTTATTGCAGATGAATTCCAGATTACTTCGCACACGCAAg                                          >  2:2677400/1‑150 (MQ=255)
                                                                                                           aaaaaCAGGGGCGGTCAAACAAGGCAATGACGTATTTCGTCTGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCTGGATATCGGTGTAATTGCTGGCGCACTGCCGTTTATTGCAGATGAATTCCAGATTACTTCGAACACGCAAg                                          >  2:2741357/1‑150 (MQ=255)
                                                                                                            aaaaCAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCAGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCTGGATATCGGTGTAATTGCTGGCGCACTGCCGTTTATTGCAGATGAATTCCAGATTACTTCGCACACGCAAGa                                         >  2:1038623/1‑150 (MQ=255)
                                                                                                                  ggggCGGTCAAACAAGGCAATGACGTTTTTCGTCAGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCTGGATATCGGTGTAATTGCTGGCGCACTGCCGTTTATTGCAGATGAATTCCAGATTACTTCGCACACGCAAGAATGGGt                                   >  1:4021877/1‑150 (MQ=255)
                                                                                                                   gggCGGTCAAACAAGGCAATGACGTTTTTCGTCTGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCTGGATATCGGTGTAATTGCTGGCGCACTGCCGTTTATTGCAGATGAATTCCAGATTACTTCGCACACGCAAGAATGGGTc                                  <  1:3000696/150‑1 (MQ=255)
                                                                                                                      cGGTCAAACAAGGCAATGACGTTTTTCGTCTGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCTGGATATCGGTGTAATTGCTGGCGCACTGCCGTTTATTGCAGATGAATTCCAGATTACTTCGCACACGCAAGAATGGGTCGTa                               >  2:125617/1‑150 (MQ=255)
                                                                                                                        gTCAAACAAGGCAATGACGTTTTTCGTCAGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCTGGATATCGGTGTAATTGCTGGCGCACTGCCGTTTATTGCAGATGAATTCCAGATTACTTCGCACACGCAAGAATGGGTCGTAAg                             >  2:1944939/1‑150 (MQ=255)
                                                                                                                         tCAAACAAGGCAATGACGTTTTTCGTCTGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCTGGATATCGGTGTAATTGCTGGCGCACTGCCGTTTATTGCAGATGAATTCCAGATTACTTCGCACACGCAAGAATGGGTCGTAAGc                            >  1:697326/1‑150 (MQ=255)
                                                                                                                          cAAACAAGGCAATGACGTTTTTCGTCTGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCTGGATATCGGTGTAATTGCTGGCGCACTGCCGTTTATTGCAGATGAATTCCAGATTACTTCGCACACGCAAGAATGGGTCGTAAGCt                           <  1:2266838/150‑1 (MQ=255)
                                                                                                                          cAAACAAGGCAATGACGTTTTTCGTCAGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCTGGATATCGGTGTAATTGCTGGCGCACTGCCGTTTATTGCAGATGAATTCCAGATTACTTCGCACACGCAAGAATGGGTCGTAAGCt                           <  1:566101/150‑1 (MQ=255)
                                                                                                                            aaCAAGGCTATGACGTTTTTCGTCAGCTTCCTTGCCGCTCTGGCGGGTTTACTCTTTGGCCTGGATATCGGTGTAATTGCTGGCGCACTGCCGTTTATTGCAGATGAATTCCAGATTACTTCGCACACGCAAGAATGGGTCGTAAGCTcc                         <  1:1546013/150‑1 (MQ=255)
                                                                                                                            aaCAAGGCAATGACGTTTTTCGTCTGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCTGGATATCGGTGTAATTGCTGGCGCACTGCCGTTTATTGCAGATGAATTCCAGATTACTTCGCACACGCAAGAATGGGTCGTAAGCTcc                         >  1:369921/1‑150 (MQ=255)
                                                                                                                             aCAAGGCAATGACGTTTTTCGTCTGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCTGGATATCGGTGTAATTGCTGGCGCACTGCCGTTTATTGCAGATGAATTCCAGATTACTTCGCACACGCAAGAATGGGTCGTAAGCTCCa                        <  2:2933903/150‑1 (MQ=255)
                                                                                                                             aCAAGGCAATGACGTTTTTCGTCTGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCTGGATATCGGTGTAATTGCTGGCGCACTGCCGTTTATTGCAGATGAATTCCAGATTACTTCGCACACGCAAGAATGGGTCGTAAGCTCCa                        <  2:690105/150‑1 (MQ=255)
                                                                                                                              cAAGGCAATGACGTTTTTCGTCTGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCTGGATATCGGTGTAATTGCTGGCGCACTGCCGTTTATTGCAGATGAATTCCAGATTACTTCGCACACGCAAGAATGGGTCGTAAGCTCCat                       >  1:4605691/1‑150 (MQ=255)
                                                                                                                                   cAATGACGTTTTTCGTCTGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCTGGATATCGGTGTAATTGCTGGCGCACTGCCGTTTATTGCAGATGAATTCCAGATTACTTCGCACACGCAAGAATGGGTCGTAAGCTCCATGATGt                  <  2:3611661/150‑1 (MQ=255)
                                                                                                                                        aCGTTTTTCGTCTGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCTGGATATCGGTGTAATTGCTGGCGCACTGCCGTTTATTGCAGATGAATTCCAGATTACTTCGCACACGCAAGAATGGGTCGTAAGCTCCATGATGTTCGGt             >  1:1120835/1‑150 (MQ=255)
                                                                                                                                          gTTTTTCGTCTGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCTGGATATCGGTGTAATTGCTGGCGCACTGCCGTTTATTGCAGATGAATTCCAGATTACTTCGCACACGCAAGAATGGGTCGTAAGCTCCATGATGTTCGGTGc           <  1:4510658/150‑1 (MQ=255)
                                                                                                                                            ttttCGTCTGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCTGGATATCGGTGTAATTGCTGGCGCACTGCCGTTTATTGCAGATGAATTCCAGATTACTTCGCACACGCAAGAATGGGTCGTAAGCTCCATGATGTTCGGTGCgg         >  2:4203829/1‑150 (MQ=255)
                                                                                                                                              ttCGTCTGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCTGGATATCGGTGTAATTGCTGGCGCACTGCCGTTTATTGCAGATGAATTCCAGATTACTTCGCACACGCAAGAATGGGTCGTAAGCTCCATGATGTTCGGTGCGGCa       <  1:4328750/150‑1 (MQ=255)
                                                                                                                                                cGTCTGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCTGGATATCGGTGTAATTGCTGGCGCACTGCCGTTTATTGCAGATGAATTCCAGATTACTTCGCACACGCAAGAATGGGTCGTAAGCTCCATGATGTTCGGTGCGGCAGt     <  1:4071690/150‑1 (MQ=255)
                                                                                                                                                 gTCTGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCTGGATATCGGTGTAATTGCTGGCGCACTGCCGTTTATTGCAGATGAATTCCAGATTACTTCGCACACGCAAGAATGGGTCGTAAGCTCCATGATGTTCGGTGCGGCAGTc    >  2:226249/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
AAGCGATTACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAATAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCTGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCTGGATATCGGTGTAATTGCTGGCGCACTGCCGTTTATTGCAGATGAATTCCAGATTACTTCGCACACGCAAGAATGGGTCGTAAGCTCCATGATGTTCGGTGCGGCAGTCGG  >  CP009273/3081549‑3081845

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: ‑