Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,532,268 T→C 15.8% 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.8% 385.0 / 57.4 163intergenic (+44/‑141)cysZ/cysKsulfate transporter, sulfite inhibited/cysteine synthase A, O‑acetylserine sulfhydrolase A subunit
Reads supporting (aligned to +/- strand):  ref base T (56/81);  new base C (12/14);  total (68/95)
Fisher's exact test for biased strand distribution p-value = 6.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 165 total aligned reads displayed.

GCGCGACGGCGATGTGGGTCGATTGCTATCGCGATAAACACGCGATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATCCTTACTTCCGCATATTCTCTGAGCGGGTATG  >  NC_000913/2532169‑2532367
                                                                                                   |                                                                                                   
gcgcGACGGCGATGTGGGTCGATTGCTATCGCGATAAACACGCGATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGt                                                                                                     >  1:4031863/1‑100 (MQ=255)
   ctcttccGATCTGGGTCGATTGCTATCGCGATAAACACGCGATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGc                                                                                                  <  2:2096335/94‑1 (MQ=255)
          gATGTGGGTCGATTGCTATCGCGATAAACACGCGATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATgcc                                                                                           <  2:3158204/100‑1 (MQ=255)
          gATGTGGGTCGATTGCTATCGCGATAAACACGCGATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATgcc                                                                                           >  2:486768/1‑100 (MQ=255)
                ggTGGAGTGGTAGCGCGATAAACACGCGGTGTGGCGGGAACGATCTGGCGGTGATTGTGGAAACCGTTTGTGTGAAACAGGGGCGGCATAAGCCACCCCt                                                                                     <  2:2477970/100‑1 (MQ=255)
                ggTCGATTGCTATCGCGATAAACACGCGATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCt                                                                                     >  2:1542144/1‑100 (MQ=255)
                 gTCGATTGCTATCGCGATAAACACGCGATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGCGGCATAAGCCACCCCtg                                                                                    <  2:4075189/100‑2 (MQ=255)
                 gTCGATTGCTATCGCGATAAACACGCGATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGCGGCATAAGCCACCCCtg                                                                                    <  1:2412711/100‑2 (MQ=255)
                  tCGATTGCTATCGCGATAAACACGCGATGTGGTGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGCGGCATAAGCCACCCCt                                                                                     <  1:2276782/98‑1 (MQ=255)
                  tCGATTGCTATCGCGATAAACACGCGATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGCGGCATAAGCCACCCCtgt                                                                                   <  1:443363/100‑3 (MQ=255)
                   cGATTGCTATCGCGATAAAGACGCGATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGCGGCATAAGCCACCCCtgtt                                                                                  <  1:3536634/100‑4 (MQ=255)
                   cGATTGCTATCGCGATAAACACGCGATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTAt                                                                                  >  2:3765135/1‑100 (MQ=255)
                    gATTGCTATCGCGATAAACACGCGATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTAtt                                                                                 <  2:576022/100‑1 (MQ=255)
                     aTTGCTATCGCGATAAACACGCGATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTc                                                                                <  1:2313026/100‑1 (MQ=255)
                      ttGCTATCGCGATAAGCACGCGATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTcc                                                                               <  1:4129233/100‑1 (MQ=255)
                       tGCTATCGCGATAAACACGCGATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGCGGCATAAGCCACCCCtgtttcac                                                                              <  1:2163545/100‑8 (MQ=37)
                          tATCGCGATAAACACGCGATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCt                                                                           <  2:2793942/100‑1 (MQ=255)
                          tATCGCGATAAACACGCGATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCt                                                                           <  2:1217356/100‑1 (MQ=255)
                            tCGCGATAAACACGCGATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTg                                                                         <  1:3425515/100‑1 (MQ=255)
                             cgcgATAAACACGCGATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGc                                                                        <  2:464809/100‑1 (MQ=255)
                              gcgATAAACACGCGATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCa                                                                       <  1:4175332/100‑1 (MQ=255)
                              gcgATAAACACGCGATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCa                                                                       <  2:2000832/100‑1 (MQ=255)
                                gATAACCACGCGTTGTGGCTGTAATAGTTTAGCGGTTATTTTGTAACCCGTTTGTGTGACACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCa                                                                       <  2:2580511/98‑1 (MQ=255)
                                gATAAACACGCGATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATg                                                                     <  1:1459801/100‑1 (MQ=255)
                                gATAAACACGCGATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATg                                                                     >  2:1021402/1‑100 (MQ=255)
                                 aTAAACACGCGATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGt                                                                    <  2:261249/100‑1 (MQ=255)
                                    aaCACGCGATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCat                                                                 <  2:2816897/100‑1 (MQ=255)
                                     acacGCGATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCatt                                                                <  1:679133/100‑1 (MQ=255)
                                      cacGCGATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCatt                                                                <  1:530012/99‑1 (MQ=255)
                                       acGCGATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCattat                                                              <  2:3169077/100‑1 (MQ=255)
                                         gcgATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTAttt                                                            >  2:1146693/1‑100 (MQ=255)
                                         gcgATGGGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCattat                                                              <  2:63500/98‑1 (MQ=255)
                                          cgATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTc                                                           <  2:3665992/100‑1 (MQ=255)
                                          cgATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTc                                                           <  1:3636378/100‑1 (MQ=255)
                                           gATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTcc                                                          <  1:769540/100‑1 (MQ=255)
                                              gtgGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCtt                                                       <  1:723813/100‑1 (MQ=255)
                                               tgGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTc                                                      <  1:3165607/100‑1 (MQ=255)
                                                ggCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCt                                                     >  2:1820330/1‑100 (MQ=255)
                                                 gCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTg                                                    <  1:4010298/100‑1 (MQ=255)
                                                 gCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTg                                                    <  1:1217044/100‑1 (MQ=255)
                                                  cGGTAACAATCTACCGGTTATTTTGTAAGCCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGt                                                   <  2:3144899/100‑1 (MQ=255)
                                                   ggTAACAATCTACCGGTTATTTTGTACACCGTTTGTGTGAAACAGGGGTGGTTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGta                                                  <  2:2359616/100‑1 (MQ=255)
                                                         aaTCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGa                                            <  1:2778652/100‑1 (MQ=255)
                                                         aaTCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGa                                            >  2:1112642/1‑100 (MQ=255)
                                                          aTCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGat                                           >  2:1260791/1‑100 (MQ=255)
                                                           tCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGata                                          <  1:1760132/100‑1 (MQ=255)
                                                           tCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGata                                          <  1:3756054/100‑1 (MQ=255)
                                                            cTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGatat                                         >  1:3053364/1‑100 (MQ=255)
                                                             tACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATg                                        >  1:4060014/1‑100 (MQ=255)
                                                             tACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATg                                        >  2:4126359/1‑100 (MQ=255)
                                                             tACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATg                                        >  2:4126294/1‑100 (MQ=255)
                                                             tACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGATTATGCCGCCCCTTATTCCATCTTGCATGTCATAATTTCCCTTCTGTATATAGATATg                                        <  2:3981549/100‑1 (MQ=255)
                                                              aCCGGTTATTTTGTAGACCGTTTGTGTGAAACAGGGGCGGCTTATGCCGCCCCTTATTCCATCTGGCATGGCATTATTTCCCTTCTTTATATAGatattc                                       >  2:769546/1‑98 (MQ=255)
                                                                cGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGc                                       >  2:3425458/1‑98 (MQ=255)
                                                                cGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTa                                     >  1:3250692/1‑100 (MQ=255)
                                                                  gTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTaaa                                   >  2:3036110/1‑100 (MQ=255)
                                                                   ttATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAAt                                  >  2:1473711/1‑100 (MQ=255)
                                                                    tcttttGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATc                                 <  2:2727513/98‑1 (MQ=255)
                                                                     aTTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATcc                                <  1:1704734/100‑1 (MQ=255)
                                                                     aTTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATcc                                <  1:2066380/100‑1 (MQ=255)
                                                                       tttGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATcaga                              >  1:2727520/1‑97 (MQ=255)
                                                                       tttGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATCCtt                              <  2:1809834/100‑1 (MQ=255)
                                                                       tttGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATCCtt                              <  1:330236/100‑1 (MQ=255)
                                                                        ttGTAAACCGTTTGTGTGAAACAGGGGTGGCTCATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTCTATATAGCTATGCTAAATCCTTa                             >  2:3361087/1‑100 (MQ=255)
                                                                         tGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATCCTTAc                            <  2:4089178/100‑1 (MQ=255)
                                                                         tGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATCCTTAc                            <  2:2408055/100‑1 (MQ=255)
                                                                         tGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATCCTTAc                            <  1:1621288/100‑1 (MQ=255)
                                                                         tGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATCCTTAc                            >  2:2843740/1‑100 (MQ=255)
                                                                              aCCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATCCTTACTTcc                        >  2:490484/1‑99 (MQ=255)
                                                                                cGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATCCTTACTTCCGCa                     <  1:3575523/100‑1 (MQ=255)
                                                                                cGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATCCTTACTTCCGCa                     <  1:1038537/100‑1 (MQ=255)
                                                                                 gTTTGGGTGAAACAGGGGTGGCTGATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCGGTATATAGATATGCTAAATCCTTACTTcctc                      >  2:1209868/1‑96 (MQ=255)
                                                                                  tttGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATCCTTACTTCCGCata                   >  1:1469506/1‑100 (MQ=255)
                                                                                      tgtgAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATCCTTACTTCCGCATATtct               >  2:3449586/1‑100 (MQ=255)
                                                                                       gtgAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATCCTTACTTCCGCATATtctc              >  1:105488/1‑100 (MQ=255)
                                                                                       gtgAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATCCTTACTTCCGCATATtctc              >  1:3038910/1‑100 (MQ=255)
                                                                                       gtgAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATCCTTACTTCCGCATATtctc              <  1:117151/100‑1 (MQ=255)
                                                                                        tgAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATCCTTACTTCCGCATATTctct             >  2:678451/1‑100 (MQ=255)
                                                                                        tgAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATCCTTACTTCCGCATATTctct             >  1:1336101/1‑100 (MQ=255)
                                                                                        tgAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATCCTTACTTCCGCATATTctct             <  1:3350346/100‑1 (MQ=255)
                                                                                        ggaataaGGGGCGGCATAAGCCACCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATCCTTACTTCCGCATATTctct             >  1:2415142/7‑100 (MQ=37)
                                                                                         gAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATCCTTACTTCCGCATATTCTCTg            >  2:1903290/1‑100 (MQ=255)
                                                                                         gAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATCCTTACTTCCGCATATTCTCTg            <  2:3370052/100‑1 (MQ=255)
                                                                                         gAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATCCTTACTTCCGCATATTCTCTg            <  1:1920648/100‑1 (MQ=255)
                                                                                          aaaCAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATCCTTACTTCCGCATATTCTCTGa           <  1:97062/100‑1 (MQ=255)
                                                                                          aaaCAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATCCTTACTTCCGCATATTCTCTGa           <  1:2445164/100‑1 (MQ=255)
                                                                                           ataaGGGGCGGCATAAGCCACCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATCCTTACTTCCGCATATTCTCTGAg          >  1:3468529/4‑100 (MQ=255)
                                                                                           aaCAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATCCTTACTTCCGCATATTCTCTGAg          <  2:4070047/100‑1 (MQ=255)
                                                                                           aaCAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATCCTTACTTCCGCATATTCTCTGAg          >  1:285687/1‑100 (MQ=255)
                                                                                           aaCAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATCCTTACTTCCGCATATTCTCTGAg          >  1:549744/1‑100 (MQ=255)
                                                                                              aGGGGCGGCATAAGCCACCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATCCTTACTTCCGCATATTCTCTGAGCgg       >  1:20911/1‑100 (MQ=255)
                                                                                              aGGGGCGGCATAAGCCACCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATACGCTAAATCCTTACTTCCGCATATTCTCTGAGCgg       >  2:375670/1‑100 (MQ=255)
                                                                                               ggggTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATCCTTACTTCCGCATATTCTCTGAGCggg      <  2:3046116/100‑1 (MQ=255)
                                                                                               ggggTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATCCTTACTTCCGCATATTCTCTGAGCggg      <  1:3765159/100‑1 (MQ=255)
                                                                                                gggTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATCCTTACTTCCGCATATTCTCTGAGCGGGt     >  1:2468380/1‑100 (MQ=255)
                                                                                                gggTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATCCTTACTTCCGCATATTCTCTGAGCGGGt     >  2:1334650/1‑100 (MQ=255)
                                                                                                 ggTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATCCTTACTTCCGCATATTCTCTGAGCGGGTa    >  1:246661/1‑100 (MQ=255)
                                                                                                 ggTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATCCTTACTTCCGCATATTCTCTGAGCGGGTa    >  1:374376/1‑100 (MQ=255)
                                                                                                 ggTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATCCTTACTTCCGCATATTCTCTGAGCGGGTa    <  1:2005863/100‑1 (MQ=255)
                                                                                                   tGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATCCTTACTTCCGCATATTCTCTGAGCGGGTATg  <  2:332077/100‑1 (MQ=255)
                                                                                                   |                                                                                                   
GCGCGACGGCGATGTGGGTCGATTGCTATCGCGATAAACACGCGATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATCCTTACTTCCGCATATTCTCTGAGCGGGTATG  >  NC_000913/2532169‑2532367

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