Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,191,887 C→T 6.5% G250E (GGA→GAA)  thiG ← 1‑deoxy‑D‑xylulose 5‑phosphate:thiol sulfurtransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,191,8870CT 6.5% 484.9 / 13.9 156G250E (GGA→GAA) thiG1‑deoxy‑D‑xylulose 5‑phosphate:thiol sulfurtransferase
Reads supporting (aligned to +/- strand):  ref base C (75/69);  new base T (8/2);  total (84/72)
Fisher's exact test for biased strand distribution p-value = 1.09e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.68e-01

Only 100 of 158 total aligned reads displayed.

GCCATCATGTCGTCGCGGGTGAGTTGCGAGGCATTTAGCGCCCGCTCTACGTCAGCAGCCGTTTTGCCGTTGATACGCAGGCGGATGTCGTCCCAGTCCAGTTGTCGCCAGCGATCGCTGAAGGTTTTCATGCCGATGCCTCCAGAAATCCGGTCAGCGGGCTGGTGGCATGAGCAAAATAACTGCGGCTGCCCGGTCCGGACTGACGTGCCAGTAGGCCTGCTTCTACCGCCAGACGAAATGCCTTCGCCATGTTGACGGGATCGTCCGCGACGGCAATCGCCGTATTCACTAACA  >  NC_000913/4191738‑4192034
                                                                                                                                                     |                                                                                                                                                   
gcCATCATGTCGTCGCGGGTGAGTTGCGAGGCATTTAGCGCCCGCTCTACGTCAGCAGCCGTTTTGCCGTTGATACGCAGGCGGATGTCGTCCCAGTCCAGTTGTCGCCAGCGATCGCTGAAGGTTTTCATGCCGATGCCTCCAGAAATc                                                                                                                                                     <  1:1173751/150‑1 (MQ=255)
   atcatGTCGTCGCGGGTGAGTTGCGAGGCATTTAGCGCCCGCTCTACGTCAGCAGCCGTTTTGCCGTTGATACGCAGGCGGATGTCGTCCCAGTCCAGTCGTCGCCAGCGATCGCTGAAGGTTTTCATGCCGATGCCTCCAGAAATCCgg                                                                                                                                                  <  2:1194030/150‑1 (MQ=255)
          gggcgcgGGTGAGTTGCGAGGCATTTAGCGCCCGCTCTACGTCAGCAGCCGTTTTGCTGTTGATACGCAGGCGAATGTAGTCCCAGTCCAGTTGTCGCCAGCGATCGCAGAAGCTTTTCATGCCGATGCCTCCAGAAATTCGCTCAGCgg                                                                                                                                           >  2:82977/4‑150 (MQ=255)
          cgtcgCGGGTGAGTTGCGAGGCATTTAGCGCCCGCTCTACGTCAGCAGCCGTTTTGCCGTTGATACGCAGGCGGATGTCGTCCCAGTCCAGTTGTCGCCAGCGATCGCTGAAGGTTTTCATGCCGATGCCTCCAGAAATCCGGTCAGCgg                                                                                                                                           <  2:954984/150‑1 (MQ=255)
            tcgCGGGTGAGTTGCGAGGCATTTAGCGCCCGCTCTACGTCAGCAGCCGTTTTGCCGTTGATACGCAGGCGGATGTCGTCCCAGTCCAGTTGTCGCCAGCGATCGCTGAAGGTTTTCATGCCGATGCCTCCAGAAATCCGGTCAGCGGGc                                                                                                                                         <  1:723673/150‑1 (MQ=255)
            gggcgGGTGAGTTGCGAGGCATTTAGCGCCCGCTCTACGTCAGCAGCCGTTTTGCTGTTGATACGCAGGCGAATGTCGTCCCAGTCCAGTTGTCGCCAGCGATCGCTGAAGGTTTTCATGCCGATGCCTCCAGAAATTCGGTCAGCGGGc                                                                                                                                         >  2:695447/3‑150 (MQ=255)
            cggggggTGAGTTGCGAGGCATTTAGCGCCCGCTCTACGTCAGCAGCCGTTTTGCTGTTGATACGCAGGCGAATGTCGTCCCAGTCCAGTTGTCGCCAGCGATCGCTGAAGGTTTTCATGCCGATGCCTCCAGAAATTCGGTCAGCgggg                                                                                                                                         >  2:2206344/5‑149 (MQ=255)
               cgGGTGAGTTGCGAGGCATTTAGCGCCCGCTCTACGTCAGCAGCCGTTTTGCCGTTGATACGCAGGCGGATGTCGTCCCAGTCCAGTTGTCGCCAGCGATCGCTGAAGGTTTTCATGCCGATGCCTCCAGAAATCCGGTCAGCGGGCtgg                                                                                                                                      <  2:1026295/150‑1 (MQ=255)
               cgGGTGAGGTGCGAGGCATTTAGCGCCCGCTCTACGTCAGCAGCCGTTTTGCCGTTGATACGCAGGCGGATGTCGTCCCAGTCCAGTTGTCGCCAGCGATCGCTGAAGGTTTTCATGCCGATGCCTCCAGAAATCCGGTCAGCGGGCtgg                                                                                                                                      <  1:2078967/150‑1 (MQ=255)
                gggTGAGTTGCGAGGCATTTAGCGCCCGCTCTACGTCAGCAGCCGTTTTGCCGTTGATACGCAGGCGGATGTCGTCCCAGTCCAGTTGTCGCCAGCGATCGCTGAAGGTTTTCATGCCGATGCCTCCAGAAATCCGGTCAGCGGGCtggt                                                                                                                                     >  2:41162/1‑150 (MQ=255)
                 ggTGAGTTGCGAGGCATTTAGCGCCCGCTCTACGTCAGCAGCCGTTTTGCCGTTGATACGCAGGCGGATGTCGTCCCAGTCCAGTTGTCGCCAGCGATCGCTGAAGGTTTTCATGCCGATGCCTCCAGAAATCCGGTCAGCGGGCtggtg                                                                                                                                    <  2:1467101/150‑1 (MQ=255)
                   tGAGTTGCGAGGCATTTAGCGCCCGCTCTACGTCAGCAGCCGTTTTGCCGTTGATACGCAGGCGGATGTCGTCCCAGTCCAGTTGTCGCCAGCGATCGCTGAAGGTTTTCATGCCGATGCCTCCAGAAATCCGGTCAGCGGGCTGGTGGc                                                                                                                                  <  1:772291/150‑1 (MQ=255)
                     aGTTGCGAGGCATTTAGCGCCCGCTCTACGTCAGCAGCCGTTTTGCCGTTGATACGCAGGCGGATGTCGTCCCAGTCCAGTTGTCGCCAGCGATCGCTGAAGGTTTTCATGCCGATGCCTCCAGAAATCCGGTCAGCGGGCTGGTGGCAt                                                                                                                                >  2:1548934/1‑150 (MQ=255)
                     aGTTGCGAGGCATTTAGCGCCCGCGCTACGTCAGCAGCCGTTTTGCCGTTGATACGCAGGCGGATGTCGTCCCAGTCCAGTTGTCGCCAGCGATCGCTGAAGGTTTTCATGCCGATGCCTCCAGAAATCCGGGCAGCGGGCTGGTGGCAt                                                                                                                                <  2:154663/150‑1 (MQ=255)
                      gTTGTGAGGCATTTAGCGCCCGCTCTACGTCAGCAGCCGTTTTGCCGTTGATACGCAGGCGGATGTCGTCCCAGTCCAGTTGTCGCCAGCGATCGCTGAAGGTTTTCATGCCGATGCCTCCAGAAATCCGGTCAGCGGGCTGGTGGCATg                                                                                                                               <  1:2484227/150‑1 (MQ=255)
                      gTTGCGAGGCATTTAGCGCCCGCTCTACGTCAGCAGCCGTTTTGCCGTTGATACGCAGGCGGATGTCGTCCCAGTCCAGTTGTCGCCAGCGATCGCTGAAGGTTTTCATGCCGATGCCTCCAGAAATCCGGTCAGCGGGCTGGTGGCATg                                                                                                                               >  1:2011562/1‑150 (MQ=255)
                        tGCGAGGCATTTAGCGCCCGCTCTACGTCAGCAGCCGTTTTGCCGTTGATACGCAGGCGGATGTCGTCCCAGTCCAGTTGTCGCCAGCGATCGCTGAAGGTTTTCATGCCGATGCCTCCAGAAATCCGGTCAGCGGGCTGGTGGCATGAg                                                                                                                             >  2:1482596/1‑150 (MQ=255)
                        tGCGAGGCATTTAGCGCCCGCTCTACGTCAGCAGCCGTTTTGCCGTTGATACGCAGGCGGATGTCGTCCCAGTCCAGTTGTCGCCAGCGATCGCTGAAGGTTTTCATGCCGATGCCTCCAGAAATCCGGTCAGCGGGCTGGTGGCATGAg                                                                                                                             <  2:379938/150‑1 (MQ=255)
                         gCGAGGCATTTAGCGCCCGCTCTACGTCAGCAGCCGTTTTGCCGTTGATACGCAGGCGGATGTCGTCCCAGTCCAGTTGTCGCCAGCGATCGCTGAAGGTTTTCATGCCGATGCCTCCAGAAATCCGGTCAGCGGGCTGGTGGCATGAGc                                                                                                                            >  2:1686230/1‑150 (MQ=255)
                           gAGGCATTTAGCGCCCGCTCTACGTCAGCAGCCGTTTTGCCGTTGATACGCAGGCGGATGTCGTCCCAGTCCAGTTGTCGCCAGCGATCGCTGAAGGTTTTCATGCCGATGCCTCCAGAAATCCGGTCAGCGGGCTGGTGGCATGAGCaa                                                                                                                          <  2:459982/150‑1 (MQ=255)
                            aGGCATTTAGCGCCCGCTCTACGTCAGCAGCCGTTTTGCTGTTGATACGCAGGCGAATGTCGTCCCAGTCCAGTTGTCGCCAGCGATCGCTGAAGGTTTTCATGCCGATGCCTCCAGAAATTCGGTCAGCGGGCTGGTGGCATGAGCaaa                                                                                                                         <  1:695447/150‑1 (MQ=255)
                             ggCATTTAGCGCCCGCTCTACGTCAGCAGCCGTTTTGCCGTTGATACGCAGGCGGATGTCGTCCCAGTCCAGTTGTCGCCAGCGATCGCTGAAGGTTTTCATGCCGATGCCTCCGGAAATCCGGTCAGCGGGCTGGTGGCATGAGCaaaa                                                                                                                        <  2:1959079/150‑1 (MQ=255)
                              gCATTTAGCGCCCGCTCTACGTCAGCAGCCGTTTTGCCGTTGATACGCAGGCGGATGTCGTCCCAGTCCAGTTGTCGCCAGCGATCGCTGAAGGTTTTCATGCCGATGCCTCCAGAAATCCGGTCAGCGGGCTGGTGGCATGAGCAAAAt                                                                                                                       >  2:2001800/1‑150 (MQ=255)
                                 tttAGCGCCCGCTCTACGTCAGCAGCCGTTTTGCCGTTGATACGCAGGCGGATGTCGTCCCAGTCCAGTTGTCGCCAGCGATCGCTGAAGGTTTTCATGCAGATGCCTCCAGAAATCCGGTCAGCGGGCTGGTGGCATGAGCAAAATAAc                                                                                                                    <  1:2411854/150‑1 (MQ=255)
                                    aGCGCCCGCTCTACGTCAGCAGCCGTTTTGCCGTTGATACGCAGGCGGATGTCGTCCCAGTCCAGTTGTCGCCAGCGATCGCTGAAGGTTTTCATGCCGATGCCTCCAGAAATCCGGTCAGCGGGCTGGTGGCATGAGCAAAATAACTGc                                                                                                                 >  1:421669/1‑150 (MQ=255)
                                    aGCGCCCGCTCTACGTCAGCAGCCGTTTTGCCGTTGATACGCAGGCGGATGTCGTCCCAGTCCAGTTGTCGCCAGCGATCGCTGAAGGTTTTCATGCCGATGCCTCCAGAAATCCGGTCAGCGGGCTGGTGGCATGAGCAAAATAACTGc                                                                                                                 >  1:2259278/1‑150 (MQ=255)
                                       gcCCGCTCTACGTCAGCAGCCGTTTTGCCGTTGATACGCAGGCGGATGTCGTCCCAGTCCAGTTGTCGCCAGCGATCGCTGAAGGTTTTCATGCCGATGCCTCCAGAAATCCGGTCAGCGGGCTGGTGGCATGAGCAAAATAACTGCGGc                                                                                                              <  1:1928910/150‑1 (MQ=255)
                                       gcCCGCTCTACGTCAGCAGCCGTTTTGCCGTTGATACGCAGGCGGATGCCGTCCCAGTCCAGTTGTCGCCAGCGATCGCTGAAGGTTTTCATGCCGATGCCTCCAGAAATCCGGTCAGCGGGCTGGTGGCATGAGCAAAATAACTGCGGc                                                                                                              <  2:1905266/150‑1 (MQ=255)
                                       gcCCGCTCTACGTCAGCAGCCGGTTTGCCGTTGATACGCAGGCGGATGTCGGCCCAGTCCAGTTGTCGCCAGCGATCGCTGAAGGTTTTCATGCCGATGCCTCCAGAAATCCGGTCAGCGGGCTGGTGGCATGAGCAAAATAACTGCGGc                                                                                                              <  1:2427057/150‑1 (MQ=255)
                                          cGCTCTACGTCAGCAGCCGTTTTGCCGTTGATACGCAGGCGGATGTCGTCCCAGTCCAGTTGTCGCCAGCGATCGCTGAAGGTTTTCATGCCGATGCCTCCAGAAATCCGGTCAGCGGGCTGGTGGCATGAGCAAAATAACTGCGGCTGc                                                                                                           >  2:807728/1‑150 (MQ=255)
                                            ctctACGTCAGCAGCCGTTTTGCCGTTGATACGCAGGCGGATGTCGTCCCAGTCCAGTTGTCGCCAGCGATCGCTGAAGGTTTTCATGCCGATGCCTCCAGAAATCCGGTCAGCGGGCTGGTGGCATGAGCAAAATAACTGCGGCTGccc                                                                                                         >  1:1578917/1‑150 (MQ=255)
                                             tctACGTCAGCAGCCGTTTTGCCGTTGATACGCAGGCGGATGTCGTCCCAGTCCAGTTGTCGCCAGCGATCGCTGAAGGTTTTCATGCCGATGCCTCCAGAAATCCGGTCAGCGGGCTGGTGGCATGAGCAAAATAACTGCGGCTGCCCg                                                                                                        <  2:900691/150‑1 (MQ=255)
                                              ctACGTCAGCAGCCGTTTTGCCGTTGATACGCAGGCGGATGTCGTCCCAGTCCAGTTGTCGCCAGCGATCGCTGAAGGTTTTCATGCCGATGCCTCCAGAAATCCGGTCAGCGGGCTGGTGGCATGAGCAAAATAACTGCGGCTGCCCgg                                                                                                       <  1:234048/150‑1 (MQ=255)
                                                aCGTCAGCAGCCGTTTTGCCGTTGATACGCAGGCGGATGTCGTCCCAGTCCAGTTGTCGCCAGCGATCGCTGAAGGTTTTCATGCCGATGCCTCCAGAAATCCGGTCAGCGGGCTGGTGGCATGAGCAAAATAACTGCGGCTGCCCGGTc                                                                                                     <  1:150884/150‑1 (MQ=255)
                                                 cGTCAGCAGCCGTTTTGCCGTTGATACGCAGGCGGATCTCGTCCCAGTCCAGTTGTCGCCAGCGATCGCTGAAGGTTTTCATGCCGATGCCTCCAGAAATCCGGTCAGCGGGCTGGTGGCATGAGCAAAATAACTGCGGCTGCCCGGTcc                                                                                                    >  2:687508/1‑150 (MQ=255)
                                                     agcagcCGTTTTGCCGTTGATACGCAGGCGGATGTCGTCCCAGTCCAGTTGTCGCCAGCGATCGCTGAAGGTTTTCATGCCGATGCCTCCAGAAATCCGGTCAGCGGGCTGGTGGCATGAGCAAAATAACTGCGGCTGCCCGGTCCGGAc                                                                                                <  2:346407/150‑1 (MQ=255)
                                                      gcagcCGTTTTGCCGTTGATACGCAGGCGGATGTCGTCCCAGTCCAGTTGTCGCCAGCGATCGCTGAAGGTTTTCATGCCGATGCCTCCGGAAATCCGGTCAGCGGGCTGGTGGCATGAGCAAAATAACTGCGGCTGCCCGGTCCGGACt                                                                                               >  1:1512736/1‑150 (MQ=255)
                                                        agcCGTTTTGCCGTTGATACGCAGGCGGATGTCGTCCCAGTCCAGTTGTCGCCAGCGATCGCTGAAGGTTTTCATGCCGATGCCTCCAGAAATCCGGTCAGCGGGCTGGCGGCATGAGCAAAATAACTGCGGCTGCCCGGTCCGGACTGa                                                                                             <  2:2224630/150‑1 (MQ=255)
                                                          ccGTTTTGCCGTTGATACGCAGGCGGATGTCGTCCCAGTCCAGTTGTCGCCAGCGATCGCTGAAGGTTTTCATGCCGATGCCTCCAGAAATCCGGTCAGCGGGCTGGTGGCATGAGCAAAATAACTGCGGCTGCCCGGTCCGGACTGACg                                                                                           <  2:1952519/150‑1 (MQ=255)
                                                            gTTTTGCCGTTGATACGCAGGCGGATGTCGTCCCAGTCCAGTTGTCGCCAGCGATCGCTGAAGGTTTTCATGCCGATGCCTCCAGAAATCCGGTCAGCGGGCTGGTGGCATGAGCAAAATAACTGCGGCTGCCCGGTCCGGACTGACGTg                                                                                         >  1:1067485/1‑150 (MQ=255)
                                                            gTTTTGCCGTTGATACGCAGGCGGATGTCGTCCCAGTCCAGTTGTCGCCAGCGATCGCTGAAGGTTTTCATGCCGATGCCTCCAGAAATCCGGTCAGCGGGCTGGTGGCATGAGCAAAATAACTGCGGCTGCCCGGTCCGGACTGACGTg                                                                                         >  1:1359265/1‑150 (MQ=255)
                                                                tGCCGTTGATACGCAGGCGGATGTCGTCCCAGTCCAGTTGTCGCCAGCGATCGCTGAAGGTTTTCATGCCGATGCCTCCAGAAATCCGGTCAGCGGGCTGGTGGCATGAGCAAAATAACTGCGGCTGCCCGGTCCGGACTGACGTGCCAg                                                                                     >  1:678291/1‑150 (MQ=255)
                                                                tGCCGTTGATACGCAGGCGGATGTCGTCCCAGTCCAGTTGTCGCCAGCGATCGCTGAAGGTTTTCATGCCGATGCCTCCAGAAATCCGGTCAGCGGGCTGGTGGCATGAGCAAAATAACTGCGGCTGCCCGGTCCGGACTGACGTGCCAg                                                                                     >  2:1855722/1‑150 (MQ=255)
                                                                tGCCGTTGATACGCAGGCGGATGTCGTCCCAGTCCAGTTGTCGCCAGCGATCGCTGAAGGTTTTCATGCCGATGCCTCCAGAAATCCGGTCAGCGGGCTGGTGGCATGAGCAAAATAACTGCGGCTGCCCGGTCCGGACTGACGTGCCAg                                                                                     <  1:643705/150‑1 (MQ=255)
                                                                tGCCGTTGATACGCAGGCGGATGTCGTCCCAGTCCAGTTGTCGCCAGCGATCGCTGAAGGTTTTCATGCCGATGCCTCCAGAAATCCGGTCAGCGGGCTGGTGGCATGAGCAAAATAACTGCGGCTGCCCGGTCCGGACTGACGTGCCAg                                                                                     >  1:1671088/1‑150 (MQ=255)
                                                                 gCCGTTGATACGCAGGCGGATGTCGTCCCAGTCCAGTTGTCGCCAGCGATCGCTGAAGGTTTTCATGCCGATGCCTCCAGAAATTCGGTCAGCGGGCTGGTGGCATGAGCAAAATAACTGCGGCTGCCCGGTCCGGACTGACGTGCCAGt                                                                                    >  2:1353156/1‑150 (MQ=255)
                                                                 gCCGTTGATACGCAGGCGGATGTCGTCCCAGTCCAGTTGTCGCCAGCGATCGCTGAAGGTTTTCATGCCGATGCCTCCAGAAATCCGGTCAGCGGGCTGGTGGCATGAGCAAAATAACTGCGGCTGCCCGGTCCGGACTGACGTGCCAGt                                                                                    <  2:2257161/150‑1 (MQ=255)
                                                                    gTTGATACGCAGGCGGATGTCGTCCCAGTCCAGTTGTCGCCAGCGATCGCTGAAGGTTTTCATGCCGATGCCTCCAGAAATCCGGTCAGCGGGCTGGTGGCATGAGCAAAATAACTGCGGCTGCCCGGTCCGGACTGACGTGCCAGTAgg                                                                                 >  2:638927/1‑150 (MQ=255)
                                                                     ttGATACGCAGGCGGATGTCGTCCCAGTCCAGTTGTCGCCAGCGATCGCTGAAGGTTTTCATGCCGATGCCTCCAGAAATCCGGTCAGCGGGCTGGTGGCATGAGCAAAATAACTGCGGCTGCCCGGTCCGGACTGACGTGCCAGTAGGc                                                                                <  1:2001800/150‑1 (MQ=255)
                                                                      tGATACGCAGGCGGATGTCGTCCCAGTCCAGTTGTCGCCAGCGATCGCTGAAGGTTTTCATGCCGATGCCTCCAGAAATCCGGTCAGCGGGCTGGTGGCATGAGCAAAATAACTGCGGCTGCCCGGTCCGGACTGACGTGCCAGTAGGcc                                                                               <  2:1449272/150‑1 (MQ=255)
                                                                           cGCAGGCGGATGTCGTCCCAGTCCAGTTGTCGCCAGCGATCGCTGAAGGTTTTCATGCCGATGCCTCCAGAAATCCGGTCAGCGGGCTGGTGGCATGAGCAAAATAACTGCGGCTGCCCGGTCCGGACTGACGTGCCAGTAGGCCTGCtt                                                                          >  1:1115252/1‑150 (MQ=255)
                                                                           cGCAGGCGGATGTCGTCCCAGTCCAGTTGTCGCCAGCGATCGCTGAAGGTTTTCATGCCGATGCCTCCAGAAATCCGGTCAGCGGGCTGGTGGCATGAGCAAAATAACTGCGGCTGCCCGGTCCGGACTGACGTGCCAGTAGGCCTGCtt                                                                          >  1:1029973/1‑150 (MQ=255)
                                                                                   gATGTCGTCCCAGTCCAGTTGTCGCCAGCGATCGCTGAAGGTTTTCATGCCGATGCCTCCAGAAATCCGGTCAGCGGGCTGGTGGCATGAGCAAAATAACTGCGGCTGCCCGGTCCGGACTGACGTGCCAGTAGGCCTGCTTCTACCGcc                                                                  >  2:1221628/1‑150 (MQ=255)
                                                                                     tggcgtcCCAGTCCAGTTGTCGCCAGCGATCGCTGAAGGTTTTCATGCCGATGCCTCCAGAAATTCGGTCAGCGGGCTGGTGGCATGAGCAAAATGACTGCGGCTGCCCGGTCCGGCCTGACGTGCCAGCAGGCCTGCTTCTACCGcccg                                                                >  2:1989556/4‑148 (MQ=255)
                                                                                     tGTCGTCCCAGTCCACTTGTCGCCAGCGATCGCTGAAGGTTTTCATGCCGATGCCTCCAGAAATCCGGTCAGCGGGCTGGTGGCATGAGCAAAATAACTGCGGCTGCCCGGTCCGGACTGACGTGCCAGTAGGCCTGCTTCTACCGCCAg                                                                >  2:2274215/1‑150 (MQ=255)
                                                                                       tcgtcCCAGTCCAGTTGTCGCCAGCGATCGCTGAAGGTTTTCATGCCGATGCCTCCAGAAATCCGGTCAGCGGGCTGGTGGCATGAGCAAAATAACTGCGGCTGCCCGGTCCGGACTGACGTGCCAGTAGGCCTGCTTCTACCGCCAGAc                                                              <  2:539457/150‑1 (MQ=255)
                                                                                        cgtcCCAGTCCAGTTGTCGCCAGCGATCGCTGAAGGTTTTCATGCCGATGCCTCCAGAAATCCGGTCAGCGGGCTGGTGGCATGAGCAAAATAACTGCGGCTGCCCGGTCCGGACTGACGTGCCAGTAGGCCTGCTTCTACCGCCAGACg                                                             >  1:900503/1‑150 (MQ=255)
                                                                                        cgtcCCAGTCCAGTTGTCGCCAGCGATCGCTGAAGGTTTTCATGCCGATGCCTCCAGAAATCCGGTCAGCGGGCTGGTGGCATGAGCAAAATAACTGCGGCTGCCCGGTCCGGACTGACGTGCCAGTAGGCCTGCTTCTACCGCCAGACg                                                             >  1:889409/1‑150 (MQ=255)
                                                                                          gcccAGTCCAGTTGTCGCCAGCGATCGCTGAAGGTTTTCATGCCGATGCCTCCAGAAATCCGGTCAGCGGGCTGGTGGCATGAGCAAAATAACTGCGGCTGCCCGGTCCGGACTGACGTGCCAGTAGGCCTGCTTCTACCGCCAGACGaa                                                           <  1:512286/149‑1 (MQ=255)
                                                                                              agtccagtTGTCGCCAGCGATCGCTGAAGGTTTTCATGCCGATGCCTCCAGAAATCCGGTCAGCGGGCTGGTGGCATGAGCAAAATAACTGCGGCTGCCCGGTCCGGACTGACGTGCCAGTAGGCCTGCTTCTACCGCCAGACGAAATGc                                                       <  2:421669/150‑1 (MQ=255)
                                                                                                 ccagtTGTCGCCAGCGATCGCTGAAGGTTTTCATGCCGATGCCTCCAGAAATCCGGTCAGCGGGCTGGTGGCATGAGCAAAATAACTGCGGCTGCCCGGTCCGGACTGACGTGCCAGTAGGcc                                                                               <  2:139255/123‑1 (MQ=255)
                                                                                                 ccagtTGTCGCCAGCGATCGCTGAAGGTTTTCATGCCGATGCCTCCAGAAATCCGGTCAGCGGGCTGGTGGCATGAGCAAAATAACTGCGGCTGCCCGGTCCGGACTGACGTGCCAGTAGGcc                                                                               >  1:139255/1‑123 (MQ=255)
                                                                                                    gtTGTCGCCAGCGATCGCTGAAGGTTTTCATGCCGATGCCTCCAGAAATCCGGTCAGCGGGCTGGTGGCATGAGCAAAATAACTGCGGCTGCCCGGTCCGGACTGACGTGCCAGTAGGCCTGCTTCTACCGCCAGACGAAATGCCTTCGc                                                 >  2:1108088/1‑150 (MQ=255)
                                                                                                    gtTGTCGCCAGCGATCGCTGAAGGTTTTCATGCCGATGCCTCCAGAAATCCGGTCAGCGGGCTGGTGGCATGAGCAAAATAACTGCGGCTGCCCGGTCCGGACTGACGTGCCAGTAGGCCTGCTTCTACCGCCAGACGAAATGCCTTCGc                                                 <  2:1902206/150‑1 (MQ=255)
                                                                                                       gTCGCCAGCGATCGCTGAAGGTTTTCATGCCGATGCCTCCAGAAATCCGGTCAGCGGGCTGGTGGCATGAGCAAAATAACTGCGGCTGCCCGGTCCGGACTGACGTGCCAGTAGGCCTGCTTCTCCCGCCAGACGAAATGCCTTCGCCag                                              >  2:2481326/1‑149 (MQ=255)
                                                                                                        tCGCCAGCGATCGCTGAAGGTTTTCATGCCGATGCCTCCAGAAATCCGGTCAGCGGGCTGGTGGCATGAGCAAAATAACTGCGGCTGCCCGGTCCGGACTGACGTGCCAGTAGGCCTGCTTCTACCGCCAGACGAAATGCCTTCGCCATg                                             <  1:574630/150‑1 (MQ=255)
                                                                                                        gcGCCAGCGATCGCTGAAGGTTTTCATGCCGATGCCTCCAGAAATCCGGTCAGCGGGCTGGTGGCATGAGCAAAATAACTGCGGCTGCCCGGTCCGGACTGACGTGCCAGTAGGCCTGCTTCTACCGCCAGACGAAATGCCTTCGCCATg                                             <  2:906330/149‑1 (MQ=255)
                                                                                                         cGCCAGCGATCGCTGAAGGTTTTCATGCCGATGCCTCCAGAAATCCGGTCAGCGGGCTGGTGGCATGAGCAAAATAACTGCGGCTGCCCGGTCCGGACTGACGTGCCAGTAGGCCTGCTTCTACCGCCAGACGAAATGCCTTCGCCATGt                                            >  1:2224858/1‑150 (MQ=255)
                                                                                                         cGCCAGCGATCGCTGAAGGTTTTCATGCCGATGCATCCAGAAATCCGGTCAGCGGGCTGGTGGCATGAGCAAAATAACTGCGGCTGCCCGGTCCGGACTGACGTGCCAGTAGGCCTGCTTCTACCGCCAGACGAAATGCCTTCGCCATGt                                            >  2:1144781/1‑150 (MQ=255)
                                                                                                          gCCAGCGATCGCTGAAGGTTTTCATGCCGATGCCTCCAGAAATCCGGTCAGCGGGCTGGTGGCATGAGCAAAATAACTGCGGCTGCCCGGTCCGGACTGACGTGCCAGTAGGCCTGCTTCTACCGCCAGACGAAATGCCTTCGCCATGtt                                           <  1:807728/150‑1 (MQ=255)
                                                                                                           ccAGCGATCGCTGAAGGTTTTCATGCCGATGCCTCCAGAAATCCGGTCAGCGGGCTGGTGGCATGAGCAAAATAACTGCGGCTGCCCGGTCCGGACTGACGTGCCAGTAGGCCTGCTTCTACCGCCAGACGAAATGCCTTCGCCATGTTg                                          >  2:1097823/1‑150 (MQ=255)
                                                                                                            cAGCGATCGCTGAAGGTTTTCATGCCGATGCCTCCAGAAATCCGGTCAGCGGGCTGGTGGCATGAGCAAAATAACTGCGGCTGCCCGGTCCGGACTGACGTGCCAGTAGGCCTGCTTCTACCGCCAg                                                                >  2:930765/1‑127 (MQ=255)
                                                                                                            cAGCGATCGCTGAAGGTTTTCATGCCGATGCCTCCAGAAATCCGGTCAGCGGGCTGGTGGCATGAGCAAAATAACTGCGGCTGCCCGGTCCGGACTGACGTGCCAGTAGGCCTGCTTCTACCGCCAg                                                                <  1:930765/127‑1 (MQ=255)
                                                                                                              gCGATCGCTGAAGGTTTTCATGCCGATGCCTCCAGAAATCCGGTCAGCGGGCTGGTGGCATGAGCAAAATAACTGCGGCTGCCCGGTCCGGACTGACGTGCCAGTAGGCCTGCTTCTACCGCCAGACGAAATGCCTTCGCCATGTTGACg                                       <  2:1521110/150‑1 (MQ=255)
                                                                                                              gCGATCGCTGAAGGTTTTCATGCCGATGCCTCCAGAAATCCGGTCAGCGGGCTGGTGGCATGAGCAAAATAACTGCGGCTGCCCGGTCCGGACTGACGTGCCAGTAGGCCTGCTTCTACCGCCAGACGAAATGCCTTCGCCATGTTGACg                                       <  1:1469569/150‑1 (MQ=255)
                                                                                                                  tCGCTGAAGGTTTTCATGCCGATGCCTCCAGAAATCCGGTCAGCGGGCTGGTGGCATGAGCAAAATAACTGCGGCTGCCCGGTCCGGACTGACGTGCCAGTAGGCCTGCTTCTACCGCCAGACGAAATGCCTTCGCCATGTTGACGGGAt                                   <  1:1920188/150‑1 (MQ=255)
                                                                                                                  tCGCTGAAGGTTTTCATGCCGATGCCTCCAGAAATCCGGTCAGCGGGCTGGTGGCATGAGCAAAATAACTGCGGCTGCCCGGTCCGGACTGACGTGCCAGTAGGCCTGCTTCTACCGCCAGACGAAATGCCTTCGCCATGTTGACGGGAt                                   <  1:1974471/150‑1 (MQ=255)
                                                                                                                   cGCTGAAGGTTTTCATGCCGATGCCTCCAGAAATCCGGTCAGCGGGCTGGTGGCATGAGCAAAATAACTGCGGCTGCCCGGTCCGGACTGACGTGCCAGTAGGCCTGCTTCTACCGCCAGACGAAATGCCTTCGCCATGTTGACGGGATc                                  >  1:901544/1‑150 (MQ=255)
                                                                                                                    gCTGAAGGTTTTCATGCCGATGCCTCCAGAAATCCGGTCAGCGGGCTGGTGGCATGAGCAAAATAACTGCGGCTGCCCGGTCCGGACTGACGTGCCAGTAGGCCTGCTTCTACCGCCAGACGAAATGCCTTCGCCATGTTGACGGGATCg                                 >  2:932853/1‑150 (MQ=255)
                                                                                                                      tGAAGGTTTTCATGCCGATGCCTCCAGAAATCCGGTCAGCGGGCTGGTGGCATGAGCAAAATAACTGCGGCTGCCCGGTCCGGACTGACGTGCCAGTAGGCCTGCTTCTACCGCCAGACGAAATGCCTTCGCCATGTTGACGGGATCGTc                               >  1:1736797/1‑150 (MQ=255)
                                                                                                                       gAAGGTTTTCATGCCGATGCCTCCAGAAATCCGGTCAGCGGGCTGGTGGCATGAGCAAAATAACTGCGGCTGCCCGGTCCGGACTGACGTGCCAGTAGGCCTGCTTCTACCGCCAGACGAAATGCCTTCGCCATGTTGACGGGATCGTcc                              >  1:2250820/1‑150 (MQ=255)
                                                                                                                       gAAGGTTTTCATGCCGATGCCTCCAGAAATCCGGTCAGCGGGCTGGTGGCATGAGCAAAATAACTGCGGCTGCCCGGTCCGGACTGACGTGCCAGTAGGCCTGCTTCTACCGCCAGACGAAATGCCTTCGCCATGTTGACGGGATCGTcc                              >  1:524446/1‑150 (MQ=255)
                                                                                                                          ggTTTTCATGCCGATGCCTCCAGAAATCCGGTCAGCGGGCTGGTGGCATGAGCAAAATAACTGCGGCTGCCCGGTCCGGACTGACGTGCCAGTAGGCCTGCTTCTACCGCCAGACGAAATGCCTTCGCCATGTTGACGGGATCGTCcgcg                           >  1:1043618/1‑150 (MQ=255)
                                                                                                                          ggTTTTCATGCCGATGCCTCCAGAAATCCGGTCAGCGGGCTGGTGGCATGAGCAAAATAACTGCGGCTGCCCGGTCCGGACTGACGTGCCAGTAGGCCTGCTTCTACCGCCAGACGAAATGCCTTCGCCATGTTGACGGGATCGTCcgcg                           >  1:766192/1‑150 (MQ=255)
                                                                                                                          ggTTTTCATGCCGATGCCTCCAGAAATCCGGTCAGCGGGCTGGTGGCATGAGCAAAATAACTGCGGCTGCCCGGTCCGGACTGACGTGCCAGTAGGCCTGCTTCTACCGCCAGACGAAATGCCTTCGCCATGTTGACGGGATCGTCcgcg                           >  1:1772495/1‑150 (MQ=255)
                                                                                                                            ttttCATGCCGATGCCTCCAGAAATCCGGTCAGCGGGCTGGTGGCATGAGCAAAATAACTGCGGCTGCCCGGTCCGGACTGACGTGCCAGTAGGCCTGCTTCTACCGCCAGACGAAATGCCTTCGCCATGTTGACGGGATCGTCCGCGAc                         >  1:128099/1‑150 (MQ=255)
                                                                                                                            ttttCATGCCGATGCCTCCAGAAATCCGGTCAGCGGGCTGGTGGCATGAGCAAAATAACTGCGGCTGCCCGGTCCGGACTGACGTGCCAGTAGGCCTGCTTCTACCGCCAGACGAAATGCCTTCGCCATGTTGACGGGATCGTCCGCGAc                         <  2:2012167/150‑1 (MQ=255)
                                                                                                                            ttttCATGCCGATGCCTCCAGAAATCCGGTCAGCGGGCTGGTGGCATGAGCAAAATAACTGCGGCTGCCCGGTCCGGACTGACGTGCCAGTAGGCCTGCTTCTACCGCCAGACGAAATGCCTTCGCCATGTTGACGGGATCGTCCGCGAc                         >  2:1230074/1‑150 (MQ=255)
                                                                                                                             tttCATGCCGATGCCTCCAGAAATCCGGTCAGCGGGCTGGTGGCATGAGCAAAATAACTGCGGCTGCCCGGTCCGGACTGACGTGCCAGTAGGCCTGCTTCTACCGCCAGACGAAATGCCTTCGCCATGTTGACGGGATCGTCCGCGACg                        >  2:2394768/1‑150 (MQ=255)
                                                                                                                                    ccGATGCCTCCAGAAATCCGGTCAGCGGGCTGGTGGCATGAGCAAAATAACTGCGGCTGCCCGGTCCGGACTGACGTGCCAGTAGGCCTGCTTCTACCGCCAGACGAAATGCCTTCGCCATGTTGACGGGATCGTCCGCGACGGCAATCg                 >  1:1700876/1‑150 (MQ=255)
                                                                                                                                        tGCCTCCAGAAATCCGGTCAGCGGGCTGGTGGCATGAGCAAAATAACTGCGGCTGCCCGGTCCGGACTGACGTGCCAGTAGGCCTGCTTCTACCGCCAGACGAAATGCCTTCGCCATGTTGACGGGATCGTCCGCGACGGCAATCGCCGt             <  2:2352904/150‑1 (MQ=255)
                                                                                                                                           cTCCAGAAATCCGGTCAGCGGGCTGGTGGCATGAGCAAAATAACTGCGGCTGCCCGGTCCGGACTGACGTGCCAGTAGGCCTGCTTCTACCGCCAGACGAAATGCCTTCGCCATGTTGACGGGATCGTCCGCGACGGCAATCGCCGTAtt          >  1:1625649/1‑150 (MQ=255)
                                                                                                                                            tCCAGAAATTCGGTCAGCGGGCTGGTGGCATGAGCAAAATAACTGCGGCTGCCCGGTCCGGACTGACGTGCCAGTAGCCCTGCTTCTACCGCCAGACGAAATGCCTTCGCCATGTTGACGGGATCGTCCGCGACGGCAATCGCCGTATTc         <  1:1353156/150‑1 (MQ=255)
                                                                                                                                            tCCAGAAATCCGTTCAGCGGGCTGGTGGCATGAGCAAAATAACTGCGGCTGCCCGGTCCGGACTGACGTGCCAGTAGGCCTGCTTCTACCGCCAGACGAAATGCCTTCGCCATGTTCCCGGGATCGTCCGCGACGGCAATCGCCGTAtta         >  1:2343438/1‑149 (MQ=255)
                                                                                                                                            tCCAGAAATCCGGTCAGCGGGCTGGTGGCATGAGCAAAATAACTGCGGCTGCCCGGTCCGGACTGACGTGCCAGTAGGCCTGCTTCTACCGCCAGACGAAATGCCTTCGCCATGTTGACGGGATCGTCCGCGACGGCAATCGCCGTATTc         >  2:883544/1‑150 (MQ=255)
                                                                                                                                            tCCAGAAATCCGGTCAGCGGGCTGGTGGCATGAGCAAAATAACTGCGGCTGCCCGGTCCGGACTGACGTGCCAGTAGGCCTGCTTCTACCGCCAGACGAAATGCCTTCGCCATGTTGACGGGATCGTCCGCGACGGCAATCGCCGTATTc         >  1:1073540/1‑150 (MQ=255)
                                                                                                                                               aGAAATCCGGTCAGCGGGCTGGTGGCATGAGCAAAATAACTGCGGCTGCCCGGTCCGGACTGACGTGCCAGTAGGCCTGCTTCTACCGCCAGACGAAATGCCTTCGCCATGTTGACGGGATc                                  <  2:1384419/122‑1 (MQ=255)
                                                                                                                                                gAAATCCGGTCAGCGGGCTGGTGGCATGAGCAAAATAACTGCGGCTGCCCGGTCCGGCCTGACGTGCCAGTAGGCCTGCTTCTCCCGCCAGACGAAATGCCTTCGCCATGTTGACGGGATCGTCCGCGACGGCAATCGCCGTATTCACTa     >  1:427155/1‑150 (MQ=255)
                                                                                                                                                  aaTCCGGTCAGCGGGCTGGTGGCATGAGCAAAATAACTGCGGCTGCCCGGTCCGGACTGACGTGCCAGTAGGCCTGCTTCTACCGCCAGACGAAATGCCTTCGCCATGTTGACGGGATCGTCCGCGACGGCAATCGCCGTATTCACTAac   >  2:1421506/1‑150 (MQ=255)
                                                                                                                                                   aTCCGGTCAGCGGGCTGGTGGCATGAGCAAAATAACTGCGGCTGCCCGGTCCGGACTGACGTGCCAGTAGGCCTGCTTCTACCGCCAGACGAAATGCCTTCGCCATGTTGACGGGATCGTCCGCGACGGCAATCGCCGTATTCACTAaca  <  2:1671088/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
GCCATCATGTCGTCGCGGGTGAGTTGCGAGGCATTTAGCGCCCGCTCTACGTCAGCAGCCGTTTTGCCGTTGATACGCAGGCGGATGTCGTCCCAGTCCAGTTGTCGCCAGCGATCGCTGAAGGTTTTCATGCCGATGCCTCCAGAAATCCGGTCAGCGGGCTGGTGGCATGAGCAAAATAACTGCGGCTGCCCGGTCCGGACTGACGTGCCAGTAGGCCTGCTTCTACCGCCAGACGAAATGCCTTCGCCATGTTGACGGGATCGTCCGCGACGGCAATCGCCGTATTCACTAACA  >  NC_000913/4191738‑4192034

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