Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,532,279 G→A 5.7% intergenic (+55/‑130) cysZ → / → cysK sulfate transporter, sulfite inhibited/cysteine synthase A, O‑acetylserine sulfhydrolase A subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,532,2790GA 5.7% 781.8 / 24.3 227intergenic (+55/‑130)cysZ/cysKsulfate transporter, sulfite inhibited/cysteine synthase A, O‑acetylserine sulfhydrolase A subunit
Reads supporting (aligned to +/- strand):  ref base G (81/133);  new base A (3/10);  total (84/143)
Fisher's exact test for biased strand distribution p-value = 3.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.57e-01

Only 100 of 243 total aligned reads displayed.

GCTGCTTAATCTGTTCATCATGCCCGTTGCCGTTTGTGGCGCGACGGCGATGTGGGTCGATTGCTATCGCGATAAACACGCGATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATCCTTACTTCCGCATATTCTCTGAGCGGGTATGCTACCTGTTGTATCCCAATTTCATACAGTTAAGGACAGGCCATGAGTAAGATTTTTGAAG  >  NC_000913/2532131‑2532427
                                                                                                                                                    |                                                                                                                                                    
gctgctTATTCTGTTCATCATGCCCGTTGGCGTTTGTGGCGCGACGGCGATGTGGGTCGATTGCTATCGCGATAAACACGCGATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATgccgc                                                                                                                                                     <  2:3782957/150‑1 (MQ=255)
 ctgctTAATCTGTTCATCATGCCCGTTGCCGTTTGTGGCGCGACGGCGATGTGGGTCGATTGCTATCGCGATAAACACGCGATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATgccgcc                                                                                                                                                    <  2:674457/150‑1 (MQ=255)
 ctgctTAATCTGTTCATCATGCCCGTTGCCGTTTGTGGCGCGACGGCGATGTGGGTCGATTGCTATCGCGATAAACACGCGATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATgccgcc                                                                                                                                                    <  1:2059258/150‑1 (MQ=255)
  tgctTAATCTGTTCATCATGCCCGTTGCCGTTTGTGGCGCGACGGCGATGTGGGTCGATTGCTATCGCGATAAACACGCGATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGccc                                                                                                                                                   >  1:1382546/1‑150 (MQ=255)
  tgctTAATCTGTTCATCATGCCCGTTGCCGTTTGTGGCGCGACGGCGATGTGGGTCGATTGCTATCGCGATAAACACGCGATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGccc                                                                                                                                                   <  2:3449945/150‑1 (MQ=255)
  tgctTAATCTGTTCATCATGCCCGTTGCCGTTTGTGGCGCGACGGCGATGTGGGTCGATTGCTATCGCGATAAACACGCGATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGCGGCATAAGCCAccc                                                                                                                                                   <  1:3281075/150‑1 (MQ=255)
     ttAATCTGTTCATCATGCCCGTTGCCGTTTGTGGCGCGACGGCGATGTGGGTCGATTGCTATCGCGATAAACACGCGATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCtt                                                                                                                                                >  1:1229001/1‑150 (MQ=255)
     ttAATCTGTTCATCATGCCCGTTGCCGTTTGTGGCGCGACGGCGATGTGGGTCGATTGCTATCGCGATAAACACGCGATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCtt                                                                                                                                                >  1:1947663/1‑150 (MQ=255)
     ttAATCTGTTCATCATGCCCGTTGCCGTTTGTGGCGCGACGGCGATGTGGGTCGATTGCTATCGCGATAAACACGCGATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGCGGCATAAGCCACCCCtg                                                                                                                                                <  1:3102128/150‑2 (MQ=255)
     ttAATCTGTTCATCATGCCCGTTGCCGTTTGTGGCGCGACGGCGATGTGGGTCGATTGCTATCGCGATAAACACGCGATGTGGCGGTAACAATCTACCGGTGATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCtt                                                                                                                                                >  2:2471748/1‑150 (MQ=255)
        aTCTGTTCATCATGCCCGTTGCCGTTTGTGGCGCGACGGCGATGTGGGTCGATTGCTATCGCGATAAACACGCGATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTAtt                                                                                                                                             <  2:322961/150‑1 (MQ=255)
             ttCATCATGCCCGTTGCCGTTTGTGGCGCGACGGCGATGTGGGTCGATTGCTATCGCGATAAACACGCGATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATc                                                                                                                                        >  1:347123/1‑150 (MQ=255)
               catcatGCCCGTTGCCGTTTGTGGCGCGACGGCGATGTGGGTCGATTGCTATCGCGATAAACACGCGATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCtt                                                                                                                                      >  2:3598255/1‑150 (MQ=255)
                   atGCCCGTTGCCGTTTGTGGCGCGACGGCGATGTGGGTCGATTGCTATCGCGATAAACACGCGATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCAt                                                                                                                                  <  2:2344467/150‑1 (MQ=255)
                       ccGTTGCCGTTTGTGGCGCGACGGCGATGTGGGTCGATTGCTATCGCGATAAACACGCGATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCa                                                                                                                              <  2:978801/150‑1 (MQ=255)
                       ccGTTGCCGTTTGTGGCGCGACGGCGATGTGGGTCGATTGCTATCGCGATAAACACGCGATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCa                                                                                                                              <  2:2512451/150‑1 (MQ=255)
                       ccGTAGCCGTTTGTGGCGCGACGGCGATGTGGGTCGATTGCTATCGCGATAAACACGCGATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCa                                                                                                                              <  2:1166212/150‑1 (MQ=255)
                        cGTTTCAGTTTGTGGCGCTACGGCGATGTGTGTCTATTGCTATCGCGATATTCACGCGATGTGGCGGTATCAATCTTACGGTTTTTTTGTAATCCTTTTGTGTGAATCAGTTTTGGCTTATGCCGCGCCTTATTCCATCTTGCATGTCat                                                                                                                             <  1:2518273/150‑1 (MQ=255)
                        cGTTGCCGTTTGTGGCGCGACGGCGATGTGGGTCGATTGCTATCGCGATAAACACGCGATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCat                                                                                                                             <  1:526320/150‑1 (MQ=255)
                         gTTGCCGTTTGTGGCGCGACGGCGATGTGGGTCGATTGCTATCGCGATAAACACGCGATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCatt                                                                                                                            <  2:19163/150‑1 (MQ=255)
                          ttGCCGTTTGTGGCGCGACGGCGATGTGGGTCGATTGCTATCGCGATAAACACGCGATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCatta                                                                                                                           >  2:1788980/1‑150 (MQ=255)
                              cGTTTGTGGCGCGACGGCGATGTGGGTCGATTGCTATCGCGATAAACACGCGATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTc                                                                                                                       <  2:2042603/150‑1 (MQ=255)
                              cGTTTGTGGCGCGACGGCGATGTGGGTCGATTGCTATCGCGATAAACACGCGATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTc                                                                                                                       <  1:1203534/150‑1 (MQ=255)
                               gTTTGTGGCGCGACGGCGATGTGGGTCGATTGCTATCGCGATAAACACGCGATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTcc                                                                                                                      <  1:3734575/150‑1 (MQ=255)
                                tttGTGGCGCGACGGCGATGTGGGTCGATTGCTATCGCGATAAACACGCGATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTTTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATTTCAATATTTacc                                                                                                                     >  1:441103/1‑150 (MQ=255)
                                    tgGCGCGACGGCGATGTGGGTCGATTGCTATCGCGATAAACACGCGATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCt                                                                                                                 <  2:775898/150‑1 (MQ=255)
                                     ggCGCGACGGCGATGTGGGTCGATTGCTATCGCGATAAACACGCGATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGTGTGGCTTATGCCGCCCCTTATTCCTTCTTGCATGTCATTATTTCCCTTCTg                                                                                                                <  1:718725/150‑1 (MQ=255)
                                      gcgcGACGGCGATGTGGGTCGATTGCTATCGCGATAAACACGCGATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGt                                                                                                               <  1:1605880/150‑1 (MQ=255)
                                        gcgACGGCGATGTGGGTCGATTGCTATCGCGATAAACACGCGATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGtat                                                                                                             >  1:3275603/1‑150 (MQ=255)
                                          gACGGCGATGTGGGTCGATTGCTATCGCGATAAACACGCGATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGtatat                                                                                                           <  2:2950039/150‑1 (MQ=255)
                                           aCGGCGATGTGGGTCGATTGCTATCGCGATAAACACGCGATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGtatata                                                                                                          <  1:2670088/150‑1 (MQ=255)
                                             ggCGATGTGGGTCGATTGCTATCGCGATAAACACGCGATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCAACTTGCATGTCATTATTTCCCTTCTGTATATAGa                                                                                                        >  1:2255898/1‑150 (MQ=255)
                                                gATGTGGGTCGATTGCTATCGCGATAAACACGCGATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCGTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGatat                                                                                                     >  1:2448693/1‑150 (MQ=255)
                                                  tgtgGGTCGATTGCTATCGCGATAAACACGCGATGTGGCGGTAATAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGc                                                                                                   <  1:2756377/150‑1 (MQ=255)
                                                  tgtgGGTCGATTGCTATCGCGATAAACACGCGATGTGGCGGTAATAATCTAACGGTTATTTTGTAAACAGTTTGTGTGAAACAGTGGTGGCTTATGCCGCCACTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGc                                                                                                   <  1:764564/150‑1 (MQ=255)
                                                   gtgGGTCGATTGCTATCGCGATAAACACGCGATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCt                                                                                                  <  2:2915722/150‑1 (MQ=255)
                                                        tCGATTGCTATCGCGATAAACACGCGATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATc                                                                                             <  2:1061958/150‑1 (MQ=255)
                                                        tCGATTGCTATCGCGATAAACACGCGATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATc                                                                                             <  1:3099533/150‑1 (MQ=255)
                                                          gATTGCTATCGCGATAAACACGCGATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATCCt                                                                                           <  1:1784895/150‑1 (MQ=255)
                                                              gCTATCGCGATAAACACGCGATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTACTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATCCTTACt                                                                                       <  2:1163232/150‑1 (MQ=255)
                                                                 aTCGCGATAAACACGCGATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATCCTTACTTcc                                                                                    <  2:2618378/150‑1 (MQ=255)
                                                                 aTCGCGATAAACACGCGATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATCCTTACTTcc                                                                                    <  1:3541140/150‑1 (MQ=255)
                                                                  tCGCGATAAACACGCGATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATCCTTACTTCCg                                                                                   >  1:2873241/1‑150 (MQ=255)
                                                                   cgcgATAAACACGCGATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATCCTTACTTCCGc                                                                                  <  2:1881190/150‑1 (MQ=255)
                                                                    gcgATAAACACGCGATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATCCTTACTTCCGCa                                                                                 >  2:258289/1‑150 (MQ=255)
                                                                      gATAAACACGCGATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATCCTTACTTCCGCata                                                                               <  1:3088662/150‑1 (MQ=255)
                                                                       aTAAACACGCGATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATCCTTACTTCCGCatat                                                                              <  1:1978202/150‑1 (MQ=255)
                                                                        tAAACACGCGATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATCCTTACTTCCGCATAtt                                                                             <  2:3310079/150‑1 (MQ=255)
                                                                        tAAACACGCGATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATCCTTACTTCCGCATAtt                                                                             <  2:1590892/150‑1 (MQ=255)
                                                                            cacGCGATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATCCTTACTTCCGCATATTctct                                                                         <  1:2194140/150‑1 (MQ=255)
                                                                             acGCGATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATCCTTACTTCCGCATATTCTCTg                                                                        >  1:1173790/1‑150 (MQ=255)
                                                                              cgcgATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAGTCCTTACTTCCGCATATTCTCTGa                                                                       <  2:1771031/150‑1 (MQ=255)
                                                                              cgcgATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATCCTTACTTCCGCATATTCTCTGa                                                                       <  1:1502126/150‑1 (MQ=255)
                                                                              cgcgATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATCCTTACTTCCGCATATTCTCTGa                                                                       <  2:973619/150‑1 (MQ=255)
                                                                              cgcgATGTCGCGCTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTCAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATCCTTACTTCCGCATATTCTCTGa                                                                       <  2:3516030/150‑1 (MQ=255)
                                                                               gcgATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATCCTTACTTCCGCATATTCTCTGAg                                                                      >  1:2588568/1‑150 (MQ=255)
                                                                               gcgATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATCCTTACTTCCGCATATTCTCTGAg                                                                      >  1:3322963/1‑150 (MQ=255)
                                                                                   tgtgGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATCCTTACTTCCGCATATTCTCTGAGCggg                                                                  <  2:2883874/150‑1 (MQ=255)
                                                                                           tAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATCCTTACTTCCGCATATTCTCTGAGCGGGTATGCTAc                                                          >  1:2860574/1‑150 (MQ=255)
                                                                                                aTCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATCCTTACTTCCGCATATTCTCTGAGCGGGTATGCTACCtgtt                                                     <  1:1480023/150‑1 (MQ=255)
                                                                                                 tCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATCCTTACTTCCGCATATTCTCTGAGCGGGTATGCTACCtgttg                                                    <  1:1379143/150‑1 (MQ=255)
                                                                                                   tACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATCCTTACTTCCGCATATTCTCTGAGCGGGTATGCTACCTGTTGTa                                                  >  1:229122/1‑150 (MQ=255)
                                                                                                   tACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATCCTTACTTCCGCATATTCTCTGAGCGGGTATGCTACCTGTTGTa                                                  >  1:3360841/1‑150 (MQ=255)
                                                                                                     ccGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATCCTTACTTCCGCATATTCTCTGAGCGGGTATGCTACCTGTTGTATc                                                <  2:435016/150‑1 (MQ=255)
                                                                                                      cGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATCCTTACTTCCGCATATTCTCTGAGCGGGTATGCTACCTGTTGTATcc                                               <  1:3545009/150‑1 (MQ=255)
                                                                                                      aggTTATTTTTTAAACCGTTTGTGTGAATCAGGGGTGGCTTATGACGATCGTTTGTGCTTGTTGCTTGTCATTATTTCCCTTCTGTATATATATATGCTATATCCTTAATTCCGCATATTCTGTTAGCGGGTATGCTTGCtgttgtttgc                                               <  2:1824286/149‑5 (MQ=255)
                                                                                                              ttGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATCCTTACTTCCGCATATTCTCTGAGCGGGTATGCTACCTGTTGTATCCCAATTTCa                                       <  2:3632551/150‑1 (MQ=255)
                                                                                                                 tAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATCCTTACTTCCGCATATTCTCTGAGCGGGTATGCTACCTGTTGTATCCCAATTTCATAc                                    >  2:2196275/1‑150 (MQ=255)
                                                                                                                  aaaCCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATCCTTACTTCCGCATATTCTCTGAGCGGGTATGCTACCTGTTGTATCCCAATTTCATACa                                   >  2:748126/1‑150 (MQ=255)
                                                                                                                  aaaCCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATCCTTACTTCCGCATATTCTCTGAGCGGGTATGCTACCTGTTGTATCCCAATTTCATACa                                   <  1:3014799/150‑1 (MQ=255)
                                                                                                                   aaCCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATCCTTACTTCCGCATATTCTCTGAGCGGGTATGCTACCTGTTGTATCCCAATTTCATACAg                                  >  2:3437614/1‑150 (MQ=255)
                                                                                                                       gTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATCCTTACTTCCGCATATTCTCTGAGCGGGTATGCTACCTGTTGTATCCCAATTTCATACAGTTaa                              <  1:849006/150‑1 (MQ=255)
                                                                                                                        tttGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATCCTTACTTCCGCATATTCTCTGAGCGGGTATGCTACCTGTTGTATCCCAATTTCATACAGTTAAg                             >  1:1597939/1‑150 (MQ=255)
                                                                                                                        tttGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATCCTTACTTCCGCATATTCTCTGAGCGGGTATGCTACCTGTTGTATCCCAATTTCATACAGTTAAg                             >  1:1666267/1‑150 (MQ=255)
                                                                                                                           gtgtgAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATCCTTACTTCCGCATATTCTCTGAGCGGGTATGCTACCTGTTGTATCCCAATTTCATACAGTTAAGGAc                          >  2:3806281/1‑150 (MQ=255)
                                                                                                                            tgtgAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATCCTTACTTCCGCATATTCTCTGAGCGGGTATGCTACCTGTTGTATCCCAATTTCATACAGTTAAGGAcc                         >  2:2764718/1‑149 (MQ=255)
                                                                                                                             gtgAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATCCTTACTTCCGCATATTCTCTGAGCGGGTATGCTACCTGTTGTATCCCAATTTCATACAGTTAAGGACAg                        <  2:3124173/150‑1 (MQ=255)
                                                                                                                              tgAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATCCTTACTTCCGCATATTCTCTGAGCGGGTATGCTACCTGTTGTATCCCAATTTCATACAGTTAAGGACAgg                       >  1:1416372/1‑150 (MQ=255)
                                                                                                                              tgAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATCCTTACTTCCGCATATTCTCTGAGCGGGTATGCTACCTGTTGTATCCCAATTTCATACAGTTAAGGACAgg                       <  1:2392187/150‑1 (MQ=255)
                                                                                                                                aaaCAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATCCTTACTTCCGCATAGTCTCTGAGCGGGTATGCTACCTGTTGTATCCCAATTTCATACAGTTAAGGACAGGcc                     >  2:1228394/1‑150 (MQ=255)
                                                                                                                                  taaGGGGCGGCATAAGCCACCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATCCTTACTTCCGCATATTCTCTGAGCGGGTATGCTACCTGTTGTATCCCAATTTCATACAGTTAAGGACAGGCCAt                   >  1:3808253/3‑150 (MQ=255)
                                                                                                                                  taaGGGGCGGCATAAGCCACCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATCCTTACTTCCGCATATTCTCTGAGCGGGTATGCTACCTGTTGTATCCCAATTTCATACAGTTAAGGACAGGCCAt                   >  1:889378/3‑150 (MQ=255)
                                                                                                                                    aGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATCCTTACTTCCGCATATTCTCTGAGCGGGTATGCTACCTGTTGTATCCCAATTTCATACAGTTAAGGACAGGCCATGa                 >  1:7323/1‑150 (MQ=255)
                                                                                                                                    aGGGGCGGCATAAGCCACCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATCCTTACTTCCGCATATTCTCTGAGCGGGTATGCTACCTGTTGTATCCCAATTTCATACAGTTAAGGACAGGCCATGa                 >  1:2957929/1‑150 (MQ=255)
                                                                                                                                      gggTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATTCTAAATACTTACTTCCGCATATTCTCTGAGCGGGTATGCTACCTGTTGTATCCCAATTTCATACAGTTAAGGACAGGCCATGAGt               >  1:761927/1‑150 (MQ=255)
                                                                                                                                      gggTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTACCTTCTGTATATAGATATTCTAAATCCTTACTTCCGCATATTCTCTGAGCGGGTATGCTACCTGTTGTATCCCAATTTCATACAGTTAATGACAGGCCATGAGt               >  1:1121620/1‑150 (MQ=255)
                                                                                                                                       ggTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATCCTTACTTCCGCATATTCTCTGAGCGGGTATGCTACCTGTTGTATCCCAATTTCATACAGTTAAGGACAGGCCATGAGTa              >  2:1305360/1‑150 (MQ=255)
                                                                                                                                           gCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATCCTTACTTCCGCATATTCTCTGAGCGGGTATGCTACCTGTTGTATCCCAATTTCATACAGTTAAGGACAGGCCATGAGTAAGAt          >  2:3690076/1‑150 (MQ=255)
                                                                                                                                             ttATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATCCTTACTTCCGCATATTCTCTGAGCGGGTATGCTACCTGTTGTATCCCAATTTCATACAGTTAAGGACAGGCCATGAGTAAGAttt        >  2:346347/1‑150 (MQ=255)
                                                                                                                                             ttATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATCCTTACTTCCGCATATTCTCTGAGCGGGTATGCTACCTGTTGTATCCCAATTTCATACAGTTAAGGACAGGCCATGAGTAAGAttt        >  1:2028582/1‑150 (MQ=255)
                                                                                                                                              tATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATCCTTACTTCCGCATATTCTCTGAGCGGGTATGCTACCTGTTGTATCCCAATTTCATACAGTTAAGGACAGGCCATGAGTAAGAtttt       <  1:2024369/150‑1 (MQ=255)
                                                                                                                                              tATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATCCTTACTTCCGCATATTCTCTGAGCGGGTATGCTACCTGTTGTATCCCAATTTCATACAGTTAAGGACAGGCCATGAGTAAGAtttt       <  1:1884486/150‑1 (MQ=255)
                                                                                                                                               aTGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATCCTTACTTCCGCATATTCTCTGAGCGGGTATGCTACCTGTTGTATCCCAATTTCATACAGTTAAGGACAGGCCATGAGTAAGAttttt      <  2:3008337/150‑1 (MQ=255)
                                                                                                                                                 gccgccCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATCCTTACTTCCGCATATTCTCTGAGCGGGTATGCTACCTGTTGTATCCCAATTTCATACAGTTAAGGACAGGCCATGAGTAAGATTTTTGa    <  2:2000653/150‑1 (MQ=255)
                                                                                                                                                  ccgccCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATCGATATGCTAAATCCTTACTTCCGCATATTCTCTGAGCGGGTATGCTACCTGTTGTATCCCAATTTCATACAGTTAAGGACAGGCCATGAGTAAGATTTTTGaa   <  1:3917194/150‑1 (MQ=255)
                                                                                                                                                  ccgccCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATCCTTACTTCCGCATATTCTCTGAGCGGGTATGCTACCTGTTGTATCCCAATTTCATACAGTTAAGGACAGGCCATGAGTAAGATTTTTGaa   >  2:1876844/1‑150 (MQ=255)
                                                                                                                                                  ccgccCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATCCTTACTTACGCATATTCTATGAGCGGGTATGCTACCTGTTGTATCCCAAATTCATACAGTTAAGGACAGGCCATGAGTAAGATTTTTTaa   >  1:1052110/1‑150 (MQ=255)
                                                                                                                                                   cgccCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATCCTTACTTCCGCATATTCTCTGAGCGGGTATGCTACCTGTTGTATCCCAATTTCATACAGTTAAGGACAGGCCATGAGTAAGATTTTTGAAg  >  1:1011885/1‑150 (MQ=255)
                                                                                                                                                   cgccCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATCCTTACTTCCGCATATTCTCTGAGCGGGTATGCTACCTGTTGTATCCCAATTTCATACAGTTAAGGACAGGCCATGAGTAAGATTTTTGAAg  >  2:1391435/1‑150 (MQ=255)
                                                                                                                                                   cgccCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATCCTTACTTCCGCATATTCTCTGAGCGGGTATGCTACCTGTTGTATCCCAATTTCATACAGTTAAGGACAGGCCATGAGTAAGATTTTTGAAg  >  1:1009548/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GCTGCTTAATCTGTTCATCATGCCCGTTGCCGTTTGTGGCGCGACGGCGATGTGGGTCGATTGCTATCGCGATAAACACGCGATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATCCTTACTTCCGCATATTCTCTGAGCGGGTATGCTACCTGTTGTATCCCAATTTCATACAGTTAAGGACAGGCCATGAGTAAGATTTTTGAAG  >  NC_000913/2532131‑2532427

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