Predicted mutation
evidence seq id position mutation annotation gene description
RA AP012306 3,616,653 G→A M697I (ATG→ATA rpoC → RNA polymerase, beta prime subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AP0123063,616,6530GA95.0% 1096.2 / 31.9 297M697I (ATG→ATArpoCRNA polymerase, beta prime subunit
Reads supporting (aligned to +/- strand):  ref base G (9/6);  new base A (128/154);  total (137/160)
Fisher's exact test for biased strand distribution p-value = 2.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.21e-01

Only 100 of 313 total aligned reads displayed.

GAAGAAACACGAAATCATCTCCGAGGCAGAAGCAGAAGTTGCTGAAATTCAGGAGCAGTTCCAGTCTGGTCTGGTAACTGCGGGCGAACGCTACAACAAAGTTATCGATATCTGGGCTGCGGCGAACGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCT  >  AP012306/3616506‑3616799
                                                                                                                                                   |                                                                                                                                                  
gaagaaACACGAAATCATCTCCGAGGCAGAAGCAGAAGTTGCTGAAATTCAGGAGCAGTTCCAGTCTGGTCTGGTAACTGCGGGCGAACGCTACAACAAAGTTATCGATATCTGGGCTGCGGCGAACGATCGTGTATCCAAAGCgatgat                                                                                                                                                  >  1:2806903/1‑150 (MQ=255)
gaagaaACACGAAATCATCTCCGAGGCAGAAGCAGAAGTTGCTGAAATTCAGGAGCAGTTCCAGTCTGGTCTGGTAACTGCGGGCGAACGCTACAACAAAGTTATCGATATCTGGGCTGCGGCGAACGATCGTGTATCCAAAGCgataat                                                                                                                                                  <  1:987969/150‑1 (MQ=255)
gaagaaACACGAAATCATCTCCGAGGCAGAAGCAGAAGTTGCTGAAATTCAGGAGCAGTTCCAGTCTGGTCTGGTAACTGCGGGCGAACGCTACAACAAAGTTATCGATATCTGGGCTGCGGCGAACGATCGTGTATCCAAAGCgataat                                                                                                                                                  <  2:554642/150‑1 (MQ=255)
gaagaaACACGAAATCATCTCCGAGGCAGAAGCAGAAGTTGCTGAAATTCAGGAGCAGTTCCAGTCTGGTCTGGTAACTGCGGGCGAACGCTACAACAAAGTTATCGATATCTGGGCTGCGGCGAACGATCGTGTATCCAAAGCgataat                                                                                                                                                  <  1:1935818/150‑1 (MQ=255)
 aagaaACACGAAATCATCTCCGAGGCAGAAGCAGAAGTTGCTGAAATTCAGGAGCAGTTCCAGTCTGGTCTGGTAACTGCGGGCGAACGCTACAACAAAGTTATCGATATCTGGGCTGCGGCGAACGATCGTGTATCCAAAGCGataatg                                                                                                                                                 >  2:1444239/1‑150 (MQ=255)
 aagaaACACGAAATCATCTCCGAGGCAGAAGCAGAAGTTGCTGAAATTCAGGAGCAGTTCCAGTCTGGTCTGGTAACTGCGGGCGAACGCTACAACAAAGTTATCGATATCTGGGCTGCGGCGAACGATCGTGTATCCAAAGCGataatg                                                                                                                                                 >  2:2106973/1‑150 (MQ=255)
  agaaACACGAAATCATCTCCGAGGCAGAAGCAGAAGTTGCTGAAATTCAGGAGCAGTTCCAGTCTGGTCTGGTAACTGCGGGCGAACGCTACAACAAAGTTATCGATATCTGGGCTGCGGCGAACGATCGTGTATCCAAAGCGATAATgg                                                                                                                                                <  1:4049382/150‑1 (MQ=255)
   gaaACACGAAATCATCTCCGAGGCAGAAGCAGAAGTTGCTGAAATTCAGGAGCAGTTCCAGTCTGGTCTGGTAACTGCGGGCGAACGCTACAACAAAGTTATCGATATCTGGGCTGCGGCGAACGATCGTGTATCCAAAGCGATAATGGa                                                                                                                                               >  1:1525381/1‑150 (MQ=255)
   gaaACACGAAATCATCTCCGAGGCAGAAGCAGAAGTTGCTGAAATTCAGGAGCAGTTCCAGTCTGGTCTGGTAACTGCGGGCGAACGCTACAACAAAGTTATCGATATCTGGGCTGCGGCGAACGATCGTGTATCCAAAGCGATAATGGa                                                                                                                                               >  1:3979320/1‑150 (MQ=255)
     aaCACGAAATCATCTCCGAGGCAGAAGCAGAAGTTGCTGAAATTCAGGAGCAGTTCCAGTCTGGTCTGGTAACTGCGGGCGAACGCTACAACAAAGTTATCGATATCTGGGCTGCGGCGAACGATCGTGTATCCAAAGCGATAATGGATa                                                                                                                                             >  2:2340852/1‑150 (MQ=255)
      acacGAAATCATCTCCGAGGCAGAAGCAGAAGTTGCTGAAATTCAGGAGCAGTTCCAGTCTGGTCTGGTAACTGCGGGCGAACGCTACAACAAAGTTATCGATATCTGGGCTGCGGCGAACGATCGTGTATCCAAAGCGATAATGGATaa                                                                                                                                            <  1:2682978/150‑1 (MQ=255)
       cacGAAATCATCTCCGAGGCAGAAGCAGAAGTTGCTGAAATTCAGGAGCAGTTCCAGTCTGGTCTGGTAACTGCGGGCGAACGCTACAACAAAGTTATCGATATCTGGGCTGCGGCGAACGATCGTGTATCCAAAGCGATAATGGATAAc                                                                                                                                           <  2:704005/150‑1 (MQ=255)
        acGAAATCATCTCCGAGGCAGAAGCAGAAGTTGCTGAAATTCAGGAGCAGTTCCAGTCTGGTCTGGTAACTGCGGGCGAACGCTACAACAAAGTTATCGATATCTGGGCTGCGGCGAACGATCGTGTATCCAAAGCGATAATGGATAAcc                                                                                                                                          >  2:3438088/1‑150 (MQ=255)
         cGAAATCATCTCCGAGGCAGAAGCAGAAGTTGCTGAAATTCAGGAGCAGTTCCAGTCTGGTCTGGTAACTGCGGGCGAACGCTACAACAAAGTTATCGATATCTGGGCTGCGGCGAACGATCGTGTATCCAAAGCGATAATGGATAACCt                                                                                                                                         <  2:850951/150‑1 (MQ=255)
         cGAAATCATCTCCGAGGCAGAAGCAGAAGTTGCTGAAATTCAGGAGCAGTTCCAGTCTGGTCTGGTAACTGCGGGCGAACGCTACAACAAAGTTATCGATATCTGGGCTGCGGCGAACGATCGTGTATCCAAAGCGATAATGGATAACCt                                                                                                                                         <  1:3113253/150‑1 (MQ=255)
            aaTCATCTCCGAGGCAGAAGCAGAAGTTGCTGAAATTCAGGAGCAGTTCCAGTCTGGTCTGGTAACTGCGGGCGAACGCTACAACAAAGTTATCGATATCTGGGCTGCGGCGAACGATCATGTATCCAAAGCGATAATGGATAACCTGCa                                                                                                                                      >  1:1640314/1‑150 (MQ=255)
             atcatcTCCGAGGCAGAAGCAGAAGTTGCTGAAATTCAGGAGCAGTTCCAGTCTGGTCTGGTAACTGCGGGCGAACGCTACAACAAAGTTATCGATATCTGGGCTGCGGCGAACGATCGTGTATCCAAAGCGATAATGGATAACCTGCaa                                                                                                                                     <  2:2576081/150‑1 (MQ=255)
             atcatcTCCGAGGCAGAAGCAGAAGTTGCTGAAATTCAGGAGCAGTTCCAGTCTGGTCTGGTAACTGCGGGCGAACGCTACAACAAAGTTATCGATATCTGGGCTGCGGCGAACGATCGTGTATCCAAAGCGATAATGGATAACCTGCaa                                                                                                                                     <  1:4258443/150‑1 (MQ=255)
              tcatcTCCGAGGCAGAAGCAGAAGTTGCTGAAATTCAGGAGCAGTTCCAGTCTGGTCTGGTAACTGCGGGCGAACGCTACAACAAAGTTATCGATATCTGGGCTGCGGCGAACGATCGTGTATCCAAAGCGATAATGGATAACCTGCaaa                                                                                                                                    >  2:485415/1‑150 (MQ=255)
               catcTCCGAGGCAGAAGCAGAAGTTGCTGAAATTCAGGAGCAGTTCCAGTCTGGTCTGGTAACTGCGGGCGAACGCTACAACAAAGTTATCGATATCTGGGCTGCGGCGAACGATCGTGTATCCAAAGCGATAATGGATAACCTGCAAAc                                                                                                                                   <  2:1830926/150‑1 (MQ=255)
               catcTCCGAGGCAGAAGCAGAAGTTGCTGAAATTCAGGAGCAGTTCCAGTCTGGTCTGGTAACTGCGGGCGAACGCTACAACAAAGTTATCGATATCTGGGCTGCGGCGAACGATCGTGTATCCAAAGCGATAATGGATAACCTGCAAAc                                                                                                                                   <  1:2101167/150‑1 (MQ=255)
                atcTCCGAGGCAGAAGCAGAAGTTGCTGAAATTCAGGAGCAGTTCCAGTCTGGTCTGGTAACTGCGGGCGAACGCTACAACAAAGTTATCGATATCTGGGCTGCGGCGAACGATCGTGTATCCAAAGCGATAATGGATAACCTGCAAACt                                                                                                                                  <  2:3300335/150‑1 (MQ=255)
                 tctcCGAGGCAGAAGCAGAAGTTGCTGAAATTCAGGAGCAGTTCCAGTCTGGTCTGGTAACTGCGGGCGAACGCTACAACAAAGTTATCGATATCTGGGCTGCGGCGAACGATCGTGTATCCAAAGCGATAATGGATAACCTGCAAACTg                                                                                                                                 <  2:596037/150‑1 (MQ=255)
                   tcCGAGGCAGAAGCAGAAGTTGCTGAAATTCAGGAGCAGTTCCAGTCTGGTCTGGTAACTGCGGGCGAACGCTACAACAAAGTTATCGATATCTGGGCTGCGGCGAACGATCGTGTATCCAAAGCGATAATGGATAACCTGCAAACTGaa                                                                                                                               >  2:2129819/1‑150 (MQ=255)
                      gAGGCAGAAGCAGAAGTTGCTGAAATTCAGGAGCAGTTCCAGTCTGGTCTGGTAACTGCGGGCGAACGCTACAACAAAGTTATCGATATCTGGGCTGCGGCGAACGATCGTGTATCCAAAGCGATAATGGATAACCTGCAAACTGAAAcc                                                                                                                            <  2:3483191/150‑1 (MQ=255)
                      gAGGCAGAAGCAGAAGTTGCTGAAATTCAGGAGCAGTTCCAGTCTGGTCTGGTAACTGCGGGCGAACGCTACAACAAAGTTATCGATATCTGGGCTGCGGCGAACGATCGTGTATCCAAAGCGATAATGGATAACCTGCAAACTGAAAcc                                                                                                                            <  2:3385053/150‑1 (MQ=255)
                         gcagaagcagaaGTTGCTGAAATTCAGGAGCAGTTCCAGTCTGGTCTGGTAACTGCGGGCGAACGCTACAACAAAGTTATCGATATCTGGGCTGCGGCGAACGATCGTGTATCCAAAGCGATAATGGATAACCTGCAAACTGAAACCGTg                                                                                                                         <  1:316602/150‑1 (MQ=255)
                           agaagcagaagTTGCTGAAATTCAGGAGCAGTTCCAGTCTGGTCTGGTAACTGCGGGCGAACGCTACAACAAAGTTATCGATATCTGGGCTGCGGCGAACGATCGTGTATCCAAAGCGATAATGGATAACCTGCAAACTGAAACCGTGAt                                                                                                                       >  1:1763517/1‑150 (MQ=255)
                           agaagcagaagTTGCTGAAATTCAGGAGCAGTTCCAGTCTGGTCTGGTAACTGCGGGCGAACGCTACAACAAAGTTATCGATATCTGGGCTGCGGCGAACGATCGTGTATCCAAAGCGATAATGGATAACCTGCAAACTGAAACCGTGAt                                                                                                                       <  2:1563681/150‑1 (MQ=255)
                           agaagcagaagTTGCTGAAATTCAGGAGCAGTTCCAGTCTGGTCTGGTAACTGCGGGCGAACGCTACAACAAAGTTATCGATATCTGGGCTGCGGCGAACGATCGTGTATCCAAAGCGATAATGGATAACCTGCAAACTGAAACCGTGAt                                                                                                                       <  1:1483996/150‑1 (MQ=255)
                           agaagcagaagTTGCTGAAATTCAGGAGCAGTTCCAGTCTGGTCTGGTAACTGCGGGCGAACGCTACAACAAAGTTATCGATATCTGGGCTGCGGCGAACGATCGTGTATCCAAAGCGATAATGGATAACCTGCAAACTGAAACCGTGAt                                                                                                                       >  2:2492856/1‑150 (MQ=255)
                            gaagcagaagTTGCTGAAATTCAGGAGCAGTTCCAGTCTGGTCTGGTAACTGCGGGCGAACGCTACAACAAAGTTATCGATATCTGGGCTGCGGCGAACGATCGTGTATCCAAAGCGATAATGGATAACCTGCAAACTGAAACCGTGAtt                                                                                                                      <  2:2524555/150‑1 (MQ=255)
                                 agaagTTGCTGAAATTCAGGAGCAGTTCCAGTCTGGTCTGGTAACTGCGGGCGAACGCTACAACAAAGTTATCGATATCTGGGCTGCGGCGAACGATCGTGTATCCAAAGCGATAATGGATAACCTGCAAACTGAAACCGTGATTAACCg                                                                                                                 <  1:3452189/150‑1 (MQ=255)
                                 agaagTTGCTGAAATTCAGGAGCAGTTCCAGTCTGGTCTGGTAACTGCGGGCGAACGCTACAACAAAGTTATCGATATCTGGGCTGCGGCGAACGATCGTGTATCCAAAGCGATAATGGATAACCTGCAAACTGAAACCGTGATTAACCg                                                                                                                 <  1:3559649/150‑1 (MQ=255)
                                 agaagTTGCTGAAATTCAGGAGCAGTTCCAGTCTGGTCTGGTAACTGCGGGCGAACGCTACAACAAAGTTATCGATATCTGGGCTGCGGCGAACGATCGTGTATCCAAAGCGATAATGGATAACCTGCAAACTGAAACCGTGATTAACCg                                                                                                                 <  1:2408961/150‑1 (MQ=255)
                                   aagTTGCTGAAATTCAGGAGCAGTTCCAGTCTGGTCTGGTAACTGCGGGCGAACGCTACAACAAAGTTATCGATATCTGGGCTGCGGCGAACGATCGTGTATCCAAAGCGATAATGGATAACCTGCAAACTGAAACCGTGATTAACCGTg                                                                                                               >  2:3086467/1‑150 (MQ=255)
                                   aagTTGCTGAAATTCAGGAGCAGTTCCAGTCTGGTCTGGTAACTGCGGGCGAACGCTACAACAAAGTTATCGATATCTGGGCTGCGGCGAACGATCGTGTATCCAAAGCGATAATGGATAACCTGCAAACTGAAACCGTGATTAACCGTg                                                                                                               <  1:1983878/150‑1 (MQ=255)
                                    agTTGCTGAAATTCAGGAGCAGTTCCAGTCTGGTCTGGTAACTGCGGGCGAACGCTACAACAAAGTTATCGATATCTGGGCTGCGGCGAACGATCGTGTATCCAAAGCGATAATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGa                                                                                                              <  1:123705/150‑1 (MQ=255)
                                     gTTGCTGAAATTCAGGAGCAGTTCCAGTCTGGTCTGGTAACTGCGGGCGAACGCTACAACAAAGTTATCGATATCTGGGCTGCGGCGAACGATCGTGTATCCAAAGCGATAATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGAc                                                                                                             <  1:1371143/150‑1 (MQ=255)
                                      ttGCTGAAATTCAGGAGCAGTTCCAGTCTGGTCTGGTAACTGCGGGCGAACGCTACAACAAAGTTATCGATATCTGGGCTGCGGCGAACGATCGTGTATCCAAAGCGATAATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGACg                                                                                                            >  2:1242181/1‑150 (MQ=255)
                                        gCTGAAATTCAGGAGCAGTTCCAGTCTGGTCTGGTAACTGCGGGCGAACGCTACAACAAAGTTATCGATATCTGGGCTGCGGCGAACGATCGTGTATCCAAAGCGATAATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGACGGt                                                                                                          >  1:2370422/1‑150 (MQ=255)
                                           gAAATTCAGGAGCAGTTCCAGTCTGGTCTGGTAACTGCGGGCGAACGCTACAACAAAGTTATCGATATCTGGGCTGCGGCGAACGATCGTGTATCCAAAGCGATAATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGACGGTCAg                                                                                                       <  2:1520625/150‑1 (MQ=255)
                                            aaaTTCAGGAGCAGTTCCAGTCTGGTCTGGTAACTGCGGGCGAACGCTACAACATAGTTATCGATATCTGGGCTGCGGCGAACGATCGTGTATCCAAAGCGATAATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGACGGTCAgg                                                                                                      <  2:3755671/150‑1 (MQ=255)
                                            aaaTTCAGGAGCAGTTCCAGTCTGGTCTGGTAACTGCGGGCGAACGCTACAACAAAGTTATCGATATCTGGGCTGCGGCGAACGATCGTGTATCCAAAGCGATAATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGACGGTCAgg                                                                                                      <  2:950535/150‑1 (MQ=255)
                                            aaaTTCAGGAGCAGTTCCAGTCTGGTCTGGTAACTGCGGGCGAACGCTACAACAAAGTTATCGATATCTGGGCTGCGGCGAACGATCGTGTATCCAAAGCGATAATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGACGGTCAgg                                                                                                      >  1:738440/1‑150 (MQ=255)
                                              aTTCAGGAGCAGTTCCAGTCTGGTCTGGTAACTGCGGGCGAACGCTACAACAAAGTTATCGATATCTGGGCTGCGGCGAACGATCGTGTATCCAAAGCGATAATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGACGGTCAGGaa                                                                                                    >  1:4143406/1‑150 (MQ=255)
                                                  aGGAGCAGTTCCAGTCTGGTCTGGTAACTGCGGGCGAACGCTACAACAAAGTTATCGATATCTGGGCTGCGGCGAACGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACAGTGATTAACCGTGACGGTCAGGAAgaga                                                                                                <  1:336800/150‑1 (MQ=255)
                                                     aGCAGTTCCAGTCTGGTCTGGTAACTGCGGGCGAACGCTACAACAAAGTTATCGATATCTGGGCTGCGGCGAACGATCGTGTATCCAAAGCGATAATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGACGGTCAGGAAGAGAAGc                                                                                             <  2:3729388/150‑1 (MQ=255)
                                                         gTTCCAGTCTGGTCTGGTAACTGCGGGCGAACGCTACAACAAAGTTATCGATATCTGGGCTGCGGCGAACGATCGTGTATCCAAAGCGATAATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGt                                                                                         <  2:2787060/150‑1 (MQ=255)
                                                            ccAGTCTGGTCTGGTAACTGCGGGCGAACGCTACAACAAAGTTATCGATATCTGGGCTGCGGCGAACGATCGTGTATCCAAAGCGATAATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTc                                                                                      <  1:3268579/150‑1 (MQ=255)
                                                                    gtctggtAACTGCGGGCGAACGCTACAACAAAGTTATCGATATCTGGGCTGCGGCGAACGATCGTGTATCCAAAGCGATAATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCAACa                                                                              <  1:2895721/150‑1 (MQ=255)
                                                                        ggtAACTGCGGGCGAACGCTACAACAAAGTTATCGATATCTGGGCTGCGGCGAACGATCGTGTATCCAAAGCGATAATCGATAACCTGCAAACTGAAACCGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCAACAGCAt                                                                          <  1:5079/150‑1 (MQ=255)
                                                                         gtAACTGCGGGCGAACGCTACAACAAAGTTATCGATATCTGGGCTGCGGCGAACGATCGTGTATCCAAAGCGATAATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCAACAGCATc                                                                         >  2:1600443/1‑150 (MQ=255)
                                                                          tAACTGCGGGCGAACGCTACAACAAAGTTATCGATATCTGGGCTGCGGCGAACGATCGTGTATCCAAAGCGATAATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCAACAGCATCt                                                                        <  2:2921657/150‑1 (MQ=255)
                                                                           aaCTGCGGGCGAACGCTACAACAAAGTTATCGATATCTGGGCTGCGGCGAACGATCGTGTATCCAAAGCGATAATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCAACAGCATCTa                                                                       <  2:2634272/150‑1 (MQ=255)
                                                                            aCTGCGGGCGAACGCTACAACAAAGTTATCGATATCTGGGCTGCGGCGAACGATCGTGTATCCAAAGCGATAATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCAACAGCATCTAc                                                                      <  1:4114303/150‑1 (MQ=255)
                                                                              tGCGGGCGAACGCTACAACAAAGTTATCGATATCTGGGCTGCGGCGAACGATCGTGTATCCAAAGCGATAATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCAACAGCATCTACat                                                                    <  2:1922262/150‑1 (MQ=255)
                                                                              tGCGGGCGAACGCTACAACAAAGTTATCGATATCTGGGCTGCGGCGAACGATCGTGTATCCAAAGCGATAATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCAACAGCATCTACat                                                                    <  2:4228125/150‑1 (MQ=255)
                                                                               gCGGGCGAACGCTACAACAAAGTTATCGATATCTGGGCTGCGGCGAACGATCGTGTATCCAAAGCGATAATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCAACAGCATCTACatg                                                                   <  1:3748867/150‑1 (MQ=255)
                                                                                cGGGCGAACGCTACAACAAAGTTATCGATATCTGGGCTGCGGCGAACGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCAACAGCATCTACatga                                                                  >  2:2327380/1‑150 (MQ=255)
                                                                                 gggCGAACGCTACAACAAAGTTATCGATATCTGGGCTGCGGCGAACGATCGTGTATCCAAAGCGATAATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCAACAGCATCTACatgat                                                                 <  2:3149844/150‑1 (MQ=255)
                                                                                    cGAACGCTACAACAAAGTTATCGATATCTGGGCTGCGGCGAACGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCAACAGCATCTACATGATGGc                                                              >  2:2610417/1‑150 (MQ=255)
                                                                                       aCGCTACAACAAAGTTATCGATATCTGGGCTGCGGCGAACGATCGTGTATCCAAAGCGATAATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCAACAGCATCTACATGATGGCCGa                                                           <  2:809997/150‑1 (MQ=255)
                                                                                        cGCTACAACAAAGTTATCGATATCTGGGCTGCGGCGAACGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCAACAGCATCTACATGATGGCCGAc                                                          <  1:1055962/150‑1 (MQ=255)
                                                                                          cTACAACAAAGTTATCGATATCTGGGCTGCGGCGAACGATCGTGTATCCAAAGCGATAATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCAACAGCATCTACATGATGGCCGACTc                                                        >  1:1151977/1‑150 (MQ=255)
                                                                                           tACAACAAAGTTATCGATATCTGGGCTGCGGCGAACGATCGTGTATCCAAAGCGATAATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCAACAGCATCTACATGATGGCCGACTcc                                                       >  1:2595845/1‑150 (MQ=255)
                                                                                           tACAACAAAGTTATCGATATCTGGGCTGCGGCGAACGATCGTGTATCCAAAGCGATAATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCAACAGCATCTACATGATGGCCGACTcc                                                       >  2:1343851/1‑150 (MQ=255)
                                                                                            acaacaAAGTTATCGATATCTGGGCTGCGGCGAACGATCGTGTATCCAAAGCGATAATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCg                                                      <  2:105213/150‑1 (MQ=255)
                                                                                             caacaaAGTTATCGATATCTGGGCTGCGGCGAACGATCGTGTATCCAAAGCGATAATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCgg                                                     <  2:1119584/150‑1 (MQ=255)
                                                                                                caaAGTTATCGATATCTGGGCTGCGGCGAACGATCGTGTATCCAAAGCGATAATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTgc                                                  >  2:3521553/1‑150 (MQ=255)
                                                                                                caaAGTTATCGATATCTGGGCTGCGGCGAACGATCGTGTATCCAAAGCGATAATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTgc                                                  >  2:4065819/1‑150 (MQ=255)
                                                                                                caaAGTTATCGATATCTGGGCTGCGGCGAACGATCGTGTATCCAAAGCGATAATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTgc                                                  <  1:3201620/150‑1 (MQ=255)
                                                                                                 aaaGTTATCGATATCTGGGCTGCGGCGAACGATCGTGTATCCAAAGCGATAATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTgcg                                                 >  2:1986941/1‑150 (MQ=255)
                                                                                                     ttATCGATATCTGGGCTGCGGCGAACGATCGTGTATCCAAAGCGATAATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTg                                             >  2:3864253/1‑150 (MQ=255)
                                                                                                      tATCGATATCTGGGCTGCGGCGAACGATCGTGTATCCAAAGCGATAATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTgg                                            >  2:2236196/1‑150 (MQ=255)
                                                                                                         cGATATCTGGGCTGCGGCGAACGATCGTGTATCCAAAGCGATAATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTc                                         >  2:2014774/1‑150 (MQ=255)
                                                                                                         cGATATCTGGGCTGCGGCGAACGATCGTGTATCCAAAGCGATAATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTc                                         <  2:2839570/150‑1 (MQ=255)
                                                                                                          gATATCTGGGCTGCGGCGAACGATCGTGTATCCAAAGCGATAATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCt                                        >  2:1095425/1‑150 (MQ=255)
                                                                                                                tGGGCTGCGGCGAACGATCGTGTATCCAAAGCGATAATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGca                                  >  1:468812/1‑150 (MQ=255)
                                                                                                                  ggCTGCGGCGAACGATCGTGTATCCAAAGCGATAATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGcaca                                <  2:3550041/150‑1 (MQ=255)
                                                                                                                   gCTGCGGCGAACGATCGTGTATCCAAAGCGATAATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAg                               >  1:686083/1‑150 (MQ=255)
                                                                                                                     tGCGGCGAACGATCGTGTATCCAAAGCGATAATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGAt                             >  2:2038359/1‑150 (MQ=255)
                                                                                                                      gcggcgAACGATCGTGTATCCAAAGCGATAATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGAtt                            <  1:2518485/150‑1 (MQ=255)
                                                                                                                           gAACGATCGTGTATCCAAAGCGATAATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCa                       >  1:937512/1‑150 (MQ=255)
                                                                                                                             aCGATCGTGTATCCAAAGCGATAATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGc                     >  2:2422008/1‑150 (MQ=255)
                                                                                                                               gATCGTGTATCCAAAGCGATAATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCtt                   <  1:1013619/150‑1 (MQ=255)
                                                                                                                               gATCGTGTATCCAAAGCGATAATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCtt                   >  1:2969288/1‑150 (MQ=255)
                                                                                                                                aTCGTGTATCCAAAGCGATAATGGATAACCTGCAAACTGAAACCATGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTg                  >  2:2668291/1‑150 (MQ=255)
                                                                                                                                   gtgtATCCAAAGCGATAATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTg               >  2:4006817/1‑150 (MQ=255)
                                                                                                                                    tgtATCCAAAGCGATAATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTgg              <  1:2875073/150‑1 (MQ=255)
                                                                                                                                    tgtATCCAAAGCGATAATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTgg              >  1:850603/1‑150 (MQ=255)
                                                                                                                                    tgtATCCAAAGCGATAATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTAAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTgg              >  2:2356780/1‑150 (MQ=255)
                                                                                                                                     gtATCCAAAGCGATAATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGt             >  2:4053718/1‑150 (MQ=255)
                                                                                                                                      tATCCAAAGCGATAATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCAACAGCATCTACATAATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTa            >  2:39596/1‑150 (MQ=255)
                                                                                                                                        tCCAAAGCGATAATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATg          <  2:281169/150‑1 (MQ=255)
                                                                                                                                           aaaGCGATAATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGt       >  1:1402651/1‑150 (MQ=255)
                                                                                                                                            aaGCGATAATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTg      >  2:2187980/1‑150 (MQ=255)
                                                                                                                                               cGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTc   <  1:2576867/150‑1 (MQ=255)
                                                                                                                                                gataatGGATAACCTGCAAACTGAAACCGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCt  >  1:1881346/1‑150 (MQ=255)
                                                                                                                                                gataatGGATAACCTGCAAACTGAAACCGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCt  >  2:1538502/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
GAAGAAACACGAAATCATCTCCGAGGCAGAAGCAGAAGTTGCTGAAATTCAGGAGCAGTTCCAGTCTGGTCTGGTAACTGCGGGCGAACGCTACAACAAAGTTATCGATATCTGGGCTGCGGCGAACGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCT  >  AP012306/3616506‑3616799

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 28 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: