Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,532,268 T→C 15.1% intergenic (+44/‑141) 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,2680TC15.1% 242.0 / 32.7 105intergenic (+44/‑141)cysZ/cysKsulfate transporter, sulfite inhibited/cysteine synthase A, O‑acetylserine sulfhydrolase A subunit
Reads supporting (aligned to +/- strand):  ref base T (40/49);  new base C (4/12);  total (44/61)
Fisher's exact test for biased strand distribution p-value = 1.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.95e-01

Only 100 of 106 total aligned reads displayed.

CGCGACGGCGATGTGGGTCGATTGCTATCGCGATAAACACGCGATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATCCTTACTTCCGCATATTCTCTGAGCGGGTAT  >  NC_000913/2532170‑2532366
                                                                                                  |                                                                                                  
cgcgACGGCGATGTGGGTCGATTGCTATCGCGATAAACACGCGATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTg                                                                                                   <  1:521993/100‑1 (MQ=255)
cgcgACGGCGATGTGGGTCGATTGCTATCGCGATAAACACGCGATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTg                                                                                                   >  2:405273/1‑100 (MQ=255)
   gACGGCGATGTGGGTCGATTGCTATCGCGATAAACACGCGATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCt                                                                                                >  1:291507/1‑100 (MQ=255)
   gACGGCGATGTGGGTCGATTGCTATCGCGATAAACACGCGATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCt                                                                                                <  1:1217375/100‑1 (MQ=255)
   gACGGCGATGTGGGTCGATTGCTATCGCGATAAACACGCGATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCt                                                                                                <  1:1848252/100‑1 (MQ=255)
     cGGCGATGTGGGTCGATTGCTATCGCGATAAACACGCGATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTa                                                                                              >  2:2123577/1‑100 (MQ=255)
        cGATGTGGGTCGATTGCTATCGCGATAAACACGCGATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATgn                                                                                           <  1:2253596/100‑2 (MQ=255)
        cGATGTGGGTCGATTGCTATCGCGATAAACACGCGATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATgc                                                                                           >  2:111694/1‑100 (MQ=255)
         gATGTGGGTCGATTGCTATCGCGATAAACACGCGATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATgcc                                                                                          >  1:1368138/1‑100 (MQ=255)
            gtgGGTCGATTGCTATCGCGATAAACACGCGATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATgccgcc                                                                                       <  1:2061214/100‑1 (MQ=255)
             tgGGTCGATTGCTATCGCGATAAACACGCGATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGccc                                                                                      >  2:754539/1‑100 (MQ=255)
             tgGGTCGATTGCTATCGCGATAAACACGCGATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGccc                                                                                      <  2:1350054/100‑1 (MQ=255)
              gggTCGATTGCTATCGCGATGAACACGCGATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGcccc                                                                                     <  1:140946/100‑1 (MQ=255)
              gggTCGATTGCTATCGCGATAAACACGCGATGTGGCGGTAACAATCTAGCGGTTATTTTGTAAACCGGTTGTGTGAAACAGGGGTGGCTTATGCCGcccc                                                                                     <  2:2446406/100‑1 (MQ=255)
              gggTCGATTGCTATCGCGATAAACACGCGATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGcccc                                                                                     <  1:1490233/100‑1 (MQ=255)
              gggTCGATTGCTATCGCGATAAACACGCGATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGcccc                                                                                     >  2:1331987/1‑100 (MQ=255)
              gggTCGATTGCTATCGCGATAAACACGCGATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGcccc                                                                                     >  1:1252573/1‑100 (MQ=255)
              gggTCGATTGCTATCGCGATAAACACGCGATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGCGGCATAAGCCAcccc                                                                                     <  1:239281/100‑1 (MQ=255)
              gggTCGATTGCTATCGCGATAAACACGCGATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGCGGCATAAGCCAcccc                                                                                     <  1:912811/100‑1 (MQ=255)
              gggTCGATTGCTATCGCGATAAACACGCGATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGCGGCATAAGCCAcccc                                                                                     <  2:303512/100‑1 (MQ=255)
              gggTCGATTGCTATCGCGATAAACACGCGATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGCGGCATAAGCCAcccc                                                                                     <  1:2234866/100‑1 (MQ=255)
               ggTCGATTGCTATCGCGATAAACACGCGATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCt                                                                                    >  2:2399823/1‑100 (MQ=255)
               ggTCGATTGCTATCGCGATAAACACGCGATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGCGGCATAAGCCACCCCt                                                                                    <  1:1482449/100‑1 (MQ=255)
               ggTCGATGGCTATCGCGATAAACACGCGATGTGGCGGTAAGAATCTACCGGTTATTTTGGAGACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCctct                                                                                    <  2:320595/100‑1 (MQ=255)
                 tCGATTGCTATCGCGATAAACACGCGATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGCGGCATAAGCCACCCCtgt                                                                                  <  2:961301/100‑3 (MQ=255)
                 tCGATTGCTATCGCGATAAACACGCGATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGCGGCATAAGCCACCCCtgt                                                                                  <  2:1500876/100‑3 (MQ=255)
                  cGATTGCTATCGCGATAAACACGCGATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGCGGCATAAGCCACCCCtgtt                                                                                 <  1:739040/100‑4 (MQ=255)
                   gATTGCTATCGCGATAAACACGCGATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTAtt                                                                                <  1:2161676/100‑1 (MQ=255)
                   gATTGCTATCGCGATAAACACGCGATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTAtt                                                                                >  1:1955536/1‑100 (MQ=255)
                   gATTGCTATCGCGATAAACACGCGATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTAtt                                                                                >  2:610128/1‑100 (MQ=255)
                   gATTGCTATCGCGATAAACACGCGATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTAt                                                                                 <  2:1672287/99‑1 (MQ=255)
                    aTTGCTATCGCGATAAACACGCGATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTc                                                                               <  2:1938065/100‑1 (MQ=255)
                    aTTGCTATCGCGATAAACACGCGATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTc                                                                               <  1:2123591/100‑1 (MQ=255)
                     ttGCTATCGCGATAAACACGCGATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGCGGCATAAGCCACCCCtgtttc                                                                               <  2:1941571/99‑6 (MQ=37)
                      tGCTATCGCGATAAACACGCGATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGCGGCATAAGCCACCCCtgtttc                                                                               <  2:1957089/98‑6 (MQ=37)
                       gCTATCGCGATAAACACGCGATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCAt                                                                            >  2:2303119/1‑100 (MQ=255)
                       gCTATCGCGATAAACACGCGATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCAt                                                                            <  1:1159552/100‑1 (MQ=255)
                        cTATCGCGATAAACACGCGATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCTATc                                                                           >  2:259373/1‑100 (MQ=255)
                        cTATCGCGATAAACACGCGATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATc                                                                           >  1:432123/1‑100 (MQ=255)
                        cTATCGCGATAAACACGCGATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGCGGCATAAGCCACCCCtgtttcacac                                                                           <  1:1328967/100‑10 (MQ=37)
                         tATCGCGATAAACACGCGATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCt                                                                          >  1:432190/1‑100 (MQ=255)
                          aTCGCGATAAACACGCGATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCt                                                                          <  1:1118894/99‑1 (MQ=255)
                            cgcgATAAACACGCGATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGc                                                                       <  1:2218685/100‑1 (MQ=255)
                              cgATAAACACGCGATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCAt                                                                     >  1:2299811/1‑100 (MQ=255)
                                   aaCACGCGATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCat                                                                <  2:1846900/100‑1 (MQ=255)
                                   aaCACGCGATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCat                                                                <  2:1088102/100‑1 (MQ=255)
                                     cacGCGATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCatta                                                              >  1:2227943/1‑100 (MQ=255)
                                        gcgATGTGGCGGTAACAATCTACCGGTTATTTTGTAACCCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTAttt                                                           >  2:1369112/1‑100 (MQ=255)
                                         cgATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTAtttn                                                          <  2:1124243/100‑2 (MQ=255)
                                         cgATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTc                                                          <  2:2162775/100‑1 (MQ=255)
                                         cgATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTc                                                          >  1:2325137/1‑100 (MQ=255)
                                           aTGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATTCCGCCCCTTATTCCATCTTGCATGTCATTATTTccc                                                        <  1:1809705/100‑1 (MQ=255)
                                           aTGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTccc                                                        <  1:2454403/100‑1 (MQ=255)
                                            tgtgGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCt                                                       >  1:1242770/1‑100 (MQ=255)
                                               ggCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCt                                                    >  2:1913973/1‑100 (MQ=255)
                                               ggCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCt                                                    >  2:135375/1‑100 (MQ=255)
                                                gcttgaaCAATCTACCGGTTATTTTGTACACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTg                                                   <  2:1620239/95‑1 (MQ=255)
                                                  ggTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGta                                                 <  2:2076880/100‑1 (MQ=255)
                                                     aaCAATCTACCGGTTATTTTGTAAACCGCTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGtatat                                              <  2:631839/100‑1 (MQ=255)
                                                      aCAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGtatata                                             <  2:2441021/100‑1 (MQ=255)
                                                      aCAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGtatata                                             <  1:2352804/100‑1 (MQ=255)
                                                       cAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAg                                            >  1:223301/1‑100 (MQ=255)
                                                       cAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAg                                            >  1:629548/1‑100 (MQ=255)
                                                        aaTCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGa                                           <  2:1158328/100‑1 (MQ=255)
                                                        aaTCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGa                                           >  2:1608202/1‑100 (MQ=255)
                                                        aaTCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGa                                           <  1:135377/100‑1 (MQ=255)
                                                              ccGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCt                                     >  1:213552/1‑100 (MQ=255)
                                                               cGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTa                                    >  1:765356/1‑100 (MQ=255)
                                                                ggTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTaa                                   >  2:1513431/1‑100 (MQ=255)
                                                                 gTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTaaa                                  <  1:2470521/100‑1 (MQ=255)
                                                                  ttATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAAt                                 <  2:1984743/100‑1 (MQ=255)
                                                                    aTTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATcc                               >  1:1382119/1‑100 (MQ=255)
                                                                       ttGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATCCTTa                            >  1:1328038/1‑100 (MQ=255)
                                                                       ttGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATCCTTa                            >  2:2081860/1‑100 (MQ=255)
                                                                       ttGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATCCTTa                            <  2:961900/100‑1 (MQ=255)
                                                                        tGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATCCTTAc                           <  2:616738/100‑1 (MQ=255)
                                                                         gTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATCCTTACt                          >  1:1907747/1‑100 (MQ=255)
                                                                          tAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATCCTTACtt                         >  2:2061580/1‑100 (MQ=255)
                                                                           aaaCCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATCCTTACTTc                        <  1:2514098/100‑1 (MQ=255)
                                                                              ccGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATCCTTACTTCCGc                     <  2:2028184/100‑1 (MQ=255)
                                                                              ccGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATCCTTACTTCCGc                     <  1:1608213/100‑1 (MQ=255)
                                                                               cGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATCCTTACTTCCGCa                    <  2:914983/100‑1 (MQ=255)
                                                                                      gtgAAACAGGGGTGGCTTATGCCGCCCCTTATTGCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATCCTTACTTCCGCATATtctc             >  2:1694769/1‑100 (MQ=255)
                                                                                       tgAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATCCTTACTTCCGCATATTctct            <  1:1821489/100‑1 (MQ=255)
                                                                                        gaataaGGGGCGGCATAAGCCACCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATCCTTACTTCCGCATATTCTCTg           >  2:1639485/6‑100 (MQ=38)
                                                                                        gaataaGGGGCGGCATAAGCCACCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATCCTTACTTCCGCATATTCTCTg           >  1:1334589/6‑100 (MQ=38)
                                                                                        gAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATCCTTACTTCCGCATATTCTCTg           <  1:610135/100‑1 (MQ=255)
                                                                                         aaaCAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATCCTTACTTCCGCATATTCTCTGa          <  1:675623/100‑1 (MQ=255)
                                                                                         aaaCAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATCCTTACTTCCGCATATTCTCTGa          >  1:349052/1‑100 (MQ=255)
                                                                                           aCAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATCCTTACTTCCGCATATTCTCTGAGc        <  1:311898/100‑1 (MQ=255)
                                                                                            cAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATCCTTACTTCCGCATATTCTCTGAGCg       <  1:622304/100‑1 (MQ=255)
                                                                                              ggggTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATCCTTACTTCCGCATATTCTCTGAGCggg     <  1:2249676/100‑1 (MQ=255)
                                                                                              ggggTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATCCTTACTTCCGCATATTCTCTGAGCggg     >  1:270067/1‑100 (MQ=255)
                                                                                              ggggTGGCTTATGCCGCCCCTGATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATCCTTACTTCCGCATATTCTCTGAGCggg     <  2:2162345/100‑1 (MQ=255)
                                                                                              ggggCGGCATAGGCCACCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATCCTTACTTCCGCATATTCTCTGAGCggg     >  2:2041310/1‑100 (MQ=255)
                                                                                               gggTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATCCTTACTTCCGCATATTCTCTGAGCGGGt    <  2:2238077/100‑1 (MQ=255)
                                                                                                ggTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTGCTGTATATAGATATGCTAAATCCTTACTTCCGCATATTCTCTGAGCGGGTa   <  2:2518059/100‑1 (MQ=255)
                                                                                                 gTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATCCTTACTTCCGCATATTCTCTGAGCGGGTAt  <  2:2227929/100‑1 (MQ=255)
                                                                                                 gTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCCTCTGTATATAGATATGCCAAATCCTTACTTCCGCACATTCTCTGAGCGCGTa   <  2:2040726/99‑1 (MQ=255)
                                                                                                  tGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATCCTTACTTCCGCATATTCTCTGAGCGGGTa   <  2:1728923/98‑1 (MQ=255)
                                                                                                  |                                                                                                  
CGCGACGGCGATGTGGGTCGATTGCTATCGCGATAAACACGCGATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATCCTTACTTCCGCATATTCTCTGAGCGGGTAT  >  NC_000913/2532170‑2532366

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