Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,252,556 A→T 5.3% L320F (TTA→TTT dusA → tRNA‑dihydrouridine synthase A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,252,5560AT 5.3% 929.1 / 28.6 301L320F (TTA→TTTdusAtRNA‑dihydrouridine synthase A
Reads supporting (aligned to +/- strand):  ref base A (161/123);  new base T (7/9);  total (168/133)
Fisher's exact test for biased strand distribution p-value = 4.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.76e-01

Only 100 of 309 total aligned reads displayed.

ACCGATGCCGATCCGGTGGCGGTAGTGCGCGCCATGTATCCGTACATTGAGCGTGAACTCAGCCAGGGGACGTATCTCGGCCATATTACCCGGCATATGTTGGGCTTGTTCCAGGGTATTCCTGGCGCGCGGCAGTGGCGGCGTTATTTAAGTGAAAATGCCCATAAAGCGGGTGCAGACATTAATGTGCTGGAACACGCGCTCAAACTGGTGGCGGATAAGCGTTAACTTTTCACCAAAAAGTAGTCAAATTCACCACGCCCTGCGCACCGTCGCGGGGCGTTTTGCTGTTAAATCAA  >  CP009273/4252407‑4252705
                                                                                                                                                     |                                                                                                                                                     
aCCGATGCCGATCCGGTGGCGGTAGTGCGCGCCATGTATCCGTACATTGAGCGTGAACTCAGCCAGGGGACGTATCTCGGCCATATTACCCGGCATATGTTGGGCTTGTTCCAGGGTATTCCTGGCGCGCGGCAGTGGCGGCGTTATTTa                                                                                                                                                       <  1:5189159/150‑1 (MQ=255)
aCCGATGCCGATCCGGTGGCGGTAGTGCGCGCCATGTATCCGTACATTGAGCGTGAACTCAGCCAGGGGACGTATCTCGGCCATATTACCCGGCATATGTTGGGCTTGTTCCAGGGTATTCCTGGCGCGCGGCAGTGGCGGCGTTATTTa                                                                                                                                                       <  2:4056565/150‑1 (MQ=255)
   gATGCCGATCCGGTGGCGGTAGTGCGCGCCATGTATCCGTACATTGAGCGTGAACTCAGCCAGGGGACGTATCTCGGCCATATTACCCGGCATATGTTGGGCTTGTTCCAGGGTATTCCTGGCGCGCGGCAGTGGCGGCGTTATTTAAGt                                                                                                                                                    <  2:3734377/150‑1 (MQ=255)
       ccGATCCGGTGGCGGTAGTGCGCGCCATGTATCCGTACATTGAGCGTGAACTCAGCCAGGGGACGTATCTCGGCCATATTACCCGGCATATGTTGGGCTTGTCCCAGGGTATTCCTGGCGCGCGGCAGTGGCGGCGTTATTTAAGTGaaa                                                                                                                                                <  2:1956761/150‑1 (MQ=255)
        cGATCCGGTGGCGGTAGTGCGCGCCATGTATCCGTACATTGAGCGTGAACTCAGCCAGGGGACGTATCTCGGCCATATTACCCGGCATATGTTGGGCTTGTTCCAGGGTATTCCTGGCGCGCGGCAGTGGCGGCGTTATTTAAGTGaaaa                                                                                                                                               >  1:224450/1‑150 (MQ=255)
         gATCCGGTGGCGGTAGTGCGCGCCATGTATCCGTACATTGAGCGTGAACTCAGCCAGGGGACGTATCTCGGCCATATTACCCGGCATATGTTGGGCTTGTTCCAGGGTATTCCTGGCGCGCGGCAGTGGCGGCGTTATTTAAGTGAAAAt                                                                                                                                              >  1:4382033/1‑150 (MQ=255)
         gATCCGGTGGCGGTAGTGCGCGCCATGTATCCGTACATTGAGCGTGAACTCAGCCAGGGGACGTATCTCGGCCATATTACCCGGCATATGTTGGGCTTGTTCCAGGGTATTCCTGGCGCGCGGCAGTGGCGGCGTTATTTAAGTGAAAAt                                                                                                                                              <  2:4991327/150‑1 (MQ=255)
         gATCCGGTGGCGGTAGTGCGCGCCATGTATCCGTACATTGAGCGTGAACTCAGCCAGGGGACGTATCTCGGCCATATTACCCGGCATATGTTGGGCTTGTTCCAGGGTATTCCTGGCGCGCGGCAGTGGCGGCGTTATTTAAGTGAAAAt                                                                                                                                              >  2:4869722/1‑150 (MQ=255)
         gATCCGGTGGCGGTAGTGCGCGCCATGTATCCGTACATTGAGCGTGAACTCAGCCAGGGGACGTATCTCGGCCATATTACCCGGCATATGTTGGGCTTGTTCCAGGGTATTCCTGGCGCGCGGCAGTGGCGGCGTTATTTAAGTGAAAAt                                                                                                                                              >  2:4785460/1‑150 (MQ=255)
          aTCCGGTGGCGGTAGTGCGCGCCATGTATCCGTACATTGAGCGTGAACTCAGCCAGGGGACGTATCTCGGCCATATTACCCGGCATATGTTGGGCTTGTTCCAGGGTATTCCTGGCGCGCGGCAGTGGCGGCGTTATTTAAGTGAAAATg                                                                                                                                             >  2:298035/1‑150 (MQ=255)
               gTGGCGGTAGTGCGCGCCATGTATCCGTACATTGAGCGTGAACTCAGCCAGGGGACGTATCTCGGCCATATTACCCGGCATATGTTGGGCTTGTTCCAGGGTATTCCTGGCGCGCGGCAGTGGCGGCGTTATTTAAGTGAAAATGCCCAt                                                                                                                                        >  2:3112516/1‑150 (MQ=255)
                 ggCGGTAGTGCGCGCCATGTATCCGTACATTGAGCGTGAACTCAGCCAGGGGACGTATCTCGGCCATATTACCCGGCATATGTTGGGCTTGTTCCAGGGTATTCCTGGCGCGCGGCAGTGGCGGCGTTATTTAAGTGAAAATGCCCATaa                                                                                                                                      >  1:3087526/1‑150 (MQ=255)
                   cGGTAGTGCGCGCCATGTATCCGTACATTGAGCGTGAACTCAGCCAGGGGACGTATCTCGGCCATATTACCCGGCATATGTTGGGCTTGTTCCAGGGTATTCCTGGCGCGCGGCAGTGGCGGCGTTATTTAAGTGAAAATGCCCATAAAg                                                                                                                                    <  2:1530085/150‑1 (MQ=255)
                    ggTAGTGCGCGCCATGTATCCGTACATTGAGCGTGAACTCAGCCAGGGGACGTATCTCGGCCATATTACCCGGCATATGTTGGGCTTGTTCCAGGGTATTCCTGGCGCGCGGCAGTGGCGGCGTTATTTTAGTGAAAATGCCCATAAAGc                                                                                                                                   <  1:153600/150‑1 (MQ=255)
                       aGTGCGCGCCATGTATCCGTACATTGAGCGTGAACTCAGCCAGGGGACGTATCTCGGCCATATTACCCGGCATATGTTGGGCTTGTTCCAGGGTATTCCTGGCGCGCGGCAGTGGCGGCGTTATTTTAGTGAAAATGCCCATAAAGCggg                                                                                                                                >  1:2992270/1‑150 (MQ=255)
                       aGTGCGCGCCATGTATCCGTACATTGAGCGTGAACTCAGCCAGGGGACGTATCTCGGCCATATTACCCGGCATATGTTGGGCTTGTTCCAGGGTATTCCTGGCGCGCGGCAGTGGCGGCGTTATTTAAGTGAAAATGCCCATAAAGCggg                                                                                                                                >  1:2024910/1‑150 (MQ=255)
                       aGTGCGCGCCATGTATCCGTACATTGAGCGTGAACTCAGCCAGGGGACGTATCTCGGCCATATTACCCGGCATATGTTGGGCTTGTTCCAGGGTATTCCTGGCGCGCGGCAGTGGCGGCGTTATTTAAGTGAAAATGCCCATAAAGCggg                                                                                                                                >  1:4327963/1‑150 (MQ=255)
                        gTGCGCGCCATGTATCCGTACATTGAGCGTGAACTCAGCCAGGGGACGTATCTCGGCCATATTACCCGGCATATGTTGGGCTTGTTCCAGGGTATTCCTGGCGCGCGGCAGTGGCGGCGTTATTTAAGTGAAAATGCCCATAAAGCGGGt                                                                                                                               <  2:1021363/150‑1 (MQ=255)
                          gcgcgcCATGTATCCGTACATTGAGCGTGAACTCAGCCAGGGGACGTATCTCGGCCATATTACCCGGCATATGTTGGGCTTGTTCCAGGGTATTCCTGGCGCGCGGCAGTGGCGGCGTTATTTAAGTGAAAATGCCCATAAAGCGGGTGc                                                                                                                             <  1:2844328/150‑1 (MQ=255)
                          gcgcgcCATGTATCCGTACATTGAGCGTGAACTCAGCCAGGGGACGTATCTCGGCCATATTACCCGGCATATGTTGGGCTTGTTCCAGGGTATTCCTGGCGCGCGGCAGTGGCGGCGTTATTTAAGTGAAAATGCCCATAAAGCGGGTGc                                                                                                                             >  2:4918536/1‑150 (MQ=255)
                            gcgcCATGTATCCGTACATTGAGCGTGAACTCAGCCAGGGGACGTATCTCGGCCATATTACCCGGCATATGTTGGGCTTGTTCCAGGGTATTCCTGGCGCGCGGCAGTGGCGGCGTTATTTAAGTGAAAATGCCCATAAAGCGGGTGCAg                                                                                                                           <  2:4708204/150‑1 (MQ=255)
                             cgcCATGTATCCGTACATTGAGCGTGAACTCAGCCAGGGGACGTATCTCGGCCATATTACCCGGCATATGTTGGGCTTGTTCCAGGGTATTCCTGGCGCGCGGCAGTGGCGGCGTTATTTAAGTGAAAATGCCCATAAAGCGGGTGCAGa                                                                                                                          >  2:3240091/1‑150 (MQ=255)
                             cgcCATGTATCCGTACATTGAGCGTGAACTCAGCCAGGGGACGTATCTCGGCCATATTACCCGGCATATGTTGGGCTAGTTCCAGGGTAATCCTGTCGCGCGGCAGTGGCGTCGTTATTTAAGTTAAAATTCCCATAAATCGGGTTCAGa                                                                                                                          >  1:1390029/1‑150 (MQ=255)
                                  tGTATCCGTACATTGAGCGTGAACTCAGCCAGGGGACGTATCTCGGCCATATTACCCGGCATATGTTGGGCTTGTTCCAGGGTATTCCTGGCGCGCGGCAGTGGCGGCGTTATTTAAGTGAAAATGCCCATAAAGCGGGTGCAGACagtc                                                                                                                     >  2:1159703/1‑147 (MQ=255)
                                     aTCCGTACATTGAGCGTGAACTCAGCCAGGGGACGTATCTCGGCCATATTACCCGGCATATGTTGGGCTTGTTCCAGGGTATTCCTGGCGCGCGGCAGTGGCGGCGTTATTTTAGTGAAAATGCCCATAAAGCGGGTGCAGACATTAAtg                                                                                                                  <  1:2555464/150‑1 (MQ=255)
                                     aTCCGTACATTGAGCGTGAACTCAGCCAGGGGACGTATCTAGGCCATATTACCCGGCATATGTTGGGCTTGTTCCAGGGTATTCCTGGCGCGCGGCAGTGGCGGCGTTATTTAAGTGAAAATGCCCATAAAGCGGGTGCAGACATTAAtg                                                                                                                  >  2:2586997/1‑150 (MQ=255)
                                     aTCAGTACATTGAGAGTGAACTCAGCCAGGGGACGTATCTCGTCCATATTACCCGGCATATGTTTGTCTTGTTCCAGGGTATTCCTGGCGCGCGGCAGTGGCGGCGTTATTTAATTGAAAATGCCCATAAAGCGGGTGCAGACATTAAtg                                                                                                                  >  2:3378009/1‑150 (MQ=255)
                                      tCCGTACATTGAGCGTGAACTCAGCCAGGGGACGTATCTCGGCCATATTACCCGGCATATGTTGGGCTTGTTCCAGGGTATTCCTGGCGCGCGGCAGTGGCGGCGTTATTTTAGTGAAAATGCCCATAAAGCGGGTGCAGACATTAAtgt                                                                                                                 <  1:1038719/150‑1 (MQ=255)
                                        cGTACATTGAGCGTGAACTCAGCCAGGGGACGTATCTCGGCCATATTACCCGGCATATGTTGGGCTTGTTCCAGGGTATTCCTGGCGCGCGGCAGTGGCGGCGTTATTTAAGTGAAAATGCCCATAAAGCGGGTGCAGACATTAATGTGc                                                                                                               >  2:772537/1‑150 (MQ=255)
                                         gTACATTGAGCGTGAACTCAGCCAGGGGACGTATCTCGGCCATATTACCCGGCATATGTTGGGCTTGTTCCAGGGTATTCCTGGCGCGCGGCAGTGGCGGCGTTATTTAAGTGAAAATGCCCATAAAGCGGGTGCAGACATTAATGTGCt                                                                                                              >  2:3169069/1‑150 (MQ=255)
                                         gTACATTGAGCGTGAACTCAGCCAGGGGACGTATCTCGGCCATATTACCCGGCATATGTTGGGCTTGTTCCAGGGTATTCCTGGCGCGCGGCAGTGGCGGCGTTATTTAAGTGAAAATGCCCATAAAGCGGGTGCAGACATTAATGTGCt                                                                                                              >  2:4390985/1‑150 (MQ=255)
                                              ttGAGCGTGAACTCAGCCAGGGGACGTATCTCGGCCATATTACCCGGCATATGTTGGGCTTGTTCCAGGGTATTCCTGGCGCGCGGCAGTGGCGGCGTTATTTAAGTGAAAATGCCCATAAAGCGGGTGCAGACATTAATGTGCTGGAac                                                                                                         <  1:2172763/150‑1 (MQ=255)
                                                    gTGAACTCAGCCAGGGGACGTATCTCGGCCATATTACCCGGCATATGTTGGGCTTGTTCCAGGGTATTCCTGGCGCGCGGCAGTGGCGGCGTTATTTAAGTGAAAATGCCCATAAAGCGGGTGCAGACATTAATGTGCTGGAACACgcgc                                                                                                   >  2:801209/1‑150 (MQ=255)
                                                      gAACTCAGCCAGGGGACGTATCTCGGCCATATTACCCGGCATATGTTGGGCTTGTTCCAGGGTATTCCTGGCGCGCGGCAGTGGCGGCGTTATTTAAGTGAAAATGCCCATAAAGCGGGTGCAGACATTAATGTGCTGGAACACGCGCTc                                                                                                 >  1:2400405/1‑150 (MQ=255)
                                                      gAACTCAGCCAGGGGACGTATCTCGGCCATATTACCCGGCATATGTTGGGCTTGTTCCAGGGTATTCCTGGCGCGCGGCAGTGGCGGCGTTATTTAAGTGAAAATGCCCATAAAGCGGGTGCAGACATTAATGTGCTGGAACACGCGCTc                                                                                                 >  2:5197732/1‑150 (MQ=255)
                                                      gAACTCAGCCAGGGGACGTATCTCGGCCATATTACCCGGCATATGTTGGGCTTGTTCCAGGGTATTCCTGGCGCGCGGCAGTGGCGGCGTTATTTAAGTGAAAATGCCCATAAAGCGGGTGCAGACATTAATGTGCTGGAACACGCGCTc                                                                                                 >  1:5238913/1‑150 (MQ=255)
                                                        aCTCAGCCAGGGGACGTATCTCGGCCATATTACCCGGCATATGTTGGGCTTGTTCCAGGGTATTCCTGGCGCGCGGCAGTGGCGGCGTTATTTAAGTGAAAATGCCCATAAAGCGGGTGCAGACATTAATGTGCTGGAACACGCGCTCaa                                                                                               >  2:3612298/1‑150 (MQ=255)
                                                          tCAGCCAGGGGACGTATCTCGGCCATATTACCCGGCATATGTTGGGCTTGTTCCAGGGTCTTCCAGGCGCGCGGCCGTGGCGGCGTTATTTAACTGAAAATGCCCATAAAGCGGGTGCAGCCCTTAATGTGCTCGAACACGCGCTCAAAc                                                                                             >  2:2000206/1‑150 (MQ=255)
                                                          tCAGCCAGGGGACGTATCTCGGCCATATTACCCGGCATATGTTGGGCTTGTTCCAGGGTATTCCTGGCGCGCGGCAGTGGCGGCGTTATTTAAGTGAAAATGCCCATAAAGCGGGTGCAGACAGTAATGTGCTGGAACACGCGCTCAAAc                                                                                             >  2:4502894/1‑150 (MQ=255)
                                                           cAGCCAGGGGACGTATCTCGGCCATATTACCCGGCATATGTTGGGCTTGTTCCAGGGTATTCCTGGCGCGCGGCAGTGGCGGCGTTATTTAAGTGAAAATGCCCATAAAGCGGGTGCAGACATTAATGTGCTGGAACACGCGCTCAAACt                                                                                            <  2:4554918/150‑1 (MQ=255)
                                                           cAGCCAGGGGACGTATCTCGGCCATATTACCCGGCATATGTTGGGCTTGTTCCAGGGTATTCCTGGCGCGCGGCAGTGGCGGCGTTATTTAAGTGAAAATGCCCATAAAGCGGGTGCAGACATTAATGTGCTGGAACACGCGCTCAAACt                                                                                            <  1:3633058/150‑1 (MQ=255)
                                                              ccAGGGGACGTATCTCGGCCATATTACCCGGCATATGTTGGGCTTGTTCCAGGGTATTCCTGGCGCGCGGCAGTGGCGGCGTTATTTAAGTGAAAATGCCCATAAAGCGGGTGCAGACATTAATGTGCTGGAACACGCGCTCAAACtggt                                                                                         <  2:4070572/150‑1 (MQ=255)
                                                              ccAGGGGACGTATCTCGGCCATATTACCCGGCATATGTTGGGCTTGTTCCAGGGTATTCCTGGCGCGCGGCAGTGGCGGCGTTATTTAAGTGAAAATGCCCATAAAGCGGGTGCAGACATTAATGTGCTGGAACACGCGCTCAAACtggt                                                                                         <  1:1089966/150‑1 (MQ=255)
                                                                aGGGGACGTATCTCGGCCCTATTACCCGGCATATGTTGGGCTTGTTCCAGGGTATTCCTGGCGCGCGGCAGTGGCGGCGTTATTTAAGTGAAAATGCCCATAAAGCGGGTGCAGACATTAATGTGCTGGAACACGCGCTCAAACtggtgg                                                                                       >  2:3502521/1‑150 (MQ=255)
                                                                aGGGGACGTATCTCGGCCATATTACCCGGCATATGTTGGGCTTGTTCCAGGGTATTCCTGGCGCGCGGCAGTGGCGGCGTTATTTAAGTGAAAATGCCCATAAAGCGGGTGCAGACATTAATGTGCTGGAACACGCGCTCAAACtggtgg                                                                                       >  2:2225424/1‑150 (MQ=255)
                                                                    gACGTATCTCGGCCATATTACCCGGCATATGTTGGGCTTGTTCCAGGGTATTCCTGGCGCGCGGCAGTGGCGGCGTTATTTAAGTGAAAATGCCCATAAAGCGGGTGCAGACATTAATGTGCTGGAACACGCGCTCAAACTGGTGGCGGa                                                                                   <  2:1591485/150‑1 (MQ=255)
                                                                       gTATCTCGGCCATATTACCCGGCATATGTTGGGCTTGTTCCAGGGTATTCCTGGCGCGCGGCAGTGGCGGCGTTATTTAAGTGAAAATGCCCATAAAGCGGGTGCAGACATTAATGTGCTGGAACACGCGCTCAAACTGGTGGCGGATaa                                                                                <  2:2569039/150‑1 (MQ=255)
                                                                       gTATCTCGGCCATATTACCCGGCATATGTTGGGCTTGTTCCAGGGTATTCCTGGCGCGCGGCAGTGGCGGCGTTATTTAAGTGAAAATGCCCATAAAGCGGGTGCAGACATTAATGTGCTGGAACACGCGCTCAAACTGGTGGCGGATaa                                                                                >  1:265813/1‑150 (MQ=255)
                                                                        tATCTCGGCCATATTACCCGGCATATGTTGGGCTTGTTCCAGGGTATTCCTGGCGCGCGGCAGTGGCGGCGTTATTTAAGTGAAAATGCCCATAAAGCGGGTGCAGACATTAATGTGCTGGAACACGCGCTCAAACTGGTGGCGGATAAg                                                                               >  2:1865266/1‑150 (MQ=255)
                                                                        tATCTCGGCCATATTACCCGGCATATGTTGGGCTTGGTCCAGGGTTTTCCTGGCGCGCGGCAGTGGCGGCGTTATTTAAGTGAAAATGCCCATAAAGCGGGTGCAGACATTAATGTGCTGGAACACGCGCTCAAACTGGTGGCGGATAAg                                                                               <  1:4504963/150‑1 (MQ=255)
                                                                             cGGCCATATTACCCGGCATATGTTGGGCTTGTTCCAGGGTATTCCTGGCGCGCGGCAGTGGCGGCGTTATTTAAGTGAAAATGCCCATAAAGCGGGTGCAGACATTAATGTGCTGGAACACGCGCTCAAACTGGTGGCGGATAAGCGTTa                                                                          <  1:2503135/150‑1 (MQ=255)
                                                                              ggCCATATTACCCGGCATATGTTGGGCTTGTTCCAGGGTATTCCTGGCGCGCGGCAGTGGCGGCGTTATTTAAGTGAAAATGCCCATAAAGCGGGTGCAGACATTAATGTGCTGGAACACGCGCTCAAACTGGTGGCGGATAAGCGTTaa                                                                         <  2:5007632/150‑1 (MQ=255)
                                                                                ccATATTACCCGGCATATGTTGGGCTTGTTCCAGGGTATTCCTGGCGCGCGGCAGTGGCGGCGTTATTTAAGTGAAAATGCCCATAAAGCGGGTGCAGACATTAATGTGCTGGAACACGCGCTCAAACTGGTGGCGGATAAGCGTTAACt                                                                       >  1:3744613/1‑150 (MQ=255)
                                                                                 cATATTACCCGGCATATGTTGGGCTTGTTCCAGGGTATTCCTGGCGCGCGGCAGTGGCGGCGTTATTTAAGTGAAAATGCCCATAAAGCGGGTGCAGACATTAATGTGCTGGAACACGCGCTCACACTGGTGGCGGATAAGCGTTAACtt                                                                      >  2:1270403/1‑150 (MQ=255)
                                                                                 cATATTACCCGGCATATGTTGGGCTTGTTCCAGGGTATTCCTGGCGCGCGGCAGTGGCGGCGTTATTTAAGTGAAAATGCCCATAAAGCGGGTGCAGACATTAATGTGCTGGAACACGCGCTCAAACTGGTGGCGGATAAGCGTTAACtt                                                                      >  2:1478872/1‑150 (MQ=255)
                                                                                   tatTACCCGGCATATGTTGGGCTTGTTCCAGGGTATTCCTGGCGCGCGGCAGTGGCGGCGTTATTTAAGTGAAAATGCCCATAAAGCGGGTGCAGACATTAATGTGCTGGAACACGCGCTCAAACTGGTGGCGGATAAGCGTTAACtttt                                                                    >  1:1712577/1‑150 (MQ=255)
                                                                                    atTACCCGGCATATGTTGGGCTTGGTCCAGGGTATTCCTGGCGCGCGGCAGTGGCGGCGTTATTTAAGTGAAAATGCCCATAAAGCGGGTGCAGACATTAATGTGCTGGAACACGCGCTCAAACTGGTGGCGGATAAGCGTTAACTTTTc                                                                   <  1:4733916/150‑1 (MQ=255)
                                                                                             cATATGTTGGGCTTGTTCCAGGGTATTCCTGGCGCGCGGCAGTGGCGGCGTTATTTAAGTGAAAATGCCCATAAAGCGGGTGCAGACATTAATGTGCTGGAACACGCGCTCAAACTGGTGGCGGATAAGCGTTAACTTTTCACCAAAAag                                                          <  2:3998056/150‑1 (MQ=255)
                                                                                             cATATGTTGGGCTTGTTCCAGGGTATTCCTGGCGCGCGGCAGTGGCGGCGTTATTTAAGTGAAAATGCCCATAAAGCGGGTGCAGACATTAATGTGCTGGAACACGCGCTCAAACTGGTGGCGGATAAGCGTTAACTTTTCACCAAAAag                                                          >  1:4824746/1‑150 (MQ=255)
                                                                                              atatGTTGGGCTTGTTCCAGGGTATTCCTGGCGCGCGGCAGTGGCGGCGTTATTTAAGTGAAAATGCCCATAAAGCGGGTGCAGACATTAATGTGCTGGAACACGCGCTCAAACTGGTGGCGGATAAGCGTTAACTTTTCACCAAAAagt                                                         >  1:1031819/1‑150 (MQ=255)
                                                                                               tatGTTGGGCTTGTTCCAGGGTATTCCTGGCGCGCGGCAGTGGCGGCGTTATTTAAGTGAAAATGCCCATAAAGCGGGTGCAGACATTAATGTGCTGGAACACGCGCTCAAACTGGTGGCGGATAAGCGTTAACTTTTCACCAAAAagta                                                        >  1:1555937/1‑150 (MQ=255)
                                                                                               tatGTTGGGCTTGTTCCAGGGTATTCCTGGCGCGCGGCAGTGGCGGCGTTATTTAAGTGAAAATGCCCATAAAGCGGGTGCAGACATTAATGTGCTGGAACACGCGCTCAAACTGGTGGCGGATAAGCGTTAACTTTTCACCAAAAagta                                                        <  2:2497051/150‑1 (MQ=255)
                                                                                                atGTTGGGCTTGTTCCAGGGTATTCCTGGCGCGCGGCAGTGGCGGCGTTATTTAAGTGAAAATGCCCATAAAGCGGGTGCAGACATTAATGTGCTGGAACACGCGCTCAAACTGGTGGCGGATAAGCGTTAACTTTTCACCAAAAagtag                                                       <  2:2706839/150‑1 (MQ=255)
                                                                                                   ttGGGCTTGTTCCAGGGTATTCCTGGCGCGCGGCAGTGGCGGCGTTATTTAAGTGAAAATGCCCATAAAGCGGGTGCAGACATTAATGTGCTGGAACACGCGCTCAAACTGGTGGCGGATAAGCGTTAACTTTTCACCAAAAAGTAGTCa                                                    <  1:3845840/150‑1 (MQ=255)
                                                                                                     gggCTTGTTCCAGGGTATTCCTGGCGCGCGGCAGTGGCGGCGTTATTTAAGTGAAAATGCCCATAAAGCGGGTGCAGACATTAATGTGCTGGAACACGCGCTCAAACTGGTGGCGGATAAGCGTTAACTTTTCACCAAAAAGTAGTCaaa                                                  <  2:2450334/150‑1 (MQ=255)
                                                                                                      ggCTTGTTCCAGGGTATTCCTGGCGCGCGGCAGTGGCGGCGTTATTTAAGTGAAAATGCCCATAAAGCGGGTGCAGACATTAATGTGCTGGAACACGCGCTCAAACTGGTGGCGGATAAGCGTTAACTTTTCACCAAAAAGTAGTCAAAt                                                 <  1:5042038/150‑1 (MQ=255)
                                                                                                       gCTTGTTCCAGGGTATTCCTGGCGCGCGGCAGTGGCGGCGTTATTTAAGTGAAAATGCCCATAAAGCGGGTGCAGACATTAATGTGCTGGAACACGCGCTCAAACTGGTGGCGGATAAGCGTTAACTTTTCACCAAAAAGTAGTCAAAtt                                                <  1:2484393/150‑1 (MQ=255)
                                                                                                        cTTGTTCCAGGGTATTCCTGGCGCGCGGCAGTGGCGGCGTTATTTAAGTGAAAATGCCCATAAAGCGGGTGCAGACATTAATGTGCTGGAACACGCGCTCAAACTGGTGGCGGATAAGCGTTAACTTTTCACCAAAAAGTAGTCAAATTc                                               >  1:1770014/1‑150 (MQ=255)
                                                                                                           gTTCCAGGGTATTCCTGGCGCGCGGCAGTGGCGGCGTTATTTAAGTGAAAATGCCCATAAAGCGGGTGCAGACATTAATGTGCTGGAACACGCGCTCAAACTGGTGGCGGATAAGCGTTAACTTTTCACCAAAAAGTAGTCAAATTcacc                                            <  2:2963845/150‑1 (MQ=255)
                                                                                                               cAGGGTATTCCTGGCGCGCGGCAGTGGCGGCGTTATTTAAGTGAAAATGCCCATAAAGCGGGTGCAGACATTAATGTGCTGGAACACGCGCTCAAACTGGTGGCGGATAAGCGTTGACTTTTCACCAAAAAGTAGTCAAATTCACCACGc                                        >  1:867337/1‑150 (MQ=255)
                                                                                                               cAGGGTATTCCTGGCGCGCGGCAGTGGCGGCGTTATTTAAGTGAAAATGCCCATAAAGCGGGTGCAGACATTAATGTGCTGGAACACGCGCTCAAACTGGTGGCGGATAAGCGTTAACTTTTCACCAAAAAGTAGTCAAATTCACCACGc                                        <  1:4869795/150‑1 (MQ=255)
                                                                                                                aGGGTATTCCTGGCGCGCGGCAGTGGCGGCGTTATTTAAGTGAAAATGCCCATAAAGCGGGTGCAGACATTAATGTGCTGGAACACGCGCTCAAACTGGTGGCGGATAAGCGTTAACTTTTCACCAAAAAGTAGTCAAATTCACCACGcc                                       >  1:5131141/1‑150 (MQ=255)
                                                                                                                 gggTATTCCTGGCGCGCGGCAGTGGCGGCGTTATTTAAGTGAAAATGCCCATAAAGCGGGTGCAGACATTAATGTGCTGGAACACGCGCTCAAACTGGTGGCGGATAAGCGTTAACTTTTCACCAAAAAGTAGTCAAATTCACCACGccc                                      >  1:4440326/1‑150 (MQ=255)
                                                                                                                  ggTATTCCTGGCGCGCGGCAGTGGCGGCGTTATTTAAGTGAAAATGCCCATAAAGCGGGTGCAGACATTAATGTGCTGGAACACGCGCTCAAACTGGTGGCGGATAAGCGTTAACTTTTCACCAAAAAGTAGTCAAATTCACCACGCCCt                                     <  1:1670546/150‑1 (MQ=255)
                                                                                                                   gTATTCCTGGCGCGCGGCAGTGGCGGCGTTATTTTAGTGAAAATGCCCATAAAGCGGGTGCAGACATTAATGTGCTGGAACACGCGCTCAAACTGGTGGCGGATAAGCGTTAACTTTTCACCAAAAAGTAGTCAAATTCACCACGCCCTg                                    <  1:526739/150‑1 (MQ=255)
                                                                                                                   gTATTCCTGGCGCGCGGCAGTGGCGGCGTTATTTAAGTGAAAATGCCCATAAAGCGGGTGCAGACATTAATGTGCTGGAACACGCGCTCAAACTGGTGGCGGATAAGCGTTAACTTTTCACCAAAAAGTAGTCAAATTCACCACGCCCTg                                    <  1:3391151/150‑1 (MQ=255)
                                                                                                                     aTTCCTGGCGCGCGGCAGTGGCGGCGTTATTTAAGTGAAAATGCCCATAAAGCGGGTGCAGACATTAATGTGCTGGAACACGCGCTCAAACTGGTGGCGGATAAGCGTTAACTTTTCACCAAAAAGTAGTCAAATTCACCACGCCCTgcg                                  >  2:2270108/1‑150 (MQ=255)
                                                                                                                       tCCTGGCGCGCGGCAGTGGCGGCGTTATTTAAGTGAAAATGCCCATAAAGCGGGTGCAGACATTAATGTGCTGGAACACGCGCTCAAACTGGTGGCGGATAAGCGTTAACTTTTCACCAAAAAGTAGTCAAATTCACCACGCCCTGCGCa                                <  2:1335143/150‑1 (MQ=255)
                                                                                                                        ccTGGCGCGCGGCAGTGGCGGCGTTATTTAAGTGAAAATGCCCATAAATCGGGTGCAGACATTAATGTGCTGGAACACGCGCTCAAACTGGTGGCGGATAAGCGTTAACTTTTCACCAAAAAGTAGTCAAATTCACCACGCCCTGCGCAc                               >  1:4028150/1‑150 (MQ=255)
                                                                                                                            gcgcgcGGCAGTGGCGGCGTTATTTTAGTGAAAATGCCCATAAAGCGGGTGCAGACATTAATGTGCTGGAACACGCGCTCAAACTGGTGGCGGATAAGCGTTAACTTTTCACCAAAAAGTAGTCAAATTCACCACGCCCTGCGCACCGTc                           <  2:4337992/150‑1 (MQ=255)
                                                                                                                             cgcgcgGCAGTGGCGGCGTTATTTAAGTGAAAATGCCCATAAAGCGGGTGCAGACATTAATGTGCTGGAACACGCGCTCAAACTGGTGGCGGATAAGCGTTAACTTTTCACCAAAAAGTAGTCAAATTCACCACGCCCTGCGCACCGTcg                          >  1:3532166/1‑150 (MQ=255)
                                                                                                                             cgcgcgGCAGTGGCGGCGTTATTTAAGTGAAAATGCCCATAAAGCGGGTGCAGACATTAATGTGCTGGAACACGCGCTCAAACTGGTGGCGGATAAGCGTTAACTTTTCACCAAAAAGTAGTCAAATTCACCACGCCCTGCGCACCGTcg                          <  1:3186327/150‑1 (MQ=255)
                                                                                                                                 cgGCAGTGGCGGCGTTATTTAAGTGAAAATGCCCATAAAGCGGGTGCAGACATTAATGTGCTGGAACACGCGCTCAAACTGGTGGCGGATAAGCGTTAACTTTTCACCAAAAAGTAGTCAAATTCACCACGCCCTGCGCACCGTCGCggg                      <  2:1713624/150‑1 (MQ=255)
                                                                                                                                    cAGTGGCGGCGTTATTTAAGTGAAAATGCCCATAAAGCGGGTGCAGACATTAATGTGCTGGAACACGCGCTCAAACTGGTGGCGGATAAGCGTTAACTTTTCACCAAAAAGTAGTCAAATTCACCACGCCCTGCGCACCGTCGCGGGGCg                   <  1:186668/150‑1 (MQ=255)
                                                                                                                                    cAGTGGCGGCGTTATTTAAGTGAAAATGCCCATAAAGCGGGTGCAGACATTAATGTGCTGGAACACGCGCTCAAACTGGTGGCGGATAAGCGTTAACTTTTCACCAAAAAGTAGTCAAATTCACCACGCCCTGCGCACCGTCGCGGGGCg                   >  2:4951359/1‑150 (MQ=255)
                                                                                                                                    cAGTGGCGGCGTTATTTAAGTGAAAATGCCCATAAAGCGGGTGCAGACATTAATGTGCTGGAACACGCGCTCAAACTGGTGGCGGATAAGCGTTAACTTTTCACCAAAAAGTAGTCAAATTCACCACGCCCTGCGCACCGTCGCGGGGCg                   <  1:1339771/150‑1 (MQ=255)
                                                                                                                                    cAGTGGCGGCGTTATTTAAGTGAAAATGCCCATAAAGCGGGTGCAGACATTAATGTGCTGGAACACGCGCTCAAACTGGTGGCGGATAAGCGTTAACTTTTCACCAAAAAGTAGTCAAATTCACCACGCCCTGCGCACCGTCGCGGGGCg                   <  1:669436/150‑1 (MQ=255)
                                                                                                                                      gTGGCGGCGTTATTTAAGTGAAAATGCCCATAAAGCGGGTGCAGACATTAATGTGCTGGAACACGCGCTCAAACTGGTGGCGGATAAGCGTTAACTTTTCACCAAAAAGTAGTCAAATTCACCACGCCCTGCGCACCGTCGCGGGGCGtt                 >  1:1967547/1‑150 (MQ=255)
                                                                                                                                       tGGCGGCGTTATTTAAGTGACAATGCCCATAAAGCGGGTGCAGACATTAATGTGCTGGAACACGCGCTCAAACTGGTGGCGGATAAGCGTTACCTTTTCACCAAAAAGGAGTCAAATTCACCACGCCCTGCGCACCGTCGCGGGGCGttt                <  2:443748/150‑1 (MQ=255)
                                                                                                                                       tGGCGGCGTTATTTAAGTGAAAATGCCCATAAAGCGGGTGCAGACATTAATGTGCTGGAACACGCGCTCAAACTGGTGGCGGATAAGCGTTAACTTTTCACCAAAAAGTAGTCAAATTCACCACGCCCTGCGCACCGTCGCGGGGCGttt                >  2:307193/1‑150 (MQ=255)
                                                                                                                                        ggcggcGTTATTTAAGTGAAAATGCCCATAAAGCGGGTGCAGACATTAATGTGCTGGAACACGCGCTCAAACTGGTGGCGGATAAGCGTTAACTTTTCACCAAAAAGTAGTCAAATTCACCACGCCCTGCGCACCGTCGCGGGGCGtttt               <  2:4328473/150‑1 (MQ=255)
                                                                                                                                         gcggcgTTATTTAAGTGAAAATGCCCATACAGCGGGTGCAGACATTAATGTGCTGGAACACGCGCTCAACCTGGTGGCGGATAAGCGTTAACTTTTCACCAAAAAGTAGTCAAATTCACCACGCCCTGCGCACCGTCGCGGGGCGTTTTg              <  2:5326786/150‑1 (MQ=255)
                                                                                                                                          cggcgTTATTTAAGTGAAAATGCCCATAAAGCGGGTGCAGACATTAATGTGCTGGAACACGCGCTCAAACTGGTGGCGGATAAGCGTTAACTTTTCACCAAAAAGTAGTCAAATTCACCACGCCCTGCGCACCGTCGCGGGGCGTTTTGc             >  1:2814707/1‑150 (MQ=255)
                                                                                                                                            gcgTTATTTAAGTGAAAATGCCCATAAAGCGGGTGCAGACATTAATGTGCTGGAACACGCGCTCAAACTGGTGGCGGATAAGCGTTAACTTTTCACCAAAAAGTAGTCAAATTCACCACGCCCTGCGCACCGTCGCGGGGCGTTTTGCTg           >  1:3048386/1‑150 (MQ=255)
                                                                                                                                              gTTATTTAAGTGAAAATGCCCATAAAGCGGGTTCAGACATTAATGTGCTGGAACACGCGCTCAAACTGGTGGCGGATAAGCGTTAACTATTCACCAAAAAGTAGTCATATTCACCACGCCCTTCGCACCGTCGCGGGGCGTTTTGCTGtt         <  2:1390029/150‑1 (MQ=255)
                                                                                                                                              gTTATTTAAGTGAAAATGCCCATAAAGCGGGTGCAGACATTAATGTGCTGGAACACGCGCTCAAACTGGTGGCGGATAAGCGTTAACTTTTCACCAAAAAGTAGTCAAATTCACCACGCCCTGCGCACCGTCGCGGGGCGTTTTGCTGtt         >  1:1140012/1‑150 (MQ=255)
                                                                                                                                              gTTATTTAAGTGAAAATGCCCATAAAGCGGGTGCAGACATTAATGTGCTGGAACACGCGCTCAAACTGGTGGCGGATAAGCGTTAACTTTTCACCAAAAAGTAGTCAAATTCACCACGCCCTGCGCACCGTCGCGGGGCGTTTTGCTGtt         >  2:2877557/1‑150 (MQ=255)
                                                                                                                                                tATTTAAGTGAAAATGCCCATAAAGCGGGTGCAGACATTAATGTGCTGGAACACGCGCTCAAACTGGTGGCGGATAAGCGTTAACTTTTCACCAAAAAGTAGTCAAATTCACCACGCCCTGCGCACCGTCGCGGGGCGTTTTGCTGTTaa       <  1:3917407/150‑1 (MQ=255)
                                                                                                                                                tATTTAAGTGAAAATGCCCATAAAGCGGGTGCAGACATTAATGTGCTGGAACACGCGCTCAAACTGGTGGCGGATAAGCGTTAACTTTTCACCAAAAAGTAGTCAAATTCACCACGCCCTGCGCACCGTCGCGGGGCGTTTTGCTGTTaa       >  2:3899254/1‑150 (MQ=255)
                                                                                                                                                     aaGTGAAAATGCCCATAAAGCGGGTGCAGACATTAATGTGCTGGAACACGCGCTCAAACTGGTGGCGGATAAGCGTTAACTTTTCACCAAAAAGTAGTCAAATTCACCACGCCCTGCGCACCGTCGCGGGGCGTTTTGCTGTTAAATCaa  <  2:4195661/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
ACCGATGCCGATCCGGTGGCGGTAGTGCGCGCCATGTATCCGTACATTGAGCGTGAACTCAGCCAGGGGACGTATCTCGGCCATATTACCCGGCATATGTTGGGCTTGTTCCAGGGTATTCCTGGCGCGCGGCAGTGGCGGCGTTATTTAAGTGAAAATGCCCATAAAGCGGGTGCAGACATTAATGTGCTGGAACACGCGCTCAAACTGGTGGCGGATAAGCGTTAACTTTTCACCAAAAAGTAGTCAAATTCACCACGCCCTGCGCACCGTCGCGGGGCGTTTTGCTGTTAAATCAA  >  CP009273/4252407‑4252705

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: