Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,099,851 C→A 6.9% S260* (TCG→TAG)  nupG → nucleoside transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,099,8510CA 6.9% 918.4 / 28.8 247S260* (TCG→TAG) nupGnucleoside transporter
Reads supporting (aligned to +/- strand):  ref base C (95/135);  new base A (9/8);  total (104/143)
Fisher's exact test for biased strand distribution p-value = 4.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.58e-01

Only 100 of 247 total aligned reads displayed.

GCAATCTTCTTTATCTTCTCAATGCTGCTGGGCGCGGAACTGCAGATTACCAACATGTTCGGTAATACCTTCCTGCACAGCTTCGACAAAGATCCGATGTTTGCCAGCAGCTTTATTGTGCAGCATGCGTCAATCATCATGTCGATTTCGCAGATCTCTGAAACCCTGTTCATTCTGACCATCCCGTTCTTCTTAAGCCGCTACGGTATTAAGAACGTAATGATGATCAGTATTGTGGCGTGGATCCTGCGTTTTGCGCTGTTTGCTTACGGCGACCCGACTCCGTTCGGTACTGT  >  CP009273/3099703‑3099998
                                                                                                                                                    |                                                                                                                                                   
gCAATCTTCTTTATCTTCTCAATGCTGCTGGGCGCGGAACTGCAGATTACCAACATGTTCGGTAATACCTTCCTGCACAGCTTCGACAAAGATCCGATGTTTGCCAGCAGCTTTATTGTGCAGCATGCGTCAATCATCATGTCGATTTCg                                                                                                                                                    <  1:4399611/150‑1 (MQ=255)
gCAATCTTCTTTATCTTCTCAATGCTGCTGGGCGCGGAACTGCAGATTACCAACATGTTCGGTAATACCTTCCTGCACAGCTTCGACAAAGATCCGATGTTTGCCAGCAGCTTTATTGTGCAGCATGCGTCAATCATCATGTCGATTTCg                                                                                                                                                    >  1:1294374/1‑150 (MQ=255)
gCAATCTTCTTTATCTTCTCAATGCTGCTGGGCGCGGAACTGCAGATTACCAACATGTTCGGTAATACCTTCCTGCACAGCTTCGACAAAGATCCGATGTTTGCCAGCAGCTTTATTGTGCAGCATGCGTCAATCATCATGTCGATTTCg                                                                                                                                                    <  1:719041/150‑1 (MQ=255)
   aTCTTCTTTATCTTCTCAATGCTGCTGGGCGCGGAACTGCAGATTACCAACATGTTCGGTAATACCTTCCTGCACAGCTTCGACAAAGATCCGATGTTTGCCAGCAGCTTTATTGTGCAGCATGCGTCAATCATCATGTCGATTTCGCAg                                                                                                                                                 <  1:1888736/150‑1 (MQ=255)
    tcttctTTATCTTCTCAATGCTGCTGGGCGCGGAACTGCAGATTACCAACATGTTCGGTAATACCTTCCTGCACAGCTTCGACAAAGATCCGATGTTTGCCAGCAGCTTTATTGTGCAGCATGCGTCAATCATCATGTCGATTTCGCAGa                                                                                                                                                <  2:1549005/150‑1 (MQ=255)
        cttTATCTTCTCAATGCTGCTGGGCGCGGAACTGCAGATTACCAACATGTTCGGTAATACCTTCCTGCACAGCTTCGACAAAGATCCGATGTTTGCCAGCAGCTTTATTGTGCAGCATGCGTCAATCATCATGTCGATTTCGCAGAtctc                                                                                                                                            <  2:3219567/150‑1 (MQ=255)
        cttTATCTTCTCAATGCTGCTGGGCGCGGAACTGCAGATTACCAACATGTTCGGTAATACCTTCCTGCACAGCTTCGACAAAGATCCGATGTTTGCCAGCAGCTTTATTGTGCAGCATGCGTCAATCATCATGTCGATTTCGCAGAtctc                                                                                                                                            <  2:2775310/150‑1 (MQ=255)
            aTCTTCTCAATGCTGCTGGGCGCGGAACTGCAGATTACCAACATGTTCGGTAATACCTTCCTGCACAGCTTCGACAAAGATCCGATGTTTGCCAGCAGCTTTATTGTGCAGCATGCGTCAATCATCATGTCGATTTCGCAGATCTCTGaa                                                                                                                                        <  2:273163/150‑1 (MQ=255)
               ttctCAATGCTGCTGGGCGCGGAACTGCAGATTACCAACATGTTCGGTAATACCTTCCTGCACAGCTTCGACAAAGATCCGATGTTTGCCAGCAGCTTTATTGTGCAGCATGCGTCAATCATCATGTCGATTTAGCAGATCTCTGAAAcc                                                                                                                                     <  1:4260450/150‑1 (MQ=255)
               ttctCAACGCTGCTGGGCGCGGAACTGCAGATTACCAACATGTTCGGTAATACCTTCCTGCACAGCTTCGACAAAGATCCGATGTTTGCCAGCAGCTTTATTGTGCAGCATGCGTCAATCATCATGTCGATTTAGCAGATCTCTGAAAcc                                                                                                                                     <  1:4260508/150‑1 (MQ=255)
                 ctcAATGCTGCTGGGCGCGGAACTGCAGATTACCAACATGTTCGGTAATACCTTCCTGCACAGCTTCGACAAAGATCCGATGTTTGCCAGCAGCTTTATTGTGCAGCATGCGTCAATCATCATGTCGATTTCGCAGATCTCTGAAACCCt                                                                                                                                   <  1:3389895/150‑1 (MQ=255)
                 ctcAATGCTGCTGGGCGCGGAACTGCAGATTACCAACATGTTCGGTAATACCTTCCTGCACAGCTTCGACAAAGATCCGATGTTTGCCAGCAGCTTTATTGTGCAGCATGCGTCAATCATCATGTCGATTTCGCAGATCTCTGAAACCCt                                                                                                                                   <  1:570684/150‑1 (MQ=255)
                  tcAATGCTGCTGGGCGCGGAACTGCAGATTACCAACATGTTCGGTAATACCTTCCTGCACAGCTTCGACAAAGATCCGATGTTTGCCAGCAGCTTTATTGTGCAGCATGCGTCAATCATCATGTCGATTTCGCAGATCTCTGAAACCCTg                                                                                                                                  <  1:771961/150‑1 (MQ=255)
                    aaTGCTGCTGGGCGCGGAACTGCAGATTACCAACATGTTCGGTAATACCTTCCTGCACAGCTTCGACAAAGATCCGATGTTTGCCAGCAGCTTTATTGTGCAGCATGCGTCAATCATCATGTCGATTTCGCAGATCTCTGAAACCCTGtt                                                                                                                                >  1:3039914/1‑150 (MQ=255)
                    aaTGCTGCTGGGCGCGGAACTGCAGATTACCAACATGTTCGGTAATACCTTCCTGCACAGCTTCGACAAAGATCCGATGTTTGCCAGCAGCTTTATTGTGCAGCATGCGTCAATCATCATGTCGATTTCGCAGATCTCTGAAACCCTGtt                                                                                                                                <  2:2428879/150‑1 (MQ=255)
                    aaTGCTGCTGGGCGCGGAACTGCAGATTACCAACATGTTCGGTAATACCTTCCTGCACAGCTTCGACAAAGATCCGATGTTTGCCAGCAGCTGTATTGTGCAGCATGCGTCAATCATCATGTCGAGTTCGCAGATCTCTGAAACCCGGtt                                                                                                                                >  2:3403534/1‑150 (MQ=255)
                     aTGCTGCTGTGAGCTTAACTTCAGATTACCAACATGTTCGGTAATACCTTCCTTCACAGATTCGACAAAGATCAGATGTTTGCCAGCAGCTTTATTGTGCAGCATGAGTCAATCATCATGTCTATTTCGCAGATCTCTGAAACCCTGTTc                                                                                                                               <  2:1165156/150‑1 (MQ=255)
                     aTGCTGCTGGGCGCGGAACTGCAGATTACCAACATGTTCGGTAATACCTTCCTGCACAGCTTCGACAAAGATCCGATGTTTGCCAGCAGCTTTATTGTGCAGCATGCGTCAATCATCATGTCGATTTCGCAGATCTCTGAAACCCTGTTc                                                                                                                               <  2:4488980/150‑1 (MQ=255)
                     aTGCTGCTGGGCGCGGAACTGCAGATTACCAACATGTTCGGTAATACCTTCCTGCACAGCTTCGACAAAGATCCGATGTTTGCCAGCAGCTTTATTGTGCAGCATGCGTCAATCATCATGTCGATTTCGCAGATCTCTGAAACCCTGTTc                                                                                                                               <  1:3234578/150‑1 (MQ=255)
                     aTGCTGCTGGGCGCGGAACTGCAGATTACCAACATGTTCGGTAATACCTTCCTGCACAGCTTCGACAAAGATCCGATGTTTGCCAGCAGCTTTATTGTGCAGCATGCGTCAATCATCATGTCGATTTCGCAGATCTCTGAAACCCTGTTc                                                                                                                               <  2:4177827/150‑1 (MQ=255)
                     aTGCTGCTGGGCGCGGAACTGCAGATTACCAACATGTTCGGTAATACCTTCCTGCACAGCTTCGACAAAGATCCGATGTTTGCCAGCAGCTTTATTGTGCAGCATGCGTCAATCATCATGTCGATTTCGCAGATCTCTGAAACCCTGTTc                                                                                                                               <  1:1857792/150‑1 (MQ=255)
                      tgctgcTGGGCGCGGAACTGCAGATTACCAACATGTTCGGTAATACCTTCCTGCACAGCTTCGACAAAGATCCGATGTTTGCCAGCAGCTTTATTGTGCAGCATGCGTCAATCATCATGTCGATTTCGCAGATCTCTGAAACCCTGTTCa                                                                                                                              <  1:1389400/150‑1 (MQ=255)
                      tgctgcTGGGCGCGGAACTGCAGATTACCAACATGTTCGGTAATACCTTCCTGCACAGCTTCGACAAAGATCCGATGTTTGCCAGCAGCTTTATTGTGCAGCATGCGTCAATCATCATGTCGATTTCGCAGATCTCTGAAACCCTGTTCa                                                                                                                              <  1:1224894/150‑1 (MQ=255)
                      tgctgcTGGGCGCGGAACTGCAGATTACCAACATGTTCGGTAATACCTTCCTGCACAGCTTCGACAAAGATCCGATGTTTGCCAGCAGCTTTATTGTGCAGCATGCGTCAATCATCATGTCGATTTCGCAGATCTCTGAAACCCTGTTCa                                                                                                                              <  2:3968337/150‑1 (MQ=255)
                       gctgctGGGCGCGGAACTGCAGATTACCAACATGTTCGGTAATACCTTCCTGCACAGCTTCGACAAAGATCCGATGTTTGCCAGCAGCTTTATTGTGCAGCATGCGTCAATCATCATGTCGATTTCGCAGATCTCTGAAACCCTGTTCAt                                                                                                                             <  1:4195779/150‑1 (MQ=255)
                       gctgctGGGCGCGGAACTGCAGATTACCAACATGTTCGGTAATACCTTCCTGCACAGCTTCGACAAAGATCCGATGTTTGCCAGCAGCTTTATTGTGCAGCATGCGTCAATCATCATGTCGATTTCGCAGATCTCTGAAACCCTGTTCAt                                                                                                                             <  1:4098817/150‑1 (MQ=255)
                        ctgctgGGCGCGGAACTGCAGATTACCAACATGTTCGGTAATACCTTCCTGCACAGCTTCGACAAAGATCCGATGTTTGCCAGCAGCTTTATTGTGCAGCATGCGTCAATCATCATGTCGATTTCGCAGATCTCTGAAACCCTGTTCAtt                                                                                                                            >  1:2027516/1‑150 (MQ=255)
                         tgctgGGCGCGGAACTGCAGATTACCAACATGTTCGGTAATACCTTCCTGCACAGCTTCGACAAAGATCCGATGTTTGCCAGCAGCTTTATTGTGCAGCATGCGTCAATCATCATGTCGATTTCGCAGATCTCTGAAACCCTGTTCATTc                                                                                                                           <  1:161210/150‑1 (MQ=255)
                          gctgGGCGCGGAACTGCAGATTACCAACATGTTCGGTAATACCTTCCTGCACAGCTTCGACAAAGATCCGATGTTTGCCAGCAGCTTTATTGTGCAGCATGCGTCAATCATCATGTCGATTTCGCAGATCTCTGAAACCCTGTTCATTCt                                                                                                                          <  1:3020955/150‑1 (MQ=255)
                                  cgGAACTGCAGATTACCAACATGTTCGGTAATACCTTCCTGCACAGCTTCGACAAAGATCCGATGTTTGCCAGCAGCTTTATTGTGCAGCATGCGTCAATCATCATGTCGATTTCGCAGATCTCTGAAACCCTGTTCATTCTGACCATcc                                                                                                                  >  2:4112427/1‑150 (MQ=255)
                                    gAACTGCAGATTACCAACATGTTCGGTAATACCTTCCTGCACAGCTTCGACAAAGATCCGATGTTTGCCAGCAGCTTTATTGTGCAGCATGCGTCAATCATCATGTCGATTTCGCAGATCTCTGAAACCCTGTTCATTCTGACCATCCCg                                                                                                                <  1:3667411/150‑1 (MQ=255)
                                    gAACTGCAGATTACCAACATGTTCGGTAATACCTTCCTGCACAGCTTCGACAAAGATCCGATGTTTGCCAGCAGCTTTATTGTGCAGCATGCGTCAATCATCATGTCGATTTCGCAGATCTCTGAAACCCTGTTCATTCTGACCATCCCg                                                                                                                <  2:4352557/150‑1 (MQ=255)
                                    gAACTGCAGATTACCAACATGTTCGGTAATACCTTCCTGCACAGCTTCGACAAAGATCCGATGTTTGCCAGCAGCTTTATTGTGCAGCATGCGTCAATCATCATGTCGATTTCGCAGATCTCTGAAACCCTGTTCATTCTGACCATCCCg                                                                                                                <  2:4012584/150‑1 (MQ=255)
                                    gAACTGCAGATTACCAACATGTTCGGTAATACCTTCCTGCACAGCTTCGACAAAGATCCGATGTTTGCCAGCAGCTTTATTGTGCAGCATGCGTCAATCATCATGTCGATTTAGCAGATCTCTGAAACCCTGTTCATTCTGACCATCCCg                                                                                                                <  1:3049011/150‑1 (MQ=255)
                                     aaCTGCAGATTACCAACATGTTCGGTAATACCTTCCTGCACAGCTTCGACAAAGATCCGATGTTTGCCAGCAGCTTTATTGTGCAGCATGCGTCAATCATCATGTCGATTTCGCAGATCTCTGAAACCCTGTTCATTCTGACCATCCCGt                                                                                                               >  1:2258490/1‑150 (MQ=255)
                                          cAGATTACCAACATGTTCGGTAATACCTTCCTGCACAGCTTCGACAAAGATCCGATGTTTGCCAGCAGCTTTATTGTGCAGCATGCGTCAATCATCATGTCGATTTCGCAGATCTCTGAAACCCTGTTCATTCTGACCATCCCGttcttc                                                                                                          >  1:583216/1‑150 (MQ=255)
                                          cAGATTACCAACATGTTCGGTAATACCTTCCTGCACAGCTTCGACAAAGATCCGATGTTTGCCAGCAGCTTTATTGTGCAGCATGCGTCAATCATCATGTCGATTTCGCAGATCTCTGAAACCCTGTTCATTCTGACCATCCCGttcttc                                                                                                          <  1:1628117/150‑1 (MQ=255)
                                           aGATTACCAACATGTTCGGTAATACCTTCCTGCACAGCTTCGACAAAGATCCGATGTTTGCCAGCAGCTTTATTGTGCAGCATGCGTCAATCATCATGTCGATTTCGCAGATCTCTGAAACCCTGTTCATTCTGACCATCCCGTtcttct                                                                                                         <  1:826493/150‑1 (MQ=255)
                                            gATTACCAACATGTTCGGTAATACCTTCCTGCACAGCTTCGACAAAGATCCGATGTTTGCCAGCAGCTTTATTGTGCAGCATGCGTCAATCATCATGTCGATTTCGCAGATCTCTGAAACCCTGTTCATTCTGACCATCCCGTTcttctt                                                                                                        <  2:2558543/150‑1 (MQ=255)
                                              ttACCAACATGTTCGGTAATACCTTCCTGCACAGCTTCGACAAAGATCCGATGTTTGCCAGCAGCTTTATTGTGCAGCATGCGTCAATCATCATGTCGATTTCGCAGATCTCTGAAACCCTGTTCATTCTGACCATCCCGTTCTTCTTaa                                                                                                      >  1:1403590/1‑150 (MQ=255)
                                               tACCAACATGTTCGGTAATACCTTCCTGCACAGCTTCGACAAAGATCCGATGTTTGCCAGCAGCTTTATTGTGCAGCATGCGTCAATCATCATGTCGATTTCGCAGATCTCTGAAACCCTGTTCATTCTGACCATCCCGTTCTTCTTAAg                                                                                                     <  1:174009/150‑1 (MQ=255)
                                               tACCAACATGTTCGGTAATACCTTCATGCACAGCTTCGACAAAGATCCGATTTTAGCCATCAGCTTTAATGTGCAGCATGCGACAATCATCATGTCGATTTCGCAGATCTCTGAAACCCTGTTCATTCTGACCATCCCGTTCTTCTTAAg                                                                                                     >  2:1853383/1‑150 (MQ=255)
                                                 ccAACATGTTCGGTAATACCTTCCTGCACAGCTTCGACAAAGATCCGATGTTTGCCAGCAGCTTTATTGTGCAGCATGCGTCAATCATCATGTCGATTTCGCAGATCTCTGAAACCCTGTTCATTCTGACCATCCCGTTCTTCTTAAGcc                                                                                                   <  2:2329771/150‑1 (MQ=255)
                                                 ccAACATGTTCGGTAATACCTTCCTGCACAGCTTCGACAAAGATCCGATGTTTGCCAGCAGCTTTATTGTGCAGCATGCGTCAATCATCATGTCGATTTCGCAGATCTCTGAAACCCTGTTCATTCTGACCATCCCGTTCTTCTTAAGcc                                                                                                   <  2:2329904/150‑1 (MQ=255)
                                                              tAATACCTTCCTGCACAGCTTCGACAAAGATCCGATGTTTGCCAGCAGCTTTATTGTGCAGCATGCGTCAATCATCATGTCGATTTCGCAGATCTCTGAAACCCTGTTCATTCTGACCATCCCGTTCTTCTTAAGCCGCTACGGTATTaa                                                                                      >  2:828938/1‑150 (MQ=255)
                                                                aTACCTTCCTGCACAGCTTCGACAAAGATCCGATGTTTGCCAGCAGCTTTATTGTGCAGCATGCGTCAATCATCATGTCGATTTCGCAGATCTCTGAAACCCTGTTCATTCTGACCATCCCGTTCTTCTTAAGCCGCTACGGTCTTAAGa                                                                                    >  2:3540439/1‑150 (MQ=255)
                                                                aTACCTTCCTGCACAGCTTCGACAAAGATCCGATGTTTGCCAGCAGCTTTATTGTGCAGCATGCGTCAATCATCATGTCGATTTCGCAGATCTCTGAAACCCTGTTCATTCTGACCATCCCGTTCTTCTTAAGCCGCTACGGTATTAAGa                                                                                    >  2:1473101/1‑150 (MQ=255)
                                                                aTACCTTCCTGCACAGCTTCGACAAAGATCCGATGTTTGCCAGCAGCTTTATTGTGCAGCATGCGTCAATCATCATGTCGATTTCGCAGATCTCTGAAACCCTGTTCATTCTGACCATCCCGTTCTTCTTAAGCCGCTACGGTATTAAGa                                                                                    >  2:2557842/1‑150 (MQ=255)
                                                                 tACCTTCCTGCACAGCTTCGACAAAGATCCGATGTTTGCCAGCAGCTTTATTGTGCAGCATGCGTCAATCATCATGTCGATTTCGCAGATCTCTGAAACCCTGTTCATTCTGACCATCCCGTTCTTCTTAAGCCGCTACGGTATTAAGaa                                                                                   >  2:850381/1‑150 (MQ=255)
                                                                 tACCTTCCTGCACAGCTTCGACAAAGATCCGATGTTTGCCAGCAGCTTTATTGTGCAGCATGCGTCAATCATCATGTCGATTTCGCAGATCTCTGAAACCCTGTTCATTCTGACCATCCCGTTCTTCTTAAGCCGCTACGGTATTAAGaa                                                                                   <  2:2971450/150‑1 (MQ=255)
                                                                               gCTTCGACAAAGATCCGATGTTTGCCAGCAGCTTTATTGTGCAGCATGCGTCAATCATCATGTCGATTTCGCAGATCTCTGAAACCCTGTTCATTCTGACCATCCCGTTCTTCTTAAGCCGCTACGGTATTAAGAACGTAATGATGATCa                                                                     <  2:2559624/150‑1 (MQ=255)
                                                                                     aCAAAGATCCGATGTTTGCCAGCAGCTTTATTGTGCAGCATGCGTCAATCATCATGTCGATTTAGCAGATCTCTGAAACCCTGTTCATTCTGACCATCCCGTTCTTCTTAAGCCGCTACGGTATTAAGAACGTAATGATGATCAGTATtg                                                               >  2:1409723/1‑150 (MQ=255)
                                                                                      cAAAGATCCGATGTTTGCCAGCAGCTTTATTGTGCAGCATGCGTCAATCATCATGTCGATTTCGCAGATCTCTGAAACCCTGTTCATTCTGACCATCCCGTTCTTCTTAAGCCGCTACGGTATTAAGAACGTAATGATGATCAGTATtgt                                                              >  1:3535916/1‑150 (MQ=255)
                                                                                         aGATCCGATGTTTGCCAGCAGCTTTATTGTGCAGCATGCGTCAATCATCATGTCGATTTCGCAGATCTCTGAAACCCTGTTCATTCTGACCATCCCGTTCTTCTTAAGCCGCTACGGTATTAAGAACGTAATGATGATCAGTATTGTGGc                                                           >  1:86923/1‑150 (MQ=255)
                                                                                         aGATCCGATGTTTGCCAGCAGCTTTATTGTGCAGCATGCGTCAATCATCATGTCGATTTCGCAGATCTCTGAAACCCTGTTCATTCTGACCATCCCGTTCTTCTTAAGCCGCTACGGTATTAAGAACGTAATGATGATCAGTATTGTGGc                                                           <  2:3633020/150‑1 (MQ=255)
                                                                                           aTCCGATGTTTGCCAGCAGCTTTATTGTGCAGCATGCGTCAATCATCATGTCGATTTCGCAGATCTCTGAAACCCTGTTCATTCTGACCATCCCGTTCTTCTTAAGCCGCTACGGTATTAAGAACGTaa                                                                              <  2:1692176/129‑1 (MQ=255)
                                                                                           aTCCGATGTTTGCCAGCAGCTTTATTGTGCAGCATGCGTCAATCATCATGTCGATTTCGCAGATCTCTGAAACCCTGTTCATTCTGACCATCCCGTTCTTCTTAAGCCGCTACGGTATTAAGAACGTAATGATGATCAGTATTGTGGCGt                                                         <  2:30359/150‑1 (MQ=255)
                                                                                            tCCGATGTTTGCCAGCAGCTTTATTGTGCAGCATGCGTCAATCATCATGTCGATTTCGCAGATCTCTGAAACCCTGTTCATTCTGACCATCCCGTTCTTCTTAAGCCGCTACGGTATTAAGAACGTAATGATGATCAGTATTGTGGCGTg                                                        <  1:612540/150‑1 (MQ=255)
                                                                                            tCCGATGTTTGCCAGCAGCTTTATTGTGCAGCATGCGTCAATCATCATGTCGATTTCGCAGATCTCTGAAACCCTGTTCATTCTGACCATCCCGTTCTTCTTAAGCCGCTACGGTATTAAGAACGTAATGATGATCAGTATTGTGGCGTg                                                        <  1:2605551/150‑1 (MQ=255)
                                                                                            tCCGATGTTTGCCAGCAGCTTTATTGTGCAGCATGCGTCAATCATCATGTCGATTTCGCAGATCTCTGAAACCCTGTTCATTCTGACCATCCCGTTCTTCTTAAGCCGCTACGGTATTAAGAACGTAATGATGATCAGTATTGTGGCGTg                                                        <  1:2957132/150‑1 (MQ=255)
                                                                                             ccGATGTTTGCCAGCAGCTTTATTGTGCAGCATGCGTCAATCATCATGTCGATTTCGCAGATCTCTGAAACCCTGTTCATTCTGACCATCCCGTTCTTCTTAAGCCGCTACGGTATTAAGAACGTAATGATGATCAGTATTGTGGCGTgg                                                       >  1:2546083/1‑150 (MQ=255)
                                                                                                aTGTTTGCCAGCAGCTTTATTGTGCAGCATGCGTCAATCATCATGTCGATTTCGCAGATCTCTGAAACCCTGTTCATTCTGACCATCCCGTTCTTCTTAAGCCGCTACGGTATTAAGAACGTAATGATGATCAGTATTGTGGCGTGGATc                                                    <  2:1503491/150‑1 (MQ=255)
                                                                                                 tGTTTGCCAGCAGCTTTATTGTGCAGCATGCGTCAATCATCATGTCGATTTCGCAGATCTCTGAAACCCTGTTCATTCTGACCATCCCGTTCTTCTTAAGCCGCTACGGTATTAAGAACGTAATGATGATCAGTATTGTGGCGTGGATcc                                                   <  2:1548270/150‑1 (MQ=255)
                                                                                                 tGTTTGCCAGCAGCTTTATTGTGCAGCATGCGTCAATCATCATGTCGATTTCGCAGATCTCTGAAACCCTGTTCATTCTGACCATCCCGTTCTTCTTAAGCCGCTACGGTATTAAGAACGTAATGATGATCAGTATTGTGGCGTGGATcc                                                   <  1:2778805/150‑1 (MQ=255)
                                                                                                 tGTTTGCCAGCAGCTTTATTGTGCAGCATGCGTCAATCATCATGTCGATTTCGCAGATCTCTGAAACCCTGTTCATTCTGACCATCCCGTTCTTCTTAAGCCGCTACGGTATTAAGAACGTAATGATGATCAGTATTGTGGCGTGGATcc                                                   <  2:1548173/150‑1 (MQ=255)
                                                                                                 tGTTTGCCAGCAGCTTTATTGTGCAGCATGCGTCAATCATCATGTCGATTTCGCAGATCTCTGAAACCCTGTTCATTCTGACCATCCCGTTCTTCTTAAGCCGCTACGGTATTAAGAACGTAATGATGATCAGTATTGTGGCGTGGATcc                                                   <  2:619866/150‑1 (MQ=255)
                                                                                                       ccAGCAGCTTTATTGTGCAGCATGCGTCAATCATCATGTCGATTTCGCAGATCTCTGAAACCCTGTTCATTCTGACCATCCCGTTCTTCTTAAGCCGCTACGGTATTAAGAACGTAATGATGATCAGTATTGTGGCGTGGATCCTGCGtt                                             >  1:3579744/1‑150 (MQ=255)
                                                                                                        cagcagCTTTATTGTGCAGCATGCGTCAATCATCATGTCGATTTCGCAGATCTCTGAAACCCTGTTCATTCTGACCATCCCGTTCTTCTTAAGCCGCTACGGTATTAAGAACGTAATGATGATCAGTATTGTGGCGTGGATCCTGCGttt                                            >  1:2838685/1‑150 (MQ=255)
                                                                                                         agcagcTTTATTGTGCAGCATGCGTCAATCATCATGTCGATTTCGCAGATCTCTGAAACCCTGTTCATTCTGACCATCCCGTTCTTCTTAAGCCGCTACGGTATTAAGAACGTAATGATGATCAGTATTGTGGCGTGGATCCTGCGtttt                                           <  2:1294374/150‑1 (MQ=255)
                                                                                                         agcagcTTTATTGTGCAGCATGCGTCAATCATCATGTCGATTTCGCAGATCTCTGAAACCCTGTTCATTCTGACCATCCCGTTCTTCTTAAGCCGCTACGGTATTAAGAACGTAATGATGATCAGTATTGTGGCGTGGATCCTGCGtttt                                           >  1:1531610/1‑150 (MQ=255)
                                                                                                               tttATTGTGCAGCATGCGTCAATCATCATGTCGATTTCGCAGATCTCTGAAACCCTGTTCATTCTGACCATCCCGTTCTTCTTAAGCCGCTACGGTATTAAGAACGTAATGATGATCAGTATTGTGGCGTGGATCCTGCGTTTTGCGCTg                                     <  2:2113075/150‑1 (MQ=255)
                                                                                                               tttATTGTGCAGCATGCGTCAATCATCATGTCGATTTCGCAGATCTCTGAAACCCTGTTCATTCTGACCATCCCGTTCTTCTTAAGCCGCTACGGTATTAAGAACGTAATGATGATCAGTATTGTGGCGTGGATCCTGCGTTTTGCGCTg                                     <  2:4355362/150‑1 (MQ=255)
                                                                                                               tttATTGTGCAGCATGCGTCAATCATCATGTCGATTTAGCAGATCTCTGAAACCCTGTTCATTCTGACCATCCCGTTCTTCTTAAGCCGCTACGGTATTAAGAACGTAATGATGATCAGTATTGTGGCGTGGATCCTGCGTTTTGCGCTg                                     >  2:4115128/1‑150 (MQ=255)
                                                                                                                 tATTGTGCAGCATGCGTCAATCATCATGTCGATTTCGCAGATCTCTGAAACCCTGTTCATTCTGACCATCCCGTTCTTCTTAAGCCGCTACGGTATTAAGAACGTAATGATGATCAGTATTGTGGCGTGGATCCTGCGTTTTGCGCTGtt                                   <  2:2169196/150‑1 (MQ=255)
                                                                                                                       gcagcaTGCGTCAATCATCATGTCGATTTCGCAGATCTCTGAAACCCTGTTCATTCTGACCATCCCGTTCTTCTTAAGCCGCTACGGTATTAAGAACGTAATGATGATCAGTATTGTGGCGTGGATCCTGCGTTTTGCGCTGTTTGCTTa                             >  1:1686982/1‑150 (MQ=255)
                                                                                                                        cagcaTGCGTCAATCATCATGTCGATTTCGCAGATCTCTGAAACCCTGTTCATTCTGACCATCCCGTTCTTCTTAAGCCGCTACGGTATTAAGAACGTAATGATGATCAGTATTGTGGCGTGGATCCTGCGTTTTGCGCTGTTTGCTTAc                            <  1:1059516/150‑1 (MQ=255)
                                                                                                                        cagcaTGCGTCAATCATCATGTCGATTTCGCAGATCTCTGAAACCCTGTTCATTCTGACCATCCCGTTCTTCTTAAGCCGCTACGGTATTAAGAACGTAATGATGATCAGTATTGTGGCGTGGATCCTGCGTTTTGCGCTGTTTGCTTAc                            <  1:2634594/150‑1 (MQ=255)
                                                                                                                           caTGCGTCAATCATCATGTCGATTTCGCAGATCTCTGAAACCCTGTTCATTCTGACCATCCCGTTCTTCTTAAGCCGCTACGGTATTAAGAACGTAATGATGATCAGTATTGTGGCGTGGATCCTGCGTTTTGCGCTGTTTGCTTACGGc                         >  1:2283938/1‑150 (MQ=255)
                                                                                                                            aTGCGTCAATCATCATGTCGATTTCGCAGATCTCTGAAACCCTGTTCATTCTGACCATCCCGTTCTTCTTAAGCCGCTACGGTATTAAGAACGTAATGATGATCAGTATTGTGGCGTGGATCCTGCGTTTTGCGCTGTTTGCTTACGGCg                        >  2:1025262/1‑150 (MQ=255)
                                                                                                                               cGTCAATCATCATGTCGATTTCGCAGATCTCTGAAACCCTGTTCATTCTGACCATCCCGTTCTTCTTAAGCCGCTACGGTATTAAGAACGTAATGATGATCAGTATTGTGGCGTGGATCCTGCGTTTTGCGCTGTTTGCTTACGGCGAcc                     >  1:4459710/1‑150 (MQ=255)
                                                                                                                               cGTCAATCATCATGTCGATTTCGCAGATCTCTGAAACCCTGTTCATTCTGACCATCCCGTTCTTCTTAAGCCGCTACGGTATTAAGAACGTAATGATGATCAGTATTGTGGCGTGGATCCTGCGTTTTGCGCTGTTTGCTTACGGCGAcc                     <  2:1247881/150‑1 (MQ=255)
                                                                                                                               cGTCAATCATCATGTCGATTTCGCAGATCTCTGAAACCCTGTTCATTCTGACCATCCCGTTCTTCTTAAGCCGCTACGGTATTAAGAACGTAATGATGATCAGTATTGTGGCGTGGATCCTGCGTTTTGCGCTGTTTGCTTACGGCGAcc                     <  2:1248294/150‑1 (MQ=255)
                                                                                                                                gTCAATCATCATGTCGATTTCGCAGATCTCTGAAACCCTGTTCATTCTGACCATCCCGTTCTTCTTAAGCCGCTACGGTATTAAGAACGTAATGATGATCAGTATTGTGGCGTGGATCCTGCGTTTTGCGCTGTTTGCTTACGGCGAccc                    >  1:1837988/1‑150 (MQ=255)
                                                                                                                                  cAATCATCATGTCGATTTCGCAGATCTCTGAAACCCTGTTCATTCTGACCATCCCGTTCTTCTTAAGCCGCTACGGTATTAAGAACGTAATGATGATCAGTATTGTGGCGTGGATCCTGCGTTTTGCGCTGTTTGCTTACGGCGACCCGa                  <  2:3751814/150‑1 (MQ=255)
                                                                                                                                   aaTCATCATGTCGATTTCGCAGATCTCTGAAACCCTGTTCATTCTGACCATCCCGTTCTTCTTAAGCCGCTACGGTATTAAGAACGTAATGATGATCAGTATTGTGGCGTGGATCCTGCGTTTTGCGCTGTTTGCTTACGGCGACCCGAc                 <  2:2503927/150‑1 (MQ=255)
                                                                                                                                   aaTCATCATGTCGATTTCGCAGATCTCTGAAACCCTGTTCATTCTGACCATCCCGTTCTTCTTAAGCCGCTACGGTATTAAGAACGTAATGATGATCAGTATTGTGGCGTGGATCCTGCGTTTTGCGCTGTTTGCTTACGGCGACCCGAc                 >  2:1701233/1‑150 (MQ=255)
                                                                                                                                   aaTCATCATGTCGATTTCGCAGATCTCTGAAACCCTGTTCATTCTGACCATCCCGTTCTTCTTAAGCCGCTACGGTATTAAGAACGTAATGATGATCAGTATTGTGGCGTGGATCCTGCGTTTTGCGCTGTTTGCTTACGGCGACCCGAc                 >  1:1010364/1‑150 (MQ=255)
                                                                                                                                     tcatcaTGTCGATTTCGCAGATCTCTGAAACCCTGTTCATTCTGACCATCCCGTTCTTCTTAAGCCGCTACGGTATTAAGAACGTAATGATGATCAGTATTGTGGCGTGGATCCTGCGTTTTGCGCTGTTTGCTTACGGCGACCCGACTc               >  1:1109999/1‑150 (MQ=255)
                                                                                                                                      catcatGTCGATTTCGCAGATCTCTGAAACCCTGTTCATTCTGACCATCCCGTTCTTCTTAAGCCGCTACGGTATTAAGAACGTAATGATGATCAGTATTGTGGCGTGGATCCTGCGTTTTGCGCTGTTTGCTTACGGCGACCCGACTcc              <  1:2047662/150‑1 (MQ=255)
                                                                                                                                      catcatGTCGATTTCGCAGATCTCTGAAACCCTGTTCATTCTGACCATCCCGTTCTTCTTAAGCCGCTACGGTATTAAGAACGTAATGATGATCAGTATTGTGGCGTGGATCCTGCGTTTTGCGCTGTTTGCTTACGGCGACCCGACTcc              <  1:1438431/150‑1 (MQ=255)
                                                                                                                                         catGTCGATTTCGCAGATCTCTGAAACCCTGTTCATTCTGACCATCCCGTTCTTCTTAAGCCGCTACGGTATTAAGAACGTAATGATGATCAGTATTGTGGCGTGGATCCTGCGTTTTGCGCTGTTTGCTTACGGCGACCCGACTCCGtt           <  2:1434011/150‑1 (MQ=255)
                                                                                                                                          atGTCGATTTCGCAGATCTCTGAAACCCTGTTCATTCTGACCATCCCGTTCTTCTTAAGCCGCTACGGTATTAAGAACGTAATGATGATCAGTATTGTGGCGTGGATCCTGCGTTTTGCGCTGTTTGCTTACGGCGACCCGACTCCGTTc          >  1:2906995/1‑150 (MQ=255)
                                                                                                                                           tGTCGATTTCGCAGATCTCTGAAACCCTGTTCATTCTGACCATCCCGTTCTTCTTAAGCCGCTACGGTATTAAGAACGTAATGATGATCAGTATTGTGGCGTGGATCCTGCGTTTTGCGCTGTTTGCTTACGGCGACCCGACTCCGTTCg         >  2:2450680/1‑150 (MQ=255)
                                                                                                                                           tGTCGATTTCGCAGATCTCTGAAACCCTGTTCATTCTGACCATCCCGTTCTTCTTAAGCCGCTACGGTATTAAGAACGTAATGATGATCAGTATTGTGGCGTGGATCCTGCGTTTTGCGCTGTTTGCTTACGGCGACCCGACTCCGTTCg         >  2:2450628/1‑150 (MQ=255)
                                                                                                                                            gTCGATTTCGCAGATCTCTGAAACCCTGTTAATTCTGACCATCCCGTTCTTCTTAAGCCGCTACGGTATTAAGAACGTAATGATGATCAGTATTGTGGCGTGGATCCTGCGTTTTGCGCTGTTTGCTTACGGCGACCCGACTCCGTTCgg        >  2:1913722/1‑150 (MQ=255)
                                                                                                                                               gATTTCGCAGATCTCTGAAACCCTGTTCATTCTGACCATCCCGTTCTTCTTAAGCCGCTACGGTATTAAGAACGTAATGATGATCAGTATTGTGGCGTGGATCCTGCGTTTTGCGCTGTTTGCTTACGGCGACCCGACTCCGTTCGgtac     >  1:702577/1‑150 (MQ=255)
                                                                                                                                               gATTTCGCAGATCTCTGAAACCCTGTTCATTCTGACCATCCCGTTCTTCTTAAGCCGCTACGGTATTAAGAACGTAATGATGATCAGTATTGTGGCGTGGATCCTGCGTTTTGCGCTGTTTGCTTACGGCGACCCGACTCCGTTCGgtac     <  2:3853822/150‑1 (MQ=255)
                                                                                                                                               gATTTCGCAGATCTCTGAAACCCTGTTCATTCTGACCATCCCGTTCTTCTTAAGCCGCTACGGTATTAAGAACGTAATGATGATCAGTATTGTGGCGTGGATCCTGCGTTTTGCGCTGTTTGCTTACGGCGACCCGACTCCGTTCGgtac     >  2:2199555/1‑150 (MQ=255)
                                                                                                                                               gATTTCGCAGATCTCTGAAACCCTGTTCATTCTGACCATCCCGTTCTTCTTAAGCCGCTACGGTATTAAGAACGTAATGATGATCAGTATTGTGGCGTGGATCCTGCGTTTTGCGCTGTTTGCTTACGGCGACCCGACTCCGTTCGgtac     >  2:2308460/1‑150 (MQ=255)
                                                                                                                                                  ttCGCAGATCTCTGAAACCCTGTTCATTCTGACCATCCCGTTCTTCTTAAGCCGCTACGGTATTAAGAACGTAATGATGATCAGTATTGTGGCGTGGATCCTGCGTTTTGCGCTGTTTGCTTACGGCGACCCGACTCCGTTCGgtactgt  <  2:1074346/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
GCAATCTTCTTTATCTTCTCAATGCTGCTGGGCGCGGAACTGCAGATTACCAACATGTTCGGTAATACCTTCCTGCACAGCTTCGACAAAGATCCGATGTTTGCCAGCAGCTTTATTGTGCAGCATGCGTCAATCATCATGTCGATTTCGCAGATCTCTGAAACCCTGTTCATTCTGACCATCCCGTTCTTCTTAAGCCGCTACGGTATTAAGAACGTAATGATGATCAGTATTGTGGCGTGGATCCTGCGTTTTGCGCTGTTTGCTTACGGCGACCCGACTCCGTTCGGTACTGT  >  CP009273/3099703‑3099998

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: