Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,472,200 A→T 11.1% S485S (TCT→TCA valS ← valyl‑tRNA synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,472,2000AT11.1% 950.8 / 93.0 308S485S (TCT→TCAvalSvalyl‑tRNA synthetase
Reads supporting (aligned to +/- strand):  ref base A (138/131);  new base T (14/20);  total (157/151)
Fisher's exact test for biased strand distribution p-value = 2.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.36e-01

Only 100 of 315 total aligned reads displayed.

CGGTTTGCCATTTTCATCTTTGATGAAGTGCATGGTCATCATGATCATGCGGGCAATCCAGAAGAAAATGATGTCGAAACCAGATACCATCACGCTGGTTGGGTGGAACTGACGCAGGGCGTCGGTATTTTCCGGCCAGCCAAGGGTAGAGAAGGTCCACAGCGCAGAAGAGAACCAGGTATCGAGAACGTCTTCGTCCTGACGCAGGACAACATCAGCACCGAGGTTATTTTCTTTACGCACTTCGTCTTCGTTGCGGCCAACATAAACGTTACCCGCTTCGTCATACCATGCCGG  >  CP009273/4472053‑4472349
                                                                                                                                                   |                                                                                                                                                     
cGGTTTGCCATTTTCATCTTTGATGAAGTGCATGGTCATCATGATCATGCGGGCAATCCAGAAGAAAATGATGTCGAAACCAGATACCATCACGCTGGTTGGGTGGAACTGACGCAGGGCGTCGGTATTTTCCGGCCAGCCAAGGGTaga                                                                                                                                                     >  2:2584507/1‑150 (MQ=255)
 ggTTTGCCATTTTCATCTTTGATGAAGTGCATGGTCATCATGATCATGCGGGCAATCCAGAAGAAAATGATGTCGAAACCAGATACCATCACGCTGGTTGGGTGGAACTGACGCAGGGCGTCGGTATTTTCCGGCCAGCCAAGGGTagag                                                                                                                                                    >  1:3467333/1‑150 (MQ=255)
    ttGCCATTTTCATCTTTGATGAAGTGCATGGTCATCATGATCATGCGGGCAATCCAGAAGAAAATGATGTCGAAACCAGATACCATCACGCTGGTTGGGTGGAACTGACGCAGGGCGTCGGTATTTTCCGGCCAGCCAAGGGTTGAGAAg                                                                                                                                                 >  1:3309685/1‑150 (MQ=255)
    ttGCCATTTTCATCTTTGATGAAGTGCATGGTCATCATGATCATGCGGGCAATCCAGAAGAAAATGATGTCGAAACCAGATACCATCACGCTGGTTGGGTGGAACTGACGCAGGGCGTCGGTATTTTCCGGCCAGCCAAGGGTAGAGAAg                                                                                                                                                 <  1:2143669/150‑1 (MQ=255)
      gCCATTTTCATCTTTGATGAAGTGCATGGTCATCATGATCATGCGGGCAATCCAGAAGAAAATGATGTCGAAACCAGATACCATCACGCTGGTTGGGTGGAACTGACGCAGGGCGTCGGTATTTTCCGGCCAGCCAAGGGTCGAGAAGGt                                                                                                                                               >  2:3907176/1‑150 (MQ=255)
      gCCATTTTCATCTTTGATGAAGTGCATGGTCATCATGATCATGCGGGCAATCCAGAAGAAAATGATGTCGAAACCAGATACCATCACGCTGGTTGGGTGGAACTGACGCAGGGCGTCGGTATTTTCCGGCCAGCCAAGGGTAGAGAAGGt                                                                                                                                               <  2:3928545/150‑1 (MQ=255)
      gCCATTTTCATCTTTGATGAAGTGCATGGTCATCATGATCATGCGGGCAATCCAGAAGAAAATGATGTCGAAACCAGATACCATCACGCTGGTTGGGTGGAACTGACGCAGGGCGTCGGTATTTTCCGGCCAGCCAAGGGTAGAGAAGGt                                                                                                                                               <  2:3062735/150‑1 (MQ=255)
        cATTTTCATCTTTGATGAAGTGCATGGTCATCATGATCATGCGGGCAATCCAGAAGAAAATGATGTCGAAACCAGATACCATCACGCTGGTTGGGTGGAACTGACGCAGGGCGTCGGTATTTTCCGGCCAGCCAAGGGTAGAGAAGGTcc                                                                                                                                             >  1:1896742/1‑150 (MQ=255)
          ttttCATCTTTGATGAAGTGCATGGTCATCATGATCATGCGGGCAATCCAGAAGAAAATGATGTCGAAACCAGATACCATCACGCTGGTTGGGTGGAACTGACGCAGGGCGTCGGTATTTTCCGGCCAGCCAAGGGTAGAGAAGGTCcac                                                                                                                                           <  1:1965048/150‑1 (MQ=255)
            ttCATCTTTGATGAAGTGCATGGTCATCATGATCATGCGGGCAATCCAGAAGAAAATGATGTCGAAACCAGATACCATCACGCTGGTTGGGTGGAACTGACGCAGGGCGTCGGTATTTTCCGGCCAGCCAAGGGTAGAGAAGGTCCACAg                                                                                                                                         >  1:280191/1‑150 (MQ=255)
            ttCATCTTTGATGAAGTGCATGGTCATCATGATCATGCGGGCAATCCAGAAGAAAATGATGTCGAAACCAGATACCATCACGCTGGTTGGGTGGAACTGACGCAGGGCGTCGGTATTTTCCGGCCAGCCAAGGGTAGAGAAGGTCCACAg                                                                                                                                         >  2:3776468/1‑150 (MQ=255)
            ttCATCTTTGATGAAGTGCATGGTCATCATGATCATGCGGGCAATCCAGAAGAAAATGATGTCGAAACCAGATACCATCACGCTGGTTGGGTGGAACTGACGCAGGGCGTCGGTATTTTCCGGCCAGCCAAGGGTAGAGAAGGTCCACAg                                                                                                                                         <  1:560423/150‑1 (MQ=255)
                 cTTTGATGAAGTGCATGGTCATCATGATCATGCGGGCAATCCAGAAGAAAATGATGTCGAAACCAGATACCATCACGCTGGTTGGGTGGAACTGACGCAGGGCGTCGGTATTTTCCGGCCAGCCAAGGGTAGAGAAGGTCCACAGCGCag                                                                                                                                    >  1:537731/1‑150 (MQ=255)
                     gatgaAGTGCATGGTCATCATGATCATGCGGGCAATCCAGAAGAAAATGATGTCGAAACCAGATACCATCACGCTGGTTGGGTGGAACTGACGCAGGGCGTCGGTATTTTCCGGCCAGCCAAGGGTAGAGAAGGTCCACAGCGCagaaga                                                                                                                                >  2:3805565/1‑150 (MQ=255)
                          aGTGCATGGTCATCATGATCATGCGGGCAATCCAGAAGAAAATGATGTCGAAACCAGATACCATCACGCTGGTTGGGTGGAACTGACGCAGGGCGTCGGTATTTTCCGGCCAGCCAAGGGTAGAGAAGGTCCACAGCGCAGAAGAGAAcc                                                                                                                           >  1:3805028/1‑150 (MQ=255)
                            tGCATGGTCATCATGATCATGCGGGCAATCCAGAAGAAAATGATGTCGAAACCAGATACCATCACGCTGGTTGGGTGGAACTGACGCAGGGCGTCGGTATTTTCCGGCCAGCCAAGGGTAGAGAAGGTCCACAGCGCAGAAGAGAACCAg                                                                                                                         >  1:1995457/1‑150 (MQ=255)
                            tGCATGGTCATCATGATCATGCGGGCAATCCAGAAGAAAATGATGTCGAAACCAGATACCATCACGCTGGTTGGGTGGAACTGACGCAGGGCGTCGGTATTTTCCGGCCAGCCAAGGGTAGAGAAGGTCCACAGCGCAGAAGAGAACCAg                                                                                                                         <  1:202219/150‑1 (MQ=255)
                            tGCATGGTCATCATGATCATGCGGGCAATCCAGAAGAAAATGATGTCGAAACCAGATACCATCACGCTGGTTGGGTGGAACTGACGCAGGGCGTCGGTATTTTCCGGCCAGCCAAGGGTAGAGAAGGTCCACAGCGCAGAAGAGAACCAg                                                                                                                         >  2:2705301/1‑150 (MQ=255)
                              cATGGTCATCATGATCATGCGGGCAATCCAGAAGAAAATGATGTCGAAACCAGATACCATCACGCTGGTTGGGTGGAACTGACGCAGGGCGTCGGTATTTTCCGGCCAGCCAAGGGTAGAGAAGGTCCACAGCGCAGAAGAGAACCAGGt                                                                                                                       >  2:2397559/1‑150 (MQ=255)
                                   tcatcaTGATCATGCGGGCAATCCAGAAGAAAATGATGTCGAAACCAGATACCATCACGCTGGTTGGGTGGAACTGACGCAGGGCGTCGGTATTTTCCGGCCAGCCAAGGGTAGAGAAGGTCCACAGCGCAGAAGAGAACCAGGTATCga                                                                                                                  >  2:1526511/1‑150 (MQ=255)
                                     atcatgatcatgCGGGCAATCCAGAAGAAAATGATGTCGAAACCAGATACCATCACGCTGGTTGGGTGGAACTGACGCAGGGCGTCGGTATTTTCCGGCCAGCCAAGGGTAGAGAAGGTCCACAGCGCAGAAGAGAACCAGGTATCgaga                                                                                                                >  2:782561/1‑150 (MQ=255)
                                     atcatgatcatgCGGGCAATCCAGAAGAAAATGATGTCGAAACCAGATACCATCACGCTGGTTGGGTGGAACTGACGCAGGGCGTCGGTATTTTCCGGCCAGCCAAGGGTAGAGAAGGTCCACAGCGCAGAAGAGAACCAGGTATCgaga                                                                                                                >  2:1993230/1‑150 (MQ=255)
                                      tcatgatcatgCGGGCAATCCAGAAGAAAATGATGTCGAAACCAGATACCATCACGCTGGTTGGGTGGAACTGACGCAGGGCGTCGGTATTTTCCGGCCAGCCAAGGGTAGAGAAGGTCCACAGCGCAGAAGAGAACCAGGTATCGAGaa                                                                                                               >  2:440820/1‑150 (MQ=255)
                                       catgatcatgCGGGCAATCCAGAAGAAAATGATGTCGAAACCAGATACCATCACGCTGGTTGGGTGGAACTGACGCAGGGCGTCGGTATTTTCCGGCCAGCCAAGGGTAGAGAAGGTCCACAGCGCAGAAGAGAACCAGGTATCGAGAAc                                                                                                              <  2:3257046/150‑1 (MQ=255)
                                           atcatgCGGGCAATCCAGAAGAAAATGATGTCGAAACCAGATACCATCACGCTGGTTGGGTGGAACTGACGCAGGGCGTCGGTATTTTCCGGCCAGCCAAGGGTTGAGAAGGTCCACAGCGCAGAAGAGAACCAGGTATCGAGAACGTCt                                                                                                          <  2:666391/150‑1 (MQ=255)
                                           atcatgCGGGCAATCCAGAAGAAAATGATGTCGAAACCAGATACCATCACGCTGGTTGGGTGGAACTGACGCAGGGCGTCGGTATTTTCCGGCCAGCCAAGGGTAGAGAAGGTCCACAGCGCAGAAGAGAACCAGGTATCGAGAACGTCt                                                                                                          <  1:1706438/150‑1 (MQ=255)
                                             catgCGGGCAATCCAGAAGAAAATGATGTCGAAACCAGATACCATCACGCTGGTTGGGTGGAACTGACGCAGGGCGTCGGTATTTTCCGGCCAGCCAAGGGTAGAGAAGGTCCACAGCGCAGAAGAGAACCAGGTATCGAGAACGTCTTc                                                                                                        <  1:808182/150‑1 (MQ=255)
                                             catgCGGGCAATCCAGAAGAAAATGATGTCGAAACCAGATACCATCACGCTGGTTGGGTGGAACTGACGCAGGGCGTCGGTATTTTCCGGCCAGCCAAGGGTAGAGAAGGTCCACAGCGCAGAAGAGAACCAGGTATCGAGAACGTCTTc                                                                                                        <  2:1424385/150‑1 (MQ=255)
                                              atgCGGGCAATCCAGAAGAAAATGATGTCGAAACCAGATACCATCACGCTGGTTGGGTGGAACTGACGCAGGGCGTCGGTATTTTCCGGCCAGCCAAGGGTAGAGAAGGTCCACAGCGCAGAAGAGAACCAGGTATCGAGAACGTCTTCg                                                                                                       >  2:3470951/1‑150 (MQ=255)
                                              atgCGGGCAATCCAGAAGAAAATGATGTCGAAACCAGATACCATCACGCTGGTTGGGTGGAACTGACGCAGGGCGTCGGTATTTTCCGGCCAGCCAAGGGTAGAGAAGGTCCACAGCGCAGAAGAGAACCAGGTATCGAGAACGTCTTCg                                                                                                       >  2:53270/1‑150 (MQ=255)
                                                gCGGGCAATCCAGAAGAAAATGATGTCGAAACCAGATACCATCACGCTGGTTGGGTGGAACTGACGCAGGGCGTCGGTATTTTCCGGCCAGCCAAGGGTAGAGAAGGTCCACAGCGCAGAAGAGAACCAGGTATCGAGAACGTCTTCGTc                                                                                                     <  2:3601681/150‑1 (MQ=255)
                                                gCGGGCAATCCAGAAGAAAATGATGTCGAAACCAGATACCATCACGCTGGTTGGGTGGAACTGACGCAGGGCGTCGGTATTTTCCGGCCAGCCAAGGGTAGAGAAGGTCCACAGCGCAGAAGAGAACCAGGTATCGAGAACGTCTTCGTc                                                                                                     <  2:3165944/150‑1 (MQ=255)
                                                 cGGGCAATCCAGAAGAAAATGATGTCGAAACCAGATACCATCACGCTGGTTGGGTGGAACTGACGCAGGGCGTCGGTATTTTCCGGCCAGCCAAGGGTAGAGAAGGTCCACAGCGCAGAAGAGAACCAGGTATCGAGAACGTCTTCGTcc                                                                                                    >  2:2550642/1‑150 (MQ=255)
                                                      aaTCCAGAAGAAAATGATGTCGAAACCAGATACCATCACGCTGGTTGGGTGGAACTGACGCAGGGCGTCGGTATTTTCCGGCCAGCCAAGGGTAGAGAAGGTCCACAGCGCAGAAGAGAACCAGGTATCGAGAACGTCTTCGTCCTGACg                                                                                               >  2:163103/1‑150 (MQ=255)
                                                      aaTCCAGAAGAAAATGATGTCGAAACCAGATACCATCACGCTGGTTGGGTGGAACTGACGCAGGGCGTCGGTATTTTCCGGCCAGCCAAGGGTAGAGAAGGTCCACAGCGCAGAAGAGAACCAGGTATCGAGAACGTCTTCGTCCTGACg                                                                                               <  1:4368271/150‑1 (MQ=255)
                                                      aaTCCAGAAGAAAATGATGTCGAAACCAGATACCATCACGCTGGTTGGGTGGAACTGACGCAGGGCGTCGGTATTTTCCGGCCAGCCAAGGGTAGAGAAGGTCCACAGCGCAGAAGAGAACCAGGTATCGAGAACGTCTTCGTCCTGACg                                                                                               >  2:4290045/1‑150 (MQ=255)
                                                        tCCAGAAGAAAATGATGTCGAAACCAGATACCATCACGCTGGTTGGGTGGAACTGACGCAGGGCGTCGGTATTTTCCGGCCAGCCAAGGGTAGAGAAGGTCCACAGCGCAGAAGAGAACCAGGTATCGAGAACGTCTTCGTCCTGACGCa                                                                                             >  2:2052125/1‑150 (MQ=255)
                                                        tCCAGAAGAAAATGATGTCGAAACCAGATACCATCACGCTGGTTGGGTGGAACTGACGCAGGGCGTCGGTATTTTCCGGCCAGCCAAGGGTAGAGAAGGTCCACAGCGCAGAAGAGAACCAGGTATCGAGAACGTCTTCGTCCTGACGCa                                                                                             <  2:3543565/150‑1 (MQ=255)
                                                         ccAGAAGAAAATGATGTCGAAACCAGATACCATCACGCTGGTTGGGTGGAACTGACGCAGGGCGTCGGTATTTTCCGGCCAGCCAAGGGTAGAGAAGGTCCACAGCGCAGAAGAGAACCAGGTATCGAGAACGTCTTCGTCCTGACGCAg                                                                                            >  1:1241526/1‑150 (MQ=255)
                                                         ccAGAAGAAAATGATGTCGAAACCAGATACCATCACGCTGGTTGGGTGGAACTGACGCAGGGCGTCGGTATTTTCCGGCCAGCCAAGGGTAGAGAAGGTCCACAGCGCAGAAGAGAACCAGGTATCGAGAACGTCTTCGTCCTGACGCAg                                                                                            <  1:1092602/150‑1 (MQ=255)
                                                          cAGAAGAAAATGATGTCGAAACCAGATACCATCACGCTGGTTGGGTGGAACTGACGCAGGGCGTCGGTATTTTCCGGCCAGCCAAGGGTAGAGAAGGTCCACAGCGCAGAAGAGAACCAGGTATCGAGAACGTCTTCGTCCTGACGCAgg                                                                                           <  2:1974388/150‑1 (MQ=255)
                                                          cAGAAGAAAATGATGTCGAAACCAGATACCATCACGCTGGTTGGGTGGAACTGACGCAGGGCGTCGGTATTTTCCGGCCAGCCAAGGGTAGAGAAGGTCCACAGCGCAGAAGAGAACCAGGTATCGAGAACGTCTTCGTCCTGACGCAgg                                                                                           >  2:1316311/1‑150 (MQ=255)
                                                           agaagaAAATGATGTCGAAACCAGATACCATCACGCTGGTTGGGTGGAACTGACGCAGGGCGTCGGTATTTTCCGGCCAGCCAAGGGTAGAGAAGGTCCACAGCGCAGAAGAGAACCAGGTATCGAGAACGTCTTCGTCCTGACGCAGGa                                                                                          >  2:437880/1‑150 (MQ=255)
                                                             aagaaAATGATGTCGAAACCAGATACCATCACGCTGGTTGGGTGGAACTGACGCAGGGCGTCGGTATTTTCCGGCCAGCCAAGGGTAGAGAAGGTCCACAGCGCAGAAGAGAACCAGGTATCGAGAACGTCTTCGTCCTGACGCAGGaca                                                                                        <  1:716069/150‑1 (MQ=255)
                                                             aagaaAATGATGTCGAAACCAGATACCATCACGCTGGTTGGGTGGAACTGACGCAGGGCGTCGGTATTTTCCGGCCAGCCAAGGGTAGAGAAGGTCCACAGCGCAGAAGAGAACCAGGTATCGAGAACGTCTTCGTCCTGACGCAGGaca                                                                                        >  1:2717932/1‑150 (MQ=255)
                                                              agaaAATGATGTCGAAACCAGATACCATCACGCTGGTTGGGTGGAACTGACGCAGGGCGTCGGTATTTTCCGGCCAGCCAAGGGTAGAGAAGGTCCACAGCGCAGAAGAGAACCAGGTATCGAGAACGTCTTCGTCCTGACGCAGGacaa                                                                                       >  1:4269707/1‑150 (MQ=255)
                                                                aaaaTGATGTCGAAACCAGATACCATCACGCTGGTTGGGTGGAACTGACGCAGGGCGTCGGTATTTTCCGGCCAGCCAAGGGTAGAGAAGGTCCACAGCGCAGAAGAGAACCAGGTATCGAGAACGTCTTCGTCCTGACGCAGGacaaca                                                                                     <  2:541447/150‑1 (MQ=255)
                                                                  aaTGATGTCGAAACCAGATACCATCACGCTGGTTGGGTGGAACTGACGCAGGGCGTCGGTATTTTCCGGCCAGCCAAGGGTAGAGAAGGTCCACAGCGCAGAAGAGAACCAGGTATCGAGAACGTCTTCGTCCTGACGCAGGACAACATc                                                                                   <  2:1005287/150‑1 (MQ=255)
                                                                   atgatgTCGAAACCAGATACCATCTCGCTGGTTGGGTGGAACTGACGCAGGGCGTCGGTATTTTCCGGCCAGCCAAGGGTAGAGAAGGTCCACAGCGCAGAAGAGAACCAGGTATCGAGAACGTCTTCGTCCTGACGCAGGACAACATCa                                                                                  <  1:4232123/150‑1 (MQ=255)
                                                                     gatgTCGAAACCAGATACCATCACGCTGGTTGGGTGGAACTGACGCAGGGCGTCGGTATTTTCCGGCCAGCCAAGGGTAGAGAAGGTCCACAGCGCAGAAGAGAACCAGGTATCGAGAACGTCTTCGTCCTGACGCAGGACAACATCAGc                                                                                <  1:3869884/150‑1 (MQ=255)
                                                                        gTCGAAACCAGATACCATCACGCTGGTTGGGTGGAACTGACGCAGGGCGTCGGTATTTTCCGGCCAGCCAAGGGTAGAGAAGGTCCACAGCGCAGAAGAGAACCAGGTATCGAGAACGTCTTCGTCCTGACGCAGGACAACATCAGCAcc                                                                             <  1:3041847/150‑1 (MQ=255)
                                                                         tCGAAACCAGATACCATCACGCTGGTTGGGTGGAACTGACGCAGGGCGTCGGTATTTTCCGGCCAGCCAAGGGTTGAGAAGGTCCACAGCGCAGAAGAGAACCAGGTATCGAGAACGTCTTCGTCCTGACGCAGGACAACATCAGCACCg                                                                            <  2:3331105/150‑1 (MQ=255)
                                                                         tCGAAACCAGATACCATCACGCTGGTTGGGTGGAACTGACGCAGGGCGTCGGTATTTTCCGGCCAGCCAAGGGTAGAGAAGGTCCACAGCGCAGAAGAGAACCAGGTATCGAGAACGTCTTCGTCCTGACGCAGGACAACATCAGCACCg                                                                            <  1:874221/150‑1 (MQ=255)
                                                                         tCGAAACCAGATACCATCACGCTGGTTGGGTGGAACTGACGCAGGGCGTCGGTATTTTCCGGCCAGCCAAGGGTAGAGAAGGTCCACAGCGCAGAAGAGAACCAGGTATCGAGAACGTCTTCGTCCTGACGCAGGACAACATCAGCACCg                                                                            >  1:2268451/1‑150 (MQ=255)
                                                                          cGAAACCAGATACCATCACGCTGGTTGGGTGGAACTGACGCAGGGCGTCGGTATTTTCCGGCCAGCCAAGGGTAGAGAAGGTCCACAGCGCAGAAGAGAACCAGGTATCGAGAACGTCTTCGTCCTGACGCAGGACAACATCAGCACCGa                                                                           >  1:2504386/1‑150 (MQ=255)
                                                                           gAAACCAGATACCATCACGCTGGTTGGGTGGAACTGACGCAGGGCGTCGGTATTTTCCGGCCAGCCAAGGGTTGAGAAGGTCCACAGCGCAGAAGAGAACCAGGTATCGAGAACGTCTTCGTCCTGACGCAGGACAACATCAGCACCGAg                                                                          <  2:323542/150‑1 (MQ=255)
                                                                           gAAACCAGATACCATCACGCTGGTTGGGTGGAACTGACGCAGGGCGTCGGTATTTTCCGGCCAGCCAAGGGTAGAGAAGGTCCACAGCGCAGAAGAGAACCAGGTATCGAGAACGTCTTCGTCCTGACGCAGGACAACATCAGCACCGAg                                                                          >  1:803216/1‑150 (MQ=255)
                                                                              aCCAGATACCATCACGCTGGTTGGGTGGAACTGACGCAGGGCGTCGGTATTTTCCGGCCAGCCAAGGGTAGAGAAGGTCCACAGCGCAGAAGAGAACCAGGTATCGAGAACGTCTTCGTCCTGACGCAGGACAACATCAGCACCGAGGtt                                                                       <  1:254581/150‑1 (MQ=255)
                                                                              aCCAGATACCATCACGCTGGTTGGGTGGAACTGACGCAGGGCGTCGGTATTTTCCGGCCAGCCAAGGGTAGAGAAGGTCCACAGCGCAGAAGAGAACCAGGTATCGAGAACGTCTTCGTCCTGACGCAGGACAACATCAGCACCGAGGtt                                                                       <  1:253532/150‑1 (MQ=255)
                                                                                 aGATACCATCACGCTGGTTGGGTGGAACTGACGCAGGGCGTCGGTATTTTCCGGCCAGCCAAGGGTAGAGAAGGTCCACAGCGCAGAAGAGAACCAGGTATCGAGAACGTCTTCGTCCTGACGCAGGACAACATCAGCACCGAGGTTAtt                                                                    >  1:3500383/1‑150 (MQ=255)
                                                                                   aTACCATCACGCTGGTTGGGTGGAACTGACGCAGGGCGTCGGTATTTTCCGGCCAGCCAAGGGTAGAGAAGGTCCACAGCGCAGAAGAGAACCAGGTATCGAGAACGTCTTCGTCCTGACGCAGGACAACATCAGCACCGAGGTTAtttt                                                                  <  1:4190912/150‑1 (MQ=255)
                                                                                     aCCATCACGCTGGTTGGGTGGAACTGACGCAGGGCGTCGGTATTTTCCGGCCAGCCAAGGGTAGAGAAGGTCCACAGCGCAGAAGAGAACCAGGTATCGAGAACGTCTTCGTCCTGACGCAGGACAACATCAGCACCGAGGTTATTTTCt                                                                <  1:2221973/150‑1 (MQ=255)
                                                                                     aCCATCACGCTGGTTGGGTGGAACTGACGCAGGGCGTCGGTATTTTCCGGCCAGCCAAGGGTAGAGAAGGTCCACAGCGCAGAAGAGAACCAGGTATCGAGAACGTCTTCGTCCTGACGCAGGACAACATCAGCACCGAGGTTATTTTCt                                                                >  1:4151317/1‑150 (MQ=255)
                                                                                         tCACGCTGGTTGGGTGGAACTGACGCAGGGCGTCGGTATTTTCCGGCCAGCCAAGGGTAGAGAAGGTCCACAGCGCAGAAGAGAACCAGGTATCGAGAACGTCTTCGTCCTGACGCAGGACAACATCAGCACCGAGGTTATTTTCTTTAc                                                            <  1:4521244/150‑1 (MQ=255)
                                                                                         tCACGCTGGTTGGGTGGAACTGACGCAGGGCGTCGGTATTTTCCGGCCAGCCAAGGGTAGAGAAGGTCCACAGCGCAGAAGAGAACCAGGTATCGAGAACGTCTTCGTCCTGACGCAGGACAACATCAGCACCGAGGTTATTTTCTTTAc                                                            <  2:363529/150‑1 (MQ=255)
                                                                                         tCACGCTGGTTGGGTGGAACTGACGCAGGGCGTCGGTATTTTCCGGCCAGCCAAGGGTAGAGAAGGTCCACAGCGCAGAAGAGAACCAGGTATCGAGAACGTCTTCGTCCTGACGCAGGACAACATCAGCACCGAGGTTATTTTCTTTAc                                                            <  2:2626798/150‑1 (MQ=255)
                                                                                            cGCTGGTTGGGTGGAACTGACGCAGGGCGTCGGTATTTTCCGGCCAGCCAAGGGTAGAGAAGGTCCACAGCGCAGAAGAGAACCAGGTATCGAGAACGTCTTCGTCCTGACGCAGGACAACATCAGCACCGAGGTGATTTTCTTTACGCa                                                         <  1:482551/150‑1 (MQ=255)
                                                                                               tGGTTGGGTGGAACTGACGCAGGGCGTCGGTATTTTCCGGCCAGCCAAGGGTAGAGAAGGTCCACAGCGCAGAAGAGAACCAGGTATCGAGAACGTCTTCGTCCTGACGCAGGACAACATCAGCACCGAGGTTATTTTCTTTACGCActt                                                      >  1:1136125/1‑150 (MQ=255)
                                                                                               tGGTTGGGTGGAACTGACGCAGGGCGTCGGTATTTTCCGGCCAGCCAAGGGTAGAGAAGGTCCACAGCGCAGAAGAGAACCAGGTATCGAGAACGTCTTCGTCCTGACGCAGGACAACATCAGCACCGAGGTTATTTTCTTTACGCActt                                                      <  2:2757554/150‑1 (MQ=255)
                                                                                                 gTTGGGTGGAACTGACGCAGGGCGTCGGTATTTTCCGGCCAGCCAAGGGTAGAGAAGGTCCACAGCGCAGAAGAGAACCAGGTATCGAGAACGTCTTCGTCCTGACGCAGGACAACATCAGCACCGAGGTTATTTTCTTTACGCActtcg                                                    <  2:4156857/150‑1 (MQ=255)
                                                                                                 gTTGGGTGGAACTGACGCAGGGCGTCGGTATTTTCCGGCCAGCCAAGGGTAGAGAAGGTCCACAGCGCAGAAGAGAACCAGGTATCGAGAACGTCTTCGTCCTGACGCAGGACAACATCAGCACCGAGGTTATTTTCTTTACGCActtcg                                                    >  2:1701308/1‑150 (MQ=255)
                                                                                                 gTTGGGTGGAACTGACGCAGGGCGTCGGTATTTTCCGGCCAGCCAAGGGTAGAGAAGGTCCACAGCGCAGAAGAGAACCAGGTATCGAGAACGTCTTCGTCCTGACGCAGGACAACATCAGCACCGAGGTTATTTTCTTTACGCActtcg                                                    <  2:1258623/150‑1 (MQ=255)
                                                                                                   tGGGTGGAACTGACGCAGGGCGTCGGTATTTTCCGGCCAGCCAAGGGTAGAGAAGGTCCACAGCGCAGAAGAGAACCAGGTATCGAGAACGTCTTCGTCCTGACGCAGGACAACATCAGCACCGAGGTTATTTTCTTTACGCActtcgtc                                                  <  2:337597/150‑1 (MQ=255)
                                                                                                    gggTGGAACTGACGCAGGGCGTCGGTATTTTCCGGCCAGCCAAGGGTAGAGAAGGTCCACAGCGCAGAAGAGAACCAGGTATCGAGAACGTCTTCGTCCTGACGCAGGACAACATCAGCACCGAGGTTATTTTCTTTACGCActtcgtct                                                 <  2:841421/150‑1 (MQ=255)
                                                                                                    gggTGGAACTGACGCAGGGCGTCGGTATTTTCCGGCCAGCCAAGGGTAGAGAAGGTCCACAGCGCAGAAGAGAACCAGGTATCGAGAACGTCTTCGTCCTGACGCAGGACAACATCAGCACCGAGGTTATTTTCTTTACGCActtcgtct                                                 <  2:2228643/150‑1 (MQ=255)
                                                                                                       tGGAACTGACGCAGGGCGTCGGTATTTTCCGGCCAGCCAAGGGTTGAGAAGGTCCACAGCGCAGAAGAGAACCAGGTATCGAGAACGTCTTCGTCCTGACGCAGGACAACATCAGCACCGAGGTTATTTTCTTTACGCActtcgtcttcg                                              <  2:2097354/150‑1 (MQ=255)
                                                                                                       tGGAACTGACGCAGGGCGTCGGTATTTTCCGGCCAGCCAAGGGTAGAGAAGGTCCACAGCGCAGAAGAGAACCAGGTATCGAGAACGTCTTCGTCCTGACGCAGGACAACATCAGCACCGAGGTTATTTTCTTTACGCActtcgtcttcg                                              <  1:650112/150‑1 (MQ=255)
                                                                                                                cGCAGGGCGTCGGTATTTTCCGGCCAGCCAAGGGTAGAGAAGGTCCACAGCGCAGAAGAGAACCAGGTATCGAGAACGTCTTCGTCCTGACGCAGGACAACATCAGCACCGAGGTTATTTTCTTTACGCACTTCGTCTTCGTTGCGGCCa                                     >  1:1021900/1‑150 (MQ=255)
                                                                                                                  cAGGGCGTCGGTATTTTCCGGCCAGCCAAGGGTTGAGAAGGTCCACAGCGCAGAAGAGAACCAGGTATCGAGAACGTCTTCGTCCTGACGCAGGACAACATCAGCACCGAGGTTATTTTCTTTACGCACTTCGTCTTCGTTGCGGCCAAc                                   <  1:2132398/150‑1 (MQ=255)
                                                                                                                   aGGGCGTCGGTATTTTCCGGCCAGCCAAGGGTAGAGAAGGTCCACAGCGCAGAAGAGAACCAGGTATCGAGAACGTCTTCGTCCTGACGCAGGACAACATCAGCACCGAGGTTATTTTCTTTACGCACTTCGTCTTCGTTGCGGCCAACa                                  >  1:4447125/1‑150 (MQ=255)
                                                                                                                       cGTCGGTATTTTCCGGCCAGCCAAGGGTAGAGAAGGTCCACAGCGCAGAAGAGAACCAGGTATCGAGAACGTCTTCGTCCTGACGCAGGACAACATCAGCACCGAGGTTATTTTCTTTACGCACTTCGTCTTCGTTGCGGCCAACATaaa                              >  1:1280558/1‑150 (MQ=255)
                                                                                                                         tCGGTATTTTCCGGCCAGCCAAGGGTAGAGAAGGTCCACAGCGCAGAAGAGAACCAGGTATCGAGAACGTCTTCGTCCTGACGCAGGACACCATCAGCACCGAGGTTATTTTCTTTACGCACTTCGTCTTCGTTGCGGCCAACATAAACg                            >  1:154788/1‑150 (MQ=255)
                                                                                                                          cGGTATTTTCCGGCCAGCCAAGGGTAGAGAAGGTCCACAGCGCAGAAGAGAACCAGGTCTCGAGAACGTCTTCGTCCTGACGCAGGACAACATCAGCACCGAGGTTATTTTCTTTACGCACTTCGTCTTCGTTGCGGCCAACATAAACGt                           >  2:4041244/1‑150 (MQ=255)
                                                                                                                          cGGTATTTTCCGGCCAGCCAAGGGTAGAGAAGGTCCACAGCGCAGAAGAGAACCAGGTATCGAGAACGTCTTCGTCCTGACGCAGGACAACATCAGCACCGAGGTTATTTTCTTTACGCACTTCGTCTTCGTTGCGGCCAACATAAACGt                           >  2:1347577/1‑150 (MQ=255)
                                                                                                                           ggTATTTTCCGGCCAGCCAAGGGTAGAGAAGGTCCACAGCGCAGAAGAGAACCAGGTATCGAGAACGTCTTCGTCCTGACGCAGGACAACATCAGCACCGAGGTTATTTTCTTTACGCACTTCGTCTTCGTTGCGGCCAACATCAACGtt                          >  1:2332216/1‑150 (MQ=255)
                                                                                                                            gTATTTTCCGGCCAGCCAAGGGTTGAGAAGGTCCACAGCGCAGAAGAGAACCAGGTATCGAGAACGTCTTCGTCCTGACGCAGGACAACATCAGCACCGAGGTTATTTTCTTTACGCACTTCGTCTTCGTTGCGGCCAACATAAACGTTa                         <  2:1113002/150‑1 (MQ=255)
                                                                                                                            gTATTTTCCGGCCAGCCAAGGGTAGAGAAGGTCCACAGCGCAGAAGAGAACCAGGTATCGAGAACGTCTTCGTCCTGACGCAGGACAACATCAGCACCGAGGTTATTTTCTTTACGCACTTCGTCTTCGTTGCGGCCAACATAAACGTTa                         >  2:2426622/1‑150 (MQ=255)
                                                                                                                            gTATTTTCCGGCCAGCCAAGGGTAGAGAAGGTCCACAGCGCAGAAGAGAACCAGGTATCGAGAACGTCTTCGTCCTGACGCAGGACAACATCAGCACCGAGGTTATTTTCTTTACGCACTTCGTCTTCGTTGCGGCCAACATAAACGTTa                         <  1:627398/150‑1 (MQ=255)
                                                                                                                            gTATTTTCCGGCCAGCCAAGGGTAGAGAAGGTCCACAGCGCAGAAGAGAACCAGGTATCGAGAACGTCTTCGTCCTGACGCAGGACAACATCAGCACCGAGGTTATTTTCTTTACGCACTTCGTCTTCGTTGCGGCCAACATAAACGTTa                         <  1:2891274/150‑1 (MQ=255)
                                                                                                                             tATTTTCCGGCCAGCCAAGGGTAGAGAAGGTCCACAGCGCAGAAGAGAACCAGGTATCGAGAACGTCTTCGTCCTGACGCAGGACAACATCAGCACCGAGGTTATTTTCTTTACGCACTTCGTCTTCGTTGCGGCCAACATAAACGTTAc                        >  1:4134726/1‑150 (MQ=255)
                                                                                                                               ttttCCGGCCAGCCAAGGGTAGAGAAGGTCCACAGCGCAGAAGAGAACCAGGTATCGAGAACGTCTTCGTCCTGACGCAGGACAACATCAGCACCGAGGTTATTTTCTTTACGCACTTCGTCTTCGTTGCGGCCAACATAAACGTTAccc                      <  2:2878861/150‑1 (MQ=255)
                                                                                                                                  tCCGGCCAGCCAAGGGTTGAGAAGGTCCACAGCGCAGAAGAGAACCAGGTATCGAGAACGTCTTCGTCCTGACGCAGGACAACATCAGCACCGAGGTTATTTTCTTTACGCACTTCGTCTTCGTTGCGGCCAACATAAACGTTACCCGCt                   <  2:566105/150‑1 (MQ=255)
                                                                                                                                    cGGCCAGCCAAGGGTTGAGAAGGTCCACAGCGCAGAAGAGAACCAGGTATCGAGAACGTCTTCGTCCTGACGCAGGACAACATCAGCACCGAGGTTATTTTCTTTACGCACTTCGTCTTCGTTGCGGCCAACATAAACGTTACCCGCTTc                 >  1:3103859/1‑150 (MQ=255)
                                                                                                                                    cGGCCAGCCAAGGGTAGAGAAGGTCCACAGCGCAGAAGAGAACCAGGTATCGAGAACGTCTTCGTCCTGACGCAGGACAACATCAGCACCGAGGTGATTTTCTTTACGCACTTCGTCTTCGTTGCGGCCAACATAAACGTTACCCGCTTc                 <  1:2750163/150‑1 (MQ=255)
                                                                                                                                       ccagccaAGGGTAGAGAAGGTCCACAGCGCAGAAGAGAACCAGGTATCGAGAACGTCTTCGTCCTGACGCAGGACAACATCAGCACCGAGGTTATTTTCTTTACGCACTTCGTCTTCGTTGCGGCCAACATAAACGTTACCCGCTTCGTc              >  1:3653257/1‑150 (MQ=255)
                                                                                                                                          gccaAGGGTAGAGAAGGTCCACAGCGCAGAAGAGAACCAGGTATCGAGAACGTCTTCGTCCTGACGCAGGACAACATCAGCACCGAGGTTATTTTCTTTACGCACTTCGTCTTCGTTGCGGCCAACATAAACGTTACCCGCTTCGTCATa           <  2:2991143/150‑1 (MQ=255)
                                                                                                                                           ccaAGGGTAGAGAAGGTCCACAGCGCAGAAGAGAACCAGGTATCGAGAACGTCTTCGTCCTGACGCAGGACAACATCAGCACCGAGGTTATTTTCTTTACGCACTTCGTCTTCGTTGCGGCCAACATAAACGTTACCCGCTTCGTCATAc          <  1:2630653/150‑1 (MQ=255)
                                                                                                                                           ccaAGGGTAGAGAAGGTCCACAGCGCAGAAGAGAACCAGGTATCGAGAACGTCTTCGTCCTGACGCAGGACAACATCAGCACCGAGGTTATTTTCTTTACGCACTTCGTCTTCGTTGCGGCCAACATAAACGTTACCCGCTTCGTCATAc          >  1:1451242/1‑150 (MQ=255)
                                                                                                                                                 gTAGAGAAGGTCCACAGCGCAGAAGAGAACCAGGTATCGAGAACGTCTTCGTCCTGACGCAGGACAACATCAGCACCGAGGTTATTTTCTTTACGCACTTCGTCTTCGTTGCGGCCAACATAAACGTTACCCGCTTCGTCATACCATGcc    <  2:3978503/150‑1 (MQ=255)
                                                                                                                                                 gTAGAGAAGGTCCACAGCGCAGAAGAGAACCAGGTATCGAGAACGTCTTCGTCCTGACGCAGGACAACATCAGCACCGAGGTTATTTTCTTTACGCACTTCGTCTTCGTTGCGGCCAACATAAACGTTACCCGCTTCGTCATACCATGcc    <  1:3397567/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
CGGTTTGCCATTTTCATCTTTGATGAAGTGCATGGTCATCATGATCATGCGGGCAATCCAGAAGAAAATGATGTCGAAACCAGATACCATCACGCTGGTTGGGTGGAACTGACGCAGGGCGTCGGTATTTTCCGGCCAGCCAAGGGTAGAGAAGGTCCACAGCGCAGAAGAGAACCAGGTATCGAGAACGTCTTCGTCCTGACGCAGGACAACATCAGCACCGAGGTTATTTTCTTTACGCACTTCGTCTTCGTTGCGGCCAACATAAACGTTACCCGCTTCGTCATACCATGCCGG  >  CP009273/4472053‑4472349

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: