Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,531,596 (G)6→5 100% coding (134/762 nt) cysZ → sulfate transporter, sulfite inhibited

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,531,5910G.98.3% 524.1 / 1.3 116coding (129/762 nt)cysZsulfate transporter, sulfite inhibited
Reads supporting (aligned to +/- strand):  ref base G (2/0);  new base . (51/63);  total (53/63)
Fisher's exact test for biased strand distribution p-value = 2.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.28e-01

Only 100 of 130 total aligned reads displayed.

CAGCGGTTTTTACTATTTTGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGATGGGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATCCCGACTCTCATGAGTTACGTTCCGGACTGGCTACAATGGCTGAGTTATC  >  NC_000913/2531492‑2531685
                                                                                                   |                                                                                              
cAGCGGTTTTTACTATTTTGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGATg                                                                                                >  1:3408216/1‑100 (MQ=255)
 aGCGGTTTTTACTATTTTGCGCAAGGCTGGAAGCTCGTCTCGCTACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGATgg                                                                                               <  2:595121/100‑1 (MQ=255)
 aGCGGTTTTTACTATTTTGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGATgg                                                                                               >  1:2240752/1‑100 (MQ=255)
 aGCGGTTTTTACTATTTTGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGATgg                                                                                               <  1:1385923/100‑1 (MQ=255)
  gCGGTTTTTACTATTTTGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGATggg                                                                                              >  2:135289/1‑100 (MQ=255)
  gCGGTTTTTACTATTTTGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGATggg                                                                                              >  1:3032048/1‑100 (MQ=255)
   cGGTTTTTACTATTTTGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGATgggg                                                                                             >  1:2431526/1‑100 (MQ=255)
    ggTTTTTACTATTTTGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGATgggg                                                                                             >  1:3796280/1‑99 (MQ=255)
     gTTTTTACTATTTTGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGATgggggc                                                                                           <  1:514665/100‑2 (MQ=255)
     gTTTTTACTATTTTGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGATgggggc                                                                                           >  1:2885754/1‑99 (MQ=255)
       ttttACTATTTTGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGATgggggcgc                                                                                         >  1:1264501/1‑97 (MQ=255)
         ttACTATTTTGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGAT‑GGGGGCGCAt                                                                                      >  1:2019518/1‑100 (MQ=255)
          tACTATTTTGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGAT‑GGGGGCGCAtt                                                                                     >  2:2844675/1‑100 (MQ=255)
           aCTATTTTGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGAT‑GGGGGCGCATTc                                                                                    >  1:1263929/1‑100 (MQ=255)
           aCTATTTTGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGAT‑GGGGGCGCATTc                                                                                    >  1:1658017/1‑100 (MQ=255)
           aCTATTTTGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGAT‑GGGGGCGCATTc                                                                                    <  1:2731757/100‑1 (MQ=255)
              tccgatctcAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGAT‑GGGGGCGCATTCtgg                                                                                 <  2:2618258/92‑1 (MQ=255)
              aTTTTGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGAT‑GGGGGCGCATTCtgg                                                                                 >  1:977163/1‑100 (MQ=255)
              aTTTTGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGAT‑GGGGGCGCATTCtgg                                                                                 <  2:3748743/100‑1 (MQ=255)
              aTTTTGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGAT‑GGGGGCGCATTCtgg                                                                                 >  1:1036289/1‑100 (MQ=255)
               ttttGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGAT‑GGGGGCGCATTCtggt                                                                                <  1:1967013/100‑1 (MQ=255)
                tttGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGAT‑GGGGGCGCATTCtggtg                                                                               <  1:2983585/100‑1 (MQ=255)
                tttCCGCAATGCTGGAAGCTCGGCTCGCAACCGGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGAT‑GGGGGCGCATTCtggtg                                                                               <  1:3408856/100‑1 (MQ=255)
                 ttGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGAT‑GGGGGCGCATTCtggtgg                                                                              <  1:1375454/100‑1 (MQ=255)
                  tGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGAT‑GGGGGCGCATTCTGGTGGc                                                                             <  1:3474501/100‑1 (MQ=255)
                  tGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGAT‑GGGGGCGCATTCTGGTGGc                                                                             <  2:1366361/100‑1 (MQ=255)
                  tGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGAT‑GGGGGCGCATTCTGGTGGc                                                                             >  2:2934130/1‑100 (MQ=255)
                  tGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGAT‑GGGGGCGCATTCTGGTGGc                                                                             >  2:1036825/1‑100 (MQ=255)
                    cgcAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGATGGGGGGCGCATTCTGGTGGct                                                                            >  1:1282846/1‑100 (MQ=255)
                    cgcAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGAT‑GGGGGCGCATTCTGGTGGctc                                                                           >  1:748559/1‑100 (MQ=255)
                    cgcAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGAT‑GGGGGCGCATTCTGGTGGctc                                                                           >  1:265564/1‑100 (MQ=255)
                     gcAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGATGGGGGGCGCATGCTGGTGGctc                                                                           >  1:1524250/1‑100 (MQ=255)
                      cAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGAT‑GGGGGCGCATTCTGGTGGCTCtt                                                                         <  1:2374483/100‑1 (MQ=255)
                       aaGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGAT‑GGGGGCGCATTCTGGTGGCTCttt                                                                        <  2:390284/100‑1 (MQ=255)
                         ggCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGAT‑GGGGGCGCATTCTGGTGGCTCTTTac                                                                      >  1:2250149/1‑100 (MQ=255)
                         ggCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGAT‑GGGGGCGCATTCTGGTGGCTCTTTac                                                                      <  1:856172/100‑1 (MQ=255)
                         ggCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGAT‑GGGGGCGCATTCTGGTGGCTCTTTac                                                                      <  1:2267140/100‑1 (MQ=255)
                          gCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGAT‑GGGGGCGCATTCTGGTGGCTCTTTaca                                                                     >  1:1374087/1‑100 (MQ=255)
                          gCTGCAAGCTTGTCTCGCATCTTGGGATTCGGCGTTTCGTTAGTTTACCGCTGCTGGTCAATATTTTGTTGAT‑GGGGGCGCATTCTGGTGGCTCtt                                                                         <  2:1019639/96‑1 (MQ=255)
                           cTGGAAGCTCGTCTCGCAACCTGGGATTCGGAGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGAT‑GGGGGCGCATTCTGGTGGCTCTTTac                                                                      <  2:148070/98‑1 (MQ=255)
                                aGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGAT‑GGGGGCGCATTCTGGTGGCTCTTTACACAGCTc                                                               <  1:2666715/100‑1 (MQ=255)
                                    cGTCTCGCAACATGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGAT‑GGAGGCGCATTCTGGTGGCTCTTTACACAGCTCGATg                                                           >  2:2004878/1‑100 (MQ=255)
                                     gTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGAT‑GGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATGt                                                          >  2:3106099/1‑100 (MQ=255)
                                      tctcGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGAT‑GGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATGTc                                                         >  2:2269026/1‑100 (MQ=255)
                                       ctcGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGAT‑GGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATGTCt                                                        >  1:3779961/1‑100 (MQ=255)
                                        tcGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGAT‑GGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTg                                                       <  2:1666247/100‑1 (MQ=255)
                                        tcGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGAT‑GGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTg                                                       <  2:2387840/100‑1 (MQ=255)
                                        tcGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGAT‑GGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTg                                                       <  2:737542/100‑1 (MQ=255)
                                          gCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGAT‑GGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTg                                                       <  2:2532868/98‑1 (MQ=255)
                                              ccTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGAT‑GGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATccc                                                 >  1:1178146/1‑100 (MQ=255)
                                               cTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGAT‑GGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATccc                                                 >  2:3609255/1‑99 (MQ=255)
                                               cTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGAT‑GGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATccc                                                 <  2:897134/99‑1 (MQ=255)
                                               cTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGAT‑GGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATCCCg                                                <  1:1936080/100‑1 (MQ=255)
                                               cTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGAT‑GGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATCCCg                                                >  2:2041645/1‑100 (MQ=255)
                                               cTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGAT‑GGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATCCCg                                                >  1:500473/1‑100 (MQ=255)
                                                tGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGAT‑GGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATCCCGa                                               <  2:3232980/100‑1 (MQ=255)
                                                tGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGAT‑GGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATCCCGa                                               >  1:168305/1‑100 (MQ=255)
                                                 gggATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGAT‑GGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATCCCGAc                                              <  2:1262525/100‑1 (MQ=255)
                                                      tCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGAT‑GGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATCCCGACTCTCa                                         <  2:945808/100‑1 (MQ=255)
                                                        ggCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGAT‑GGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATCCCGACTCTCATg                                       <  1:3681350/100‑1 (MQ=255)
                                                          cGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGAT‑GGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATCCCGACTCTCATGAg                                     >  2:924448/1‑100 (MQ=255)
                                                           gTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGAT‑GGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATCCCGACTCTCATGAGt                                    <  2:3779933/100‑1 (MQ=255)
                                                           gTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGAT‑GGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATCCCGACTCTCATGAGt                                    >  1:2781192/1‑100 (MQ=255)
                                                           gTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGAT‑GGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATCCCGACTCTCATGAGt                                    >  2:351530/1‑100 (MQ=255)
                                                            tttCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGAT‑GGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATCCCGACTCTCATGAGtt                                   >  2:1039101/1‑100 (MQ=255)
                                                             ttCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGAT‑GGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATCCCGACTCTCATGAGTTa                                  <  1:135289/100‑1 (MQ=255)
                                                               cGTTATTTTACCGCTGCTGGTCAATATTTTGTTGAT‑GGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATCCCGACTCTCATGAGTTACg                                >  2:1396267/1‑100 (MQ=255)
                                                               cGTTATTTTACCGCTGCTGGTCAATATTTTGTTGAT‑GGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATCCCGACTCTCATGAGTTACg                                >  2:2425570/1‑100 (MQ=255)
                                                                gTTATTTTACCGCTGCTGGTCAATATTTTGTTGAT‑GGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATCCCGACTCTCATGAGTTACGt                               >  1:2855076/1‑100 (MQ=255)
                                                                    ttttACCGCTGCTGGTCAATATTTTGTTGAT‑GGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATCCCGACTCTCATGAGTTACGTTCCg                           <  1:1092228/100‑1 (MQ=255)
                                                                    ttttACCGCTGCTGGTCAATATTTTGTTGAT‑GGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATCCCGACTCTCATGAGTTACGTTCCg                           >  2:3264206/1‑100 (MQ=255)
                                                                    ttttACCGCTGCTGGTCAATATTTTGTTGAT‑GGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATCCCGACTCTCATGAGTTACGTTCCg                           >  1:763065/1‑100 (MQ=255)
                                                                    ttttACCGCTGCTGGTCAATATTTTGTTGAT‑GGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATCCCGACTCTCACGAGTTACGTTcc                            <  1:2512614/99‑1 (MQ=255)
                                                                          cGCTGCTGGTCAATATTTTGTTGAT‑GGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATCCCGACTCTCATGAGTTACGTTCCGGACTgg                     <  2:117106/100‑1 (MQ=255)
                                                                           gctgctGGTCAATATTTTGTTGAT‑GGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATCCCGACTCTCATGAGTTACGTTCCGGACTGGc                    <  1:3088660/100‑1 (MQ=255)
                                                                            ctgctgGTCAATATTTTGTTGAT‑GGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATCCCGACTCTCATGAGTTACGTTCCGGACTGGCt                   <  2:168305/100‑1 (MQ=255)
                                                                            ctgctgGTCAATATTTTGTTGAT‑GGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATCCCGACTCTCATGAGTTACGTTCCGGACTGGCt                   <  1:2627636/100‑1 (MQ=255)
                                                                              gctgGTCAATATTTTGTTGAT‑GGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATCCCGACTCTCATGAGTTACGTTCCGGACTGGCTAc                 <  1:1441256/100‑1 (MQ=255)
                                                                              gctgGTCAATATTCTGTTGAT‑GGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATCCCGACTCTCATGAGTTACGTTCCGGACTGGCTAc                 >  2:2357660/1‑100 (MQ=255)
                                                                                tgGTCAATATTTTGTTGAT‑GGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATCCCGACTCTCATGAGTTACGTTCCGGACTGGCTACa                <  1:506013/99‑1 (MQ=255)
                                                                                    cAATATTTTGTTGAT‑GGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATCCCGACTCTCATGAGTTACGTTCCGGACTGGCTACAATGGc           <  2:2251896/100‑1 (MQ=255)
                                                                                    cAATATTTTGTTGAT‑GGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATCCCGACTCTCATGAGTTACGTTCCGGACTGGCTACAATGGc           <  1:1129793/100‑1 (MQ=255)
                                                                                    cAATATTTTGTTGAT‑GGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATCCCGACTCTCATGAGTTACGTTCCGGACTGGCTACAATGGc           >  1:1761796/1‑100 (MQ=255)
                                                                                    cAATATTTTGTTGAT‑GGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATCCCGACTCTCATGAGTTACGTTCCGGACTGGCTACAATGGc           >  1:3738905/1‑100 (MQ=255)
                                                                                    cAATATTTTGTTGAT‑GGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATCCCGACTCTCATGAGTTACGTTCCGGACTGGCTACAATGGc           >  2:1055931/1‑100 (MQ=255)
                                                                                    cAATATTTTGTTGAT‑GGGGGCGCATGCTGGTGGCTCTTTACACAGCTCGATGTCTGGATCCCGACTCTCATGAGTTACGTTCCGGACTGGCTACAATGGc           <  2:3387372/100‑1 (MQ=255)
                                                                                     aaTATTTTGTTGAT‑GGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATCCCGACTCTCATGAGTTACGTTCCGGACTGGCTACAATGGCt          >  2:2748190/1‑100 (MQ=255)
                                                                                      atatTTTGTTGAT‑GGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATCCCGACTCTCATGAGTTACGTTCCGGACTGGCTACAATGGCTg         <  2:2855056/100‑1 (MQ=255)
                                                                                      atatTTTGTTGAT‑GGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATCCCGACTCTCATGAGTTACGTTCCGGACTGGCTACAATGGCTg         >  1:2542129/1‑100 (MQ=255)
                                                                                      atatTTTGTTCAT‑GGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATCCCGACTCTTATCAGTTACGTTCCGGACTGGCTACAATGGCTg         <  2:2514179/100‑1 (MQ=255)
                                                                                        atTTTGTTGAT‑GGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATCCCGACTCTCATGAGTTACGTTCCGGACTGGCTACAATGGCTGAg       <  2:1016924/100‑1 (MQ=255)
                                                                                        atTTTGTTGAT‑GGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATCCCGACTCTCATGAGTTACGTTCCGGACTGGCTACAATGGCTGAg       <  1:934736/100‑1 (MQ=255)
                                                                                        atTTTGTTGAT‑GGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATCCCGACTCTCATGAGTTACGTTCCGGACTGGCTACAATGGCTGAg       <  2:1204715/100‑1 (MQ=255)
                                                                                        atTTTGTTGAT‑GGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATCCCGACTCTCATGAGTTACGTTCCGGACTGGCTACAATGGCTGAg       >  2:384421/1‑100 (MQ=255)
                                                                                        atTTTGTTGAT‑GGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATCCCGACTCTCATGAGTTACGTTCCGGACTGGCTACAATGGCTGAg       <  1:1549154/100‑1 (MQ=255)
                                                                                         ttttGTTGAT‑GGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATCCCGACTCTCATGAGTTACGTTCCGGACTGGCTACAATGGCTGAGt      <  2:748552/100‑1 (MQ=255)
                                                                                          tttGTTGAT‑GGGGGCGCATTCTGGTAGCTCTTTACACTGCTCGATGTCTGTATCCCGACTCGCATGAGTTACGTTCCGGACTGGCTACAATGGCTGAGtt     >  2:1577186/1‑100 (MQ=255)
                                                                                           ttgttgAT‑GGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATCCCGACTCTCATGAGTTACGTTCCGGACTGGCTACAATGGCTGCGTTa    >  2:1627890/1‑100 (MQ=255)
                                                                                            tgttgAT‑GGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATCCCGACTCTCATGAGTTACGTTCCGGACTGGCTACAATGGCTGAGTTAt   <  2:1374078/100‑1 (MQ=255)
                                                                                             gttgAT‑GGGGGCGCATTCTGGTGGCTCTTTACACAGGTCGATGTCTGGATCCCGACTCTCATGAGTTACGTTGCGGACTGGCTACAAATGCTGAGTTATc  >  2:1485597/1‑100 (MQ=255)
                                                                                                   |                                                                                              
CAGCGGTTTTTACTATTTTGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGATGGGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATCCCGACTCTCATGAGTTACGTTCCGGACTGGCTACAATGGCTGAGTTATC  >  NC_000913/2531492‑2531685

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: ‑