Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,275,040 C→A 13.2% intergenic (+56/‑295) agaA → / → agaS pseudogene, N‑acetylgalactosamine‑6‑phosphate deacetylase fragment; putative N‑acetylgalactosamine‑6‑phosphate deacetylase/tagatose‑6‑phosphate ketose/aldose isomerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,275,0400CA13.2% 1199.1 / 113.5 368intergenic (+56/‑295)agaA/agaSpseudogene, N‑acetylgalactosamine‑6‑phosphate deacetylase fragment; putative N‑acetylgalactosamine‑6‑phosphate deacetylase/tagatose‑6‑phosphate ketose/aldose isomerase
Reads supporting (aligned to +/- strand):  ref base C (137/182);  new base A (18/31);  total (155/213)
Fisher's exact test for biased strand distribution p-value = 4.41e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.28e-01

Only 100 of 404 total aligned reads displayed.

GCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGATAGTTTGCTCCTTTATTGGGCCTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTCGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTAACGCGTCTTATCCGACCTACAGTTGGTGACCGCAAGGCCGGATAAAGCGTTTGCGCCGCATCCAGCAATCCCTTTTGCTTCCTTTATCTTTTCTTTCAACGATCACAAATTTC  >  CP009273/3274891‑3275189
                                                                                                                                                     |                                                                                                                                                     
gCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGATAGTTTGCTCCTTTATTGGGCCTTCACTTCCCCCGTAAGGCCTTTCTTtttctttc                                                                                                                                                       <  1:1402164/150‑1 (MQ=255)
gCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGATAGTTTGCTCCTTTATTGGGCCTTCACTTCCCCCGTAAGGCCTTTCTTtttctttc                                                                                                                                                       >  1:285716/1‑150 (MQ=255)
 cTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGATAGTTTGCTCCTTTATTGGGCCTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTCg                                                                                                                                                      >  2:354910/1‑150 (MQ=255)
 cTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGATAGTTTGCTCCTTTATTGGGCCTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTCg                                                                                                                                                      >  2:2920477/1‑150 (MQ=255)
   tAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGATAGTTTGCTCCTTTATTGGGCCTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTAGtt                                                                                                                                                    >  2:501428/1‑150 (MQ=255)
       ccGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGATAGTTTGCTCCTTTATTGGGCCTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTCGTTTTGa                                                                                                                                                <  1:4414902/150‑1 (MQ=255)
       ccGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGATAGTTTGCTCCTTTATTGGGCCTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTCGTTTTGa                                                                                                                                                <  1:1861977/150‑1 (MQ=255)
        cGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGATAGTTTGCTCCTTTATTGGGGGTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTCGTTTTGAt                                                                                                                                               >  1:3574474/1‑150 (MQ=255)
          ggCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGATAGTTTGCTCCTTTATTGGGGGTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTCGTTTTGATCt                                                                                                                                             >  1:5025468/1‑150 (MQ=255)
           gCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGATAGTTTGCTCCTTTATTGGGCCTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTCGTTTTGATCtg                                                                                                                                            <  2:2577940/150‑1 (MQ=255)
            cAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGATAGTTTGCTCCTTTATTGGGCCTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTCGTTTTGATCtgt                                                                                                                                           <  2:4935534/150‑1 (MQ=255)
            cAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGATAGTTTGCTCCTTTATTGGGCCTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTCGTTTTGATCtgt                                                                                                                                           <  1:546478/150‑1 (MQ=255)
             aaaCGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGATAGTTTGCTCCTTTATTGGGCCTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTCGTTTTGATCtgtg                                                                                                                                          <  2:5132467/150‑1 (MQ=255)
              aaCGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGATAGTTTGCTCCTTTATTGGGCCTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTCGTTTTGATCTGTGc                                                                                                                                         <  2:2005880/150‑1 (MQ=255)
               aCGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGATAGTTTGCTCCTTTATTGGGCCTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTCGTTTTGATCTGTGCa                                                                                                                                        <  2:4637005/150‑1 (MQ=255)
                 gcgcCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGATAGTTTGCTCCTTTATTGGGCCTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTCGTTTTGATCTGTGCAGc                                                                                                                                      >  2:1243583/1‑150 (MQ=255)
                   gcCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGATAGTTTGCTCCTTTATTGGGCCTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTCGTTTTGATCTGTGCAGCgg                                                                                                                                    <  2:1404836/150‑1 (MQ=255)
                   gcCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGATAGTTTGCTCCTTTATTGGGCCTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTCGTTTTGATCTGTGCAGCgg                                                                                                                                    <  2:3379311/150‑1 (MQ=255)
                     cAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGATAGTTTGCTCCTTTATTGGGGCTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTCGTTTTGATCTGTGCAGCGgtg                                                                                                                                  >  1:2638884/1‑150 (MQ=255)
                     cAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGATAGTTTGCTCCTTTATTGGGCGTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTCGTTTTGATCTGTGCAGCGgtg                                                                                                                                  >  1:4119930/1‑150 (MQ=255)
                      aGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGATAGTTTGCTCCTTTATTGGGGCTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTCGTTTTGATCTGTGCAGCGgtgt                                                                                                                                 >  1:3686094/1‑150 (MQ=255)
                      aGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCAGGATTCAGGGTCAATTAGCTTAGTTTTGATAGTTTGCTCCTTTATTGGGGGTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTCGTTTTGATCTGTGCAGCGgtgt                                                                                                                                 >  1:4544188/1‑150 (MQ=255)
                        cgtcgtTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGATAGTTTGCTCCTTTATTGGGCCTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTCGTTTTGATCTGTGCAGCGGTGTCg                                                                                                                               <  2:1937427/150‑1 (MQ=255)
                         gtcgtTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGATAGTTTGCTCCTTTATTGGGGGTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTCGTTTTGATCTGTGCAGCGGTGTCgg                                                                                                                              >  1:1131768/1‑150 (MQ=255)
                           cgtTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGATAGTTTGCTCCTTTATTGGGCCTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTCGTTTTGATCTGTGCAGCGGTGTCGGAt                                                                                                                            <  2:1762626/150‑1 (MQ=255)
                             ttGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGATAGTTTGCTCCTTTATTGGGCCTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTCGTTTTGATCTGTGCAGCGGTGTCGGATGc                                                                                                                          <  1:3309898/150‑1 (MQ=255)
                                cgcTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGATAGTTTGCTCCTTTATTGGGCCTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTCGTTTTGATCTGTGCAGCGGTGTCGGATGCGAc                                                                                                                       >  2:2786392/1‑150 (MQ=255)
                                cgcTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGATAGTTTGCTCCTTTATTGGGCCTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTCGTTTTGATCTGTGCAGCGGTGTCGGATGCGAc                                                                                                                       <  1:4969196/150‑1 (MQ=255)
                                   tGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGATAGTTTGCTCCTTTATTGGGGCTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTCGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCt                                                                                                                    >  1:2990189/1‑150 (MQ=255)
                                      aTAGCGGGCTACATGTGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGATAGTTTGCTCCTTTATTGGGCCTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTCGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTAAc                                                                                                                 <  1:504790/150‑1 (MQ=255)
                                       tAGCGGGCTACATGTGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGATAGTTTGCTCCTTTATTGGGCCTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTCGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTAAcg                                                                                                                <  2:1469879/150‑1 (MQ=255)
                                        aGCGGGCTACATGTGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGATAGTTTGCTCCTTTATTGGGCCTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTCGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTAAcgc                                                                                                               <  2:3675091/150‑1 (MQ=255)
                                         gCGGTCTACATGTGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGATAGTTTGCTCCTTTATTGGGCCTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTAGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTAAcgcg                                                                                                              <  1:3836298/150‑1 (MQ=255)
                                         gCGGGCTACATGTGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGATAGTTTGCTCCTTTATTGGGCCTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTAGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTAAcgcg                                                                                                              <  1:1264631/150‑1 (MQ=255)
                                                aCATGTGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGATAGTTTGCTCCTTTATTGGGCCTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTCGTTTTTATCTGTGCATCGGTGTCGGATGCGACGCTAACGCGTCTTAta                                                                                                       >  2:4049106/1‑149 (MQ=255)
                                                  aTGTGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGATAGTTTGCTCCTTTATTGGGCCTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTCGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTAACGCGTCTTATCCg                                                                                                     >  1:4793170/1‑150 (MQ=255)
                                                             aTCTGGATTCAGGGTCAATTAGCTTCGTTTTGATAGTTTGCTCCTTTATTGGGCCTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTCGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTAACGCGTCTTATCCGACCTACAGTTg                                                                                          >  1:1757815/1‑150 (MQ=255)
                                                                 ggATTCAGGGTCAATTAGCTTCGTTTTGATAGTTTGCTCCTTTATTGGGCCTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTCGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTAACGCGTCTTATCCGACCTACAGTTGGTGa                                                                                      >  2:4906384/1‑150 (MQ=255)
                                                                   aTTCAGGGTCAATTAGCTTCGTTTTGATAGTTTGCTCCTTTATTGGGCCTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTCGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTAACGCGTCTTATCCGACCTACAGTTGGTGAcc                                                                                    <  2:912759/150‑1 (MQ=255)
                                                                   aTTCAGGGTCAATTAGCTTCGTTTTGATAGTTTGCTCCTTTATTGGGCCTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTCGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTAACGCGTCTTATCCGACCTACAGTTGGTGAcc                                                                                    <  2:1087523/150‑1 (MQ=255)
                                                                   aTTCAGGGTCAATTAGCTTCGTTTTGATAGTTTGCTCCTTTATTGGGCCTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTAGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTAACGCGTCTTATCCGACCTACAGTTGGTGAcc                                                                                    <  1:4531253/150‑1 (MQ=255)
                                                                    ttCAGGGTCAATTAGCTTCGTTTTGATAGTTTGCTCCTTTATTGGGCCTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTCGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTAACGCGTCTTATCCGACCTACAGTTGGTGACCg                                                                                   >  2:3762115/1‑150 (MQ=255)
                                                                     tCAGGGTCAATTAGCTTCGTTTTGATAGTTTGCTCCTTTATTGGGCCTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTCGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTAACGCGTCTTATCCGACCTACAGTTGGTGACCGc                                                                                  <  2:533630/150‑1 (MQ=255)
                                                                     tCAGGGTCAATTAGCTTCGTTTTGATAGTTTGCTCCTTTATTGGGCCTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTCGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTAACGCGTCTTATCCGACCTACAGTTGGTGACCGc                                                                                  >  1:398028/1‑150 (MQ=255)
                                                                       aGGGTCAATTAGCTTCGTTTTGATAGTTTGCTCCTTTATTGGGCCTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTCGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTAACGCGTCTTATCCGACCTACAGTTGGTGACCGCaa                                                                                >  2:3839902/1‑150 (MQ=255)
                                                                        gggTCAATTAGCTTCGTTTTGATAGTTTGCTCCTTTATTGGGCCTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTCGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTAACGCGTCTTATCCGACCTACAGTTGGTGACCGCAAg                                                                               <  2:4357683/150‑1 (MQ=255)
                                                                                tAGCTTCGTTTTGATAGTTTGCTCCTTTATTGGGGGTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTCGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTAACGCGTCTTATCCGACCTACAGTTGGTGACCGCAAGGCCGGATa                                                                       >  1:1556634/1‑150 (MQ=255)
                                                                                    ttCGTTTTGATAGTTTGCTCCTTTATTGGGCCTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTAGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTAACGCGTCTTATCCGACCTACAGTTGGTGACCGCAAGGCCGGATAAAGc                                                                   <  2:1464443/150‑1 (MQ=255)
                                                                                    ttCGTTTTGATAGTTTGCTCCTTTATTGGGCCTTCACTTCCACCGTAAGGCCTTTCTTTTTCTTTCGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTAACGCGTCTTATCCGACCTACAGTTGGTGACCGCAAGGCCGGATAAAGc                                                                   <  1:4645744/150‑1 (MQ=255)
                                                                                        ttttGATAGTTTGCTCCTTTATTGGGCCTTCACGTCCCCCGTAAGGCCTTTCTTTTTCTTTCGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTAACGCGTCTTATCCGACCTACAGTTGGTGACCGCAAGGCCGGATAAAGCGttt                                                               <  1:1990447/150‑1 (MQ=255)
                                                                                          ttGATAGTTTGCTCCTTTATTGGGCCTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTCGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTAACGCGTCTTATCCGACCTACAGTTGGTGACCGCAAGGCCGGATAAAGCGTTTgc                                                             <  2:5101326/150‑1 (MQ=255)
                                                                                          ttGATAGTTTGCTCCTTTATTGGGCCTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTCGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTAACGCGTCTTATCCGACCTACAGTTGGTGACCGCAAGGCCGGATAAAGCGTTTgc                                                             <  1:4206832/150‑1 (MQ=255)
                                                                                           tGATAGTTTGCTCCTTTATTGGGCCTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTCGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTAACGCGTCTTATCCGACCTACAGTTGGTGACCGCAAGGCCGGATAAAGCGGTTgcg                                                            >  2:4482203/1‑150 (MQ=255)
                                                                                             aTAGTTTGCTCCTTTATTGGGCCTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTCGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTAACGCGTCTTATCCGACCTACAGTTGGTGACCGCAAGGCCGGATAAAGCGTTTGcgcc                                                          >  1:4485666/1‑150 (MQ=255)
                                                                                                gTTTGCTCCTTTATTGGGCCTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTCGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTAACGCGTCTTATCCGACCTACAGTTGGTGACCGCAAGGCCGGATAAAGCGTTTGCGCCGCa                                                       >  1:2735105/1‑150 (MQ=255)
                                                                                                   tGCTCCTTTATTGGGCCTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTCGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTAACGCGTCTTATCCGACCTACAGTTGGTGACCGCAAGGCCGGATAAAGCGTTTGCGCCGCATcc                                                    <  1:3535883/150‑1 (MQ=255)
                                                                                                    gCTCCTTTATTGGGCCTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTCGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTAACGCGTCTTATCCGACCTACAGTTGGTGACCGCAAGGCCGGATAAAGCGTTTGCGCCGCATCCa                                                   <  2:3222572/150‑1 (MQ=255)
                                                                                                      tCCTTTATTGGGCCTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTAGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTAACGCGTCTTATCCGACCTACAGTTGGTGACCGCAAGGCCGGATAAAGCGTTTGCGCCGCATCCAGc                                                 >  1:4952997/1‑150 (MQ=255)
                                                                                                       ccTTTATTGGGCCTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTCGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTAACGCGTCTTATCCGACCTACAGTTGGTGACCGCAAGGCCGGATAAAGCGTTTGCGCCGCATCCAGCa                                                <  2:1369870/150‑1 (MQ=255)
                                                                                                         tttATTGGGCCTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTCGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTAACGCGTCTTATCCGACCTACAGTTGGTGACCGCAAGGCCGGATAAAGCGTTTGCGCCGCATCCAGCAAt                                              <  2:1633759/150‑1 (MQ=255)
                                                                                                         tttATTGGGCCTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTAGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTAACGCGTCTTATCCGACCTACAGTTGGTGACCGCAAGGCCGGATAAAGCGTTTGCGCCGCATCCAGCAAt                                              <  1:2455406/150‑1 (MQ=255)
                                                                                                         tttATTGGGCCGTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTCGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTAACGCGTCTTATCCGACCTACAGTTGGTGACCGCAAGGCCGGATAAAGCGTTTGCGCCGCATCCAGCAAt                                              <  1:228493/150‑1 (MQ=255)
                                                                                                          ttATTGGGCCTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTCGATTTGATCTGTGCAGCGGTGTCGGATGCGACGCTAACGCGTCTTATCCGACCTACAGTTGGTGACCGCAAGGCCGGATAAAGCGTTTGCGCCGCATCCAGCAATc                                             <  1:2558138/150‑1 (MQ=255)
                                                                                                          ttATTGGGCCTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTAGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTAACGCGTCTTATCCGACCTACAGTTGGTGACCGCAAGGCCGGATAAAGCGTTTGCGCCGCATCCAGCAATc                                             >  2:2829835/1‑150 (MQ=255)
                                                                                                           tATTGGGCCTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTCGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTAACGCGTCTTATCCGACCTACAGTTGGTGACCGCAAGGCCGGATAAAGCGTTTGCGCCGCATCCAGCAATcc                                            >  1:1918625/1‑150 (MQ=255)
                                                                                                             ttGGGCCTTCACTTCCCCCGTCAGGCCTTTCTTTTTCTTTCGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTAACGCGTCTTATCCGACCTACAGTTGGTGACCGCAAGGCCGGATAAAGCGTTTGCGCCGCATCCAGCAATCCCt                                          >  2:2371879/1‑150 (MQ=255)
                                                                                                             ttGGGCCTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTCGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTAACGCGTCTTATCCGACCTACAGTTGGTGACCGCAAGGCCGGATAAAGCGTTTGCGCCGCATCCAGCAATCCCt                                          <  1:997121/150‑1 (MQ=255)
                                                                                                             ttGGGCCTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTCGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTAACGCGTCTTATCCGACCTACAGTTGGTGACCGCAAGGCCGGATAAAGCGTTTGCGCCGCATCCAGCAATCCCt                                          <  2:3558837/150‑1 (MQ=255)
                                                                                                              tGGGCCTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTCGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTAACGCGTCTTATCCGACCTACAGTTGGTGACCGCAAGGCCGGATAAAGCGTTTGCGCCGCATCCAGCAATCCCtt                                         <  1:3205818/150‑1 (MQ=255)
                                                                                                              tGGGCCTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTCGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTAACGCGTCTTATCCGACCTACAGTTGGTGACCGCAAGGCCGGATAAAGCGTTTGCGCCGCATCCAGCAATCCCtt                                         >  2:2521584/1‑150 (MQ=255)
                                                                                                                tcccTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTCGATTTGATCTGTGCAGCGTTGTCGGATGCTACGCTAACGCGTCTTATCCGACCTACAGTTGGTGACCGCAAGGCCGGATAAATCGTTTGCGCCGCATCCAGCAATCCCtttt                                       >  2:476743/3‑150 (MQ=255)
                                                                                                                ggCCTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTCGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTAACGCGTCTTATCCGACCTACAGTTGGTGACCGCAAGGCCGGATAAAGCGTTTGCGCCGCATCCAGCAATCCCtttt                                       <  1:4806022/150‑1 (MQ=255)
                                                                                                                ggCCTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTCGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTAACGCGTCTTATCCGACCTACAGTTGGTGACCGCAAGGCCGGATAAAGCGTTTGCGCCGCATCCAGCAATCCCtttt                                       <  2:4418037/150‑1 (MQ=255)
                                                                                                                  ccTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTCGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTAACGCGTCTTATCCGACCTACAGTTGGTGACCGCAAGGCCGGATAAAGCGTTTGCGCCGCATCCAGCAATCCCTTTTGc                                     >  1:1529417/1‑150 (MQ=255)
                                                                                                                    ttCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTAGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTAACGCGTCTTATCCGACCTACAGTTGGTGACCGCAAGGCCGGATAAAGCGTTTGCGCCGCATCCAGCAATCCCTTTTGCtt                                   <  1:819867/150‑1 (MQ=255)
                                                                                                                        cTTCCCCCGTAAGGCCTTTCTTTTTCTTTCGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTAACGCGTCTTATCCGACCTACAGTTGGTGACCGCAAGGCCGGATAAAGCGTTTGCGCCGCATCCAGCAATCCCTTTTGCTTCCtt                               <  1:676144/150‑1 (MQ=255)
                                                                                                                         ttCCCCCGTAAGGCCTTTCTTTTTCTTTCGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTAACGCGTCTTATCCGACCTACAGTTGGTGACCGCAAGGCCGGATAAAGCGTTTGCGCCGCATCCAGCAATCCCTTTTGCTTCCttt                              >  2:5044305/1‑150 (MQ=255)
                                                                                                                          tcacccGTAAGGCCTTTCTTTTTCTTTCGTTTTTATCTGTGCATCGGTGTCGGATGAGACGCTAACGCGTCTTATCCGACATAAAGTTGGTGACCGCAAGGCCGGATAAAGCGTTTGCGCCGCATCCAGCAATCCCTTTTGCTTCCTTTa                             >  2:3458120/4‑150 (MQ=255)
                                                                                                                           cccccGTAAGGCCTTTCTTTTTCTTTCGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTAACGCGTCTTATCCGACCTACAGTTGGTGACCGCAAGGCCGGATAAAGCGTTTGCGCCGCATCCAGCAATCCCTTTTGCTTCCTTTAt                            <  2:2638884/150‑1 (MQ=255)
                                                                                                                                gTAAGGCCTTTCTTTTTCTTTCGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTAACGCGTCTTATCCGACCTACAGTTGGTGACCGCAAGGCCGGATAAAGCGTTTGCGCCGCATCCAGCAATCCCTTTTGCTTCCTTTAtctttt                       >  1:2431898/1‑150 (MQ=255)
                                                                                                                                gTAAGGCCTTTCTTTTTCTTTCGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTAACGCGTCTTATCCGACCTACAGTTGGTGACCGCAAGGCCGGATAAAGCGTTTGCGCCGCATCCAGCAATCCCTTTTGCTTCCTTTAtctttt                       <  1:1631327/150‑1 (MQ=255)
                                                                                                                                  aaGGCCTTTCTTTTTCTTTCGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTAACGCGTCTTATCCGACCTACAGTTGGTGACCGCAAGGCCGGATAAAGCGTTTGCGCCGCATCCAGCAATCCCTTTTGCTTCCTTTAtcttttct                     <  2:1543762/150‑1 (MQ=255)
                                                                                                                                   aGGCCTTTCTTTTTCTTTCGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTAACGCGTCTTATCCGACCTACAGTTGGTGACCGCAAGGCCGGATAAAGCGTTTGCGCCGCATCCAGCAATCCCTTTTGCTTCCTTTAtcttttctt                    >  2:3256005/1‑150 (MQ=255)
                                                                                                                                    ggCCTTTCTTTTTCTTTCGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTAACGCGTCTTATCCGACCTACAGTTGGTGACCGCAAGGCCGGATAAAGCGTTTGCGCCGCATCCAGCAATCCCTTTTGCTTCCTTTAtcttttcttt                   >  1:1046741/1‑150 (MQ=255)
                                                                                                                                    ggCCTTTCTTTTTCTTTAGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTAACGCGTCTTATCCGACCTACAGTTGGTGACCGCAAGGCCGGATAAAGCGTTTGCGCCGCATCCAGCAATCCCTTTTGCTTCCTTTAtcttttcttt                   <  1:1316928/150‑1 (MQ=255)
                                                                                                                                      ccTTTCTTTTTCTTTCGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTAACGCGTCTTATCCGACCTACAGTTGGTGACCGCAAGGCCGGATAAAGCGTTTGCGCCGCATCCAGCAATCCCTTTTGCTTCCTTTATCTTTTCTTTCa                 <  2:3968043/150‑1 (MQ=255)
                                                                                                                                       ctttctttTTCTTTCGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTAACGCGTCTTATCCGACCTACAGTTGGTGACCGCAAGGCCGGATAAAGCGTTTGCGCCGCATCCAGCAATCCCTTTTGCTTCCTTTATCTTTTCTTTCaa                <  1:3954662/150‑1 (MQ=255)
                                                                                                                                        tttctttttcttTCGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTAACGCGTCTTATCCGACCTACAGTTGGTGACCGCAAGGCCGGATAAAGCGTTTGCGCCGCATCCAGCAATCCCTTTTGCTTCCTTTATCTTTTCTTTCAAc               >  1:4190572/1‑150 (MQ=255)
                                                                                                                                           ctttttctttCGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTAACGCGTCTTATCCGACCTACAGTTGGTGACCGCAAGGCCGGATAAAGCGTTTGCGCCGCATCCAGCAATCCCTTTTGCTTCCTTTATCTTTTCTTTCAACGAt            <  2:3188880/150‑1 (MQ=255)
                                                                                                                                           ctttttctttCGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTAACGCGTCTTATCCGACCTACAGTTGGTGACCGCAAGGCCAGATAAAGCGTTTGCGCCGCATCCAGCAATCCCTTTTGCTTCCTTTATCTTTTCTTTCAACGAt            <  2:808748/150‑1 (MQ=255)
                                                                                                                                           ctttttctttCGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTAACGCGTCTTATCCGACCTACAGTTGGTGACCGCAAGGCCAGATAAAGCGTTTGCGCCGCATCCAGCAATCCCTTTTGCTTCCTTTATCTTTTCTTTCAACGAt            <  2:3386131/150‑1 (MQ=255)
                                                                                                                                            tttttctttCGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTAACGCGTCTTATCCGACCTACAGTTGGTGACCGCAAGGCCGGATAAAGCGTTTGCGCCGCATCCAGCAATCCCTTTTGCTTCCTTTATCTTTTCTTTCAACGATc           >  2:2987994/1‑150 (MQ=255)
                                                                                                                                             ttttctttCGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTAACGCGTCTTATCCGACCTACAGTTGGTGACCGCAAGGCCGGATAAAGCGTTTGCGCCGCATCCAGCAATCCCTTTTGCTTCCTTTATCTTTTCTTTCAACGATca          >  2:2866992/1‑150 (MQ=255)
                                                                                                                                              tttctttcGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTAACGCGTCTTATCCGACCTACAGTTGGTGACCGCAAGGCCGGATAAAGCGTTTGCGCCGCATCCAGCAATCCCTTTTGCTTCCTTTATCTTTTCTTTCAACGATcac         >  2:4886312/1‑150 (MQ=255)
                                                                                                                                                 ctttcGTTTTGATCTTTGCAGCGGTGTCGGATGCGTCGCTAACTCGTATTATCCGACCTACAGTTGGTGACCGCAAGGACGGATAAAGCGTTTGCGACGCATCCAGCAATCCCTTTTGCTTCCTTTATCTTTTCTTTCAACGATCACaaa      <  2:648631/150‑1 (MQ=255)
                                                                                                                                                 ctttcGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTAACGCGTCTTATCCGACCTACAGTTGGTGACCGCAAGGCCGGATAAAGCGTTTGCGCCGCATCCAGCAATCCCTTTTGCTTCCTTTATCTTTTCTTTCAACGATCACaaa      >  1:3946301/1‑150 (MQ=255)
                                                                                                                                                  tttcGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTAACGCGTCTTATCCGACCTACAGTTGGTGACCGCAAGGCCGGATAAAGCGTTTGCGCCGCATCCAGCAATCCCTTTTGCTTCCTTTATCTTTTCTTTCAACGATCACAAAt     <  2:4197570/150‑1 (MQ=255)
                                                                                                                                                  tttcGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTAACGCGTCTTATCCGACCTACAGTTGGTGACCGCAAGGCCGGATAAAGCGTTTGCGCCGCATCCAGCAATCCCTTTTGCTTCCTTTATCTTTTCTTTCAACGATCACAAAt     <  1:2137152/150‑1 (MQ=255)
                                                                                                                                                  tttaGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTAACGCGTCTTATCCGACCTACAGTTGGTGACCGCAAGGCCGGATAAAGCGTTTGCGCCGCATCCAGCAATCCCTTTTGCTTCCTTTATCTTTTCTTTCAACGATCACAAAt     <  1:2931345/150‑1 (MQ=255)
                                                                                                                                                   ttcGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTAACGCGTCTTATCCGACCTACAGTTGGTGACCGCAAGGCCGGATAAAGCGTTTGCGCCGCATCCAGCAATCCCTTTTGCTTCCTTTATCTTTTCTTTCAACGATCACAAAtt    <  2:4870256/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGATAGTTTGCTCCTTTATTGGGCCTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTCGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTAACGCGTCTTATCCGACCTACAGTTGGTGACCGCAAGGCCGGATAAAGCGTTTGCGCCGCATCCAGCAATCCCTTTTGCTTCCTTTATCTTTTCTTTCAACGATCACAAATTTC  >  CP009273/3274891‑3275189

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: