Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,353,069 A→T 6.4% I19I (ATT→ATA)  yjdC ← putative transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,353,0690AT 6.4% 991.9 / 41.3 297I19I (ATT→ATA) yjdCputative transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base A (118/160);  new base T (11/8);  total (129/168)
Fisher's exact test for biased strand distribution p-value = 2.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

Only 100 of 304 total aligned reads displayed.

AGCATCAGCTGACGACGCCAGACATCTATCTGTTGGCTAAGATAACGCAGCGCATCGTAGAGGATTGCCTCTTTGTCTGGCCAGAAGCGGCGTAGCTCGTCCAGTGGATAATCCACACGTTCAGCAACCATCTCCAGCGTGGTGTTGGCAATCCCTTGTAATTCTAATAATTTCAGGGCTTCTCCCAGTACATCTTCACGTTGCACGCTATTTTCCTCCGTCTTTCCCACTGCAAGTGTCGTTCACGGTTGGCGATCGCGCAAATGTGCGCTGAAGGTTTCAGCATCCATAAAGCCCGT  >  CP009273/4352920‑4353218
                                                                                                                                                     |                                                                                                                                                     
aGCATCAGCTGACGACGCCAGACATCTATCTGTTGGCTAAGATAACGCAGCGCATCGTAGAGGATTGCCTCTTTGTCTGGCCAGAAGCGGCGTAGCTCGTCCAGTGGATAATCCACACGTTCAGCAACCATCTCCAGCGTGGTGTTGGCa                                                                                                                                                       >  2:1212857/1‑150 (MQ=255)
  cATCAGCTGACGACGCCAGACATCTATCTGTTGGCTAAGATAACGCAGCGCATCGTAGAGGATTGCCTCTTTGTCTGGCCAGAAGCGGCGTAGCTCGTCCAGTGGATAATCCACACGTTCAGCAACCATCTCCAGCGTGGTGTTGGCact                                                                                                                                                     >  2:2481111/1‑148 (MQ=255)
   aTCAGCTGACGACGCCAGACATCTATCTGTTGGCTAAGATAACGCAGCGCATCGTAGAGGATTGCCTCTTTGTCTGGCCAGAAGCGGCGTAGCTCGTCCAGTGGATAATCCACACGTTCAGCAACCATCTCCAGCGTGGTGTTGGCAATc                                                                                                                                                    >  2:2744836/1‑150 (MQ=255)
       gCTGACGACGCCAGACATCTATCTGTTGGCTAAGATAACGCAGCGCATCGTAGAGGATTGCCTCTTTGTCTGGCCAGAAGCGGCGTAGCTCGTCCAGTGGATAATCCACACGTTCAGCAACCATCTCCAGCGTGGTGTTGGCTATCCCtt                                                                                                                                                <  1:305011/150‑1 (MQ=255)
       gCTGACGACGCCAGACATCTATCTGTTGGCTAAGATAACGCAGCGCATCGTAGAGGATTGCCTCTTTGTCTGGCCAGAAGCGGCGTAGCTCGTCCAGTGGATAATCCACACGTTCAGCAACCATCTCCAGCGTGGTGTTGGCAATCCCtt                                                                                                                                                <  2:360262/150‑1 (MQ=255)
       gCTGACGACGCCAGACATCTATCTGTTGGCTAAGATAACGCAGCGCATCGTAGAGGATTGCCTCTTTGTCTGGCCAGAAGCGGCGTAGCTCGTCCAGTGGATAATCCACACGTTCAGCAACCATCTCCAGCGTGGTGTTGGCAATCCCtt                                                                                                                                                <  1:1376229/150‑1 (MQ=255)
         tGACGACGCCAGACATCTATCTGTTGGCTAAGATAACGCAGCGCATCGTAGAGGATTGCCTCTTTGTCTGGCCAGAAGCGGCGTAGCTCGTCCAGTGGATAATCCACACGTTCAGCAACCATCTCCAGCGTGGTGTTGGCAATCCCTTGt                                                                                                                                              <  1:2577498/150‑1 (MQ=255)
         tGACGACGCCAGACATCTATCTGTTGGCTAAGATAACGCAGCGCATCGTAGAGGATTGCCTCTTTGTCTGGCCAGAAGCGGCGTAGCTCGTCCAGTGGATAATCCACACGTTCAGCAACCATCTCCAGCGTGGTGTTGGCAATCCCTTGt                                                                                                                                              <  1:2952449/150‑1 (MQ=255)
               cgCCAGACATCTATCTGTTGGCTAAGATAACGCAGCGCATCGTAGAGGATTGCCTCTTTGTCTGGCCAGAAGCGGCGTAGCTCGTCCAGTGGATAATCCACACGTTCAGCAACCATCTCCAGCGTGGTGTTGGCAATCCCTTGTAATTCt                                                                                                                                        <  1:1610433/150‑1 (MQ=255)
                  cAGACATCTATCTGTTGGCTAAGATAACGCAGCGCATCGTAGAGGATTGCCTCTTTGTCTGGCCAGAAGCGGCGTAGCTCGTCCAGTGGATAATCCACACGTTCAGCAACCATCTCCAGCGTGGTGTTGGCAATCCCTTGTAATTCtaat                                                                                                                                     <  1:2129222/150‑1 (MQ=255)
                    gACATCTATCTGTTGGCTAAGATAACGCAGCGCATCGTAGAGGATTGCCTCTTTGTCTGGCCAGAAGCGGCGTAGCTCGTCCAGTGGATAATCCACACGTTCAGCAACCATCTCCAGCGTGGTGTTGGCAATCCCTTGTAATTCtaataa                                                                                                                                   >  1:1912593/1‑150 (MQ=255)
                         ctatctGTTGGCTAAGATAACGCAGCGCATCGTAGAGGATTGCCTCTTTGTCTGGCCAGAAGCGGCGTAGCTCGTCCAGTGGATAATCCACACGTTCAGCAACCATCTCCAGCGTGGTGTTGGCAATCCCTTGTAATTCTAATAATTTCa                                                                                                                              <  1:2586325/150‑1 (MQ=255)
                          tatctGTTGGCTAAGATAACGCAGCGCATCGTAGAGGATTGCCTCTTTGTCTGGCCAGAAGCGGCGTAGCTCGTCCAGTGGATAATCCACACGTTCAGCAACCATCTCCAGCGTGGTGTTGGCAATCCCTTGTAATTCTAATAATTTCAg                                                                                                                             <  2:2124511/150‑1 (MQ=255)
                            tctGTTGGCTAAGATAACGCAGCGCATCGTAGAGGATTGCCTCTTTGTCTGGCCAGAAGCGGCGTAGCTCGTCCAGTGGATAATCCACACGTTCAGCAACCATCTCCAGCGTGGTGTTGGCAATCCCTTGTAATTCTAATAATTTCAggg                                                                                                                           <  2:1271795/150‑1 (MQ=255)
                             ctGTTGGCTAAGATAACGCAGCGCATCGTAGAGGATTGCCTCTTTGTCTGGCCAGAAGCGGCGTAGCTCGTCCAGTGGATAATCCACACGTTCAGCAACCATCTCCAGCGTGGTGTTGGCAATCCCTTGTAATTCTAATAATTTCAGGGc                                                                                                                          >  2:2865988/1‑150 (MQ=255)
                              tGTTGGCTAAGATAACGCAGCGCATCGTAGAGGATTGCCTCTTTGTCTGGCCAGAAGCGGCGTAGCTCGTCCAGTGGATAATCCACACGTTCAGCAACCATCTCCAGCGTGGTGTTGGCAATCCCTTGTAATTCTAATAATTTCAGGGCt                                                                                                                         <  1:87556/150‑1 (MQ=255)
                               gTTGGCTAAGATAACGCAGCGCATCGTAGAGGATTGCCTCTTTGTCTGGCCAGAAGCGGCGTAGCTCGTCCAGTGGATAATCCACACGTTCAGCAACCATCTCCAGCGTGGTGTTGGCAATCCCTTGTAATTCTAATAATTTCAGGGCtt                                                                                                                        <  1:3225015/150‑1 (MQ=255)
                                 tGGCTAAGATAACGCAGCGCAACGTAGAGGATTGACTCTTTGTCTGGCCAGAAGCGGCGTAGCTCGTCCAGTGGATAATCCACACGTTCAGCAACCAACTCCAGCGTGGTGTTCGCAATCCATTGTAATTCTAATAATATCATGGCTtca                                                                                                                      >  1:1895574/1‑149 (MQ=255)
                                    cTAAGATAACGCAGCGCATCGTAGAGGATTGCCTCTTTGTCTGGCCAGAAGCGGCGTAGCTCGTCCAGTGGATAATCCACACGTTCAGCAACCATCTCCAGCGTGGTGTTGGCTATCCCTTGTAATTCTAATAATTTCAGGGCTTCTccc                                                                                                                   <  1:869017/150‑1 (MQ=255)
                                    cTAAGATAACGCAGCGCATCGTAGAGGATTGCCTCTTTGTCTGGCCAGAAGCGGCGTAGCTCGTCCAGTGGATAATCCACACGTTCAGCAACCATCTCCAGCGTGGTGTTGGCAATCCCTTGTAATTCTAATAATTTCAGGGCTTCTccc                                                                                                                   >  1:2213034/1‑150 (MQ=255)
                                    cTAAGATAACGCAGCGCATCGTAGAGGATTGCCTCTTTGTCTGGCCAGAAGCGGCGTAGCTCGTCCAGTGGATAATCCACACGTTCAGCAACCATCTCCAGCGTGGTGTTGGCAATCCCTTGTAATTCTAATAATTTCAGGGCTTCTccc                                                                                                                   >  1:2834607/1‑150 (MQ=255)
                                    cTAAGATAACGCAGCGCATCGTAGAGGATTGCCTCTTTGTCTGGCCAGAAGCGGCGTAGCTCGTCCAGTGGATAATCCACACGTTCAGCAACCATCTCCAGCGTGGTGTTGGCAATCCCTTGTAATTCTAATAATTTCAGGGCTTCTccc                                                                                                                   <  2:109241/150‑1 (MQ=255)
                                    cTAAGATAACGCAGCGCATCGTAGAGGATTGCCTCTTTGTCTGGCCAGAAGCGGCGTAGCTCGTCCAGTGGATAATCCACACGTTCAGCAACCATCTCCAGCGTGGTGTTGGCAATCCCTTGTAATTCTAATAATTTCAGGGCTTCTccc                                                                                                                   <  1:3393301/150‑1 (MQ=255)
                                     tAAGATAACGCAGCGCATCGTAGAGGATTGCCTCTTTGTCTGGCCAGAAGCGGCGTAGCTCGTCCAGTGGATAATCCACACGTTCAGCAACCATCTCCAGCGTGGTGTTGGCAATCCCTTGTAATTCTAATAATTTCAGGGCTTCTCCCa                                                                                                                  >  1:2750528/1‑150 (MQ=255)
                                       aGATAACGCAGCGCATCGTAGAGGATTGCCTCTTTGTCTGGCCAGAAGCGGCGTAGCTCGTCCAGTGGATAATCCACACGTTCAGCAACCATCTCCAGCGTGGTGATGGCAATCACTTGTAATTCTAATAATTTCAGGGCTTCTCCCAGt                                                                                                                >  1:4006696/1‑150 (MQ=255)
                                           aaCGCAGCGCATCGTAGAGGATTGCCTCTTTGTCTGGCCAGAAGCGGCGTAGCTCGTCCAGTGGATAATCCACACGTTCAGCAACCATCTCCAGCGTGGTGTTGGCAATCCCTTGTAATTCTAATAATTTCAGGGCTTCTCCCAGTACAt                                                                                                            >  2:465382/1‑150 (MQ=255)
                                              gCAGCGCATCGTAGAGGATTGCCTCTTTGTCTGGCCAGAAGCGGCGTAGCTCGTCCAGTGGATAATCCACACGTTCAGCAACCATCTCCAGCGTGGTGTTGGCAATCCCTTGTAATTCTAATAATTTCAGGGCTTCTCCCAGTACATCtt                                                                                                         <  2:705431/150‑1 (MQ=255)
                                               cAGCGCATCGTAGAGGATTGCCTCTTTGTCTGGCCAGAAGCGGCGTAGCTCGTCCAGTGGATAATCCACACGTTCAGCAACCATCTCCAGCGTGGTGTTGGCAATCCCTTGTAATTCTAATAATTTCAGGGCTTCTCCCAGTACATCTTc                                                                                                        <  1:1464472/150‑1 (MQ=255)
                                                aGCGCATCGTAGAGGATTGCCTCTTTGTCTGGCCAGAAGCGGCGTAGCTCGTCCAGTGGATAATCCACACGTTCAGCAACCATCTCCAGCGTGGTGTTGGCAATCCCTTGTAATTCTAATAATTTCAGGGCTTCTCCCAGTACATCTTCa                                                                                                       <  2:1901130/150‑1 (MQ=255)
                                                aGCGCATCGTAGAGGATTGCCTCTTTGTCTGGCCAGAAGCGGCGTAGCTCGTCCAGTGGATAATCCACACGTTCAGCAACCATCTCCAGCGTGGTGTTGGCAATCCCTTGTAATTCTAATAATTTCAGGGCTTCTCCCAGTACATCTTCa                                                                                                       <  1:3615500/150‑1 (MQ=255)
                                                aGCGCATCGTAGAGGATTGCCTCTTTGTCTAGCCAGAAGCGGCGTAGCTCGTCCAGTGGATAATCCACACGTTCAGCAACCATCTCCAGCGTGGTGTTGGCAATCCCTTGTAATTCTAATAATTTCAGCGCTTCTCCCAGTACATCTTCa                                                                                                       >  2:1327388/1‑150 (MQ=255)
                                                  cgcATCGTAGAGGATTGCCTCTTTGTCTGGCCAGAAGCGGCGTAGCTCGTCCAGTGGATAATCCACACGTTCAGCAACCATCTCCAGCGTGGTGTTGGCAATCCCTTGTAATTCTAATAATTTCAGGGCTTCTCCCAGTACATCTTCANg                                                                                                     <  2:1470771/150‑1 (MQ=255)
                                                  cgcATCGTAGAGGATTGCCTCTTTGTCTGGCCAGAAGCGGCGTAGCTCGTCCAGTGGATAATCCACACGTTCAGCAACCATCTCCAGCGTGGTGTTGGCAATCCCTTGTAATTCTAATAATTTCAGGGCTTCTCCCAGTACATCTTCACg                                                                                                     >  2:4430410/1‑150 (MQ=255)
                                                     aTCGTAGAGGATTGCCTCTTTGTCTGGCCAGAAGCGGCGTAGCTCGACCAGTGGATAATCCACACGTTCAGCAACCATCTCCAGCGTGGTGTTGGCAATCCCTTGTAATTCTAATAATTTCAGGGCTTCTCCCAGTACATCTTCACGTTg                                                                                                  >  2:2547150/1‑150 (MQ=255)
                                                       cGTAGAGGATTGCCTCTTTGTCTGGCCAGAAGCGGCGTAGCTCGTCCAGTGGATAATCCACACGTTCAGCAACCATCTCCAGCGTGGTGTTGGCAATCCCTTGTAATTCTAATAATTTCAGGGCTTCTCCCAGTACATCTTCACGTTGCa                                                                                                >  1:3854728/1‑150 (MQ=255)
                                                        gTAGAGGATTGCCTCTTTGTCTGGCCAGAAGCTGCGTAGCTCGTCCAGTGGATAATCCACACGTTCAGCAACCATCTCCAGCGTGGTGTTGGCAATCCCTTGTAATTCTAATAATTTCAGGGCTTCTCCCAGTACATCTTCACGTTGCAc                                                                                               <  2:2271371/150‑1 (MQ=255)
                                                           gagGATTGCCTCTTTGTCTGGCCAGAAGCGGCGTAGCTCGTCCAGTGGATAATCCACACGTTCAGCAACCATCTCCAGCGTGGTGTTGGCAATCCCTTGTAATTCTAATAATTTCAGGGCTTCTCCCAGTACATCTTCACGTTGCACGCt                                                                                            <  1:480149/150‑1 (MQ=255)
                                                           gagGATTGCCTCTTTGTCTGGCCAGAAGCGGCGTAGCTCGTCCAGTGGATAATCCACACGTTCAGCAACCATCTCCAGCGTGGTGTTGGCAATCCCTTGTAATTCTAATAATTTCAGGGCTTCTCCCAGTACATCTTCACGTTGCACGCt                                                                                            <  2:1682071/150‑1 (MQ=255)
                                                              gATTGCCTCTTTGTCTGGCCAGAAGCGGCGTAGCTCGTCCAGTGGATAATCCACACGTTCAGCAACCATCTCCAGCGTGGTGTTGGCAATCCCTTGTAATTCTAATAATTTCAGGGCTTCTCCCAGTACATCTTCACGTTGCACGCTAtt                                                                                         >  2:4539518/1‑150 (MQ=255)
                                                              gATTGCCTCTTTGTCTGGCCAGAAGCGGCGTAGCTCGTCCAGTGGATAATCCACACGTTCAGCAACCATCTCCAGCGTGGTGTTGGCAATCCCTTGTAATTCTAATAATTTCAGGGCTTCTCCCAGTACATCTTCACGTTGCACGCTAtt                                                                                         >  2:3878539/1‑150 (MQ=255)
                                                               aTTGCCTCTTTGTCTGGCCAGAAGCGGCGTAGCTCGTCCAGTGGATAATCCACACGTTCAGCAACCATCTCCAGCGTGGTGTTGGCAATCCCTTGTAATTCTAATAATTTCAGGGCTTCTCCCAGTACATCTTCACGTTGCACGCTAttt                                                                                        <  2:1623292/150‑1 (MQ=255)
                                                                 tGCCTCTTTGTCTGGCCAGAAGCGGCGTAGCTCGTCCAGTGGATAATCCACACGTTCAGCAACCATCTCCAGCGTGGTGTTGGCAATCCCTTGTAATTCTAATAATTTCAGGGCTTCTCCCAGTACATCTTCACGTTGCACGCTATTTtc                                                                                      >  1:999283/1‑150 (MQ=255)
                                                                  gCCTCTTTGTCTGGCCAGAAGCGGCGTAGCTCGTCCAGTGGATAATCCACACGTTCAGCAACCATCTCCAGCGTGGTGTTGGCAATCCCTTGTAATTCTAATAATTTCAGGGCTTCTCCCAGTACATCTTCACGTTGCACGCTATTTtcc                                                                                     <  1:2958691/150‑1 (MQ=255)
                                                                  gCCTCTTTGTCTGGCCAGAAGCGGCGTAGCTCGTCCAGTGGATAATCCACACGTTCAGCAACCATCTCCAGCGTGGTGTTGGCAATCCCTTGTAATTCTAATAATTTCAGGGCTTCTCCCAGTACATCTTCACGTTGCACGCTATTTtcc                                                                                     <  1:1331286/150‑1 (MQ=255)
                                                                     tctTTGTCTGGCCAGAAGCGGCGTAGCTCGTCCAGTGGATAATCCACACGTTCAGCAACCATCTCCAGCGTGGTGTTGGCAATCCCTTGTAATTCTAATAATTTCAGGGCTTCTCCCAGTACATCTTCACGTTGCACGCTATTTtcctcc                                                                                  <  1:2460381/150‑1 (MQ=255)
                                                                      ctTTGTCTGGCCAGAAGCGGCGTAGCTCGTCCAGTGGATAATCCACACGTTCAGCAACCATCTCCAGCGTGGTGTTGGCAATCCCTTGTAATTCTAATAATTTCAGGGCTTCTCCCAGTACATCTTCACGTTGCACGCTATTTTCCTCCg                                                                                 >  1:33275/1‑150 (MQ=255)
                                                                       tttGTCTGGCCAGAAGCGGCGTAGCTCGTCCAGTGGATAATCCACACGTTCAGCAACCATCTCCAGCGTGGTGTTGGCAATCCCTTGTAATTCTAATAATTTCAGGGCTTCTCCCAGTACATCTTCACGTTGCACGCTATTTTCCTCCGt                                                                                >  1:1231158/1‑150 (MQ=255)
                                                                         tGTCTGGCCAGAAGCGGCGTAGCTCGTCCAGTGGATAATCCACACGTTCAGCAACCATCTCCAGCGTGGTGTTGGCAATCCCTTGTAATTCTAATAATTTCAGGGCTTCTCCCAGTACATCTTCACGTTGCACGCTATTTTCCTCCGTCt                                                                              <  1:3150186/150‑1 (MQ=255)
                                                                          gTCTGGCCAGAAGCGGCGTAGCTCGTCCAGTGGATAATCCACACGTTCAGCAACCATCTCCAGCGTGGTGTTGGCAATCCCTTGTAATTCTAATAATTTCAGGGCTTCTCCCAGTACATCTTCACGTTGCACGCTATTTTCCTCCGTCtt                                                                             <  1:917508/150‑1 (MQ=255)
                                                                           tCTGGCCAGAAGCGGCGTAGCTCGTCCAGTGGATAATCCACACGTTCAGCAACCATCTCCAGCGTGGTGTTGGCAATCCCTTGTAATTCTAATAATTTCAGGGCTTCTCCCAGTACATCTTCACGTTGCTCGCTATTTTCCTCCGTCttt                                                                            >  2:3208458/1‑150 (MQ=255)
                                                                            cTGGCCAGAAGCGGCGTAGCTCGTCCAGTGGATAATCCACACGTTCAGCAACCATCTCCAGCGTGGTGTTGGCAATCCCTTGTAATTCTAATAATTTCAGGGCTTCTCCCAGTACATCTTCACGTTGCACGCTATTTTCCTCCGTCTTTc                                                                           <  2:4316992/150‑1 (MQ=255)
                                                                            cTGGCCAGAAGCGGCGTAGCTCGTCCAGTGGATAATCCACACGTTCAGCAACCATCTCCAGCGTGGTGTTGGCAATCCCTTGTAATTCTAATAATTTCAGGGCTTCTCCCAGTACATCTTCACGTTGCACGCTATTTTCCTCCGTCTTTc                                                                           >  1:3733924/1‑150 (MQ=255)
                                                                              ggCCAGAAGCGGCGTAGCTCGTCCAGTGGATAATCCACACGTTCAGCAACCATCTCCAGCGTGGTGTTGGCAATCCCTTGTAATTCTAATAATTTCAGGGCTTCTCCCAGTACATCTTCACGTTGCACGCTATTTTCCTCCGTCTTTccc                                                                         <  2:1705893/150‑1 (MQ=255)
                                                                                 cAGTAGCGGCGTAGCTCGTCCAGTGGATAATCCACACGTTCAGCAACCATCTCCAGCGTGGTGTTGGCAATCCCTTGTAATTCTAATAATTTCAGGGCTTCTCCCAGTACATCTTCACGTTGCACGCTATTTTCCTCCGTCTTTCCCACt                                                                      >  2:3293598/1‑150 (MQ=255)
                                                                                 cAGAAGCGGCGTAGCTCGTCCAGTGGATAATCCACACGTTCAGCAACCATCTCCAGCGTGGTGTTGGCTATCCCTTGTAATTCTAATAATTTCAGGGCTTCTCCCAGTACATCTTCACGTTGCACGCTATTTTCCTCCGTCTTTCCCACt                                                                      <  1:65037/150‑1 (MQ=255)
                                                                                 cAGAAGCGGCGTAGCTCGTCCAGTGGATAATCCACACGTTCAGCAACCATCTCCAGCGTGGTGTTGGCAATCCCTTGTAATTCTAATAATTTCAGGGCTTCTCCCAGTACATCTTCACGTTGCACGCTATTTTCCTCCGTCTTTCCCACt                                                                      >  1:3500623/1‑150 (MQ=255)
                                                                                   gAAGCGGCGTAGCTCGTCCAGTGGATAATCCACACGTTCAGCAACCATCTCCAGCGTGGTGTTGGCAATCCCTTGTAATTCTAATAATTTCAGGGCTTCTCCCAGTACATCTTCACGTTGCACGCTATTTTCCTCCGTCTTTCCCACTGc                                                                    <  2:4521930/150‑1 (MQ=255)
                                                                                    aaGCGGCGTAGCTCGTCCAGTGGATAATCCACACGTTCAGCAACCATCTCCAGCGTGGTGTTGGCAATCCCTTGTAATTCTAATAATTTCAGGGCTTCTCCCAGTACATCTTCACGTTGCACGCTATTTTCCTCCGTCTTTCCCACTGCa                                                                   <  1:116106/150‑1 (MQ=255)
                                                                                    aaGCGGCGTAGCTCGTCCAGTGGATAATCCACACGTTCAGCAACCATCTCCAGCGTGGTGTTGGCAATCCCTTGTAATTCTAATAATTTCAGGGCTTCTCCCAGTACATCTTCACGTTGCACGCTATTTTCCTCCGTCTTTCCCACTGCa                                                                   >  1:2291648/1‑150 (MQ=255)
                                                                                    aaGCGGCGTAGCTCGTCCAGTGGATAATCCACACGTTCAGCAACCATCTCCAGCGTGGTGTTGGCAATCCCTTGTAATTCTAATAATTTCAGGGCTTCTCCCAGTACATCTTCACGTTGCACGCTATTTTCCTCCGTCTTTCCCACTGCa                                                                   >  1:2780621/1‑150 (MQ=255)
                                                                                        ggcgTAGCTCGTCCAGTGGATAATCCACACGTTCAGCAACCATCTCCAGCGTGGTGTTGGCAATCCCTTGTAATTCTAATAATTTCAGGGCTTCTCCCAGTACATCTTCACGTTGCACGCTATTTTCCTCCGTCTTTCCCACTGCAAgtg                                                               >  1:439810/1‑150 (MQ=255)
                                                                                         gcgTAGCTCGTCCAGTGGATAATCCACACGTTCAGCAACCATCTCCAGCGTGGTGTTGGCAATCCCTTGTAATTCTATGAATTTCAGGGCTTCTCCCAGTACATCTTCACGTTGCACGCTATTTTCCTCCGTCTTTCCCACTGCAAgtgt                                                              <  2:289475/150‑1 (MQ=255)
                                                                                             aGCTCGTCCAGTGGATAATCCACACGTTCAGCAACCATCTCCAGCGTGGTGTTGGCAATCCCTTGTAATTCTAATAATTTCAGGGCTTCTCCCAGTACATCTTCACGTTGCACGCTATTTTCCTCCGTCTTTCCCACTGCAAGTGTCGtt                                                          >  1:1054432/1‑150 (MQ=255)
                                                                                              gCTCGTCCAGTGGATAATCCACACGTTCAGCAACCATCTCCAGCGTGGTGTTGGCAATCCCTTGTAATTCTAATAATTTCAGGGCTTCTCCCAGTACATCTTCACGTTGCACGCTATTTTCCTCCGTCTTTCCCACTGCAAGTGTCGTTc                                                         >  2:1934007/1‑150 (MQ=255)
                                                                                                tCGTCCAGTGGATAATCCACACGTTCAGCAACCATCTCCAGCGTGGTGTTGGCAATCCCTTGTAATTCTAATAATTTCAGGGCTTCTCCCAGTACATCTTCACGTTGCACGCTATTTTCCTCCGTCTTTCCCACTGCAAGTGTCGTTCAc                                                       <  2:1885293/150‑1 (MQ=255)
                                                                                                tCGTCCAGTGGATAATCCACACGTTCAGCAACCATCTCCAGCGTGGTGTTGGCAATCCCTTGTAATTCTAATAATTTCAGGGCTTCTCCCAGTACATCTTCACGTTGCACGCTATTTTCCTCCGTCTTTCCCACTGCAAGTGTCGTTCAc                                                       >  1:2727672/1‑150 (MQ=255)
                                                                                                   tCCAGTGGATAATCCACACGTTCAGCAACCATCTCCAGCGTGGTGTTGGCAATCCCTTGTAATTCTAATAATTTCAGGGCTTCTCCCAGTACATCTTCACGTTGCACGCTATTTTCCTCCGTCTTTCCCACTGCAAGTGTCGTTCACGGt                                                    >  1:3061439/1‑150 (MQ=255)
                                                                                                   tCCAGTGGATAATCCACACGTTCAGCAACCATCTCCAGCGTGGTGTTGGCAATCCCTTGTAATTCTAATAATTTCAGGGCTTCTCCCAGTACATCTTCACGTTGCACGCTATTTTCCTCCGTCTTTCCCACTGCAAGTGTCGTTCACGGt                                                    >  2:3630501/1‑150 (MQ=255)
                                                                                                    ccAGTGGATAATCCACACGTTCAGCAACCATCTCCAGCGTGGTGTTGGCAATCCCTTGTAATTCTAATAATTTCAGGGCTTCTCCCAGTACATCTTCACGTTGCACGCTATTTTCCTCCGTCTTTCCCACTGCAAGTGTCGTTCACGGtt                                                   <  1:1212857/150‑1 (MQ=255)
                                                                                                    ccAGTGGATAATCCACACGTTCAGCAACCATCTCCAGCGTGGTGTTGGCAATCCCTTGTAATTCTAATAATTTCAGGGCTTCTCCCAGTACATCTTCACGTTGCACGCTATTTTCCTCCGTCTTTCCCACTGCAAGTGTCGTTCACGGtt                                                   <  1:1250085/150‑1 (MQ=255)
                                                                                                     cAGTGGATAATCCACACGTTCAGCAACCATCTCCAGCGTGGTGTTGGCAATCCCTTGTAATTCTAATAATTTCAGGGCTTCTCCCAGTACATCTTCACGTTGCACGCTATTTTCCTCCGTCTTTCCCACTGCAAGTGTCGTTCACGGTTg                                                  <  2:191144/150‑1 (MQ=255)
                                                                                                     cAGTGGATAATCCACACGTTCAGCAACCATCTCCAGCGTGGTGTTGGCAATCCCTTGTAATTCTAATAATTTCAGGGCTTCTCCCAGTACATCTTCACGTTGCACGCTATTTTCCTCCGTCTTTCCCACTGCAAGTGTCGTTCACGGTTg                                                  <  2:4490835/150‑1 (MQ=255)
                                                                                                     cAGTGGATAATCCACACGTTCAGCAACCATCTCCAGCGTGGTGTTGGCAATCCCTTGTAATTCTAATAATTTCAGGGCTTCTCCCAGTACATCTTCACGTTGCACGCTATTTTCCTCCGTCTTTCCCACTGCAAGTGTCGTTCACGGTTg                                                  <  1:4124414/150‑1 (MQ=255)
                                                                                                      aGTGGATAATCCACACGTTCAGCAACCATCTCCAGCGTGGTGTTGGCAATCCCTTGTAATTCTAATAATTTCAGGGCTTCTCCCAGTACATCTTCACGTTGCACGCTATTTTCCTCCGTCTTTCCCACTGCAAGTGTCGTTCACGGTTgg                                                 >  2:676551/1‑150 (MQ=255)
                                                                                                           aTAATCCACACGTTCAGCAACCATCTCCAGCGTGGTGTTGGCAATCCCTTGTAATTCTAATAATTTCAGGGCTTCTCCCAGTACATCTTCACGTTGCACGCTATTTTCCTCCGTCTTTCCCACTGCAAGTGTCGTTCACGGTTGGCGATc                                            <  1:3651541/150‑1 (MQ=255)
                                                                                                           aTAATCCACACGTTCAGCAACCATCTCCAGCGTGGTGTTGGCAATCCCTTGTAATTCTAATAATTTCAGGGCTTCTCCCAGTACATCTTCACGTTGCACGCTATTTTCCTCCGTCTTTCCCACTGCAAGTGTCGTTCACGGTTGGCGATc                                            <  1:3651002/150‑1 (MQ=255)
                                                                                                            tAATCCACACGTTCAGCAACCATCTCCAGCGTGGTGTTGGCAATCCCTTGTAATTCTAATAATTTCAGGGCTTCTCCCAGTACATCTTCACGTTGCACGCTATTTTCCTCCGTCTTTCCCACTGCAAGTGTCGTTCACGGTTGGCGATcg                                           >  1:4163061/1‑150 (MQ=255)
                                                                                                            tAATCCACACGTTCAGCAACCATCTCCAGCGTGGTGTTGGCAATCCCTTGTAATTCTAATAATTTCAGGGCTTCTCCCAGTACATCTTCACGTTGCACGCTATTTTCCTCCGTCTTTCCCACTGCAAGTGTCGTTCACGGTTGGCGATcg                                           >  1:17505/1‑150 (MQ=255)
                                                                                                              aTCCACACGTTCAGCAACCATCTCCAGCGTGGTGTTGGCTATCCCTTGTAATTCTAATAATTTCAGGGCTTCTCCCAGTACATCTTCACGTTGCACGCTATTTTCCTCCGTCTTTCCCACTGCAAGTGTCGTTCACGGTTGGCGATcgcg                                         >  1:511390/1‑150 (MQ=255)
                                                                                                                   cacGTTCAGCAACCATCTCCAGCGTGGTGTTGGCAATCCCTTGTAATTCTAATAATTTCAGGGCTTCTCCCAGTACATCTTCACGTTGCACGCTATTTTCCTCCGTCTTTCCCACTGCAAGTGTCGTTCACGGTTGGCGATCGCGCAAAt                                    <  2:300492/150‑1 (MQ=255)
                                                                                                                      gTTCAGCAACCATCTCCAGCGTGGTGTTGGCAATCCCTTGTAATTCTAATAATTTCAGGGCTTCTCCCAGTACATCTTCACGTTGCACGCTATTTTCCTCCGTCTTTCCCACTGCAAGTGTCGTTCACGGTTGGCGATCGCGCAAAtgtg                                 >  1:45646/1‑150 (MQ=255)
                                                                                                                        tCAGCAACCATCTCCAGCGTGGTGTTGGCAATCCCTTGTAATTCTAATAATTTCAGGGCTTCTCCCAGTACATCTTCACGTTGCACGCTATTTTCCTCCGTCTTTCCCACTGCAAGTGTCGTTCACGGTTGGCGATCGCGCAAATGTgcg                               >  2:3944909/1‑150 (MQ=255)
                                                                                                                          aGCAACCATCTCCAGCGTGGTGTTGGCAATCCCTTGTAATTCTAATAATTTCAGGGCTTCTCCCAGTACATCTTCACGTTGCACGCTATTTTCCTCCGTCTTTCCCACTGCAAGTGTCGTTCACGGTTGGCGATCGCGCAAATGTGCGCt                             <  2:4232110/150‑1 (MQ=255)
                                                                                                                           gCAACCATCTCCAGCGTGGTGTTGGCAATCCCTTGTAATTCTAATAATTTCAGGGCTTCTCCCAGTACATCTTCACGTTGCACGCTATTTTCCTCCGTCTTTCCCACTGCAAGTGTCGTTCACGGTTGGCGATCGCGCAAATGTGCGCTg                            >  2:3554363/1‑150 (MQ=255)
                                                                                                                             aaCCATCTCCAGCGTGGTGTTGGCTATCCCTTGTAATTCTAATAATTTCAGGGCTTCTCCCAGTACATCTTCACGTTGCACGCTATTTTCCTCCGTCTTTCCCACTGCAAGTGTCGTTCACGGTTGGCGATCGCGCAAATGTGCGCTGaa                          <  1:1548377/150‑1 (MQ=255)
                                                                                                                              aCCATCTCCAGCGTGGTGTTGGCAATCCCTTGTAATTCTAATAATTTCAGGGCTTCTCCCAGTACATCTTCACGTTGCACGCTATTTTCCTCCGTCTTTCCCACTGCAAGTGTCGTTCACGGTTGGCGATCGCGCAAATGTGCGCTGAAg                         >  2:3053841/1‑150 (MQ=255)
                                                                                                                                cATCTCCCGCGTGGTGTTGGCAATCCCTTGTAATTCTAATAATTTCAGGGCTTCTCCCAGTACATCTTCACGTTGCACGCTATTTTCCTCCGTCTTTCCCACTGCAAGTGTCGTTCACGGTTGGCGATCGCGCAAATGTGCGCTGAAGGt                       <  2:2359730/150‑1 (MQ=255)
                                                                                                                                cATCTCCAGCGTGGTGTTGGCAATCCCTTGTAATTCTAATAATTTCAGGGCTTCTCCCAGTACATCTTCACGTTGCACGCTATTTTCCTCCGTCTTTCCCACTGCAAGTGTCGTTCACGGTTGGCGATCGCGCAAATGTGCGCTGAAGGt                       >  2:2730236/1‑150 (MQ=255)
                                                                                                                                 aTCTCCAGCGTGGTGTTGGCAATCCCTTGTAATTCTAATAATTTCAGGGCTTCTCCCAGTACATCTTCACGTTGCACGCTATTTTCCTCCGTCTTTCCCACTGCAAGTGTCGTTCACGGTTGGCGATCGCGCAAATGTGCGCTGAAGGtt                      <  2:4216/150‑1 (MQ=255)
                                                                                                                                 aTCTCCAGCGTGGTGTTGGCAATCCCTTGTAATTCTAATAATTTCAGGGCTTCTCCCAGTACATCTTCACGTTGCACGCTATTTTCCTCCGTCTTTCCCACTGCAAGTGTCGTTCACGGTTGGCGATCGCGCAAATGTGCGCTGAAGGtt                      <  2:4182/150‑1 (MQ=255)
                                                                                                                                 aTCTCCAGCGTGGTGTTGGCAATCCCTTGTAATTCTAATAATTTCAGGGCTTCTCCCAGTACATCTTCACGTTGCACGCTATTTTCCTCCGTCTTTCCCACTGCAAGTGTCGTTCACGGTTGGCGATCGCGCAAATGTGCGCTGAAGGtt                      <  2:3926620/150‑1 (MQ=255)
                                                                                                                                 aTCTCCAGCGTGGTGTTGGCAATCCCTTGTAATTCTAATAATTTCAGGGCTTCTCCCAGTACATCTTCACGTTGCACGCTATTTTCCTCCGTCTTTCCCACTGCAAGTGTCGTTCACGGTTGGCGATCGCGCAAATGTGCGCTGAAGGtt                      <  2:4367/150‑1 (MQ=255)
                                                                                                                                   ctcCAGCGTGGTGTTGGCTATCCCTTGTAATTCTAATAATTTCAGGGCTTCTCCCAGTACATCTTCACGTTGCACGCTATTTTCCTCCGTCTTTCCCACTGCAAGTGTCGTTCACGGTTGGCGATCGCGCAAATGTGCGCTGAAGGTTTc                    >  2:1102254/1‑150 (MQ=255)
                                                                                                                                   ctcCAGCGTGGTGTTGGCAATCCCTTGTAATTCTAATAATTTCAGGGCTTCTCCCAGTACATCTTCACGTTGCACGCTATTTTCCTCCGTCTTTCCCACTGCAAGTGTCGTTCACGGTTGGCGATCGCGCAAATGTGCGCTGAAGGTTTc                    <  2:2569902/150‑1 (MQ=255)
                                                                                                                                     ccAGCGTGGTGTTGGCAATCCCTTGTAATTCTAATAATTTCAGGGCTTCTCCCAGTACATCTTCACGTTGCACGCTATTTTCCTCCGTCTTTCCCACTGCAAGTGTCGTTCACGGTTGGCGATCGCGCAAATGTGCGCTGAAGGTTTCAg                  <  2:604482/150‑1 (MQ=255)
                                                                                                                                      cAGCGTGGTGTTGGCAATCCCTTGTAATTCTAATAATTTCAGGGCTTCTCCCAGTACATCTTCACGTTGCACGCTATTTTCCTCCGTCTTTCCCACTGCAAGTGTCGTTCACGGTTGGCGATCGCGCAAATGTGCGCTGAAGGTTTCAGc                 <  1:2855000/150‑1 (MQ=255)
                                                                                                                                          gtggtgTTGGCAATCCCTTGTAATTCTAATAATTTCAGGGCTTCTCCCAGTACATCTTCACGTTGCACGCTATTTTCCTCCGTCTTTCCCACTGCAAGTGTCGTTCACGGTTGGCGATCGCGCAAATGTGCGCTGAAGGTTTCAGCATcc             >  1:1782851/1‑150 (MQ=255)
                                                                                                                                               gtTGGCAATCCCTTGTAATTCTAATAATTTCAGGGCTTCTCCCAGTACATCTTCACGTTGCACGCTATTTTCCTCCGTCTTTCCCACTGCAAGTGTCGTTCACGGTTGGCGATCGCGCAAATGTGCGCTGAAGGTTTCAGCATCCATaaa        >  2:3383919/1‑150 (MQ=255)
                                                                                                                                                ttGGCTATCCCTTGTAATTCTAATAATTTCAGGGCTTCTCCCAGTACATCTTCACGTTGCACGCTATTTTCCTCCGTCTTTCCCACTGCAAGTGTCGTTCACGGTTGGCGATCGCGCAAATGTGCGCTGAAGGTTTCAGCATCCATAAAg       >  1:3267132/1‑150 (MQ=255)
                                                                                                                                                  ggCAATCCCTTGTAATTCTAATAATTTCAGGGCTTCTCCCAGTACATCTTCACGTTGCACGCTATTTTCCTCCGTCTTTCCCACTGCAAGAGTCGTTCACGGTTGGCGATCGCGCAAATGAGCGCTGAAGGTTTCAGCATCCATAAAGcc     >  2:3420305/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
AGCATCAGCTGACGACGCCAGACATCTATCTGTTGGCTAAGATAACGCAGCGCATCGTAGAGGATTGCCTCTTTGTCTGGCCAGAAGCGGCGTAGCTCGTCCAGTGGATAATCCACACGTTCAGCAACCATCTCCAGCGTGGTGTTGGCAATCCCTTGTAATTCTAATAATTTCAGGGCTTCTCCCAGTACATCTTCACGTTGCACGCTATTTTCCTCCGTCTTTCCCACTGCAAGTGTCGTTCACGGTTGGCGATCGCGCAAATGTGCGCTGAAGGTTTCAGCATCCATAAAGCCCGT  >  CP009273/4352920‑4353218

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: ‑