Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,902,302 C→A 100% D242Y (GAT→TAT)  pyrG ← CTP synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,902,3020CA97.7% 687.7 / 4.5 216D242Y (GAT→TAT) pyrGCTP synthetase
Reads supporting (aligned to +/- strand):  ref base C (1/4);  new base A (88/123);  total (89/127)
Fisher's exact test for biased strand distribution p-value = 6.51e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 216 total aligned reads displayed.

ACTTCACTTACCGGGTTCGCTTCTTCGAAGATAACCTGTTCCCATTCGGACAGATTCGCTTCCGGGCAGTTTAAGCTGAATCGTTTACAAATATAATCGTCCAGCCCCTGAGATTTCAACAGGCCCGGAATTTTATAGATGGAATCGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACAATGCAATCTTCGCACGTTCGTTCGCCGGAACAGCGCGATCTGAACGACAAATCAGGATGTCAGGCTGGATACCGATGGAGAGCAGCTCTTTTACAGA  >  CP009273/2902157‑2902446
                                                                                                                                                 |                                                                                                                                                
aCTTCACTTACCGGGTTCGCTTCTTCGAAGATAACCTGTTCCCATTCGGACAGATTCGCTTCCGGGCAGTTTAAGCTGAATCGTTTACAAATATAATCGTCCAGCCCCTGAGATTTCAACAGGCCCGGAATTTTATAGATGGAATAGa                                                                                                                                                <  2:68432/148‑1 (MQ=255)
 cTTCACTTACCGGGTTCGCTTCTTCGAAGATAACCTGTTCCCATTCGGACAGATTCGCTTCCGGGCAGTTTAAGCTGAATCGTTTACAAATATAATCGTCCAGCCCCTGAGATTTCAACAGGCCCGGAATTTTATAGATGGAATAGAc                                                                                                                                               <  2:1373716/148‑1 (MQ=255)
     aCTTACCGGGTTCGCTTCTTCGAAGATAACCTGTTCCCATTCGGGCAGATTCGCTTACGGGCAGTTTAAGCTGAATCGTTTACAAATATAATCGTCCAGCCCCTGAGATTTCAACAGGCCCGGAATTTTATAGATGGAATAGACGTCt                                                                                                                                           <  1:3221751/148‑1 (MQ=255)
      cTTACCGGGTTCGCTTCTTCGAAGATAACCTGTTCCCATTCGGACAGATTCGCTTCCGGGCAGTTTAAGCTGAATCGTTTACAAATATAATCGTCCAGCCCCTGAGATTTCAACAGGCCCGGAATTTTATAGATGGAATAGACGTCtt                                                                                                                                          <  2:942366/148‑1 (MQ=255)
       ttACCGGGTTCGCTTCTTCGAAGATAACCTGTTCCCATTCGGACAGATTCGCTTCCGGGCAGTTTAAGCTGAATCGTTTACAAATATAATCGTCCAGCCCCTGAGATTTCAACAGGCCCGGAATTTTATAGATGGAATAGACGTCttt                                                                                                                                         >  2:78405/1‑148 (MQ=255)
       ttACCGGGTTCGCTTCTTCGAAGATAACCTGTTCCCATTCGGACAGATTCGCTTCCGGGCAGTTTAAGCTGAATCGTTTACAAATATAATCGTCCAGCCCCTGAGATTTCAACAGGCCCGGAATTTTATAGATGGAATAGACGTCttt                                                                                                                                         >  2:78314/1‑148 (MQ=255)
        tACCGGGTTCGCTTCTTCGAAGATAACCTGTTCCCATTCGGACAGATTCGCTTCCGGGCAGTTTAAGCTGAATCGTTTACAAATATAATCGTCCAGCCCCTGAGATTTCAACAGGCCCGGAATTTTATAGATGGAATAGACGTCTTTc                                                                                                                                        <  2:364760/148‑1 (MQ=255)
        tACCGGGTTCGCTTCTTCGAAGATAACCTGTTCCCATTCGGACAGATTCGCTTCCGGGCAGTTTAAGCTGAATCGTTTACAAATATAATCGTCCAGCCCCTGAGATTTCAACAGGCCCGGAATTTTATAGATGGAATAGACGTCTTTc                                                                                                                                        >  2:1551106/1‑148 (MQ=255)
         aCCGGGTTCGCTTCTTCGAAGATAACATGTTCCCATTCGGACAGATTCGCTTCCGGGCAGTTTAAGCTGAAAAGTTTAAAAATATAATCGTCCAGCACCAGTGATTTCAACAGGCCCGGAATTTGTTAGATGGAATCGACGTCTTTCa                                                                                                                                       >  1:3446257/1‑148 (MQ=255)
            gggTTCGCTTCTTCGAAGATAACCTGTTCCCATTCGGACAGATTCGCTTCCGGGCAGTTTAAGCTGAATCGTTTACAAATATAATCGTCCAGCCCCTGAGATTTCAACAGGCCCGGAATTTTATAGATGGAATAGACGTCTTTCagag                                                                                                                                    >  1:1296481/1‑148 (MQ=255)
              gTTCGCTTCTTCGAAGATAACCTGTTCCCATTCGGACAGATTCGCTTCCGGGCAGTTTAAGCTGAATCGTTTACAAATATAATCGTCCAGCCCCTGAGATTTCAACAGGCCCGGAATTTTATAGATGGAATAGACGTCTTTCAGAGaa                                                                                                                                  <  2:1605701/148‑1 (MQ=255)
              gTTCGCTTCTTCGAAGATAACCTGTTCCCATTCGGACAGATTCGCTTCCGGGCAGTTTAAGCTGAATCGTTTACAAATATAATCGTCCAGCCCCTGAGATTTCAACAGGCCCGGAATTTTATAGATGGAATAGACGTCTTTCAGAGaa                                                                                                                                  <  2:1606714/148‑1 (MQ=255)
               ttCGCTTCTTCGAAGATAACCTGTTCCCATTCGGACAGATTCGCTTCCGGGCAGTTTAAGCTGAATCGTTTACAAATATACTCGTCCAGCCCCTGAGATTTCAACAGGCCCGGAATTTTATAGATGGAATAGACGTCTTTCAGAGaaa                                                                                                                                 <  1:337820/148‑1 (MQ=255)
                  gCTTCTTCGAAGATAACCTGTTCCCATTCGGACAGATTCGCTTCCGGGCAGTTTAAGCTGAATCGTTTACAAATATAATCGTCCAGCCCCTGAGATTTCAACAGGCCCGGAATTTTATAGATGGAATAGACGTCTTTCAGAGAAATaa                                                                                                                              <  1:2897233/148‑1 (MQ=255)
                   cttcttCGAAGATAACCTGTTCCCATTCGGACAGATTCGCTTCCGGGCAGTTTAAGCTGAATCGTTTACAAATATAATCGTCCAGCCCCTGAGATTTCAACAGGCCCGGAATTTTATAGATGGAATAGACGTCTTTCAGAGAAATAAc                                                                                                                             <  2:2746420/148‑1 (MQ=255)
                        tcGAAGATAACCTGTTCCCATTCGGACAGATTCGCTTCCGGGCAGTTTAAGCTGAATCGTTTACAAATATAATCGTCCAGCCCCTGAGATTTCAACAGGCCCGGAATTTTATAGATGGAATAGACGTCTTTCAGAGAAATAACCGCtt                                                                                                                        <  1:1475621/148‑1 (MQ=255)
                         cGAAGATAACCTGTTCCCATTCGGACAGATTCGCTTCCGGGCAGTTTAAGCTGAATCGTTTACAAATATAATCGTCCAGCCCCTGAGATTTCAACAGGCCCGGAATTTTATAGATGGAATAGACGTCTTTCAGAGAAATAACCGCttt                                                                                                                       <  2:2517640/148‑1 (MQ=255)
                         cGAAGATAACCTGTTCCCATTCGGACAGATTCGCTTCCGGGCAGTTTAAGCTGAATCGTTTACAAATATAATCGTCCAGCCCCTGAGATTTCAACAGGCCCGGAATTTTATAGATGGAATAGACGTCTTTCAGAGAAATAACCGCttt                                                                                                                       >  2:2524964/1‑148 (MQ=255)
                           aaGATAACCTGTTCCCATTCGGACAGATTCGCTTCCGGGCAGTTTAAGCTGAATCGTTTACAAATATAATCGTCCAGCCCCTGAGATTTCAACCGGCCCGGAATTTTATATATGGAATAGACGTCTTTCAGAGAAATAACCGCttttt                                                                                                                     >  1:864283/1‑148 (MQ=255)
                           aaGATAACCTGTTCCCATTCGGACAGATTCGCTTCCGGGCAGTTTAAGCTGAATCGTTTACAAATATAATCGTCCAGCCCCTGAGATTTCAACAGGCCCGGAATTTTATAGATGGAATAGACGTCTTTCAGAGAAATAACCGCttttt                                                                                                                     <  2:2959571/148‑1 (MQ=255)
                           aaGATAACCTGTTCCCATTCGGACAGATTCGCTTCCGGGCAGTTTAAGCTGAATCGTTTACAAATATAATCGTCCAGCCCCTGAGATTTCAACAGGCCCGGAATTTTATAGATGGAATAGACGTCTTTCAGAGAAATAACCGCttttt                                                                                                                     <  2:382223/148‑1 (MQ=255)
                            aGATAACCTGTTCCCATTCGGACAGATTCGCTTCCGGGCAGTTTAAGCTGAATCGTTTACAAATATAATCGTCCAGCCCCTGAGATTTCAACAGGCCCGGAATTTTATAGATGGAATAGACGTCTTTCAGAGAAATAACCGCTTTTTc                                                                                                                    <  1:1754414/148‑1 (MQ=255)
                             gATAACCTGTTCCCATTCGGACAGATTCGCTTCCGGGCAGTTTAAGCTGAATCGTTTACAAATATAATCGTCCAGCCCCTGAGATTTCAACAGGCCCGGAATTTTATAGATGGAATAGACGTCTTTCAGAGAAATAACCGCTTTTTcc                                                                                                                   >  2:2442645/1‑148 (MQ=255)
                              aTAACCTGTTCCCATTCGGACAGATTCGCTTCCGGGCAGTTTAAGCTGAATCGTTTACAAATATAATCGTCCAGCCCCTGAGATTTCAACAGGCCCGGAATTTTATAGATGGAATAGACGTCTTTCAGAGAAATAACCGCTTTTTCCg                                                                                                                  >  1:1924672/1‑148 (MQ=255)
                              aTAACCTGTTCCCATTCGGACAGATTCGCTTCCGGGCAGTTTAAGCTGAATCGTTTACAAATATAATCGTCCAGCCCCTGAGATTTCAACAGGCCCGGAATTTTATAGATGGAATAGACGTCTTTCAGAGAAATAACCGCTTTTTCCg                                                                                                                  >  2:1177950/1‑148 (MQ=255)
                               tAACCTGTTCCCATTCGGACAGATTCGCTTCCGGGCAGTTTAAGCTGAATCGTTTACAAATATAATCGTCCAGCCCCTGAGATTTCAACAGGCCCGGAATTTTATAGATGGAATAGACGTCTTTCAGAGAAATAACCGCTTTTTCCgg                                                                                                                 <  2:2129314/148‑1 (MQ=255)
                                aaCCTGTTCCCATTCGGACAGATTCGCTTCCGGGCAGTTTAAGCTGAATCGTTTACAAATATAATCGTCCAGCCCCTGAGATTTCAACAGGCCCGGAATTTTATAGATGGAATAGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGa                                                                                                                <  1:735216/148‑1 (MQ=255)
                                aaCCTGTTCCCATTCGGACAGATTCGCTTCCGGGCAGTTTAAGCTGAATCGTTTACAAATATAATCGTCCAGCCCCTGAGATTTCAACAGGCCCGGAATTTTATAGATGGAATAGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGa                                                                                                                <  1:1049045/148‑1 (MQ=255)
                                  ccTGTTCCCATTCGGACAGATTCGCTTCCGGGCAGTTTAAGCTGAATCGTTTACAAATATAATCGTCCAGCCCCTGAGATTTCAACAGGCCCGGAATTTTATAGATGGAATAGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAAc                                                                                                              <  1:2975778/148‑1 (MQ=255)
                                  ccTGTTCCCATTCGGACAGATTCGCTTCCGGGCAGTTTAAGCTGAATCGTTTACAAATATAATCGTCCAGCCCCTGAGATTTCAACAGGCCCGGAATTTTATAGATGGAATAGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAAc                                                                                                              <  1:1226221/148‑1 (MQ=255)
                                   cTGTTCCCATTCGGACAGATTCGCTTCCGGGCAGTTTAAGCTGAATCGTTTACAAATATAATCGTCCAGCCCCTGAGATTTCAACAGGCCCGGAATTTTATAGATGGAATAGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACa                                                                                                             <  1:3067226/148‑1 (MQ=255)
                                    tGTTCCCATTCGGACAGATTCGCTTCCGTGCAGTTTAAGCTGAATCGTTTACAAATATAATCGTCCAGCCCCTGAGATTTCAACAGGCCCGGAATTTTATAGATGGAATAGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACAt                                                                                                            <  1:1110503/148‑1 (MQ=255)
                                       tCCCATTCGGACAGATTCGCTTCCGGGCAGTTTAAGCTGAATCGTTTACAAATATAATCGTCCAGCCCCTGAGATTTCAACAGGCCCGGAATTTTATAGATGGAATAGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTAc                                                                                                         <  1:2092311/148‑1 (MQ=255)
                                       tCCCATTCGGACAGATTCGCTTCCGGGCAGTTTAAGCTGAATCGTTTACAAATATAATCGTCCAGCCCCTGAGATTTCAACAGGCCCGGAATTTTATAGATGGAATAGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTAc                                                                                                         >  1:572154/1‑148 (MQ=255)
                                       tCCCATTCGGACAGATTCGCTTCCGGGCAGTTTAAGCTGAATCGTTTACAAATATAATCGTCCAGCCCCTGAGATTTCAACAGGCCCGGAATTTTATAGATGGAATAGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTAc                                                                                                         <  1:2792460/148‑1 (MQ=255)
                                       tCCCATTCGGACAGATTCGCTTCCGGGCAGTTTAAGCTGAATCGTTTACAAATATAATCGTCCAGCCCCTGAGATTTCAACAGGCCCGGAATTTTATAGATGGAATAGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTAc                                                                                                         <  2:1249497/148‑1 (MQ=255)
                                         ccATTCGGACAGATTCGCTTCCGGGCAGTTTAAGCTGAATCGTTTACAAATATAATCGTCCAGCCCCTGAGATTTCAACAGGCCCGGAATTTTATAGATGGAATAGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAg                                                                                                       <  2:2633777/148‑1 (MQ=255)
                                             tCGGACAGATTCGCTTCCGGGCAGTTTAAGCTTAATCGTTTACAAATATAATCGTCCAGCCCCTGAGATTTCAACAGGCCCGGAATTTTATAGATGGAATAGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACa                                                                                                   <  2:2948212/148‑1 (MQ=255)
                                             tCGGACAGATTCGCTTCCGGGCAGTTTAAGCTGAATCGTTTACAAATATAATCGTCCAGCCCCTGAGATTTCAACAGGCCCGGAATTTTATAGATGGAATAGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACa                                                                                                   <  2:3155913/148‑1 (MQ=255)
                                             tCGGACAGATTCGCTTCCGGGCAGTTTAAGCTGAATCGTTTACAAATATAATCGTCCAGCCCCTGAGATTTCAACAGGCCCGGAATTTTATAGATGGAATAGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACa                                                                                                   <  2:2164760/148‑1 (MQ=255)
                                              cGGACAGATTCGCTTCCGGGCAGTTTAAGCTGAATCGTTTACAAATATAATCGTCCAGCCCCTGAGATTTCAACAGGCCCGGAATTTTATAGATGGAATAGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACaa                                                                                                  >  1:2491554/1‑148 (MQ=255)
                                                 aCAGATTCGCTTCCGGGCAGTTTAAGCTGAATCGTTTACAAATATAATCGTCCAGCCCCTGAGATTTCAACAGGCCCGGAATTTTATAGATGGAATAGACGTCTTTCAGAGAAATAAACGCTTTTTACGGAACATTAAAGAACAAAGc                                                                                               >  1:2074550/1‑148 (MQ=255)
                                                 aCAGATTCGCTTCCGGGCAGTTTAAGCTGAATCGTTTACAAATATAATCGTCCACCCCCTGAGATTTCAACAGGCCCGGAATTTTATAGATGGAATAGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACAAAGc                                                                                               >  1:304894/1‑148 (MQ=255)
                                                    gATTCGCTTCCGGGCAGTTTAAGCTGAATCGTTTACAAATATAATCGTCCAGCCCCTGAGATTTCAACAGGCCCGGAATTTTATAGATGGAATAGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACAATGCAAt                                                                                            <  2:310219/148‑1 (MQ=255)
                                                    gATTCGCTTCCGGGCAGTTTAAGCTGAATCGTTTACAAATATAATCGTCCAGCCCCTGAGATTTCAACAGGCCCGGAATTTTATAGATGGAATAGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACAATGCAAt                                                                                            <  2:545203/148‑1 (MQ=255)
                                                      ttCGCTTCCGGGCAGTTTAAGCTGAATCGTTTACAAATATAATCGTCCAGCCCCTGAGATTTCAACAGGCCCGGAATTTTATAGATGGAATAGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACAATGCAATCt                                                                                          <  1:2408116/148‑1 (MQ=255)
                                                       tCGCTTCCGGGCAGTTTAAGCTGAATCGTTTACAAATATAATCGTCCAGCCCCTGAGATTTCAACAGGCCCGGAATTTTATAGATGGAATAGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACAATGCAATCtt                                                                                         <  2:1486795/148‑1 (MQ=255)
                                                       tCGCTTCCGGGCAGTTTAAGCTGAATCGTTTACAAATATAATCGTCCAGCCCCTGAGATTTCAACAGGCCCGGAATTTTATAGATGGAATAGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACAATGCAATCtt                                                                                         <  1:170449/148‑1 (MQ=255)
                                                       tCGCTTCCGGGCAGTTTAAGCTGAATCGTTTACAAATATAATCGTCCAGCCCCTGAGATTTCAACAGGCCCGGAATTTTATAGATGGAATAGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACAATGCAATCtt                                                                                         <  1:2146827/148‑1 (MQ=255)
                                                       tCGCTTCCGGGCAGTTTAAGCTGAATCGTTTACAAATATAATCGTCCAGCCCCTGAGATTTCAACAGGCCCGGAATTTTATAGATGGAATAGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACAATGCAATCtt                                                                                         >  2:1058532/1‑148 (MQ=255)
                                                         gCTTCCGGGCAGTTTAAGCTGAATCGTTTACAAATATAATCGTCCAGCCCCTGAGATTTCAACAGGCCCGGAATTTTATAGATGGAATAGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACAATGCAATCTTCg                                                                                       <  2:1089571/148‑1 (MQ=255)
                                                         gCTTCCGGGCAGTTTAAGCTGAATCGTTTACAAATATAATCGTCCAGCCCCTGAGATTTCAACAGGCCCGGAATTTTATAGATGGAATAGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACAATGCAATCTTCg                                                                                       <  2:2846506/148‑1 (MQ=255)
                                                              cGGGCAGTTTAAGCTGAATCGTTTACAAATATAATCGTCCAGCCCCTGAGATTTCAACAGGCCCGGAATTTTATAGATGGAATAGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACAATGCAATCTTCGCAcgt                                                                                  <  2:2557802/148‑1 (MQ=255)
                                                                ggCAGTTTAAGCTGAATCGTTTACAAATATAATCGTCCAGCCCCTGAGATTTCAACAGGCCCGGAATTTTATAGATGGAATAGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACAATGCAATCTTCGCAcgttc                                                                                >  2:1776727/1‑148 (MQ=255)
                                                                      ttAAGCTGAATCGTTTACAAATATAATCGTCCAGCCCCTGAGATTTCAACAGGCCCGGAATTTTATAGATGGAATAGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACAATGCAATCTTCGCACGTTCGTTCGc                                                                          <  2:3500351/148‑1 (MQ=255)
                                                                      ttAAGCTGAATCGTTTACAAATATAATCGTCCAGCCCCTGAGATTTCAACAGGCCCGGAATTTTATAGATGGAATAGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACAATGCAATCTTCGCACGTTCGTTCGc                                                                          <  1:1188320/148‑1 (MQ=255)
                                                                         aGCTGAATCGTTTACAAATATAATCGTCCAGCCCCTGAGATTTCAACAGGCCCGGAATTTTATAGATGGAATAGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACAATGCAATCTTCGCACGTTCGTTCGCCgg                                                                       >  1:2169310/1‑148 (MQ=255)
                                                                          gCTGAATCGTTTACAAATATAATCGTCCAGCCCCTGAGATTTCAACAGGCCCGGAATTTTATAGATGGAATAGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACAATGCAATCTTCGCACGTTCGTTCGCCGGa                                                                      >  2:1453582/1‑148 (MQ=255)
                                                                             gAATCGTTTACAAATATAATCGTCCAGCCCCTGAGATTTCAACAGGCCCGGAATTTTATAGATGGAATAGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACAATGCAATCTTCGCACGTTCGTTCGCCGGAACa                                                                   <  1:3202129/148‑1 (MQ=255)
                                                                             gAATCGTTTACAAATATAATCGTCCAGCCCCTGAGATTTCAACAGGCCCGGAATTTTATAGATGGAATAGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACAATGCAATCTTCGCACGTTCGTTCGCCGGAACa                                                                   <  1:158985/148‑1 (MQ=255)
                                                                              aaTCGTTTACAAATATAATCGTCCAGCCCCTGAGATTTCAACAGGCCCGGAATTTTATAGATGGAATAGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACAATGCAATCTTCGCACGTTCGTTCGCCGGAACAg                                                                  >  2:644535/1‑148 (MQ=255)
                                                                                 cGTTTACAAATATAATCGTCCAGCCCCTGAGATTTCAACAGCCCCGGAATTTTATAGATGGAATAGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACAATGCAATCTTCGCACGTTCGTTCGCCGGAACAgcgc                                                               >  2:3371125/1‑148 (MQ=255)
                                                                                  gTTTACAAATATAATCGTCCAGCCCCTGAGATTTCAACAGGCCCGGAATTTTATAGATGGAATCGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACAATGCAATCTTCGCACGTTCGTTCGCCGGAACAGcgcg                                                              <  2:1720650/148‑1 (MQ=255)
                                                                                    ttACAAATATAATCGTCCAGCCCCTGAGATTTCAACAGGCCCGGAATTTTATAGATGGAATAGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACAATGCAATCTTCGCACGTTCGTTCGCCGGAACAGCGCGAt                                                            >  1:2970962/1‑148 (MQ=255)
                                                                                       cAAATATAATCGTCCAGCCCCTGAGATTTCAACAGGCCCGGAATTTTATAGATGGAATAGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACAATGCAATCTTCGCACGTTCGTTCGCCGGAACAGCGCGATCTg                                                         <  1:2074985/148‑1 (MQ=255)
                                                                                        aaaTATACTCGTCCAGCCCCTGAGATTTCAACAGGCCCGGAATTTTATAGATGGAATAGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACAATGCAATCTTCGCACGTTCGTTCGCCGGAACAGCGCGATCTGa                                                        <  1:2983582/148‑1 (MQ=255)
                                                                                        aaaTATAATCGTCCAGCCCCTGAGATTTCAACAGGCCCGGAATTTTATAGATGGAATAGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACAATGCAATCTTCGCACGTTCGTTCGCCGGAACAGCGCGATCTGa                                                        <  1:2698270/148‑1 (MQ=255)
                                                                                        aaaTATAATCGTCCAGCCCCTGAGATTTCAACAGGCCCGGAATTTTATAGATGGAATAGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACAATGCAATCTTCGCACGTTCGTTCGCCGGAACAGCGCGATCTGa                                                        >  1:100745/1‑148 (MQ=255)
                                                                                            ataATCGTCCAGCCCCTGAGATTTCAACAGGCCCGGAATTTTATAGATGGAATAGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACAATGCAATCTTCGCACGTTCGTTCGCCGGAACAGCGCGATCTGAACGa                                                    >  2:1393068/1‑148 (MQ=255)
                                                                                              aaTCGTCCAGCCCCTGAGATTTCAACAGGCCCGGAATTTTATAGATGGAATAGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACAATGCAATCTTCGCACGTTCGTTCGCCGGAACAGCGCGATCTGAACGACa                                                  >  1:690976/1‑148 (MQ=255)
                                                                                                tCGTCCAGCCCCTGAGATTTCAACAGGCCCGGAATTTTATAGATGGAATAGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACAATGCAATCTTCGCACGTTCGTTCGCCGGAACAGCGCGATCTGAACGACaaa                                                >  1:368002/1‑148 (MQ=255)
                                                                                                 cGTCCAGCCCCTGAGATTTCAACAGGCCCGGAATTTTATAGATGGAATAGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACAATGCAATCTTCGCACGTTCGTTCGCCGGAACAGCGCGATCTGAACGACAAAt                                               >  1:1422214/1‑148 (MQ=255)
                                                                                                      aGCCCCTGAGATTTCAACAGGCCCGGAATTTTATAGATGGAATAGACGTCTTTCAGAGACATAACCGCTTTTTCCGGAACATAACAGAACAATGCAATCTTCGCACGTTCGTTCGCCGGAACCGCGCGATCTGAACCACAAATCAGGa                                          >  2:644046/1‑148 (MQ=255)
                                                                                                      aGCCCCTGAGATTTCAACAGGCCCGGAATTTTATAGATGGAATAGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACAATGCAATCTTCGCACGTTCGTTCGCCGGAACAGCGCGATCTGAACGACAAATCAGGa                                          >  2:2403421/1‑148 (MQ=255)
                                                                                                        ccccTGAGATTTCAACAGGCCCGGAATTTTATAGATGGAATCGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACAATGCAATCTTCGCACGTTCGTTCGCCGGAACAGCGCGATCTGAACGACAAATCAGGATg                                        <  1:3335901/148‑1 (MQ=255)
                                                                                                         cccTGAGATTTCAACAGGCCCGGAATTTTATAGATGGAATAGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACAATGCAATCTTCGCACGTTCGTTCGCCGGAACAGCGCGATCTGAACGACAAATCAGGATGt                                       <  2:1435436/148‑1 (MQ=255)
                                                                                                          ccTGAGATTTCAACAGGCCCGGAATTTTATAGATGGAATAGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACAATGCAATCTTCGCACGTTCGTTCGCCGGAACAGCGCGATCTGAACGACAAATCAGGATGTc                                      >  2:3217296/1‑148 (MQ=255)
                                                                                                           cTGAGATTTCAACAGGCCCGGAATTTTATAGATGGAATAGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACAATGCAATCTTCGCACGTTCGTTCGCCGGAACAGCGCGATCTGAACGACAAATCAGGATGTCa                                     >  1:774040/1‑148 (MQ=255)
                                                                                                           cTGAGATTTCAACAGGCCCGGAATTTTATAGATGGAATAGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACAATGCAATCTTCGCACGTTCGTTCGCCGGAACAGCGCGATCTGAACGACAAATCAGGATGTCa                                     >  2:3413983/1‑148 (MQ=255)
                                                                                                                aTTTCAACAGGCCCGGAATTTTATAGATGGAATAGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACAATGCAATCTTCGCACGTTCGTTCGCCGGAACAGCGCGATCTGAACGACAAATCAGGATGTCAGGCTg                                <  1:2912064/148‑1 (MQ=255)
                                                                                                                  ttCAACAGGCCCGGAATTTTATAGATGGAATAGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACAATGCAATCTTCGCACGTTCGTTCGCCGGAACAGCGCGATCTGAACGACAAATCAGGATGTCAGGCTGGa                              >  2:896100/1‑148 (MQ=255)
                                                                                                                          gCCCGGAATTTTATAGATGGAATAGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACAATGCAATCTTCGCACGTTCGTTCGCCGGAACAGCGCGATCTGAACGACAAATCAGGATGTCAGGCTGGATACCGATg                      >  2:3005584/1‑148 (MQ=255)
                                                                                                                          gCCCGGAATTTTATAGATGGAATAGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACAATGCAATCTTCGCACGTTCGTTCGCCGGAACAGCGCGATCTGAACGACAAATCAGGATGTCAGGCTGGATACCGATg                      >  2:3005730/1‑148 (MQ=255)
                                                                                                                           cccGGAATTTTATAGATGGAATCGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACAATGCAATCTTCGCACGTTCGTTCGCCGGAACAGCGCGATCTGAACGACAAATCAGGATGTCAGGCTGGATACCGATgg                     <  2:184766/148‑1 (MQ=255)
                                                                                                                           cccGGAATTTTATAGATGGAATAGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACAATGCAATCTTCGCACGTTCGTTCGCCGGAACAGCGCGATCTGAACGACAAATCAGGATGTCAGGCTGGATACCGATgg                     <  1:1820612/148‑1 (MQ=255)
                                                                                                                           cccGGAATTTTATAGATGGAATAGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACAATGCAATCTTCGCACGTTCGTTCGCCGGAACAGCGCGATCTGAACGACAAATCAGGATGTCAGGCTGGATACCGATgg                     <  1:2308026/148‑1 (MQ=255)
                                                                                                                              ggAATTTTATAGATGGAATAGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACAATGCAATCTTCGCACGTTCGTTCGCCGGAACAGCGCGATCTGAACGACAAATCAGGATGTCAGGCTGGATACCGATGgaga                  >  2:1949868/1‑148 (MQ=255)
                                                                                                                                aaTTTTATAGATGGAATAGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACAATGCAATCTTCGCACGTTCGTTCGCCGGAACAGCGCGATCTGAACGACAAATCAGGATGTCAGGCTGGATACCGATGGAGagc                >  1:980400/1‑148 (MQ=255)
                                                                                                                                  ttttATAGATGGAATAGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACAATGCAATCTTCGCACGTTCGTTCGCCGGAACAGCGCGATCTGAACGACAAATCAGGATGTCAGGCTGGATACCGATGGAGagca               <  1:89952/147‑1 (MQ=255)
                                                                                                                                   tttATAGATGGAATAGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACAATGCAATCTTCGCACGTTCGTTCGCCGGAACAGCGCGATCTGAACGACAAATCAGGATGTCAGGCTGGATACCGATGGAGagcagc             <  2:2082504/148‑1 (MQ=255)
                                                                                                                                    ttATAGATGGAATAGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACAATGCAATCTTCGCACGTTCGTTCGCCGGAACAGCGCGATCTGAACGACAAATCAGGATGTCAGGCTGGATACCGATGGAGAGCAGct            <  1:2442645/148‑1 (MQ=255)
                                                                                                                                    ttATAGATGGAATAGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACAATGCAATCTTCGCACGTTCGTTCGCCGGAACAGCGCGATCTGAACGACAAATCAGGATGTCAGGCTGGATACCGATGGAGAGCAGct            >  1:426622/1‑148 (MQ=255)
                                                                                                                                     tataGATGGAATAGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACAATGCAATCTTCGCACGTTCGTTCGCCGGAACAGCGCGATCTGAACGACAAATCAGGATGTCAGGCTGGATACCGATGGAGAGCAGctc           >  2:905927/1‑148 (MQ=255)
                                                                                                                                      ataGATGGAATAGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACAATGCAATCTTCGCACGTTCGTTCGCCGGAACAGCGCGATCTGAACGACAAATCAGGATGTCAGGCTGGATACCGATGGAGAGCAGctct          >  2:3023420/1‑148 (MQ=255)
                                                                                                                                      ataGATGGAATAGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACAATGCAATCTTCGCACGTTCGTTCGCCGGAACAGCGCGATCTGAACGACAAATCAGGATGTCAGGCTGGATACCGATGGAGAGCAGctct          <  2:945949/148‑1 (MQ=255)
                                                                                                                                       taGATGGAATAGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACAATGCAATCTTCGCACGTTCGTTCGCCGGAACAGCGCGATCTGAACGACAAATCAGGATGTCAGGCTGGATACCGATGGAGAGCAGCTCtt         <  2:2289955/148‑1 (MQ=255)
                                                                                                                                        aGATGGAATAGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACAATGCAATCTTCGCACGTTCGTTCGCCGGAACAGCGCGATCTGAACGACAAATCAGGATGTCAGGCTGGATACCGATGGAGAGCAGCTCttt        >  2:404193/1‑148 (MQ=255)
                                                                                                                                          aTGGAATAGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACAATGCAATCTTCGCACGTTCGTTCGCCGGAACAGCGCGATCTGAACGACAAATCAGGATGTCAGGCTGGATACCGATGGAGAGCAGCTCTTTTa      >  1:1082964/1‑148 (MQ=255)
                                                                                                                                          aTGGAATAGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACAATGCAATCTTCGCACGTTCGTTCGCCGGAACAGCGCGATCTGAACGACAAATCAGGATGTCAGGCTGGATACCGATGGAGAGCAGCTCTTTTa      <  1:2524964/148‑1 (MQ=255)
                                                                                                                                           tGGAATAGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACAATGCAATCTTCGCACGTTCGTTCGCCGGAACAGCGCGATCTGAACGACAAATCAGGATGTCAGGCTGGATACCGATGGAGAGCAGCTCTTTTAc     >  1:2305524/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
ACTTCACTTACCGGGTTCGCTTCTTCGAAGATAACCTGTTCCCATTCGGACAGATTCGCTTCCGGGCAGTTTAAGCTGAATCGTTTACAAATATAATCGTCCAGCCCCTGAGATTTCAACAGGCCCGGAATTTTATAGATGGAATCGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACAATGCAATCTTCGCACGTTCGTTCGCCGGAACAGCGCGATCTGAACGACAAATCAGGATGTCAGGCTGGATACCGATGGAGAGCAGCTCTTTTACAGA  >  CP009273/2902157‑2902446

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: