Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,020,049 C→G 5.2% P55A (CCG→GCG)  fre → NAD(P)H‑flavin reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,020,0490CG 5.2% 727.9 / 23.4 191P55A (CCG→GCG) freNAD(P)H‑flavin reductase
Reads supporting (aligned to +/- strand):  ref base C (87/94);  new base G (4/6);  total (91/100)
Fisher's exact test for biased strand distribution p-value = 7.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

Only 100 of 191 total aligned reads displayed.

TAAAGTGACCTCGGTAGAAGCTATCACGGATACCGTATATCGTGTCCGCATCGTGCCAGACGCGGCCTTTTCTTTTCGTGCTGGTCAGTATTTGATGGTAGTGATGGATGAGCGCGACAAACGTCCGTTCTCAATGGCTTCGACGCCGGATGAAAAAGGGTTTATCGAGCTGCATATTGGCGCTTCTGAAATCAACCTTTACGCGAAAGCAGTCATGGACCGCATCCTCAAAGATCATCAAATCGTGGTCGACATTCCCCACGGAGAAGCGTGGCTGCGCGATGATGAAGAG  >  CP009273/4019904‑4020195
                                                                                                                                                 |                                                                                                                                                  
tAAAGTGACCTCGGTAGAAGCTATCACGGATACCGTATATCGTGTCCGCATCGTGCCAGACGCGGCCTTTTCTTTTCGTGCTGGTCAGTATTTGATGGTAGTGATGGATGAGCGCGACAAACGTCCGTTCTCAATGGCTTCGACGCCg                                                                                                                                                  >  1:2381155/1‑148 (MQ=255)
  aaGTGACCTCGGTAGAAGCTATCACGGATACCGTATATCGTGTCCGCATCGTGCCAGACGCGGCCTTTTCTTTTCGTGCTGGTCAGTATTTGATGGTAGTGATGGATGAGCGCGACAAACGTCCGTTCTCAATGGCTTCGACGCCGGa                                                                                                                                                <  1:113325/148‑1 (MQ=255)
   aGTGACCTCGGTAGAAGCTATCACGGATACCGTATATCGTGTCCGCATCGTGCCAGACGCGGCCTTTTCTTTTCGTGCTGGTCAGTATTTGATGGTAGTGATGGATGAGCGCGACAAACGTCCGTTCTCAATGGCTTCGACGCCGGAt                                                                                                                                               <  2:342664/148‑1 (MQ=255)
        ccTCGGTAGAAGCTATCACGGATACCGTATATCGTGTCCGCATCGTGCCAGACGCGGCCTTTTCTTTTCGTGCTGGTCAGTATTTGATGGTAGTGATGGATGAGCGCGACAAACGTCCGTTCTCAATGGCTTCGACGCCGGATGaaaa                                                                                                                                          <  2:962993/148‑1 (MQ=255)
        ccTCGGTAGAAGCTATCACGGATACCGTATATCGTGTCCGCATCGTGCCAGACGCGGCCTTTTCTTTTCGTGCTGGTCAGTATTTGATGGTAGTGATGGATGAGCGCGACAAACGTCCGTTCTCAATGGCTTCGACGCCGGATGaaaa                                                                                                                                          >  1:3236179/1‑148 (MQ=255)
             gTAGAAGCTATCACGGATACCGTATATCGTGTCCGCATCGTGCCAGACGCGGCCTTTTCTTTTCGTGCTGGTCAGTATTTGATGGTAGTGATGGATGAGCGCGACAAACGTCCGTTCTCAATGGCTTCGACGCCGGATGAAAAAGGGt                                                                                                                                     >  2:2317804/1‑148 (MQ=255)
                gAAGCTATCACGGATACCGTATATCGTGTCCGCATCGTGCCAGACGCGGCCTTTTCTTTTCGTGCTGGTCAGTATTTGATGGTAGTGATGGATGAGCGCGACAAACGTCCGTTCTCAATGGCTTCGACGCCGGATGAAAAAGGGTTTa                                                                                                                                  >  2:2875120/1‑148 (MQ=255)
                   gCTATCACGGATACCGTATATCGTGTCCGCATCGTGCCAGACGCGGCCTTTTCTTTTCGTGCTGGTCAGTATTTGATGGTAGTGATGGATGAGCGCGACAAACGTCCGTTCTCAATGGCTTCGACGCCGGATGAAAAAGGGTTTATCg                                                                                                                               >  2:2406447/1‑148 (MQ=255)
                     tATCACGGATACCGTATATCGTGTCCGCATCGTGCCAGACGCGGCCTTTTCTTTTCGTGCTGGTCAGTATTTGATGGTAGTGATGGATGAGCGCGACAAACGTCCGTTCTCAATGGCTTCGACGCCGGATGAAAAAGGGTTTATCGAg                                                                                                                             <  1:1542554/148‑1 (MQ=255)
                      aTCACGGATACCGTATATCGTGTCCGCATCGTGCCAGACGCGGCCTTTTCTTTTCGTGCTGGTCAGTATTTGATGGTAGTGATGGATGAGCGCGACAAACGTCCGTTCTCAATGGCTTCGACGCCGGATGAAAAAGGGTTTATCGAGc                                                                                                                            <  2:389192/148‑1 (MQ=255)
                        cACGGATACCGTATATCGTGTCCGCATCGTGCCAGACGCGGCCTTTTCTTTTCGTGCTGGTCAGTATTTGATGGTAGTGATGGATGAGCGCGACAAACGTCCGTTCTCAATGGCTTCGACGCCGGATGAAAAAGGGTTTATCGAGCTg                                                                                                                          >  1:2588360/1‑148 (MQ=255)
                        cACGGATACCGTATATCGTGTCCGCATCGTGCCAGACGCGGCCTTTTCTTTTCGTGCTGGTCAGTATTTGATGGTAGTGATGGATGAGCGCGACAAACGTCCGTTCTCAATGGCTTCGACGCCGGATGAAAAAGGGTTTATCGAGCTg                                                                                                                          >  1:1485130/1‑148 (MQ=255)
                           ggATACCGTATATCGTGTCCGCATCGTGCCAGACGCGGCCTTTTCTTTTCGTGCTGGTCAGTATTTGATGGTAGTGATGGATGAGCGCGACAAACGTCCGTTCTCAATGGCTTCGACGGCGGATGAAAAAGGGTTTATCGAGCTGCat                                                                                                                       >  1:1936366/1‑148 (MQ=255)
                           ggATACCGTATATCGTGTCCGCATCGTGCCAGACGCGGCCTTTTCTTTTCGTGCTGGTCAGTATTTGATGGTAGTGATGGATGAGCGCGACAAACGTCCGTTCTCAATGGCTTCGACGGCGGATGAAAAAGGGTTTATCGAGCTGCat                                                                                                                       >  1:810526/1‑148 (MQ=255)
                           ggATACCGTATATCGTGTCCGCATCGTGCCAGACGCGGCCTTTTCTTTTCGTGCTGGTCAGTATTTGATGGTAGTGATGGATGAGCGCGACAAACGTCCGTTCTCAATGGCTTCGACGCCGGATGAAAAAGGGTTTATCGAGCTGCat                                                                                                                       <  1:3293055/148‑1 (MQ=255)
                            gATACCGTATATCGTGTCCGCATCGTGCCAGACGCGGCCTTTTCTTTTCGTGCTGGTCAGTATTTGATGGTAGTGATGGATGAGCGCGACAAACGTCCGTTCTCAATGGCTTCGACGCCGGATGAAAAAGGGTTTATCGAGCTGCata                                                                                                                      <  2:1992683/148‑1 (MQ=255)
                            gATACCGTATATCGTGTCCGCATCGTGCCAGACGCGGCCTTTTCTTTTCGTGCTGGTCAGTATTTGATGGTAGTGATGGATGAGCGCGACAAACGTCCGTTCTCAATGGCTTCGACGCCGGATGAAAAAGGGTTTATCGAGCTGCata                                                                                                                      <  1:749519/148‑1 (MQ=255)
                             aTACCGTATATCGTGTCCGCATCGTGCCAGACGCGGCCTTTTCTTTTCGTGCTGGTCAGTATTTGATGGTAGTGATGGATGAGCGCGACAAACGTCCGTTCTCAATGGCTTCGACGCCGGATGAAAAAGGGTTTATCGAGCTGCatat                                                                                                                     >  2:2671176/1‑148 (MQ=255)
                               aCCGTATATCGTGTCCGCATCGTGCCAGACGCGGCCTTTTCTTTTCGTGCTGGTCAGTATTTGATGGTAGTGATGGATGAGCGCGACAAACGTCCGTTCTCAATGGCTTCGACGCCGGATGAAAAAGGGTTTATCGAGCTGCATATTg                                                                                                                   >  2:2073412/1‑148 (MQ=255)
                                  gTATATCGTGTCCGCATCGTGCCAGACGCGGCCTTTTCTTTTCGTGCTGGTCAGTATTTGATGGTAGTGATGGATGAGCGCGACAAACGTCCGTTCTCAATGGCTTCGACGCCGGATGAAAAAGGGTTTATCGAGCTGCATATTGgcg                                                                                                                <  2:1460339/148‑1 (MQ=255)
                                   tataTCGTGTCCGCATCGTGCCAGACGCGGCCTTTTCTTTTCGTGCTGGTCAGTATTTGATGGTAGTGATGGATGAGCGCGACAAACGTCCGTTCTCAATGGCTTCGACGCCGGATGAAAAAGGGTTTATCGAGCTGCATATTGgcgc                                                                                                               <  2:2043263/148‑1 (MQ=255)
                                   tataTCGTGTCCGCATCGTGCCAGACGCGGCCTTTTCTTTTCGTGCTGGTCAGTATTTGATGGTAGTGATGGATGAGCGCGACAAACGTCCGTTCTCAATGGCTTCGACGCCGGATGAAAAAGGGTTTATCGAGCTGCATATTGgcgc                                                                                                               >  2:1761671/1‑148 (MQ=255)
                                      atCGTGTCCGCATCGTGCCAGACGCGGCCTTTTCTTTTCGTGCTGGTCAGTATTTGATGGTAGTGATGGATGAGCGCGACAAACGTCCGTTCTCAATGGCTTCGACGCCGGATGAAAAAGGGTTTATCGAGCTGCATATTGGCGCTTc                                                                                                            <  1:1892784/148‑1 (MQ=255)
                                         gtgtCCGCATCGTGCCAGACGCGGCCTTTTCTTTTCGTGCTGGTCAGTATTTGATGGTAGTGATGGATGAGCGCGACAAACGTCCGTTCTCAATGGCTTCGACGCCGGATGAAAAAGGGTTTATCGAGCTGCATATTGGCGCTTCTGa                                                                                                         <  2:141371/148‑1 (MQ=255)
                                             ccGCATCGTGCCAGACGCGGCCTTTTCTTTTCGTGCTGGTCAGTATTTGATTGTAGTGATGGATGAGCGCGACAAACGTCCGTTCTCAATGGCTTCGACGCCGGATGAAAAAGGGTTTATAGAGCTGCATATTTGCGCTTATGtaaac                                                                                                     >  1:1344135/1‑146 (MQ=255)
                                             ccGCATCGTGCCAGACGCGGCCTTTTCTTTTCGTGCTGGTCAGTATTTGATGGTAGTGATGTATTAGCGCGAAAAAAGTCCGTTATCAATGGCTTCGACGCCGGATGAAAAATGGTTTATAGAGCTGCATATTTGCGCTTCTGACATc                                                                                                     >  1:2237271/1‑148 (MQ=255)
                                             ccGCATCGTGCCAGACGCGGCCTTTTCTTTTCGTGCTGGTCAGTATTTGATGGTAGTGATGGATGAGCGCGACAAACGTCCGTTCTCAATGGCTTCGACGCCGGATGAAAAAGGGTTTATCGAGCTGCATATTGGCGCTTCTGAAATc                                                                                                     >  1:442129/1‑148 (MQ=255)
                                              cGCATCGTGCCAGACGCGGCCTTTTCTTTTCGTGCTGGTCAGTATTTGATGGTAGTGATGGATGAGCGCGACAAACGTCCGTTCTCAATGGCTTCGACGCCGGATGAAAAAGGGTTTATCGAGCTGCATATTGGCGCTTCTGAAATCa                                                                                                    <  1:2662779/148‑1 (MQ=255)
                                              cGCATCGTGCCAGACGCGGCCTTTTCTTTTCGTGCTGGTCAGTATTTGATGGTAGTGATGGATGAGCGCGACAAACGTCCGTTCTCAATGGCTTCGACGCCGGATGAAAAAGGGTTTATCGAGCTGCATATTGGCGCTTCTGAAATCa                                                                                                    <  2:2862491/148‑1 (MQ=255)
                                               gCATCGTGCCAGACGCGGCCTTTTCTTTTCGTGCTGGTCAGTATTTGATGGTAGTGATGGATGAGCGCGACAAACGTCCGTTCTCAATGGCTTCGACGCCGGATGAAAAAGGGTTTATCGAGCTGCATATTGGCGCTTCTGAAATCaa                                                                                                   <  2:1337277/148‑1 (MQ=255)
                                               gCATCGTGCCAGACGCGGCCTTTTCTTTTCGGGCTGGTCAGTATTTGATGGTAGTGATGGATGAGCGCGACAAACGTCCGTTCTCAATGGCTTCGACGCCGGATGAAAAAGGGTTTATCGAGCTGCATATTGGCGCTTCTGAAATCaa                                                                                                   >  1:1029205/1‑148 (MQ=255)
                                                cATCGTGCCAGACGCGGCCTTTTCTTTTCGTGCTGGTCAGTATTTGATGGTAGTGATGGATGAGCGCGACAAACGTCCGTTCTCAATGGCTTCGACGCCGGATGAAAAAGGGTTTATCGAGCTGCATATTGGCGCTTCTGAAATCAAc                                                                                                  >  2:2478391/1‑148 (MQ=255)
                                                   cGTGCCAGACGCGGCCTTTTCTTTTCGTGCTGGTCAGTATTTGATGGTAGTGATGGATGAGCGCGACAAACGTCCGTTCTCAATGGCTTCGACGCCGGATGAAAAAGGGTTTATCGAGCTGCCTATTGGCGCTTCTGAAATCAACCtt                                                                                               >  1:3552479/1‑148 (MQ=255)
                                                   cGTGCCAGACGCGGCCTTTTCTTTTCGTGCTGGTCAGTATTTGATGGTAGTGATGGATGAGCGCGACAAACGTCCGTTCTCAATGGCTTCGACGCCGGATGAAAAAGGGTTTATCGAGCTGCATATTGGCGCTTCTGAAATCAACCtt                                                                                               >  2:1507871/1‑148 (MQ=255)
                                                   cGTGCCAGACGCGGCCTTTTCTTTTCGTGCTGGTCAGTATTTGATGGTAGTGATGGATGAGCGCGACAAACGTCCGTTCTCAATGGCTTCGACGCCGGATGAAAAAGGGTTTATCGAGCTGCATATTGGCGCTTCTGAAATCAACCtt                                                                                               <  2:1760658/148‑1 (MQ=255)
                                                    gTGCCAGACGCGGCCTTTTCTTTTCGTGCTGGTCAGTATTTGATGGTAGTGATGGATGAGCGCGACAAACGTCCGTTCTCAATGGCTTCGACGCCGGATGAAAAAGGGTTTATCGAGCTGCATATTGGCGCTTCTGAAATCAACCttt                                                                                              <  1:913165/148‑1 (MQ=255)
                                                     tGCCAGACGCGGCCTTTTCTTTTCGTGCTGGTCAGTATTTGATGGTAGTGATGGATGAGCGCGACAAACGTCCGTTCTCAATGGCTTCGACGCCGGATGAAAAAGGGTTTATCGAGCTGCATATTGGCGCTTCTGAAATCAACCTTTa                                                                                             <  2:2917929/148‑1 (MQ=255)
                                                       ccAGACGCGGCCTTTTCTTTTCGTGCTGGTCAGTATTTGATGGTAGTGATGGATGAGCGCGACAAACGTCCGTTCTCAATGGCTTCGACGCCGGATGAAAAAGGGTTTATCGAGCTGCATATTGGCGCTTCTGAAATCAACCTTTAcg                                                                                           >  2:171736/1‑148 (MQ=255)
                                                           aCGCGGCCTTTTCTTTTCGTGCTGGTCAGTATTTGATGGTAGTGATGGATGAGCGCGACAAACGTCCGTTCTCAATGGCTTCGACGCCGGATGAAAAAGGGTTTATCGAGCTGCATATTGGCGCTTCTGAAATCAACCTTTACGCGaa                                                                                       <  1:1424587/148‑1 (MQ=255)
                                                              cgGCCTTTTCTTTTCGTGCTGGTCAGTATTTGATGGTAGTGATGGATGAGCGCGACAAACGTCCGTTCTCAATGGCTTCGACGCCGGATGAAAAAGGGTTTATCGAGCTGCATATTGGCGCTTCTGAAATCAACCTTTACGCGAAAGc                                                                                    <  2:1395314/148‑1 (MQ=255)
                                                                 ccTTTTCTTTTCGTGCTGGTCAGTATTTGATGGTAGTGATGGATGAGCGCGACAAACGTCCGTTCTCAATGGCTTCGACGCCGGATGAAAAAGGGTTTATCGAGCTGCATATTGGCGCTTCTGAAATCAACCTTTACGCGAAAGCAGt                                                                                 >  2:2768890/1‑148 (MQ=255)
                                                                    tttcttttcGTGCTGGTCAGTATTTGATGGTAGTGATGGATGAGCGCGACAAACGTCCGTTCTCAATGGCTTCGACGCCGGATGAAAAAGGGTTTATCGAGCTGCATATTGGCGCTTCTGAAATCAACCTTTACGCGAAAGCAGTCAt                                                                              >  2:1634194/1‑148 (MQ=255)
                                                                      tcttttcGTGCTGGTCAGTATTTGATGGTAGTGATGGATGAGCGCGACAAACGTCCGTTCTCAATGGCTTCGACGCCGGATGAAAAAGGGTTTATCGAGCTGCATATTGGCGCTTCTGAAATCAACCTTTACGCGAAAGCAGTCATgg                                                                            <  1:1855382/148‑1 (MQ=255)
                                                                       cttttcGTGCTGGTCAGTATTTGATGGTAGTGATGGATGAGCGCGACAAACGTCCGTTCTCAATGGCTTCGACGCCGGATGAAAAAGGGTTTATCGAGCTGCATATTGGCGCTTCTGAAATCAACCTTTACGCGAAAGCAGTCATGGa                                                                           >  2:3286214/1‑148 (MQ=255)
                                                                          ttcGTGCTGGTCAGTATTTGATGGTAGTGATGGATGAGCGCGACAAACGTCCGTTCTCAATGGCTTCGACGCCGGATGAAAAAGGGTTTATCGAGCTGCATATTGGCGCTTCTGAAATCAACCTTTACGCGAAAGCAGTCATGGACCg                                                                        <  2:2573703/148‑1 (MQ=255)
                                                                           tcGTGCTGGTCAGTATTTGATGGTAGTGATGGATGAGCGCGACAAACGTCCGTTCTCAATGGCTTCGACGCCGGATGAAAAAGGGTTTATCGAGCTGCATATTGGCGCTTCTGAAATCAACCTTTACGCGAAAGCAGTCATGGACCGc                                                                       <  2:1870203/148‑1 (MQ=255)
                                                                             gTGCTGGTCAGTATTTGATGGTAGTGATGGATGAGCGCGACAAACGTCCGTTCTCAATGGCTTCGACGCCGGATGAAAAAGGGTTTATCGAGCTGCATATTGGCGCTTCTGAAATCAACCTTTACGCGAAAGCAGTCATGGACCGCAt                                                                     <  1:1028081/148‑1 (MQ=255)
                                                                                   gTCAGTATTTGATGGTAGTGATGGATGAGCGCGACAAACGTCCGTTCTCAATGGCTTCGACGCCGGATGAAAAAGGGTTTATCGAGCTGCATATTGGCGCTTCTGAAATCAACCTTTACGCGAAAGCAGTCATGGACCGCATCCTCaa                                                               >  1:1862082/1‑148 (MQ=255)
                                                                                   gTCAGTATTTGATGGTAGTGATGGATGAGCGCGACAAACGTCCGTTCTCAATGGCTTCGACGCCGGATGAAAAAGGGTTTATCGAGCTGCATATTGGCGCTTCTGAAATCAACCTTTACGCGAAAGCAGTCATGGACCGCATCCTCaa                                                               >  1:2491096/1‑148 (MQ=255)
                                                                                    tCAGTATTTGATGGTAGTGATGGATGAGCGCGACAAACGTCCGTTCTCAATGGCTTCGACGCCGGATGAAAAAGGGTTTATCGAGCTGCATATTGGCGCTTCTGAAATCAACCTTTACGCGAAAGCAGTCATGGACCGCATCCTCaaa                                                              <  1:825710/148‑1 (MQ=255)
                                                                                        tATTTGATGGTAGTGATGGATGAGCGCGACAAACGTCCGTTCTCAATGGCTTCGACGCCGGATGAAAAAGGGTTTATCGAGCTGCATATTGGCGCTTCTGAAATCAACCTTTACGCGAAAGCAGTCATGGACCGCATCCTCAAAGatc                                                          <  1:520926/148‑1 (MQ=255)
                                                                                         aTTTGATGGTAGTGATGGATGAGCGCGACAAACGTCCGTTCTCAATGGCTTCGACGCCGGATGAAAAAGGGTTTATCGAGCTGCATATTGGCGCTTCTGAAATCAACCTTTACGCGAAAGCAGTCATGGACCGCATCCTCAAAGatca                                                         <  1:1119279/148‑1 (MQ=255)
                                                                                          tttGATGGTAGTGATGGATGAGCGCGACAAACGTCCGTTCTCAATGGCTTCGACGCCGGATGAAAAAGGGTTTATCGAGCTGCATATTGGCGCTTCTGAAATCAACCTTTACGCGAAAGCAGTCATGGACCGCATCCTCAAAGatcat                                                        >  2:2236138/1‑148 (MQ=255)
                                                                                          tttGATGGTAGTGATGGATGAGCGCGACAAACGTCCGTTCTCAATGGCTTCGACGCCGGATGAAAAAGGGTTTATCGAGCTGCATATTGGCGCTTCTGAAATCAACCTTTACGCGAAAGCAGTCATGGACCGCATCCTCAAAGatcat                                                        <  2:1525206/148‑1 (MQ=255)
                                                                                              aTGGTAGTGATTGATGAGCGCGACAAACGTCCGTTCTCAATGGCTTCGACGCCGGATGAAAAAGGGTTTATCGAGCTGCATATTGGCGCTTCTGAAATCAACCTTTACGCGAAAGCAGTCATGGACCGCATCCTCAAAGATCATCaaa                                                    <  2:704547/148‑1 (MQ=255)
                                                                                              aTGGTAGTGATGGATGAGCGCGACAAACGTCCGTTCTCAATGGCTTCGACGGCGGATGAAAAAGGGTTTATCGAGCTGCATATTGGCGCTTCTGAAATCAACCTTTACGCGAAAGCAGTCATGGACCGCATCCTCAAAGATCATCaaa                                                    <  2:810526/148‑1 (MQ=255)
                                                                                              aTGGTAGTGATGGATGAGCGCGACAAACGTCCGTTCTCAATGGCTTCGACGGCGGATGAAAAAGGGTTTATCGAGCTGCATATTGGCGCTTCTGAAATCAACCTTTACGCGAAAGCAGTCATGGACCGCATCCTCAAAGATCATCaaa                                                    <  2:1936366/148‑1 (MQ=255)
                                                                                              aTGGTAGTGATGGATGAGCGCGACAAACGTCCGTTCTCAATGGCTTCGACGCCGGATGAAAAAGGGTTTATCGAGCTGCATATTGGCGCTTCTGAAATCAACCTTTACGCGAAAGCAGTCATGGACCGCATCCTCAAAGATCATCaaa                                                    <  2:1720336/148‑1 (MQ=255)
                                                                                                 gTAGTGATGGATGAGCGCGACAAACGTCCGTTCTCAATGGCTTCGACGCCGGATGAAAAAGGGTTTATCGAGCTGCATATTGGCGCTTCTGAAATCAACCTTTACGCGAAAGCAGTCATGGACCGCATCCTCAAAGATCATCAAATCg                                                 <  1:211873/148‑1 (MQ=255)
                                                                                                  tAGTGATGGATGAGCGCGACAAACGTCCGTTCTCAATGGCTTCGACGCCGGATGAAAAAGGGTTTATCGAGCTGCATATTGGCGCTTCTGAAATCAACCTTTACGCGAAAGCAGTCATGGACCGCATCCTCAAAGATCATCAAATCGt                                                >  1:1503944/1‑148 (MQ=255)
                                                                                                   aGTGATGGATGAGCGCGACAAACGTCCGTTCTCAATGGCTTCGACGCCGGATGAAAAAGGGTTTATCGAGCTGCATATTGGCGCTTCTGAAATCAACCTTTACGCGAAAGCAGTCATGGACCGCATCCTCAAAGATCATCAAATCGTg                                               <  2:1838493/148‑1 (MQ=255)
                                                                                                       atggatgAGCGCGACAAACGTCCGTTCTCAATGGCTTCGACGCCGGATGAAAAAGGGTTTATCGAGCTGCATATTGGCGCTTCTGAAATCAACCTTTACGCGAAAGCAGTCATGGACCGCATCCTCAAAGATCATCAAATCGTGGTCg                                           >  1:3031715/1‑148 (MQ=255)
                                                                                                       atggatgAGCGCGACAAACGTCCGTTCTCAATGGCTTCGACGCCGGATGAAAAAGGGTTTATCGAGCTGCATATTGGCGCTTCTGAAATCAACCTTTACGCGAAAGCAGTCATGGACCGCATCCTCAAAGATCATCAAATCGTGGTCg                                           >  2:3472583/1‑148 (MQ=255)
                                                                                                         ggatgAGCGCGACAAACGTCCGTTCTCAATGGCTTCGACGCCGGATGAAAAAGGGTTTATCGAGCTGCATATTGGCGCTTCTGAAATCAACCTTTACGCGAAAGCAGTCATGGACCGCATCCTCAAAGATCATCAAATCGTGGTCGAc                                         >  2:2994018/1‑148 (MQ=255)
                                                                                                         ggatgAGCGCGACAAACGTCCGTTCTCAATGGCTTCGACGCCGGATGAAAAAGGGTTTATCGAGCTGCATATTGGCGCTTCTGAAATCAACCTTTACGCGAAAGCAGTCATGGACCGCATCCTCAAAGATCATCAAATCGTGGTCGAc                                         >  2:3563507/1‑148 (MQ=255)
                                                                                                           atgAGCGCGACAAACGTCCGTTCTCAATGGCTTCGACGCCGGATGAAAAAGGGTTTATCGAGCTGCATATTGGCGCTTCTGAAATCAACCTTTACGCGAAAGCAGTCATGGACCGCATCCTCAAAGATCATCAAATCGTGGTCGACAt                                       >  2:680104/1‑148 (MQ=255)
                                                                                                             gAGCGCGACAAACGTCCGTTCTCAATGGCTTCGACGCCGGATGAAAAAGGGTTTATCGAGCTGCATATTGGCGCTTCTGAAATCAACCTTTACGCGAAAGCAGTCATGGACCGCATCCTCAAAGATCATCAAATCGTGGTCGACATTc                                     <  2:527061/148‑1 (MQ=255)
                                                                                                              aGCGCGACAAACGTCCGTTCTCAATGGCTTCGACGCCGGATGAAAAAGGGTTTATCGAGCTGCATATTGGCGCTTCTGAAATCAACCTTTACGCGAAAGCAGTCATGGACCGCATCCTCAAAGATCATCAAATCGTGGTCGACATTcc                                    >  2:2624268/1‑148 (MQ=255)
                                                                                                              aGCGCGACAAAAGTCCGTTCTCAATGGCTTCGACGCCTGATGAAAAAGGGTTTATCGAGCTGAATATTGGCGCTTCTGAAATCAACCATTAAGCTAAAGCAGTCATGTACCGCATCCTCAAAGATCATCAAATCGTGGTAGACATTcc                                    >  2:2110840/1‑148 (MQ=255)
                                                                                                               gcgcGACAAACGTCCGTTCTCAATGGCTTCGACGCCGGATGAAAAAGGGTTTATCGAGCTGCATATTGGCGCTTCTGAAATCAACCTTTACGCGAAAGCAGTCATGGACCGCATCCTCAAAGATCATCAAATCGTGGTCGACATTccc                                   <  2:371751/148‑1 (MQ=255)
                                                                                                               gcgcGACAAACGTCCGTTCTCAATGGCTTCGACGCCGGATGAAAAAGGGTTTATCGAGCTGCATATTGGCGCTTCTGAAATCAACCTTTACGCGAAAGCAGTCATGGACCGCATCCTCAAAGATCATCAAATCGTGGTCGACATTccc                                   >  2:274300/1‑148 (MQ=255)
                                                                                                                cgcgACAAACGTCCGTTCTCAATGGCTTCGACGCCGGATGAAAAAGGGTTTATCGAGCTGCATATTGGCGCTTCTGAAATCAACCTTTACGCGAAAGCAGTCATGGACCGCATCCTCAAAGATCATCAAATCGTGGTCGACATTcccc                                  <  1:1074949/148‑1 (MQ=255)
                                                                                                                cgcgACAAACGTCCGTTCTCAATGGCTTCGACGCCGGATGAAAAAGGGTTTATCGAGCTGCATATTGGCGCTTCTGAAATCAACCTTTACGCGAAAGCAGTCATGGACCGCATCCTCAAAGATCATCAAATCGTGGTCGACATTcccc                                  <  2:1083738/148‑1 (MQ=255)
                                                                                                                   gACAAACGTCCGTTCTCAATGGCTTCGACGGCGGATGAAAAAGGGTTTATCGAGCTGCATATTGGCGCTTCTGAAATCAACCTTTACGCGAAAGCAGTCATGGACCGCATCCTCAAAGATCATCAAATCGTGGTCGACATTCCCCACg                               >  1:1251599/1‑148 (MQ=255)
                                                                                                                    aCAAACGTCCGTTCTCAATGGCTTCGACGCCGGATGAAAAAGGGTTTATCGAGCTGCATATTGGCGCTTCTGAAATCAACCTTTACGCGAAAGCAGTCATGGACCGCATCCTCAAAGATCATCAAATCGTGGTCGACATTCCCCACgg                              <  2:3236179/148‑1 (MQ=255)
                                                                                                                    aCAAACGTCCGTTCTCAATGGCTTCGACGCCGGATGAAAAAGGGTTTATCGAGCTGCATATTGGCGCTTCTGAAATCAACCTTTACGCGAAAGCAGTCATGGACCGCATCCTCAAAGATCATCAAATCGTGGTCGACATTCCCCACgg                              >  2:1196245/1‑148 (MQ=255)
                                                                                                                     cAAACGTCCGTTCTCAATGGCTTCGACGCCGGATGAAAAAGGGTTTATCGAGCTGCATATTGGCGCTTCTGAAATCAACCTTTACGCGAAAGCAGTCATGGACCGCATCCTCAAAGATCATCAAATCGTGGTCGACATTCCCCACGga                             >  1:50/1‑148 (MQ=255)
                                                                                                                     cAAACGTCCGTTCTCAATGGCTTCGACGCCGGATGAAAAAGGGTTTATCGAGCTGCATATTGGCGCTTCTGAAATCAACCTTTACGCGAAAGCAGTCATGGACCGCATCCTCAAAGATCATCAAATCGTGGTCGACATTCCCCACGga                             >  2:1149251/1‑148 (MQ=255)
                                                                                                                       aaCGTCCGTTCTCAATGGCTTCGACGCCGGATGAAAAAGGGTTTATCGAGCTGCATATTGGCGCTTCTGAAATCAACCTTTACGCGACAGCAGTCATGGACCGCATCCTCAAAGATCATCAAATCGTGGTCGACATTCCCCACGgaga                           >  2:3217106/1‑148 (MQ=255)
                                                                                                                       aaCGTCCGTTCTCAATGGCTTCGACGCCGGATGAAAAAGGGTTTATCGAGCTGCATATTGGCGCTTCTGAAATCAACCTTTACGCGAAAGCAGTCATGGACCGCATCCTCAAAGATCATCAAATCGTGGTCGACATTCCCCACGgaga                           <  2:1029205/148‑1 (MQ=255)
                                                                                                                           tCCGTTCTCAATGGCTTCGACGCCGGATGAAAAAGGGTTTATCGAGCTGCATATTGGCGCTTCTGAAATCAACCTTTACGCGAAAGCAGTCATGGACCGCATCCTCAAAGATCATCAAATCGTGGTCGACATTCCCCACGGAGAAGCg                       >  1:1463810/1‑148 (MQ=255)
                                                                                                                            ccGTTCTCAATGGCTTCGACGCCGGATGAAAAAGGGTTTATCGAGCTGCATATTGGCGCTTCTGAAATCAACCTTTACGCGAAAGCAGTCATGGACCGCATCCTCAAAGATCATCAAATCGTGGTCGACATTCCCCACGGAGAAGCGt                      <  1:2317804/148‑1 (MQ=255)
                                                                                                                            ccGTTCTCAATGGCTTCGACGCCGGATGAAAAAGGGTTTATCGAGCTGCATATTGGCGCTTCTGAAATCAACCTTTACGCGAAAGCAGTCATGGACCGCATCCTCAAAGATCATCAAATCGTGGTCGACATTCCCCACGGAGAAGCGt                      >  1:1288806/1‑148 (MQ=255)
                                                                                                                             cGTTCTCAATGGCTTCGACGCCGGATGAAAAAGGGTTTATCGAGCTGCATATTGGCGCTTCTGAAATCAACCTTTACGCGAAAGCAGTCATGGACCGCATCCTCAGAGATCATCAAATCGTGGTCGACATTCCCCACGGAGAAGCGTg                     >  2:2335241/1‑148 (MQ=255)
                                                                                                                              gTTCTCAATGGCTTCGACGGCGGATGAAAAAGGGTTTATCGAGCTGCATATTGGCGCTTCTGAAATCAACCTTTACGCGAAAGCAGTCATGGACCGCATCCTCAAAGATCATCAAATCGTGGTCGACATTCCCCACGGAGAAGCGTgg                    <  2:2004462/148‑1 (MQ=255)
                                                                                                                               ttCTCAATGGCTTCGACGCCGGATGAAAAAGGGTTTATCGAGCTGCATATTGGCGCTTCTGAAATCAACCTTTACGCGAAAGCAGTCATGGACCGCATCCTCAAAGATCATCAAATCGTGGTCGACATTCCCCACGGAGAAGCGTGGc                   <  1:2756341/148‑1 (MQ=255)
                                                                                                                               ttCTCAAGGGCTTCGACGCCGGATGAAAAAGGGTTTATCGAGCTGCATATTGGCGCTTCTGAAATCAACCTTTACGCGAAAGCAGTCATGGACCGCATCCTCAAAGATCATCAAATCGTGGTCGACATTCCCCACGGAGAAGCGTGGc                   <  1:475310/148‑1 (MQ=255)
                                                                                                                                tctcAATGGCTTCGACGCCGGATGAAAAAGGGTTTATCGAGCTGCATATTGGCGCTTCTGAAATCAACCTTTACGCGAAAGCAGTCATGGACCGCATCCTCAAAGATCATCAAATCGTGGTCGACATTCCCCACGGAGAAGCGTGGCt                  >  2:2005431/1‑148 (MQ=255)
                                                                                                                                 ctcAATGGCTTCGACGCCGGATGAAAAAGGGTTTATCGAGCTGCATATTGGCGCTTCTGAAATCAACCTTTACGCGAAAGCAGTCATGGACCGCATCCTCAAAGATCATCAAATCGTGGTCGACATTCCCCACGGAGAAGCGTGGCTg                 >  1:3123049/1‑148 (MQ=255)
                                                                                                                                 ctcAATGGCTTCGACGCCGGATGAAAAAGGGTTTATCGAGCTGCATATTGGCGCTTCTGAAATCAACCTTTACGCGAAAGCAGTCATGGACCGCATCCTCAAAGATCATCAAATCGTGGTCGACATTCCCCACGGAGAAGCGTAGCTg                 <  1:2899167/148‑1 (MQ=255)
                                                                                                                                     aTGGCTTCGACGGCGGATGAAAAAGGGTTTATCGAGCTGCATATTGGCGCTTCTGAAATCAACCTTTACGCGAAAGCAGTCATGGACCGCATCCTCAAAGATCATCAAATCGTGGTCGACATTCCCCACGGAGAAGCGTGGCTGcgcg             <  1:2291202/148‑1 (MQ=255)
                                                                                                                                         cTTCGACGCCGGATGAAAAAGGGTTTATCGAGCTGCATATTGGCGCTTCTGAAATCAACCTTTACGCGAAAGCAGTCATGGACCGCATCCTCAAAGATCATCAAATCGTGGTCGACATTCCCCACGGAGAAGCGTGGCTGCGCgatga         >  1:2033323/1‑148 (MQ=255)
                                                                                                                                          ttCGACGCCGGATGAAAAAGGGTTTATCGAGCTGCATATTGGCGCTTCTGAAATCAACCTTTACGCGAAAGCAGTCATGGACCGCATCCTCAAAGATCATCAAATCGTGGTCGACATTCCCCACGGAGAAGCGTGGCTGCGCgatgat        <  2:2237271/148‑1 (MQ=255)
                                                                                                                                          ttCGACGCCGGATGAAAAAGGGTTTATCGAGCTGCATATTGGCGCGTCTGAAATCAACCTTTACGCGAAAGCAGTCATGGACCGCATCCTCAAAGATCATCAAATCGTGGTCGACATTCCCCACGGAGAAGCGTGGCTGCGCgatgat        <  1:2395850/148‑1 (MQ=255)
                                                                                                                                           tCGACGCCGGATGAAAAAGGGTTTATCGAGCTGCATATTGGCGCTTCTGAAATCAACCTTTACGCGAAAGCAGTCATGGACCGCATCCTCAAAGATCATCAAATCGTGGTCGACATTCCCCACGGAGAAGCGTGGCTGCGCGatgatg       >  1:2181947/1‑148 (MQ=255)
                                                                                                                                            cGACGCCGGATGAAAAAGGGTTTATCGAGCTGCATATTGGCGCTTCTGAAATCAACCTTTACGCGAAAGCAGTCATGGACCGCATCCTCAAAGATCATCAAATCGTGGTCGACATTCCCCACGGAGAAGCGTGGCTGCGCGAtgatga      <  1:1095367/148‑1 (MQ=255)
                                                                                                                                            cGACGCCGGATGAAAAAGGGTTTATCGAGCTGCATATTGGCGCTTCTGAAATCAACCTTTACGCGAAAGCAGTCATGGACCGCATCCTCAAAGATCATCAAATCGTGGTCGACATTCCCCACGGAGAAGCGTGGCTGCGCGAtgatga      <  1:2874975/148‑1 (MQ=255)
                                                                                                                                             gACGCCGGATGAAAAAGGGTTTATCGAGCTGCATATTGGCGCTTCTGAAATCAACCTTTACGCGAAAGCAGTCATGGACCGCATCCTCAAAGATCATCAAATCGTGGTCGACATTCCCCACGGAGAAGCGTGGCTGCGCGATGATGaa     <  1:3072683/148‑1 (MQ=255)
                                                                                                                                             gACGCCGGATGAAAAAGGGTTTATCGAGCTGCATATTGGCGCTTCTGAAATCAACCTTTACGCGAAAGCAGTCATGGACCGCATCCTCAAAGATCATCAAATCGTGGTCGACATTCCCCACGGAGAAGCGTGGCTGCGCGATGATGaa     >  1:2209908/1‑148 (MQ=255)
                                                                                                                                              aCGCCGGATGAAAAAGGGTTTATCGAGCTGCATATTGGCGCTTCTGAAATCAACCTTTACGCGAAAGCAGTCATGGACCGCATCCTCAAAGATCATCAAATCGTGGTCGACATTCCCCACGGAGAAGCGTGGCTGCGCGATGATTAag    >  2:3273425/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
TAAAGTGACCTCGGTAGAAGCTATCACGGATACCGTATATCGTGTCCGCATCGTGCCAGACGCGGCCTTTTCTTTTCGTGCTGGTCAGTATTTGATGGTAGTGATGGATGAGCGCGACAAACGTCCGTTCTCAATGGCTTCGACGCCGGATGAAAAAGGGTTTATCGAGCTGCATATTGGCGCTTCTGAAATCAACCTTTACGCGAAAGCAGTCATGGACCGCATCCTCAAAGATCATCAAATCGTGGTCGACATTCCCCACGGAGAAGCGTGGCTGCGCGATGATGAAGAG  >  CP009273/4019904‑4020195

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑