Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,654,545 A→G 33.6% L130P (CTG→CCG)  sseB ← rhodanase‑like enzyme, sulfur transfer from thiosulfate

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,654,5450AG33.6% 184.9 / 166.5 179L130P (CTG→CCG) sseBrhodanase‑like enzyme, sulfur transfer from thiosulfate
Reads supporting (aligned to +/- strand):  ref base A (54/64);  new base G (28/32);  total (83/96)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.65e-01

Only 100 of 182 total aligned reads displayed.

TAAATAAGGTGGTGAGTGAATCAATCATTTGTGCTGGCGGCTCTGCGACTTCCGATAATATCAGCGATTCACCGCCTTCCAGGATTTCCTGGCTGCTCAGCGGATTTCCCTCTTCACCAATCAACAAACTGATTTCACGCGGCATAAATTCTTTACCGGTTGGCAGTTTGGCATTAAGGAAGAGCGTTTCGCCAAGTG  >  NC_000913/2654447‑2654644
                                                                                                  |                                                                                                   
 aaaTAAGGTGGTGAGTGAATCAATCATTTGTGCTGGCGGCTCTGCGACTTCCGATAATATCAGCGATTCACCGCCTTCCAGGATTTCCTGGCTGCTCAGc                                                                                                   >  2:3886351/1‑100 (MQ=255)
   aTAAGGTGGTGAGTGAATCAATCATTTGTGCTGGCGGCTCTGCGACTTCCGATAATATCAGCGATTCACCGCCTTCCAGGATTTCCTGGCTGCTcggcgg                                                                                                 >  2:768271/1‑100 (MQ=255)
   aTAAGGTGGTGAGTGAATCAATCATTTGTGCTGGCGGCTCTGCGACTTCCGATAATATCAGCGATTCACCGCCTTCCAGGATTTCCTGGCTGCTCAGCgg                                                                                                 <  2:822381/100‑1 (MQ=255)
   aTAAGGTGGTGAGTGAATCAATCATTTGTGCTGGCGGCTCTGCGACTTCCGATAATATCAGCGATTCACCGCCTTCCAGGATTTCCTGGCTGCTCAGCgg                                                                                                 <  2:641741/100‑1 (MQ=255)
    tAAGGTGGTGAGTGAATCAATCATTTGTGCTGGCGGCTCTGCGACTTCCGATAATATCAGCGATTCACCGCCTTCCAGGATTTCCTGGCTGCTCAGCGGa                                                                                                >  1:3450955/1‑100 (MQ=255)
    tAAGGTGGTGAGTGAATCAATCATTTGTGCTGGCGGCTCTGCGACTTCCGATAATATCAGCGATTCACCGCCTTCCAGGATTTCCTGGCTGCTCAGCGGa                                                                                                >  1:442542/1‑100 (MQ=255)
     aaGGTGGTGAGTGAATCAATCATTTGTGCTGGCGGCTCTGCGACTTCCGATAATATCAGCGATTCACCGCCTTCCAGGATTTCCTGGCTGCTCAGCGGAt                                                                                               <  2:282249/100‑1 (MQ=255)
       ggtggtGAGTGAATCAATCATTTGTGCTGGCGGCTCTGCGACTTCCGATAATATCAGCGATTCACTGCCTTCCAGGATTTCCTGGCTGCTCAGCGGAttt                                                                                             <  1:560987/100‑1 (MQ=255)
       ggtggtGAGTGAATCAATCATTTGTGCTGGCGGCTCTGCGACTTCCGATAATATCAGCGATTCACCGCCTTCCAGGATTTCCTGGCTGCTCAGCGGAttt                                                                                             <  2:3964462/100‑1 (MQ=255)
          ggtgAGTGAATCAATCATTTGTGCTGGCGGCTCTGCGACTTCCGATAATATCAGCGATTCACCGCCTTCCAGGATTTCCTGGCTGCTCGGCGGATTTccc                                                                                          >  1:3641482/1‑100 (MQ=255)
          ggtgAGTGAATCAATCATTTGTGCTGGCGGCTCTGCGACTTCCGATAATATCAGCGATTCACCGCCTTCCAGGATTTCCTGGCTGCTCAGCGGATTTccc                                                                                          >  2:4074180/1‑100 (MQ=255)
           gtgagtgaATCAATCATTTGTGCTGGCGGCTCTGCGACTTCCGATAATATCAGCGATTCACCGCCTTCCAGGATTTCCTGGCTGCTCAGCGGATTTCCct                                                                                         >  1:2924674/1‑100 (MQ=255)
            tgagtgaATCAATCATTTGTGCTGGCGGCTCTGCGACTTCCGATAATATCAGCGATGCACCGCCTTCCAGGATTTCCTGGCTGCTCAGCGGATTTCCctc                                                                                        >  1:2825705/1‑100 (MQ=255)
             gagtgaATCAATCATTTGTGCTGGCGGCTCTGCGACTTCCGATAATATCAGCGATTCACCGCCTTCCAGGATTTCCTGGCTGCTCGGCGGATTTCCctct                                                                                       >  1:80866/1‑100 (MQ=255)
             gagtgaATCAATCATTTGTGCTGGCGGCTCTGCGACTTCCGATAATATCAGCGATTCACCGCCTTCCAGGATTTCCTGGCTGCTCAGCGGATTTCCctc                                                                                        >  2:1785784/1‑99 (MQ=255)
              agtgaATCAATCATTTGTGCTGGCGGCTCTGCGACTTCCGATAACATCAGCGATTCACCGCCTTCCATGATTCCCTGGCTGCGCAGCGGATTTCCctc                                                                                        >  1:3553734/1‑98 (MQ=255)
               gtgaATCAATCATTTGTGCTGGCGGCTCTGCGACTTCCGATAATATCAGCGATTCACCGCCTTCCAGGATTTCCTGGCTGCTCAGCGGATTTCCCTCTTc                                                                                     <  1:838922/100‑1 (MQ=255)
                tgaATCAATCATTTGTGCTGGCGGCTCTGCGACTTCCGATAATATCAGCGATTCACCGCCTTCCAGGATTTCCTGGCTGCTCGGCGGATTTCCCTCTTCa                                                                                    <  2:296291/100‑1 (MQ=255)
                 gaATCAATCATTTGTGCTGGCGGCTCTGCGACTTCCGATAATATCAGCGATTCACCGCCTTCCAGGATTTCCTGGCTGCTCAGCGGATTTCCCTCTTCAc                                                                                   >  1:4057490/1‑100 (MQ=255)
                   atcaatcaTTTGTGCTGGCGGCTCTGCGACTTCCGATAATATCAGCGATTCACCGCCTTCCAGGATTTCCTGGCTGCTCGGCGGATTTCCCTCTTCACCa                                                                                 <  1:1548783/100‑1 (MQ=255)
                      aatcaTTTGTGCTGGCGGCTCTGCGACTTCCGATAATATCAGCGATTCACCGCCTTCCAGGATTTCCTGGCTGCTCAGCGGATTTCCCTCTTCACCAATc                                                                              >  2:1413917/1‑100 (MQ=255)
                       atcaTTTGTGCTGGCGGCTCTGCGACTTCCGATAATATCAGCGATTCACCGCCTTCCAGGATTTCCTGGCTGCTCGGCGGATTTCCCTCTTCACCAATca                                                                             <  2:4104299/100‑1 (MQ=255)
                       atcaTTTGTGCTGGCGGCTCTGCGACTTCCGATAATATCAGCGATTCACCGCCTTCCAGGATTTCCTGGCTGCTCGGCGGATTTCCCTCTTCACCAATca                                                                             >  1:1854583/1‑100 (MQ=255)
                       atcaTTTGTGCTGGCGGCTCTGCGACTTCCGATAATATCAGCGATTCACCGCCTTCCAGGATTTCCTGGCTGCTCAGCGGATTTCCCTCTTCACCAAtaa                                                                             >  1:2683665/1‑98 (MQ=255)
                       atcaTTTGTGCTGGCGGCTCTGCGACTTCCGATAATATCAGCGATTCACCGCCTTCCAGGATTTCCTGGCTGCTCAGCGGATTTCCCTCTTCACCAATca                                                                             <  2:1100288/100‑1 (MQ=255)
                        tcaTTTGTGCTGGCGGCTCTGCGACTTCTGATAATATCAGTGATTCACCGCCTTCCAGGATTTCCTGGCTGCTCGGCGGATTTCCCTCTTCACCAATca                                                                             <  1:1049379/99‑1 (MQ=255)
                         caTTTGTGCTGGCGGCTCTGCGACTTCCGATAATATCAGCGATTCACCGCCTTCCAGGATTTCCTGGCTGCTCGGCGGATTTCCCTCTTCGCCAATcaac                                                                           >  2:3859503/1‑100 (MQ=255)
                         caTTTGTGCTGGCGGCTCTGCGACTTCCGATAATATCAGCGATTCACCGCCTTCCAGGATTTCCTGGCTGCTCAGCGGATTTCCCTCTTCACCAATcaac                                                                           <  2:2537820/100‑1 (MQ=255)
                         caTTTGTGCTGGCGGCTCTGCGACTTCCGATAATATCAGCGATTCACCGCCTTCCAGGATTTCCTGGCTGCTCAGCGGATTTCCCTCTTCACCAATcaac                                                                           >  1:907301/1‑100 (MQ=255)
                          aTTTGTGCTGGCGGCTCTGCGACTTCCGATAATATCAGCGATTCACCGCCTTCCAGGATTTCCTGGCTGCTCGGCGGATTTCCCTCTTCACCAATcaaca                                                                          >  1:339666/1‑100 (MQ=255)
                            ttGTGCTGGCGGCTCTGCGACTTCCGATAATATCAGCGATTCACCGCCTTCCAGGATTTCCTGGCTGCTCAGCGGATTTCCCTCTTCACCAATCAACaaa                                                                        >  2:243484/1‑100 (MQ=255)
                                gCTGGCGGCTCTGCGACTTCCGATAATATCAGCGATTCACCGCCTTCCAGGATTTCCTGGCTGCTCAGCGGATTTCCCTCTTCACCAATCAACAAACTGa                                                                    <  2:3914204/100‑1 (MQ=255)
                                  tGGCGGCTCTGCGACTTCCGATAATATCAGCGATTCACCGCCTTCCAAGATTTCCTGGCTGCTCAGCGGATTTCCCTCTTCACCAATCAACAAACTGAtt                                                                  <  1:425007/100‑1 (MQ=255)
                                   ggcggcTCTGCGACTTCCGATAATATCAGCGATTCACCGCCTTCCAGGATTTCCTGGCTGCTCAGCGGATTTCCCTCTTCACCAATCAACAAACTGAttt                                                                 <  1:1165075/100‑1 (MQ=255)
                                   ggcggcTCTGCGACTTCCGATAATACCCGCGATTCACCGCCTTCCAGGATTTCCTGGCTGCTCAGCGGATTTCCCTCTTCACCAATcaacaatcaga                                                                    <  1:4186255/97‑6 (MQ=255)
                                    gcggcTCTGCGACTTCCGATAATATCAGCGATTCACCGCCTTCCAGGATTTCCTGGCTGCTCGGCGGATTTCCCTCTTCACCAATCAACAAACTGATTTc                                                                <  2:528190/100‑1 (MQ=255)
                                     cggcTCTGCGACTTCCGATAATATCAGCGATTCACCGCCTTCCAGGATTTCCTGGCTGCTCGGCGGATTTCCCTCTTCACCAATCAACAAACTGATTTCa                                                               >  2:4383560/1‑100 (MQ=255)
                                     cggcTCTGCGACTTCCGATAATATCAGCGATTCACCGCCTTCCAGGATTTCCTGGCTGCTCAGCGGATTTCCCTCTTCACCAATCAACAAACTGATTTCa                                                               <  1:2076980/100‑1 (MQ=255)
                                       gcTCTGCGACTTCCGATAATATCAGCGATTCACCGCCTTCCAGGATTTCCTGGCTGCTCAGCGGATTTCCCTCTTCACCAATCAACAAACTGATTTCAcg                                                             >  2:2105643/1‑100 (MQ=255)
                                        ctctGCGACTTCCGATAATATCAGCGATTCACCGCCTTCCAGGATTTCCTGGCTGCTCGGCGGATTTCCCTCTTCACCAATCAACAAACTGATTTCAcgc                                                            <  2:3183556/100‑1 (MQ=255)
                                        ctctGCGACTTCCGATAATATCAGCGATTCACCGCCTTCCAGGATTTCCTGGCTGCTCGGCGGATTTCCCTCTTCACCAATCAACAAACTGATTTCAcgc                                                            >  2:2516455/1‑100 (MQ=255)
                                        ctctGCGACTTCCGATAATATCAGCGATTCACCGCCTTCCAGGATTTCCTGGCTGCTCGGCGGATTTCCCTCTTCACCAATCAACAAACTGATTTCAcgc                                                            >  1:4323902/1‑100 (MQ=255)
                                        ctctGCGACTTCCGATAATATCAGCGATTCACCGCCTTCCAGGATTTCCTGGCTGCTCAGCGGATTTCCCTCTTCACCAATCAACAAACTGATTTCAcgc                                                            <  2:3472916/100‑1 (MQ=255)
                                        ctctGCGACTTCCGATAATATCAGCGATTCACCGCCTTCCAGGATTTCCTGGCTGCTCAGCGGATTTCCCTCTTCACCAATCAACAAACTGATTTCAcgc                                                            <  1:309793/100‑1 (MQ=255)
                                         tctGCGACTTCCGATAATATCAGCGATTCACCGCCTTCCAGGATTTCCTGGCTGCTCAGCGGATTTCCCTCTTCACCAATCAACAAACTGATTTCAcgcg                                                           >  1:4352209/1‑100 (MQ=255)
                                          ctGCGACTTCCGATAATATCAGCGATTCACCGCCTTCCAGGATTTCCTGGCTGCTCAGCGGATTTCCCTCTTCACCAATCAACAAACTGATTTCACGCgg                                                          <  2:2215042/100‑1 (MQ=255)
                                            gCGACTTCCGATAATATCAGCGATTCACCGCCTTCCAGGATTTCCTGGCTGCTCGGCGGATTTCCCTCTTCACCAATCAACAAACTGATTTCACGCGGCa                                                        <  2:2192738/100‑1 (MQ=255)
                                            gCGACTTCCGATAATATCAGCGATTCACCGCCTTCCAGGATTTCCTGGCTGCTCGGCGGATTTCCCTCTTCACCAATCAACAAACTGATTTCACGCGGCa                                                        >  1:2669306/1‑100 (MQ=255)
                                            gCGACTTCCGATAATATCAGCGATTCACCGCCTTCCAGGATTTCCTGGCTGCTCAGCGGATTTCCCTCTTCACCAATCAACAAACTGATTTCACGCGGCa                                                        >  1:414792/1‑100 (MQ=255)
                                               aCTTCCGATAATATCAGCGATTCACCGCCTTCCAGGATTTCCTGGCTGCTCAGCGGATTTCCCTCTTCACCAATCAACAAACTGATTTCACGCGGCATaa                                                     >  2:4057956/1‑100 (MQ=255)
                                               aCTTCCGATAATATCAGCGATTCACCGCCTTCCAGGATTTCCTGGCTGCTCAGCGGATTTCCCTCTTCACCAATCAACAAACTGATTTCACGCGGCATa                                                      >  2:121670/1‑99 (MQ=255)
                                                   ccGATAATATCAGCGATTCACCGCCTTCCAGGATTTCCTGGCTGCTCGGCGGATTTCCCTCTTCACCAATCAACAAACTGATTTCACGCGGCATAAATTc                                                 <  1:2840676/100‑1 (MQ=255)
                                                   ccGATAATATCAGCGATTCACCGCCTTCCAGGATTTCCTGGCTGCTCAGCGGATTTCCCTCTTCACCAATCAACAAACTGATTTCACGCGGCATAAATTc                                                 >  2:688830/1‑100 (MQ=255)
                                                    cGATAATATCAGCGATTCACCGCCTTCCAGGATTTCCTGGCTGCTCAGCGGATTTCCCTCTTCACCAATCAACAAACTGATTTCACGCGGCATAAATTCt                                                >  2:4227263/1‑100 (MQ=255)
                                                    cGATAATATCAGCGATTCACCGCCTTCCAGGATTTCCTGGCTGCTCAGCGGATTTCCCTCTTCACCAATCAACAAACTGATTTCACGCGGCATAAATTCt                                                <  2:907304/100‑1 (MQ=255)
                                                    cGATAATATCAGCGATTCACCGCCTTCCAGGATTTCCTGGCTGCTCAGCGGATTTCCCTCTTCACCAATCAACAAACTGATTTCACGCGGCATAAATTCt                                                <  1:1208820/100‑1 (MQ=255)
                                                        aataTCAGCGATTCACCGCCTTCCAGGATTTCCTGGCTGCTCAGCGGATTTCCCTCTTCACCAATCAACAAACTGATTTCACGCGGCATAAATTCTTTAc                                            <  1:161967/100‑1 (MQ=255)
                                                        aataTCAGCGATTCACCGCCTTCCAGGATTTCCTGGCTGCTCAGCGGATATCCCTCTTCACCCATCAACACACTGATTTCACGCGGCATAAATTCTTTAc                                            <  1:3512761/100‑1 (MQ=255)
                                                         atatCAGCGATTCACCGCCTTCCAGGATTTCCTGGCTGCTCAGCGGATTTCCCTCTTCACCAATCAACAAACTGATTTCACGCGGCATAAATTCTTTAcc                                           <  2:1190785/100‑1 (MQ=255)
                                                          tttCAGCGATTCAACGCCTTCCAGGATTTCCAGGCTGCTCGGCGGATTTCCCTCTTCACCAATCAACAAACTGATTTCACGCGGCATAAATTCTTTACCg                                          >  2:2968422/3‑100 (MQ=255)
                                                           atCAGCGATTCACCGCCTTCCAGGATTTCCTGGCTGCTCAGCGGATTTCCCTCTTCACCAATCAACAAACTGATTTCACGCGGCATAAATTCTTTACCgg                                         >  2:2607433/1‑100 (MQ=255)
                                                           atCAGCGATTCACCGCCTTCCAGGATTTCCTGGCTGCTCAGCGGATTTCCCTCTTCACCAATCAACAAACTGATTTCACGCGGCATAAATTCTTTACCgg                                         >  2:2607137/1‑100 (MQ=255)
                                                            tCAGCGATTCACCGCCTTCCAGGATTTCCTGGCTGCTCGGCGGATTTCCCTCTTCACCAATCAACAAACTGATTTCACGCGGCATAAATTCTTTACCGGt                                        <  1:38438/100‑1 (MQ=255)
                                                            tCAGCGATTCACCGCCTTCCAGGATTTCCTGGCTGCTCAGCGGATTTCCCTCTTCACCAATCAACAAACTGATTTCACGCGGCATAAATTCTTTACCGGt                                        >  2:607791/1‑100 (MQ=255)
                                                             cAGCGATTCACCGCCTTCCAGGATTTCCTGGCTGCTCAGCGGATTTCCCTCTTCACCAATCAACAAACTGATTTCACGCGGCATAAATTCTTTACCGGtt                                       >  1:54470/1‑100 (MQ=255)
                                                              aGCGATTCACCGCCTTCCAGGATTTCCTGGCTGCTCGGCGGATTTCCCTCTTCACCAATCAACAAACTGATTTCACGCGGCATAAATTCTTTACCGGTTg                                      >  1:4029604/1‑100 (MQ=255)
                                                              aGCGATTCACCGCCTTCCAGGATTTCCTGGCTGCTCAGCGGATTTCCCTCTTCACCAATCAACAAACTGATTTCACGCGGCATAAATTCTTTACCGGTTg                                      <  1:1069673/100‑1 (MQ=255)
                                                              aGCGATTCACCGCCTTCCAGGATTTCCTGGCTGCTCAGCGGATTTCCCTCTTCACCAATCAACAAACTGATTTCACGCGGCATAAATTCTTTACCGGTTg                                      <  2:2789136/100‑1 (MQ=255)
                                                                cGATTCACCGCCTTCCAGGATTTCCTGGCTGCTCGGCGGATTTCCCTCTTCACCAATCAACAAACTGATTTCACGCGGCATAAATTCTTTACCGGTTGGc                                    <  2:1111308/100‑1 (MQ=255)
                                                                cGATTCACCGCCTTCCAGGATTTCCTGGCTGCTCAGCGGATTTCCCTCTTCACCAATCAACAAACTGATTTCACGCGGCATAAATTCTTTACCGGTTGGc                                    <  2:2730196/100‑1 (MQ=255)
                                                                cGATTCACCGCCTTCCAGGATTTCCTGGCTGCTCAGCGGATTTCCCTCTTCACCAATCAACAAACTGATTTCACGCGGCATAAATTCTTTACCGGTTGGc                                    >  2:1152881/1‑100 (MQ=255)
                                                                  aTTCACCGCCTTCCAGGATTTCCTGGCTGCTCAGCGGATTTCCCTCTTCACCAATCAACAAACTGATTTCACGCGGCATAACTTCTTTACCGGTTGGCAg                                  >  1:703079/1‑100 (MQ=255)
                                                                     cACCGCCTTCCAGGATTTCCTGGCTGCTCAGCGGATTTCCCTCTTCACCAATCAACAAACTGATTTCACGCGGCATAAATTCTTTACCGGTTGGCAGttt                               >  1:3355305/1‑100 (MQ=255)
                                                                      aCCGCCTTCCAGTAGTTCCTGGCGGCTCAGCGGATTTCCCTCTTCACCAATCAGCAAACTGATTTCACGCGGCATAAATTCTTTACCGGTTGGCAGTTTg                              <  2:31960/100‑1 (MQ=255)
                                                                      aCCGCCTTCCAGGATTTCCTGGCTGCTCGGCGGATTTCCCTCTTCACCAATCAACAAACTGATTTCACGCGGCATAAATTCTTTACCGGTTGGCAGTTTg                              >  2:1833864/1‑100 (MQ=255)
                                                                      aCCGCCTTCCAGGATTTCCTGGCTGCTCAGCGGATTTCCCTCTTCACCAATCAACAAACTGATTTCACGCGGCATAAATTCTTTACCGGTTGGCAGTTTg                              <  1:4187719/100‑1 (MQ=255)
                                                                      aCCGCCTTCCAGGATTTCCTGGCTGCTCAGCGGATTTCCCTCTTCACCAATCAACAAACTGATTTCACGCGGCATAAATTCTTTACCGGTTGGCAGTTTg                              >  2:2218436/1‑100 (MQ=255)
                                                                       ccGCCTTCCAGGATTTCCTGGCTGCTCAGCGGATTTCCCTCTTCACCAATCAACAAACTGATTTCACGCGGCATAAATTCTTTACCGGTTGGCAGttt                               >  1:1460252/1‑98 (MQ=255)
                                                                       ccGCCTTCCAGGATTTCCTGGCTGCTCAGCGGATTTCCCTCTTCACCAATCAACAAACTGATTTCACGCGGCATAAATTCTTTACCGGTTGGCAGTTTgg                             >  2:1398610/1‑100 (MQ=255)
                                                                       ccGCCTTCCAGGATTTCCTGGCTGCTCAGCGGATTTCCCTCTTCACCAATCAACAAACTGATTTCACGCGGCATAAATTCTTTACCGGTTGGCAGTTTgg                             <  2:1646028/100‑1 (MQ=255)
                                                                        cGCCTTCCAGGATTTCCTGGCTGCTCGGCGGATTTCCCTCTTCACCAATCAACAAACTGATTTCACGCGGCATAAATTCTTTACCGGTTGGCAGTTTGGc                            <  2:1825859/100‑1 (MQ=255)
                                                                        cGCCTTCCAGGATTTCCTGGCTGCTCAGCGGATTTCCCTCTTCACCAATCAACAAACTGATTTCACGCGGCATAAATTCTTTACCGGTTGGCAGTTTGGc                            <  1:1377111/100‑1 (MQ=255)
                                                                         gCCTTCCAGGATTTCCTGGCTGCTCGGCGGATTTCCCTCTTCACCAATCAACAAACTGATTTCACGCGGCATAAATTCTTTACCGGTTGGCAGTTTGGCa                           <  1:2516442/100‑1 (MQ=255)
                                                                         gCCTTCCAGGATTTCCTGGCTGCTCAGCGGATTTCCCTCTTCACCAATCAACAAACTGATTTCACGCGGCATAAATTCTTTACCGGTTGGCAGTTTGGCa                           <  2:1802758/100‑1 (MQ=255)
                                                                            ttCCAGGATTTCCTGGCTGCTCGGCGGATTTCCCTCTTCACCAATCAACAAACTGATTTCACGCGGCATAAATTCTTTACCGGTTGGCAGTTTGGCATTa                        <  1:3756080/100‑1 (MQ=255)
                                                                            ttCCAGGATTTCCTGGCTGCTCAGCGGATTTCCCTCTTCACCAATCAACAAACTGATTTCACGCGGCATAAATTCTTTACCGGTTGGCAGTTTGGCATTa                        <  2:1840033/100‑1 (MQ=255)
                                                                            ttCCAGGATTTCCTGGCTGCTCAGCGGATTTCCCTCTTCACCAATCAACAAACTGATTTCACGCGGCATAAATTCTTTACCGGTTGGCAGTTTGGAATTa                        >  1:3992248/1‑100 (MQ=255)
                                                                             gccAGGAAATCCTGGCTGCTCGGCGGATTTCCCTCTTCACCAATCAACAAACTGATTTCACGCGGCATAAATTCTTTACCGGTTGGCAGTTTGGCATTaa                       >  2:4302065/2‑100 (MQ=255)
                                                                              ccAGGATTTCCTGGCTGCTCAGCGGATTTCCCTCTTCACCAATCAACAAACTGATTTCACGCGGCATAAATTCTTTACCGGTTGGCAGTTTGGCATTAAg                      <  1:2057696/100‑1 (MQ=255)
                                                                              ccAGGAATTCCTGGCTGCGCGGCGGATTTCCCTCTCCACCAATCGACAAACTGATTTCACACGGCATAAATTCTTTACCTGTTGGCAGTTTGGCATTAAg                      <  1:3494875/100‑1 (MQ=255)
                                                                                    tttCCTGGCTGCTCGGCGGATTTCCCTCTTCACCAATCAACAAACTGATTTCACGCGGCATAAATTCTTTACCGGTTGGCAGTTTGGCATTAAGGAagag                <  2:339667/100‑1 (MQ=255)
                                                                                      tCCTGGCTGCTCAGCGGATTTCCCTCTTCACCAATCAACAAACTGATTTCACGCGGCATAAATTCTTTACCGGTTGGCAGTTTGGCATTAAGGAAGAGCg              >  1:2258095/1‑100 (MQ=255)
                                                                                      tCCTGGCTGCTCAGCGGATTTCCCTCATCACCAATCAACAGACTGATTTCACGCGGCATAAATTCTTTACCGGTTGGCAGTTTGGCATTAAGGAAGAGCg              <  1:3074365/100‑1 (MQ=255)
                                                                                        cTGGCTGCTCAGCGGATTTCCCTCTTCACCAATCAACAAACTGATTTCACGCGGCATAAATTCTTTACCGGTTGGCAGTTTGGCATTAAGGAAGAGCg              >  1:1772176/1‑98 (MQ=255)
                                                                                          ggCTGCTCGGCGGATTTCCCTCTTCACCAATCAACAAACTGATTTCACGCGGCATAAATTCTTTACCGGTTGGCAGTTTGGCATTAAGGAAGAGCGTTTc          >  1:398465/1‑100 (MQ=255)
                                                                                          ggCTGCTCGGCGGATTTCCCTCTTCACCAATCAACAAACTGATTTCACGCGGCATAAATTCTTTACCAGTTGGCAGTTTGGCATTAAGGAagagtg              >  1:2774385/1‑94 (MQ=255)
                                                                                           gctgctCGGCGGATTTCCCTCTTCACCAATCAACAAACTGATTTCACGCGGCATAAATTCTTTACCGGTTGGCAGTTTGGCATTAAGGAAGAGCGTTTCg         <  2:2272363/100‑1 (MQ=255)
                                                                                                tCAGCGGATTTCCCTCTTCACCAATCAACAAACTGATTTCACGCGGCATAAATTCTTTACCGGTTGGCAGTTTGGCAGTAAGGAAGAGCGTTTCGCCAAg    <  1:2555102/100‑1 (MQ=255)
                                                                                                 cAGCGGATTTCCCTCTTCACCAATCAACAAACTGATTTCACGCGGCATAAATTCTTTACCGGTTGGCAGTTTGGCATTAAGGAAGAGCGTTTCGCCAAgt   >  1:1272360/1‑100 (MQ=255)
                                                                                                  aGCGGATTTCCCTCTTCACCAATCAACAAACTGATTTCACGCGGCATAAATTCTTTACCGGTTGGCAGTTTGGCATTAAGGAAGAGCGTTTCGCCAAgtg  <  2:1710363/100‑1 (MQ=255)
                                                                                                  |                                                                                                   
TAAATAAGGTGGTGAGTGAATCAATCATTTGTGCTGGCGGCTCTGCGACTTCCGATAATATCAGCGATTCACCGCCTTCCAGGATTTCCTGGCTGCTCAGCGGATTTCCCTCTTCACCAATCAACAAACTGATTTCACGCGGCATAAATTCTTTACCGGTTGGCAGTTTGGCATTAAGGAAGAGCGTTTCGCCAAGTG  >  NC_000913/2654447‑2654644

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: