Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,231,349 A→C 39.2% Q227P (CAG→CCG)  ygjR → putative NAD(P)‑dependent dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,231,3490AC39.2% 156.1 / inf 299Q227P (CAG→CCG) ygjRputative NAD(P)‑dependent dehydrogenase
Reads supporting (aligned to +/- strand):  ref base A (100/82);  new base C (44/73);  total (144/155)
Fisher's exact test for biased strand distribution p-value = 4.36e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.85e-01

Only 100 of 309 total aligned reads displayed.

TATGGATATCGGTTTTTACTGTCTGGCGTCGGCGGTGGCGTTATTTGGTGAGCCGAAAAGCGTGCAGGCAACCGCCAGTTTGCTGGCAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCAGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGCTGGTGATTGAAAAACTGTCTGAATGCCAGAAAGTGTGCTTCGTGCCGCGTGGCAGCCAAATGCAGGATCTCACCCAGCCGCAGCATATTA  >  CP009273/3231200‑3231498
                                                                                                                                                     |                                                                                                                                                     
tATGGATATCGGTTTTTACTGTCTGGCGTCGGCGGTGGCGTTATTTGGTGAGCCGAAAAGCGTGCAGGCAACCGCCAGTTTGCTGGCAAGCGGCGTTGACGCCCCGGGCGTGGGGGTGATGGATTACGGGGATTTCAGCGTCACCTTgca                                                                                                                                                       >  1:1083777/1‑150 (MQ=255)
     atatCGGTTTTTACTGTCTGGCGTCGGCGGTGGCGTTATTTGGTGAGCCGAAAAGCGTGCAGGCAACCGCCAGTTTGCTGGCAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCAGCAcc                                                                                                                                                  <  2:416916/150‑2 (MQ=255)
     atatCGGTTTTTACTGTCTGGCGTCGGCGGTGGCGTTATTTGGTGAGCCGAAAAGCGTGCAGGCAACCGCCAGTTTGCTGGCAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCAGCACt                                                                                                                                                  <  2:3601717/150‑1 (MQ=255)
      tatCGGTTTTTACTGTCTGGCGTCGGCGGTGGCGTTATTTGGTGAGCCGAAAAGCGTGCAGGCAACCGCCAGTTTGCTGGCAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCAGCACTc                                                                                                                                                 <  1:146964/150‑1 (MQ=255)
      tatCGGTTTTTACTGTCTGGCGTCGGCGGTGGCGTTATTTGGTGAGCCGAAAAGCGTGCAGGCAACCGCCAGTTTGCTGGCAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCAGCACTc                                                                                                                                                 >  2:2151119/1‑150 (MQ=255)
      gatCGGTTTTTACTGTCTGGCGTCGGCGGTGGCGTTATTTGGTGAGACGAAAAGCGTGCAGGCAACCGCCAGTTTGCTTGCAAGCGGCGTTGACGCCCCGGGCGTGTTGGTGATGGTTTACGGTGATTTAAGAGGCCCCATGCCGCAcac                                                                                                                                                 >  1:1445854/2‑148 (MQ=255)
       atCGGTTTTTACTGTCTGGCGTCGGCGGTGGCGTTATTTGGTGAGCCGAAAAGCGTGCAGGCAACCGCCAGTTTGCTGGCAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCAGCACTCa                                                                                                                                                >  1:3039526/1‑150 (MQ=255)
        tCGGTTTTTACTGTCTGGCGTCGGCGGTGGCGTTATTTGGTGAGCCGAAAAGCGTGCAGGCAACCGCCAGTTTGCTGGCAAGCGGCGTTGTCGCCCAGGGCGTGGTGGTGTTGGATTACGGTGATTTAAGCGTCACCTTTCAGCACTCaa                                                                                                                                               >  1:2917854/1‑150 (MQ=255)
        tCGGTTTTTACTGTCTGGCGTCGGCGGTGGCGTTATTTGGTGAGCCGAAAAGCGTGCAGGCAACCGCCAGTTTGCTGGCAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCAGCACTCaa                                                                                                                                               >  2:3781349/1‑150 (MQ=255)
         cGGTTTTTACTGTCTGGCGTCGGCGGTGGCGTTATTTGGTGAGCCGAAAAGCGTGCAGGCAACCGCCAGTTTGCTGGCAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTAAGGTGATTTAAGCGTCACCTTGCAGCACTCaaa                                                                                                                                              >  1:1893174/1‑150 (MQ=255)
          ggTTTTTACTGTCTGGCGTCGGCGGTGGCGTTATTTGGTGAGCCGAAAAGCGTGCAGGCAACCGCCAGTTTGCTGGCAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCAGCACTCaaaa                                                                                                                                             >  2:3446338/1‑150 (MQ=255)
            tttttACTGTCTGGCGTCGGCGGTGGCGTTATTTGGTGAGCCGAAAAGCGTGCAGGCAACCGCCAGTTTGCTGGCAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGt                                                                                                                                           <  2:2059509/150‑1 (MQ=255)
            tttttACTGTCTGGCGTCGGCGGTGGCGTTATTTGGTGAGCCGAAAAGCGTGCAGGCAACCGCCAGTTTGCTGGCAAGCGGCGGTGTCGCCCAGGGCGGGGTTGGGGTTGATTACGGTGTTTTAATCGTCACCATGCCGCACGAAAAAGt                                                                                                                                           >  1:737320/1‑150 (MQ=255)
               ttACTGTCTGGCGTCGGCGGTGGCGTTATTTGGTGAGCCGAAAAGCGTGCAGGCAACCGCCAGTTTGCTGGCAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCAGCACTCAAAAGTCAg                                                                                                                                        >  1:3577632/1‑150 (MQ=255)
               ttACTGTCTGGCGTCGGCGGTGGCGTTATTTGGTGAGCCGAAAAGCGTGCAGGCAACCGCCAGTTTGCTGGCAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCAGCACTCAAAAGTCAg                                                                                                                                        <  2:884235/150‑1 (MQ=255)
                tACTGTCTGGCGTCGGCGGTGGCGTTATTTGGTGAGCCGAAAAGCGTGCAGGCAACCGCCAGTTTGCTGGCAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCAGCACTCAAAAGTCAGt                                                                                                                                       >  1:708500/1‑150 (MQ=255)
                  cTGTCTGGCGTCGGCGGTGGCGTTATTTGGTGAGCCGAAAAGCGTGCAGGCAACCGCCAGTTTGCTGGCAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCAGCACTCAAAAGTCAGTGa                                                                                                                                     >  2:748536/1‑150 (MQ=255)
                     tCTGGCGTCGGCGGTGGCGTTATTTGGTGAGCCGAAAAGCGTGCAGGCAACCGCCAGTTTGCTGGCAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCAGCACTCAAAAGTCAGTGATTc                                                                                                                                  >  2:552534/1‑150 (MQ=255)
                     tCTGGCGTCGGCGGTGGCGTTATTTGGTGAGCCGAAAAGCGTGCAGGCAACCGCCAGTTTGCTGGCAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCAGCACTCAAAAGTCAGTGATTc                                                                                                                                  <  1:3856453/150‑1 (MQ=255)
                     tCTGGCGTCGGCGGTGGCGTTATTTGGTGAGCCGAAAAGCGTGCAGGCAACCGCCAGTTTGCTGGCAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCAGCACTCAAAAGTCAGTGATTc                                                                                                                                  <  2:3080597/150‑1 (MQ=255)
                     tCTGGCGTCGGCGGTGGCGTTATTTGGTGAGCCGAAAAGCGTGCAGGCAACCGCCAGTTTGCTGGCAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCAGCACTCAAAAGTCAGTGATTc                                                                                                                                  >  2:260513/1‑150 (MQ=255)
                     tCTGGCGTCGGCGGTGGCGTTATTTGGTGAGCCGAAAAGCGTGCAGGCAACCGCCAGTTTGCTGGCAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCAGCACTCAAAAGTCAGTGATTc                                                                                                                                  >  2:3328799/1‑150 (MQ=255)
                      cTGGCGTCGGCGGTGGCGTTATTTGGTGAGCCGAAAAGCGTGCAGGCAACCGCCAGTTTGCTGGCAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCt                                                                                                                                 >  2:3363740/1‑150 (MQ=255)
                       tGGCGTCGGCGGTGGCGTTATTTGGTGAGCCGAAAAGCGTGCAGGCAACCGCCAGTTTGCTGGCAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCAGCACTCAAAAGTCAGTGATTCTg                                                                                                                                >  2:1596170/1‑150 (MQ=255)
                           gTCGGCGGTGGCGTTATTTGGTGAGCCGAAAAGCGTGCAGGCAACCGCCAGTTTGCTGGCAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCt                                                                                                                            <  2:1331033/150‑1 (MQ=255)
                           gTCGGCGGTGGCGTTATTTGGTGAGCCGAAAAGCGTGCAGGCAACCGCCAGTTTGCTGGCAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCt                                                                                                                            <  2:4087437/150‑1 (MQ=255)
                            tCGGCGGTGGCGTTATTTGGTGAGCCGAAAAGCGTGCAGGCAACCGCCAGTTTGCTGGCAAGTGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCAGCACTCAAAAGTCAGTGATTCTGTCCTg                                                                                                                           <  2:1188203/150‑1 (MQ=255)
                             cggcggTGGCGTTATTTGGTGAGCCGAAAAGCGTGCAGGCAACCGCCAGTTTGCTGGCAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTAAGCGTCACCTTTCAGCACACAAAAGTCAGTGATTCTGTCCTgg                                                                                                                          >  1:463706/1‑150 (MQ=255)
                              ggcggTGGCGTTATTTGGTGAGCCGAAAAGCGTGCAGGCAACCGCCAGTTTGCTGGCAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTAAGCGTCACCATGCAGCACACAAAATTCAGTGATTCTGTCCTGgc                                                                                                                         >  1:2228535/1‑150 (MQ=255)
                                cggTGGCGTTATTTGGTGAGCCGAAAAGCGTGCAGGCAACCGCCAGTTTGCTGGCAAGCGGCGTTGACGCCCAGGGCGTGGGGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGgcga                                                                                                                       >  2:1872442/1‑150 (MQ=255)
                                cggTGGCGTTATTTGGTGAGCCGAAAAGAGTGCAGGCAACCGACAGTTTGCTGGCAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCAGCACTCAAAAGTCAGTGATTCTGTCCTGgcta                                                                                                                       >  2:3831301/1‑148 (MQ=255)
                                        ttATTTGGTGAGCCGAAAAGCGTGCAGGCAACCGCCAGTTTGCTGGCAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCAGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCgagagt                                                                                                               >  2:505381/1‑148 (MQ=255)
                                         tATTTGGTGAGCCGAAAAGCGTGCAGGCAACCGCCAGTTTGCTGGCAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTc                                                                                                              >  2:328111/1‑150 (MQ=255)
                                         tATTTGGTGAGCCGAAAAGCGTGCAGGCAACCGCCAGTTTGCTGGCAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCAGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTc                                                                                                              <  2:1409498/150‑1 (MQ=255)
                                          aTTTGGTGAGCCGAAAAGCGTGCAGGCAACCGCCAGTTTGCTGGCAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCAGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCa                                                                                                             <  2:3009700/150‑1 (MQ=255)
                                                 gAGCCGAAAAGCGTGCAGGCAACCGCCAGTTTGCTGGCAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCAGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGaa                                                                                                      >  2:1511736/1‑150 (MQ=255)
                                                  aGCCGAAAAGCGTGCAGGCAACCGCCAGTTTGCTGGCAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAg                                                                                                     <  1:4121654/150‑1 (MQ=255)
                                                   gCCGAAAAGCGTGCAGGCAACCGCCAGTTTGCTGGCAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGc                                                                                                    <  2:1130173/150‑1 (MQ=255)
                                                     cGAAAAGCGTGCAGGCAACCGCCAGTTTGCTGGCAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCAGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAg                                                                                                  >  1:292539/1‑150 (MQ=255)
                                                         aaGCGTGCAGGCAACCGCCAGTTTGCTGGCAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCAGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATc                                                                                              <  2:2937985/150‑1 (MQ=255)
                                                              tGCAGGCAACCGCCAGTTTGCTGGCAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGCTgg                                                                                         <  2:2519389/150‑1 (MQ=255)
                                                              tGCAGGCAACCGCCAGTTTGCTGGCAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGCTgg                                                                                         <  1:883516/150‑1 (MQ=255)
                                                                 aGGCAACCGCCAGTTTGCTGGCAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCAGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGCTGGTGa                                                                                      >  1:3385549/1‑150 (MQ=255)
                                                                 aGGCAACCGCCAGTTTGCTGGCAAGCGGCGTTGACGCCCAGGGCGGGGTGGGGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGCTGGTGa                                                                                      >  2:1888610/1‑150 (MQ=255)
                                                                  ggCAACCGCCAGTTTGCTGGCAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCAGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGCTGGTGAt                                                                                     <  2:713391/150‑1 (MQ=255)
                                                                    cAACCGCCAGTTTGCTGGCAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGCTGGTGATTg                                                                                   <  2:3862620/150‑1 (MQ=255)
                                                                    cAACCGCCAGTTTGCTGGCAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGCTGGTGATTg                                                                                   <  2:4048907/150‑1 (MQ=255)
                                                                    cAACCGCCAGTTTGCTGGCAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCAGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGCTGGTGATTg                                                                                   <  1:1609758/150‑1 (MQ=255)
                                                                      aCCGCCAGTTTGCTGGCAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCAGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGCTGGTGATTGaa                                                                                 <  1:2833432/150‑1 (MQ=255)
                                                                      aCCGCCAGTTTGCTGGCAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCAGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGCTGGTGATTGaa                                                                                 >  2:2917330/1‑150 (MQ=255)
                                                                        cGCCAGTTTGCTGGCAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCAGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGCTGGTGATTGaaaa                                                                               <  2:371548/150‑1 (MQ=255)
                                                                            aGTTTGCTGGCAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCAGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGCTGGTGATTGAAAAACTg                                                                           >  1:3537433/1‑150 (MQ=255)
                                                                            aGTTTGCTGGCAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCAGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGCTGGTGATTGAAAAACTg                                                                           >  1:2596588/1‑150 (MQ=255)
                                                                              tttGCTGGCAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGCTGGTGATTGAAAAACTGTc                                                                         <  1:2264890/150‑1 (MQ=255)
                                                                               ttGCTGGCAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGCTGGTGATTGAAAAACTGTCt                                                                        <  1:3086515/150‑1 (MQ=255)
                                                                               ttGCTGGCAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGCTGGTGATTGAAAAACTGTCt                                                                        <  1:2956518/150‑1 (MQ=255)
                                                                               ttGCTGGCAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCAGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGCTGGTGATTGAAAAACTGTCt                                                                        <  2:1389738/150‑1 (MQ=255)
                                                                               ttGCTGGCAAGCGGAGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGCTGGTGATTGAAAAACTGTCt                                                                        <  1:3297466/150‑1 (MQ=255)
                                                                                tGCTGGCAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGCTGGTGATTGAAAAACTGTCTg                                                                       >  2:622717/1‑150 (MQ=255)
                                                                                    ggCAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCAGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGCTGGTGATTGAAAAACTGTCTGAATg                                                                   >  2:2879176/1‑150 (MQ=255)
                                                                                       aaGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCAGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGCTGGTGATTGAAAAACTGTCTGAATGCCa                                                                <  2:1430931/150‑1 (MQ=255)
                                                                                       aaGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCAGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGCTGGTGATTGAAAAACTGTCTGAATGCCa                                                                <  2:7970/150‑1 (MQ=255)
                                                                                        aGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGCTGGTGATTGAAAAACTGTCTGAATGCCAg                                                               >  2:903740/1‑150 (MQ=255)
                                                                                        aGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCAGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTAAGGGCGAAGCAGGATCGCTGGTGATTGAAAAACTGTCTTAATTCCAg                                                               >  2:4291719/1‑150 (MQ=255)
                                                                                           ggcgTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCAGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGCTGGTGATTGAAAAACTGTCTGAATGCCAGaaa                                                            >  2:2363978/1‑150 (MQ=255)
                                                                                            gcgTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGCTGGTGATTGAAAAACTGTCTGAATGCCAGAAAg                                                           <  2:1470571/150‑1 (MQ=255)
                                                                                            gcgTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGCTGGTGATTGAAAAACTGTCTGAATGCCAGAAAg                                                           <  1:4302472/150‑1 (MQ=255)
                                                                                             cgTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCAGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGCTGGTGATTGAAAAACTGTCTGAATGCCAGAAAgt                                                          <  1:2204407/150‑1 (MQ=255)
                                                                                               ttGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGCTGGTGATTGAAAAACTGTCTGAATGCCAGAAAgtgt                                                        >  1:4027825/1‑150 (MQ=255)
                                                                                               ttGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGCTGGTGATTGAAAAACTGTCTGAATGCCAGAAAgtgt                                                        >  2:4221868/1‑150 (MQ=255)
                                                                                                  aCGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCAGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGCTGGTGATTGAAAAACTGTCTGAATGCCAGAAAGTGTGCt                                                     <  1:600752/150‑1 (MQ=255)
                                                                                                   cGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCAGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGCTGGTGATTGAAAAACTGTCTGAATGCCAGAAAGTGTGCtt                                                    >  2:3233566/1‑150 (MQ=255)
                                                                                                   cGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCAGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGCTGGTGATTGAAAAACTGTCTGAATGCCAGAAAGTGTGCtt                                                    >  2:1743294/1‑150 (MQ=255)
                                                                                                      ccAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGCTGGTGATTGAAAAACTGTCTGAATGCCAGAAAGTGTGCTTCGt                                                 <  2:2577683/150‑1 (MQ=255)
                                                                                                       cAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCAGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGCTGGTGATTGAAAAACTGTCTGAATGCCAGAAAGTGTGCTTCGTg                                                <  1:2463138/150‑1 (MQ=255)
                                                                                                         gggCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCAGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGCTGGTGATTGAAAAACTGTCTGAATGCCAGAAAGTGTGCTTCGTGcc                                              <  1:2147631/150‑1 (MQ=255)
                                                                                                          ggCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCAGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGCTGGTGATTGAAAAACTGTCTGAATGCCAGAAAGTGTGCTTCGTGCcg                                             <  1:1547413/150‑1 (MQ=255)
                                                                                                          ggCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCAGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGCTGGTGATTGAAAAACTGTCTGAATGCCAGAAAGTGTGCTTCGTGCcg                                             <  1:3331604/150‑1 (MQ=255)
                                                                                                            cGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGCTGGTGATTGAAAAACTGTCTGAATGCCAGAAAGTGTGCTTCGTGCcgcg                                           >  2:3509383/1‑150 (MQ=255)
                                                                                                            cGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGCTGGTGATTGAAAAACTGTCTGAATGCCAGAAAGTGTGCTTCGTGCcgcg                                           <  2:1017929/150‑1 (MQ=255)
                                                                                                                gtggtgATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGCTGGTGATTGAAAAACTGTCTGAATGCCAGAAAGTGTGCTTCGTGCCGCGTGGc                                       >  2:3666714/1‑150 (MQ=255)
                                                                                                                 tggtgATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGCTGGTGATTGAAAAACTGTCTGAATGCCAGAAAGTGTGCTTCGTGCCGCGTGGCa                                      >  1:3705626/1‑150 (MQ=255)
                                                                                                                      aTGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGCTGGTGATTGAAAAACTGTCTGAATGCCAGAAAGTGTGCTTCGTGCCGCGTGGCAGCCaa                                 >  1:3965849/1‑150 (MQ=255)
                                                                                                                          aTTACGGTGATTTCAGCGTCACCTTGCAGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGCTGGTGATTGAAAAACTGTCTGAATGCCAGAAAGTGTGCTTCGTGCCGCGTGGCAGCCAAATGc                             <  1:1359236/150‑1 (MQ=255)
                                                                                                                             aCGGTGATTTCAGCGTCACCTTGCAGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGCTGGTGATTGAAAAACTGTCTGAATGCCAGAAAGTGTGCTTCGTGCCGCGTGGCAGCCAAATGCAgg                          <  1:4217645/150‑1 (MQ=255)
                                                                                                                             aCGGTGATTTCAGCGTCACCTTGCAGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGCTGGTGATTGAAAAACTGTCTGAATGCCAGAAAGTGTGCTTCGTGCCGCGTGGCAGCCAAATGCAgg                          <  1:277348/150‑1 (MQ=255)
                                                                                                                              cGGTGATTTCAGCGTCACCTTGCAGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGCTGGTGATTGAAAAACTGTCTGAATGCCAGAAAGTGTGCTTCGTGCCGCGTGGCAGCCAAATGCAGGa                         >  2:4260083/1‑150 (MQ=255)
                                                                                                                                   aTTTCAGCGTCACCTTGCAGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGCTGGTGATTGAAAAACTGTCTGAATGCCAGAAAGTGTGCTTCGTGCCGCGTGGCAGCCAAATGCAGGATCTCa                    >  1:3465352/1‑150 (MQ=255)
                                                                                                                                   aTTTCAGCGTCACCTTGCAGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGCTGGTGATTGAAAAACTGTCTGAATGCCAGAAAGTGTGCTTCGTGCCGCGTGGCAGCCAAATGCAGGATCTCa                    >  1:2707240/1‑150 (MQ=255)
                                                                                                                                     ttCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGCTGGTGATTGAAAAACTGTCTGAATGCCAGAAGGTGTGCTTCGTGCCGCGTGGCAGCCAAATGCAGGATCTCAcc                  >  2:51823/1‑150 (MQ=255)
                                                                                                                                       cAGCGTCACCTTGCAGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGCTGGTGATTGAAAAACTGTCTGAATGCCAGAAAGTGTGCTTCGTGCCGCGTGGCAGCCAAATGCAGGATCTCACCCa                <  1:2492564/150‑1 (MQ=255)
                                                                                                                                       cAGCGTCACCTTGCAGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGCTGGTGATTGAAAAACTGTCTGAATGCCAGAAAGTGTGCTTCGTGCCGCGTGGCAGCCAAATGCAGGATCTCACCCa                <  2:1956208/150‑1 (MQ=255)
                                                                                                                                          cGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGCTGGTGATTGAAAAACTGTCTGAATGCCAGAAAGTGTGCTTCGTGCCGCGTGGCAGCCAAATGCAGGATCTCACCCAGcc             >  2:2723552/1‑150 (MQ=255)
                                                                                                                                           gTCACCTTGCAGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGCTGGTGATTGAAAAACTGTCTGAATGCCAGAAAGTGTGCTTCGTGCCGCGTGGCAGCCAAATGCAGGATCTCACCCAGCCg            <  2:3874041/150‑1 (MQ=255)
                                                                                                                                           gTCACCTTGCAGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGCTGGTGATTGAAAAACTGTCTGAATGCCAGAAAGTGTGCTTCGTGCCGCGTGGCAGCCAAATGCAGGATCTCACCCAGCCg            >  2:2652083/1‑150 (MQ=255)
                                                                                                                                             cACCTTGCAGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGCAGCAGGATCGCTGGTGATTGAAAAACTGTCTGAATGCCAGAAAGTGTGCTTCGTGCCGCGTGGCAGCCAACTGCAGGATCTCACCCAGCCgca          <  2:2485358/150‑1 (MQ=255)
                                                                                                                                              aCCTTGCAGCACTCAATAGTGAGTGATTCTGTCCTGGAGAGCGAGATTCAGGGCGAAGCAGGATCGCTGTTGATTGAAAAACTGTCTGAATGCCAGAAAGTGTGCTTCGTTCCGCGTGGCAGCCAAATGCAGGATCTCACCCAGCCgcag         <  1:748408/150‑1 (MQ=255)
                                                                                                                                              aCCTTGCAGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGCTGGTGATTGAAAAACTGTCTGAATGCCAGAAAGTGTGCTTCGTGCCGCGTGGCAGCCAAATGCAGGATCTCACCCAGCCgcag         <  1:2977946/150‑1 (MQ=255)
                                                                                                                                              aCCTTGCAGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGCTGGTGATTGAAAAACTGTCTGAATGCCAGAAAGTGTGCTTCGTGCCGCGTGGCAGCCAAATGCAGGATCTCACCCAGCCgcag         <  2:305028/150‑1 (MQ=255)
                                                                                                                                                     agcaCTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGCTGGTGATTGAAAAACTGTCTGAATGCCAGAAAGTGTGCTTCGTGCCGCGTGGCAGCCAAATGCAGGATCTCACCCAGCCGCAGCATATTa  <  1:2097922/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TATGGATATCGGTTTTTACTGTCTGGCGTCGGCGGTGGCGTTATTTGGTGAGCCGAAAAGCGTGCAGGCAACCGCCAGTTTGCTGGCAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCAGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGCTGGTGATTGAAAAACTGTCTGAATGCCAGAAAGTGTGCTTCGTGCCGCGTGGCAGCCAAATGCAGGATCTCACCCAGCCGCAGCATATTA  >  CP009273/3231200‑3231498

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: