Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,946,256 A→T 10.9% Q138L (CAG→CTG)  ilvE → branched‑chain amino‑acid aminotransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,946,2560AT10.9% 1186.9 / 105.8 395Q138L (CAG→CTG) ilvEbranched‑chain amino‑acid aminotransferase
Reads supporting (aligned to +/- strand):  ref base A (161/191);  new base T (25/18);  total (186/209)
Fisher's exact test for biased strand distribution p-value = 1.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.30e-01

Only 100 of 398 total aligned reads displayed.

CAAAAACAATCTCACCAGCGCCTATATCCGTCCGCTGATCTTCGTCGGTGATGTTGGCATGGGAGTAAACCCGCCAGCGGGATACTCAACCGACGTGATTATCGCTGCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAGCGCGCCGCCACGGTTATCAGGAAGGTATCGCGCTGGATG  >  CP009273/3946107‑3946405
                                                                                                                                                     |                                                                                                                                                     
  aaaaCAATCTCACCAGCGCCTATATCCGTCCGCTGATCTTCGTCGGTGATGTTGGCATGGGAGTAAACCCGCCAGCGGGATACTCAACCGACGTGATTATCGCTGCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCTgg                                                                                                                                                     >  2:1127435/1‑150 (MQ=255)
  aaaaCAATCTCACCAGCGCCTATATCCGTCCGCTGATCTTCGTCGGTGATGTTGGCATGGGAGTAAACCCGCCAGCGGGATACTCAACCGACGTGATTATCGCTGCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAgg                                                                                                                                                     >  1:3923357/1‑150 (MQ=255)
    aaCAATCTCACCAGCGCCTATATCCGTCCGCTGATCTTCGTCGGTGATGTTGGCATGGGAGTAAACCCGCCAGCGGGATACTCAACCGACGTGATTATCGCTGCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAgggg                                                                                                                                                   >  1:339620/1‑150 (MQ=255)
     aCAATCTCACCAGCGCCTATATCCGTCCGCTGATCTTCGTCGGTGATGTTGGCATGGGAGTAAACCCGCCAGCGGGATACTCAACCGACGTGATTATCGCTGCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCTGGGGa                                                                                                                                                  >  1:3737873/1‑150 (MQ=255)
       aaTCTCACCAGCGCCTATATCCGTCCGCTGATCTTCGTCGGTGATGTTGGCATGGGAGTAAACCCGCCAGCGGGATACTCAACCGACGTGATTATCGCTGCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATc                                                                                                                                                <  2:4495487/150‑1 (MQ=255)
       aaTCTCACCAGCGCCTATATCCGTCCGCTGATCTTCGTCGGTGATGTTGGCATGGGAGTAAACCCGCCAGCGGGATACTCAACCGACGTGATTATCGCTGCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATc                                                                                                                                                <  1:548992/150‑1 (MQ=255)
        aTCTCACCAGCGCCTATATCCGTCCGCTGATCTTCGTCGGTGATGTTGGCATGGGAGTAAACCCGCCAGCGGGATACTCAACCGACGTGATTATCGCTGCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATcg                                                                                                                                               >  2:5009684/1‑150 (MQ=255)
        aTCTCACCAGCGCCTATATCCGTCCGCTGATCTTCGTCGGTGATGTTGGCATGGGAGTAAACCCGCCAGCGGGATACTCAACCGACGTGATTATCGCTGCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATcg                                                                                                                                               <  1:276423/150‑1 (MQ=255)
         tctcACCAGCGCCTATATCCGTCCGCTGATCTTCGTCGGTGATGTTGGCATGGGAGTAAACCCGCCAGCGGGATACTCAACCGACGTGATTATCGCTGCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATcga                                                                                                                                              <  2:4339583/150‑1 (MQ=255)
         tctcACCAGCGCCTATATCCGTCCGCTGATCTTCGTCGGTGATGTTGGCATGGGAGTAAACCCGCCAGCGGGATACTCAACCGACGTGATTATCGCTGCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATcga                                                                                                                                              >  1:4626472/1‑150 (MQ=255)
          ctcACCAGCGCCTATATCCGTCCGCTGATCTTCGTCGGTGATGTTGGCATGGGAGTAAACCCGCCAGCGGGATACTCAACCGACGTGATTATCGCTGCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATcgat                                                                                                                                             >  1:1081535/1‑150 (MQ=255)
           tcACCAGCGCCTATATCCGTCCGCTGATCTTCGTCGGTGATGTTGGCATGGGAGTAAACCCGCCAGCGGGATACTCAACCGACGTGATTATCGCTGCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATcgatg                                                                                                                                            <  2:576884/150‑1 (MQ=255)
           tcACCAGCGCCTATATCCGTCCGCTGATCTTCGTCGGTGATGTTGGCATGGGAGTAAACCCGCCAGCGGGATACTCAACCGACGTGATTATCGCTGCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATcgatg                                                                                                                                            <  2:2555751/150‑1 (MQ=255)
            cACCAGCGCCTATATCCGTCCGCTGATCTTCGTCGGTGATGTTGGCATGGGAGTAAACCCGCCAGCGGGATACTCAACCGACGTGATTATCGCTGCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATcgatgc                                                                                                                                           >  1:3009983/1‑150 (MQ=255)
              ccAGCGCCTATATCCGTCCGCTGATCTTCGTCGGTGATGTTGGCATGGGAGTAAACCCGCCAGCGGGATACTCAACCGACGTGATTATCGCTGCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATcgatgcga                                                                                                                                         <  1:1180622/150‑1 (MQ=255)
                   gcCTATATCCGTCCGCTGATCTTCGTCGGTGATGTTGGCATGGGAGTAAACCCGCCAGCGGGATACTCAACCGACGTGATTATCGCTGCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGtt                                                                                                                                    >  2:127479/1‑150 (MQ=255)
                   gcCTATATCCGTCCGCTGATCTTCGTCGGTGATGTTGGCATGGGAGTAAACCCGCCAGCGGGATACTCAACCGACGTGATTATCGCTGCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGtt                                                                                                                                    >  2:1479714/1‑150 (MQ=255)
                       atatCCGTCCGCTGATCTTCGTCGGTGATGTTGGCATGGGAGTAAACCCGCCAGCGGGATACTCAACCGACGTGATTATCGCTGCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTtcct                                                                                                                                <  2:4852006/150‑1 (MQ=255)
                       atatCCGTCCGCTGATCTTCGTCGGTGATGTTGGCATGGGAGTAAACCCGCCAGCGGGATACTCAACCGACGTGATTATCGCTGCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTtcct                                                                                                                                <  1:3258208/150‑1 (MQ=255)
                            cgtccgCTGATCTTCGTCGGTGATGTTGGCATGGGAGTAAACCCGCCAGCGGGATACTCAACCGACGTGATTATCGCTGCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTgg                                                                                                                           >  2:313034/1‑150 (MQ=255)
                                cgCTGATCTTCGTCGGTGATGTTGGCATGGGAGTAAACCCGCCAGCGGGATACTCAACCGACGTGATTATCGCTGCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAAcc                                                                                                                       <  1:1704282/150‑1 (MQ=255)
                                     aTCTTCGTCGGTGATGTTGGCATGGGAGTAAACCCGCCAGCGGGATACTCAACCGACGTGATTATCGCTGCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCgca                                                                                                                  >  1:4701893/1‑150 (MQ=255)
                                       cTTCGTCGGTGATGTTGGCATGGGAGTAAACCCGCCAGCGGGATACTCAACCGACGTGATTATCGCTGCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCgcagc                                                                                                                <  2:1330809/150‑1 (MQ=255)
                                         tcgtcgGTGATGTTGGCATGGGAGTAAACCCGCCAGCGGGATACTCAACCGACGTGATTATCGCTGCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCAc                                                                                                              >  1:3193831/1‑150 (MQ=255)
                                          cgtcgGTGATGTTGGCATGGGAGTAAACCCGCCAGCGGGATACTCAACCGACGTGATTATCGCTGCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCTGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCAcc                                                                                                             >  1:1962718/1‑150 (MQ=255)
                                          cgtcgGTGATGTTGGCATGGGAGTAAACCCGCCAGCGGGATACTCAAACGAAGTGATTATCGCTGCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCTGGGGATCGATGCGATGGTTTCCTCCTGGAACAGCGCaccacc                                                                                                             >  2:4671323/1‑150 (MQ=255)
                                              ggTGATGTTGGCATGGGAGTAAACCCGCCAGCGGGATACTCAACCGACGTGATTATCGCTGCTTTCCCGTGGGGAGCGTATCTGGGTGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAac                                                                                                         >  2:279337/1‑150 (MQ=255)
                                              ggTGATGTTGGCATGGGAGTAAACCCGCCAGCGGGATACTCAACCGACGTGATTATCGCTGCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAac                                                                                                         <  2:3071232/150‑1 (MQ=255)
                                               gTGATGTTGGCATGGGAGTAAACCCGCCAGCGGGATACTCAACCGACGTGATTATCGCTGCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCTGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAaca                                                                                                        >  2:1816102/1‑150 (MQ=255)
                                                tGATGTTGGCATGGGAGTAAACCCGCCAGCGGGATACTCAACCGACGTGATTATCGCTGCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAacac                                                                                                       <  2:1643667/150‑1 (MQ=255)
                                                   tGTTGGCATGGGAGTAAACCCGCCAGCGGGATACTCAACCGACGTGATTATCGCTGCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCAt                                                                                                    <  1:2704618/150‑1 (MQ=255)
                                                    gTTGGCATGGGAGTAAACCCGCCAGCGGGATACTCAACCGACGTGATTATCGCTGCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATc                                                                                                   <  1:2253429/150‑1 (MQ=255)
                                                    gTTGGCATGGGAGTAAACCCGCCAGCGGGATACTCAACCGACGTGATTATCGCTGCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATc                                                                                                   <  1:153755/150‑1 (MQ=255)
                                                     ttGGCATGGGAGTAAACCCGCCAGCGGGATACTCAACCGACGTGATTATCGCTGCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATcc                                                                                                  <  2:1812498/150‑1 (MQ=255)
                                                     ttGGCATGGGAGTAAACCCGCCAGCGGGATACTCAACCGACGTGATTATCGCTGCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATcc                                                                                                  >  1:4440505/1‑150 (MQ=255)
                                                       ggCATGGGAGTAAACCCGCCAGCTGGATACTCAACCGACGTGATTATCGCTGCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCg                                                                                                <  1:1825569/150‑1 (MQ=255)
                                                       ggCATGGGAGTAAACCCGCCAGCGGGATACTCAACCGACGTGATTATCGCTGCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCg                                                                                                <  1:4069198/150‑1 (MQ=255)
                                                       ggCATGGGAGTAAACCCGCCAGCGGGATACTCAACCGACGTGATTATCGCTGCTTTCCCGTGGGGAGCGGATCTGGGCGCAGAAGCGCTGGAGCAGGGGAGCGCTGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCCTCCCg                                                                                                >  2:3268537/1‑150 (MQ=255)
                                                          aTGGGAGTAAACCCGCCAGCGGGATACTCAACCGACGTGATTATCGCTGCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGAcg                                                                                             <  1:204469/150‑1 (MQ=255)
                                                            gggAGGAAACCCGCCAGCGGGATACTCAACCGACGTGATTATCGCTGCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGAcggc                                                                                           <  2:2284585/150‑1 (MQ=255)
                                                             ggAGTAAACCCGCCAGCGGGATACTCAACCGACGTGATTATCGCTGCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGAcggcg                                                                                          >  2:541756/1‑150 (MQ=255)
                                                                gTAAACCCGCCAGCGGGATACTCAACCGACGTGATTATCGCTGCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCa                                                                                       >  1:1609287/1‑150 (MQ=255)
                                                                  aaaCCCGCCAGCGGGATACTCAACCGACGTGATTATCGCTGCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCaaa                                                                                     <  2:4334940/150‑1 (MQ=255)
                                                                    aCCCGCCAGCGGGATACTCAACCGACGTGATTATCGCTGCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAg                                                                                   <  1:2069937/150‑1 (MQ=255)
                                                                     cccGCCAGCGGGATACTCAACCGACGTGATTATCGCTGCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGc                                                                                  <  2:919399/150‑1 (MQ=255)
                                                                       cGCCAGCGGGATACTCAACCGACGTGATTATCGCTGCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCg                                                                                >  1:897786/1‑150 (MQ=255)
                                                                        gCCAGCGGGATACTCAACCGACGTGATTATCGCTGCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCgg                                                                               >  1:2947298/1‑150 (MQ=255)
                                                                        gCCAGCGGGATACTCAACCGACGTGATTATCGCTGCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCgg                                                                               >  2:3470063/1‑150 (MQ=255)
                                                                          cAGCGGGATACTCAACCGACGTGATTATCGCTGCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCTGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCggtg                                                                             >  2:3631311/1‑150 (MQ=255)
                                                                          cAGCGGGATACTCAACCGACGTGATTATCGCTGCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCTGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCggtg                                                                             <  2:288057/150‑1 (MQ=255)
                                                                          cAGCGGGATACTCAACCGACGTGATTATCGCTGCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCggtg                                                                             <  1:28630/150‑1 (MQ=255)
                                                                          cAGCGGGATACTCAACCGACGTGATTATCGCTGCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCggtg                                                                             <  1:1012958/150‑1 (MQ=255)
                                                                           aGCGGGATACTCAACCGACGTGATTATCGCTGCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCggtgg                                                                            >  2:1077841/1‑150 (MQ=255)
                                                                              gggATACTCAACCGACGTGATTATCGCTGCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTaa                                                                         >  1:4310140/1‑150 (MQ=255)
                                                                                 aTACTCAACCGACGTGATTATCGCTGCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCTATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTa                                                                      <  1:710006/150‑1 (MQ=255)
                                                                                     tCAACCGACGTGATTATCGCTGCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCTGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACctc                                                                  >  2:2414310/1‑150 (MQ=255)
                                                                                     tCAACCGACGTGATTATCGCTGCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACctc                                                                  <  2:2187996/150‑1 (MQ=255)
                                                                                      cAACCGACGTGATTATCGCTGCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACctct                                                                 <  2:1251911/150‑1 (MQ=255)
                                                                                         ccGACGTGATTATCGCTGCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCtt                                                              <  2:288551/150‑1 (MQ=255)
                                                                                          cGACGTGATTATCGCTGCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTc                                                             >  2:4613052/1‑150 (MQ=255)
                                                                                           gACGTGATTATCGCTGCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTcc                                                            >  2:3911802/1‑150 (MQ=255)
                                                                                             cGTGATTATCGCTGCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGACAAAAGCCGGTGGTAACTACCTCTCTTCCct                                                          >  2:2088930/1‑150 (MQ=255)
                                                                                              gTGATTATCGCTGCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCctg                                                         >  1:3837217/1‑150 (MQ=255)
                                                                                              gTGATTATCGCTGCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCctg                                                         <  1:4792431/150‑1 (MQ=255)
                                                                                                gATTATCGCTGCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCctgct                                                       <  2:3559415/150‑1 (MQ=255)
                                                                                                gATTATCGCTGCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCctgct                                                       <  1:3142922/150‑1 (MQ=255)
                                                                                                        ctgctTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTa                                               <  1:1366614/150‑1 (MQ=255)
                                                                                                        ctgctTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTa                                               <  1:4236733/150‑1 (MQ=255)
                                                                                                         tgctTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCTGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAg                                              <  2:4700073/150‑1 (MQ=255)
                                                                                                            tttCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGa                                           <  1:4160307/150‑1 (MQ=255)
                                                                                                             ttCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGaa                                          <  2:1388966/150‑1 (MQ=255)
                                                                                                             ttCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGaa                                          >  2:2711659/1‑150 (MQ=255)
                                                                                                               cccGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGACgc                                        >  2:808946/1‑150 (MQ=255)
                                                                                                               acagtgtGGATCGTATCTGTGCGCAGAATCGCTTGAGCATGGGATCGATTCGATGGTTTCCTCCTGGAACCGCTCATCACCAAACACAATCCCGACGGCGGCAAAAGCCTGTGTTAACTACCTCTCTTCACTGCTGGTGGGTATCGAAgc                                        >  2:3009910/4‑150 (MQ=255)
                                                                                                                ccGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAgcg                                       >  1:1592434/1‑150 (MQ=255)
                                                                                                                 cGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAgcgc                                      >  2:4078205/1‑150 (MQ=255)
                                                                                                                 cGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAgcgc                                      >  1:595859/1‑150 (MQ=255)
                                                                                                                  gTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAgcgcg                                     <  1:5102496/150‑1 (MQ=255)
                                                                                                                  gTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAgcgcg                                     <  2:3769217/150‑1 (MQ=255)
                                                                                                                   tGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCTGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGACgcgcgc                                    >  2:3379099/1‑150 (MQ=255)
                                                                                                                    ggggAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAGCGcgcc                                   <  2:2433674/150‑1 (MQ=255)
                                                                                                                      ggAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAGCGcgccgc                                 >  1:4479977/1‑150 (MQ=255)
                                                                                                                        aGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAGCGCGCCGCCa                               <  2:414669/150‑1 (MQ=255)
                                                                                                                         gCGTATCTGGGCGCAGAAGCGCTGGAGCTGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAGCGCGCCGCCAc                              <  2:2347295/150‑1 (MQ=255)
                                                                                                                          cGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAGCGCGCCGCCACg                             <  2:1720550/150‑1 (MQ=255)
                                                                                                                          cGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAGCGCGCCGCCACg                             <  2:3026034/150‑1 (MQ=255)
                                                                                                                           gTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAGCGCGCCGCCACgg                            >  1:1681247/1‑150 (MQ=255)
                                                                                                                            tATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAGCGCGCCGCCACGGt                           <  1:5068483/150‑1 (MQ=255)
                                                                                                                             aTCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAGCGCGCCGCCACTGtt                          <  2:1916057/150‑1 (MQ=255)
                                                                                                                                tGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAGCGCGCCGCCACGGTTATc                       >  1:742322/1‑150 (MQ=255)
                                                                                                                                     gcAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAGCGCGCCGCCACGGTTATCAGGaa                  <  1:4105899/150‑1 (MQ=255)
                                                                                                                                      cAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAGCGCGCCGCCACGGTTATCAGGAAg                 >  1:4564972/1‑150 (MQ=255)
                                                                                                                                       aGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAGCGCGCCGCCACGGTTATCAGGAAgg                >  1:1335843/1‑150 (MQ=255)
                                                                                                                                         aaGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAGCGCGCCGCCACGGTTATCAGGAAGGTa              >  2:4923432/1‑150 (MQ=255)
                                                                                                                                            cgcTGGAGCTGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAGCGCGCCGCCACGGTTATCAGGAAGGTATcg           >  1:3588472/1‑150 (MQ=255)
                                                                                                                                            cgcTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAGCGCGCCGCCACGGTTATCAGGAAGGTATcg           >  1:2464793/1‑150 (MQ=255)
                                                                                                                                              cTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAGCGCGCCGCCACGGTTATCAGGAAGGTATcgcg         <  2:3716191/150‑1 (MQ=255)
                                                                                                                                               tGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAGCGCGCCGCCACGGTTATCAGGAAGGTATCgcgc        >  1:4943993/1‑150 (MQ=255)
                                                                                                                                                ggAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAGCGCGCCGCCACGGTTATCAGGAAGGTATCGCGCt       <  2:3193831/150‑1 (MQ=255)
                                                                                                                                                 gAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAGCGCGCCGCCACGGTTATCAGGAAGGTATCGCGCTg      <  1:3854405/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CAAAAACAATCTCACCAGCGCCTATATCCGTCCGCTGATCTTCGTCGGTGATGTTGGCATGGGAGTAAACCCGCCAGCGGGATACTCAACCGACGTGATTATCGCTGCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAGCGCGCCGCCACGGTTATCAGGAAGGTATCGCGCTGGATG  >  CP009273/3946107‑3946405

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: