Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,654,545 A→G 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,5450AG98.9% 645.1 / ‑4.5 174L130P (CTG→CCG) sseBrhodanase‑like enzyme, sulfur transfer from thiosulfate
Reads supporting (aligned to +/- strand):  ref base A (0/2);  new base G (80/92);  total (80/94)
Fisher's exact test for biased strand distribution p-value = 5.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 175 total aligned reads displayed.

AAATAAGGTGGTGAGTGAATCAATCATTTGTGCTGGCGGCTCTGCGACTTCCGATAATATCAGCGATTCACCGCCTTCCAGGATTTCCTGGCTGCTCAGCGGATTTCCCTCTTCACCAATCAACAAACTGATTTCACGCGGCATAAATTCTTTACCGGTTGGCAGTTTGGCATTAAGGAAGAGCGTTTCGCCAA  >  NC_000913/2654448‑2654641
                                                                                                 |                                                                                                
aaaTAAGGTGGTGAGTGAATCAATCATTTGTGCTGGCGGCTCTGCGACTTCCGATAATATCAGCGATTCACCGCCTTCCAGGATTTCCTGGCTGCTCGGc                                                                                                <  2:3305416/100‑1 (MQ=255)
 aaTAAGGTGGTGAGTGAATCAATCATTTGTGCTGGCGGCTCTGCGACTTCCGATAATATCAGCGATTCACCGCCTTCCAGGATTTCCTGGCTGCTCGGCg                                                                                               <  1:1896074/100‑1 (MQ=255)
   tAAGGTGGTGAGTGAATCAATCATTTGTGCTGGCGGCTCTGCGACTTCCGATAATATCAGCGATTCACCGCCTTCCAGGATTTCCTGGCTGCTCGGCGGa                                                                                             >  2:649905/1‑100 (MQ=255)
     aGGTGGTGAGTGAATCAATCATTTGTGCTGGCGGCTCTGCGACTTCCGATAATATCAGCGATTCACCGCCTTCCAGGATTTCCTGGCTGCTCGGCGGAtt                                                                                           <  1:2864299/100‑1 (MQ=255)
         ggtgAGTGAATCAATCATTTGTGCTGGCGGCTCTGCGACTTCCGATAATATCAGCGATTCACCGCCTTCCAGGCTTTCCTGGCTGCTCAGCGGATTTccc                                                                                       <  1:2076575/100‑1 (MQ=255)
         ggtgAGTGAATCAATCATTTGTGCTGGCGGCTCTGCGACTTCCGATAATATCAGCGATTCACCGCCTTCCAGGATTTCCTGGCTGCTCGGCGGATTTccc                                                                                       <  2:1821615/100‑1 (MQ=255)
         ggtgAGTGAATCAATCATTTGTGCTGGCGGCTCTGCGACTTCCGATAATATCAGCGATTCACCGCCTTCCAGGATTTCCTGGCTGCTCGGCGGATTTccc                                                                                       <  1:1389559/100‑1 (MQ=255)
           tgagtgaATCAATCATTTGTGCTGGCGGCTCTGCGACTTCCGATAATATCAGCGATTCACCGCCTTCCAGGATTTCCTGGCTGCTCGGCGGATTTCCctc                                                                                     <  1:3154421/100‑1 (MQ=255)
           tgagtgaATCAATCATTTGTGCTGGCGGCTCTGCGACTTCCGATAATATCAGCGATTCACCGCCTTCCAGGATTTCCTGGCTGCTCGGCGGATTTCCctc                                                                                     >  1:885844/1‑100 (MQ=255)
            gagtgaATCAATCATTTGTGCTCGCGGCTCTACGACTACCGATAATACCAGCGATTCACCGCCTTCCAGGATTTCCTGGCTGCTCGGCGGATTTCCctct                                                                                    <  2:2299529/100‑1 (MQ=255)
             agtgaATCAATCATTTGTGCTGGCGGCTCTGCGACTTCCGATAATATCAGCGATTCACCGCCTTCCAGGATTTCCTGGCTGCTCGGCGGATTTCCCTCtt                                                                                   >  2:1016116/1‑100 (MQ=255)
             agtgaATCAATCATTTGTGCGGGCGGCTCTGCGACTTCCGATAATATCAGCGATTCACCGCCTTCCAGGATTTCCTGGCTGCTCGGCGGATTTCCCTCtt                                                                                   >  2:2088554/1‑100 (MQ=255)
              gtgaATCAATCATTTGTGCTGGCGGCTCTGCGAGTTCCGATAATATCAGCGATTCACCGCCTTCCAGGATTTCCTGGCTGCTCGGCGGATTTccc                                                                                       <  2:679812/95‑1 (MQ=255)
              gtgaATCAATCATTTGTGCTGGCGGCTCTGCGACTTCCGATAATATCAGCGATTCACCGCCTTCCTGGATTTCCTGGCTGCTCGGCGGATTTCCCTCTTc                                                                                  >  2:2845337/1‑100 (MQ=255)
              gtgaATCAATCATTTGTGCTGGCGGCTCTGCGACTTCCGATAATATCAGCGATTCACCGCCTTCCAGGATTTCCTGGCTGCTCGGCGGATTTCCCTCTTc                                                                                  <  1:695148/100‑1 (MQ=255)
              gtgaATCAATCATTTGTGCTGGCGGCTCTGCGACTTCCGATAATATCAGCGATTCACCGCCTTCCAGGATTTCCTGGCTGCTCGGCGGATTTCCCTCTTc                                                                                  >  2:2506745/1‑100 (MQ=255)
              gtgaATCAATCATTTGTGCTGGCGGCTCTGCGACTTCCGATAATATCAGCGATTCACCGCCTTCCAGGATTTCCTGGCTGCTCGGCGGATTTCCCTCTTc                                                                                  >  2:210965/1‑100 (MQ=255)
              gtgaATCAATCATTTGTGCTGGCGGCTCTGCGACTTCCGATAATATCAGCGATTCACCGCCTTCCAGGATTTCCTGGCTGCTCGGCGGATTTCCCTCTTc                                                                                  >  2:836104/1‑100 (MQ=255)
              gtgaATCAATCATTTGTGCTGGCGGCTCTGCGACTTCCGATAATATCAGCGATTCACCGCCTTCCAGGATTTCCTGGCTGCTCGGCGGATTTCCCTCTTc                                                                                  <  1:1858238/100‑1 (MQ=255)
               tgaATCAATCATTTGTGCTGGCGGCTCTGCGACTTCCGATAATATCAGCGATTCACCGCCTTCCAGGATTTCCTGGCTGCTCGGCGGATTTCCCTCTTCa                                                                                 <  1:2032935/100‑1 (MQ=255)
               tgaATCAATCATTTGTGCTGGCGGCTCTGCGACTTCCGATAATATCAGCGATTCACCGCCTTCCAGGATTTCCTGGCTGCTCGGCGGATTTCCCTCTTCa                                                                                 <  1:1814376/100‑1 (MQ=255)
                gaATCAATCATTTGTGCTGGCGGCTCTGCGACTTCCGATAATATCAGCGATTCACCGCCTTCCAGGATTTCCTGGCTGCTCGGCGGATTTCCCTCTTCAc                                                                                <  2:550582/100‑1 (MQ=255)
                 aatcaatcATTTGTGCTGGCGGCTCTGCGACTTCCGATAATATCAGCGATTCACCGCCTTCCAGGATTTCCTGGCTGCTCGGCGGATTTCCCTCTTCAcc                                                                               <  1:649905/100‑1 (MQ=255)
                 aatcaatcATTTGTGCTGGCGGCTCTGCGACTTCCGATAATATCAGCGATTCACCGCCTTCCAGGATTTCCTGGCTGCTCGGCGGATTTCCCTCTTCAcc                                                                               <  1:1530286/100‑1 (MQ=255)
                  atcaatcaTTTGTGCTGGCGGCTCTGCGACTTCCGATAATATCAGCGATTCACCGCCTTCCAGGATTTCCTGGCTGCTCGGCGGATTTCCCTCTTCACCa                                                                              >  2:1941365/1‑100 (MQ=255)
                      atcaTTTGTGCTGGCGGCTCTGCGACTTCCGATAATATCAGCGATTCACCGCCTTCCAGGATTTCCTGGCTGCTCGGCGGATTTCCCTCTTCACCAATca                                                                          >  1:991739/1‑100 (MQ=255)
                      atcaTTTGTGCTGGCGGCTCTGCGACTTCCGATAATATCAGCGATTCACCGCCTTCCAGGATTCCCTGGCTGCTCGGCGGATTTCCCTCTTCACCAATca                                                                          <  1:2394729/100‑1 (MQ=255)
                         cgatctgCTGGCGGCTCTGCGACTTCCGATAATATCAGCGATTCACCGCCTTCCAGGATTTCCTGGCTGCTCGGCGGATTTCCCTCTTCACCAATcaaca                                                                       <  2:3098649/95‑1 (MQ=255)
                         aTTTGTGCTGGCGGCTCTGCGACTTCCGATAATATCAGCGATTCACCGCCTTCCAGGATTTCCTGGCTGCTCGGCGGATTTCCCTCTTCACCAATcaaca                                                                       <  1:3399654/100‑1 (MQ=255)
                             gtgCTGGCGGCTCTGCGACTTCCGATAATATCAGCGATTCACCGCCTTCCAGGATTTCCTGGCTGCTCGGCGGATTTCCCTCTTCACCAATCAACAAACt                                                                   <  2:3009017/100‑1 (MQ=255)
                             gtgCTGGCGGCTCTGCGACTTCCGATAATATCAGCGATTCACCGCCTTCCAGGATTTCCTGGCTGCTCGGCGGATTTCCCTCTTCACCAATCAACAAACt                                                                   >  1:722957/1‑100 (MQ=255)
                               gCTGGCGGCTCTGCGACTTCCGATAATATCAGCGATTCACCGCCTTCCAGGATTTCCTGGCTGCTCGGCGGATTTCCCTCTTCACCAATcaacaagatcg                                                                 >  1:3098661/1‑95 (MQ=255)
                               gCTGGCGGCTCTGCGACTTCCGATAATATCAGCGATTCACCGCCTTCCAGGATTTCCTGGCTGCTCGGCGGATTTCCCTCTTCACCAATCAACAAACTGa                                                                 >  2:2631221/1‑100 (MQ=255)
                                cTGGCGGCTCTGCGACTTCCGATAATATCAGCGATTCACCGCCTTCCAGGATTTCCTGGCTGCTCGGCGGATTTCCCTCTTCACCAATCAACAAACTGAt                                                                >  1:977723/1‑100 (MQ=255)
                                 tGGCGGCTCTGCGACTTCCGATAATATCAGCGATTCACCGCCTTCCAGGATTTCCTGGCTGCTCGGCGGATTTCCCTCTTCACCAATCAACAAACTGAtt                                                               <  1:618008/100‑1 (MQ=255)
                                     ggcTCTGCGACTTCCGATAATATCAGCGATTCACCGCCTTCCAGGATTTCCTGGCTGCTCGGCGGATTTCCCTCTTCACCAATCAACAAACTGATTTCAc                                                           <  1:1217984/100‑1 (MQ=255)
                                     ggcTCTGCGACTTCCGATAATATCAGCGATTCACCGCCTTCCAGGATTTCCTGGCTGCTCGGCGGATTTCCCTCTTCACCAATCAACAAACTGATTTCAc                                                           >  1:153052/1‑100 (MQ=255)
                                       ctctGCGACTTCCGATAATATCAGCGATTCACCGCCTTCCAGGATTTCCTGGCTGCTCGGCGGATTTCCCTCTTCACCAATCAACAAACTGATTTCAcgc                                                         <  1:1009914/100‑1 (MQ=255)
                                       ctctGCGACTTCCGATAATATCAGCGATTCACCGCCTTCCAGGATTTCCTGGCTGCTCGGCGGATTTCCCTCTTCACCAATCAACAAACTGATTTCAcgc                                                         >  1:1429028/1‑100 (MQ=255)
                                       ctctGCGACTCCCGATAATATCAGCGATTCACCGCCTTCCAGGATTTCCTGGCTGCTCGGCGGATTTCCCTCTTCACCAATCAACAAACTGATTTCAcgc                                                         >  2:357281/1‑100 (MQ=255)
                                        tctGCGACTTCCGATAATATCAGCGATTCACCGCCTTCCAGGATTTCCTGGCTGCTCGGCGGATTTCCCTCTTCACCAATCAACAAACTGATTTCAcgcg                                                        <  1:853732/100‑1 (MQ=255)
                                        tctGCGACTTCCGATAATATCAGCGATTCACCGCCTTCCAGGATTTCCTGGCTGCTCGGCGGATTTCCCTCTTCACCAATCAACAAACTGATTTCAcgcg                                                        >  1:415548/1‑100 (MQ=255)
                                         ctGCGACTTCCGATAATATCAGCGATTCACCGCCTTCCAGGATTTCCTGGCTGCTCGGCGGATTTCCCTCTTCACCAATCAACAAACTGATTTCACGCgg                                                       <  1:899366/100‑1 (MQ=255)
                                          tGCGACTTCCGATAATATCAGCGATTCACCGCCTTCCAGGATTTCCTGGCTGCTCGGCGGATTTCCCTCTTCACCAATCAACAAACTGATTTCACGCGGc                                                      >  2:3273931/1‑100 (MQ=255)
                                          tGCGACTTCCGATAATATCAGCGATTCACCGCCTTCCAGGATTTCCTGGCTGCTCGGCGGATTTCCCTCTTCACCAATCAACAAACTGATTTCACGCGGc                                                      <  2:841952/100‑1 (MQ=255)
                                           gCGACTTCCGATAATATCAGCGATTCACCGCCTTCCAGGATTTCCTGGCTGCTCGGCGGATTTCCCTCTTCACCAATCAACAAACTGATTTCACGCGGCa                                                     >  2:423231/1‑100 (MQ=255)
                                           gCGACTCCCGATAATATCAGCGATTCACCGCCTTCCAGGATTTCCTGGCTGCTCGGCGGATTTCCCTCTTCACCAATCAACAAACTGATTTCACGCGGCa                                                     <  1:357282/100‑1 (MQ=255)
                                             gACTTCCGATAATATCAGCGATTCACCGCCTTCCAGGATTTCCTGGCTGCTCGGCGGATTTCCCTCTTCACCAATCAACAAACTGATTTCACGCGGCATa                                                   <  1:2164193/100‑1 (MQ=255)
                                              aCTTCCGATAATATCAGCGATTCACCGCCTTCCAGGATTTCCTGGCTGCTCGGCGGATTTCCCTCTTCACCAATCAACAAACTGATTTCACGCGGCATaa                                                  >  2:387331/1‑100 (MQ=255)
                                                   cGATAATATCAGCGATTCACCGCCTTCCAGGATTTCCTGGCTGCTCGGCGGATTTCCCTCTTCACCAATCAACAAACTGATTTCACGCGGCATAAATTCt                                             >  2:451558/1‑100 (MQ=255)
                                                   cGATAATATCAGCGATTCACCGCCTTCCAGGATTTCCTGGCTGCTCGGCGGATTTCCCTCTTCACCAATCAACAAACTGATTTCACGCGGCATAAATTCt                                             >  2:313876/1‑100 (MQ=255)
                                                   cGATAATATCAGCGATTCACCGCCTTCCAGGATTTCCTGGCTGCTCGGCGGATTTCCCTCTTCACCAATCAACAAACTGATTTCACGCGGCATAAATTCt                                             <  1:836104/100‑1 (MQ=255)
                                                   cGATAATATCAGCGATTCACCGCCTTCCAGGATTTCCTGGCTGCTCGGCGGATTTCCCTCTTCACCAATCAACAAACTGATTTCACGCGGCATAAATTCt                                             >  1:2231387/1‑100 (MQ=255)
                                                    gATAATATCAGCGATTCACCGCCTTCCAGGATTTCCTGGCTGCTCGGCGGATTTCCCTCTTCACCAATCAACAAACTGATTTCACGCGGCATAAATTCtt                                            <  2:2245769/100‑1 (MQ=255)
                                                    gATAATATCAGCGATTCACCGCCTTCCAGGATTTCCTGGCTGCTCGGCGGATTTCCCTCTTCACCAATCAACAAACTGATTTCACGCGGCATAAATTCtt                                            <  2:977722/100‑1 (MQ=255)
                                                     ataataTCAGCGATTCACCGCCTTCCAGGATTTCCTGGCTGCTCGGCGGATTTCCCTCTTCACCAATCAACAAACTGATTTCACGCGGCATAAATTCttt                                           <  1:2631229/100‑1 (MQ=255)
                                                       aataTCAGCGATTCACCGCCTTCCAGGATTTCCTGGCTGCTCGGCGGATTTCCCTCTTCACCAATCAACAAACTGATTTCACGCGGCATAAATTCTTTAc                                         <  1:769747/100‑1 (MQ=255)
                                                       aataTCAGCGATTCACCGCCTTCCAGGATTTCCTGGCTGCTCGGCGGATTTCCCTCTTCACCAATCAACAAACTGATTTCACGCGGCATAAATTCTTTAc                                         <  1:769291/100‑1 (MQ=255)
                                                        atatCAGCGATTCACCGCCTTCCAGGATTTCCTGGCTGCTCGGCGGATTTCCCTCTTCACCAATCAACAAACTGATTTCACGCGGCATAAATTCTTTAcc                                        <  1:217054/100‑1 (MQ=255)
                                                         tatCAGCGATTCACCGCCTTCCAGGATTTCCTGGCTGCTCGGCGGATTTCCCTCTTCACCAATCAACAAACTGATTTCACGCGGCATAAATTCTTTACCg                                       <  2:153053/100‑1 (MQ=255)
                                                          atCAGCGATTCACCGCCTTCCAGGATTTCCTGGCTGCTCGGCGGATTTCCCTCTTCACCAATCAACAAACTGATTTCACGCGGCATAAATTCTTTACCgg                                      >  1:2080660/1‑100 (MQ=255)
                                                            cAGCGATTCACCGCCTTCCAGGATTTCCTGGCTGCTCGGCGGATTTCCCTCTTCACCAATCAACAAACTGATTTCACGCGGCATAAATTCTTTACCGGtt                                    >  2:1708603/1‑100 (MQ=255)
                                                            cAGCGATTCACCGCCTTCCAGGATTTCCTGGCTGCTCGGCGGATTTCCCTCTTCACCAATCAACAAACTGATTTCACGCGGCATAAATTCTTTACCGGtt                                    <  2:753306/100‑1 (MQ=255)
                                                            cAGCGATTCACCGCCTTCCAGGATTTCCTGGCTGCTCGGCGGATTTCCCTCTTCACCAATCAACAAACTGATTTCACGCGGCATAAATTCTTTACCGGtt                                    >  1:1487488/1‑100 (MQ=255)
                                                            cAGCGATTCACCGCCTTCCAGGATTTCCTGGCTGCTCGGCGGATTTCCCTCTTCACCAATCAACAAACTGATTTCACGCGGCATAAATTCTTTACCGGtt                                    >  2:49536/1‑100 (MQ=255)
                                                             aGCGATTCACCGCCTTCCAGGATTTCCTGGCTGCTCGGCGGATTTCCCTCTTCACCAATCAACAATCTGATTTCACGCGGCATAAATTCTTTACCGGTTg                                   <  1:1145544/100‑1 (MQ=255)
                                                               cGATTCACCGCCTTCCAGGATTTCCTGGCTGCTCGGCGGATTTCCCTCTTCACCAATCAACAAACTGATTTCACGCGGCATAAATTCTTTACCGGTTGGc                                 <  2:923125/100‑1 (MQ=255)
                                                                   tCACCGCCTTCCAGGATTTCCTGGCTGCTCGGCGGATTTCCCTCTTCACCAATCAACAAACTGATTTCACGCGGCATAAATTCTTTACCGGTTGGCAGtt                             >  1:651152/1‑100 (MQ=255)
                                                                   tCACCGCCTTCCAGGATTTCCTGGCTGCTCGGCGGATTTCCCTCTTCACCAATCAACAAACTGATTTCACGCGGCATAAATTCTTTACCGGTTGGCAGtt                             <  1:3253908/100‑1 (MQ=255)
                                                                    cACCGCCTTCCAGGATTTCCTGGCTGCTCGGCGGATTTCCCTCTTCACCAATCAACAAACTGATTTCACGCGGCATAAATTCTTTACCGGTTGGCAGttt                            <  2:507558/100‑1 (MQ=255)
                                                                    cACCGCCTTCCAGGATTTCCTGGCTGCTCGGCGGATTTCCCTCTTCACCAATCAACAAACTGATTTCACGCGGCATAAATTCTTTACCGGTTGGCAGttt                            >  1:2474095/1‑100 (MQ=255)
                                                                    cACCGCCTTCCAGGATTTCCTGGCTGCTCGGCGGATTTCCCTCTTCACCAATCAACAAACTGATTTCACGCGGCATAAATTCTTTACCGGTTGGCAGttt                            >  1:2349358/1‑100 (MQ=255)
                                                                    cACCGCCTTCCAGGATTTCCTGGCTGCTCGGCGGATTTCCCTCTTCACCAATCAACAAACTGATTTCACGCGGCATAAATTCTTTACCGGTTGGCAGttt                            >  2:23985/1‑100 (MQ=255)
                                                                    cACCGCCTTCCAGGATTTCCTGGCTGCTCGGCGGATTTCCCTCTTCACCAATCAACAAACTGATTTCACGCGGCATAAATTCTTTACCGGTTGGCAGttt                            <  1:173680/100‑1 (MQ=255)
                                                                    cACCGCCTTCCAGGATTTCCTGGCTGCTCGGCGGATTTCCCTCTTCACCAATCAACAAACTGATTTCACGCGGCATAAATTCTTTACCGGTTGGCAGttt                            <  2:1521053/100‑1 (MQ=255)
                                                                     aCCGCCTTCCAGGATTTCCTGGCTGCTCGGCGGATTTCCCTCTTCACCAATCAACAAACTGATTTCACGCGGCATAAATTCTTTACCGGTTGGCAGTTTg                           >  2:3164793/1‑100 (MQ=255)
                                                                      ccGCCTTCCAGGATTTCCTGGCTGCTCGGCGGATTTCCCTCTTCACCAATCAACAAACTGATTTCACGCGGCATAAATTCTTTACCGGTTGGCAGTTTgg                          >  1:1102315/1‑100 (MQ=255)
                                                                      ccGCCTTCCAGGATTTCCTGGCTGCTCGGCGGATTTCCCTCTTCACCAATCAACAAACTGATTTCACGCGGCATAAATTCTTTACCGGTTGGCAGTTTg                           <  1:495426/99‑1 (MQ=255)
                                                                      ccGCCTTCCAGGATTTCCTGGCTGCTCGGCGGATTTCCCTCTTCACCAATCAACAAACTGATTTCACGCGGCATAAATTCTTTACCGGTTGGCAGTTTg                           <  1:2993690/99‑1 (MQ=255)
                                                                       cGCCTTCCTGGATTTCCTGGCTGCTCGGCGGATTTCCCTCTTCACCAATCAACAAACTGATTTCACGCGGCATAAATTCTTTACCGGTTGGCAGTTTGGc                         <  1:2845347/100‑1 (MQ=255)
                                                                        gCCTTCCAGGATTTCCTGGCTGCTCGGCGGATTTCCCTCTTCACCAATCAACAAACTGATTTCACGCGGCATAAATTCTTTACCGGTTGGCAGTTTGGCa                        <  1:1016117/100‑1 (MQ=255)
                                                                         ccTTCCAGGATTTCCTGGCTGCTCGGCGGATTTCCCTCTTCACCAATCAACAAACTGATTTCACGCGGCATAAATTCTTTACCGGTTGGCAGTTTGGCAt                       <  1:495157/100‑1 (MQ=255)
                                                                           ttCCAGGATTTCCTGGCTGCTCGGCGGATTTCCCTCTTCACCAATCAACAAACTGATTTCACGCGGCATAAATTCTTTACCGGTTGGCAGTTTGGCATTa                     >  1:928824/1‑100 (MQ=255)
                                                                           ttCCAGGATTTCCTGGCTGCTCGGCGGATTTCCCTCTTCACCAATCAACAAACTGATTTCACGCGGCATAAATTCTTTACCGGTTGGCAGTTTGGCATTa                     >  2:601799/1‑100 (MQ=255)
                                                                           ttCCAGGATTTCCTGGCTGCTCGGCGGATTTCCCTCTTCACCAATCAACAAACTGATTTCACGCGGCATAAATTCTTTACCGGTTGGCAGTTTGGCATTa                     <  1:3273944/100‑1 (MQ=255)
                                                                             ccAGGATTTCCTGGCTGCTCGGCGGATTTCCCTCTTCACCAATCAACAAACTGATTTCACGCGGCATAAATTCTTTACCGGTTGGCAGTTTGGCATTAAg                   >  2:2906087/1‑100 (MQ=255)
                                                                             ccAGGATTTCCTGGCTGCTCGGCGGATTTCCCTCTTCACCAATCAACAAACTGATTTCACGCGGCATAAATTCTTTACCGGTTGGCAGTTTGGCATTAAg                   <  1:3359148/100‑1 (MQ=255)
                                                                              cAGGATTTCCTGGCTGCTCGGCGGATTTCCCTCTTCACCAATCAACAAACTGATTTCACGCGGCATAAATTCTTTACCGGTTGGCAGTTTGGCATTAAgg                  <  2:2891216/100‑1 (MQ=255)
                                                                               aGGATTTCCTGGCTGCTCGGCGGATTTCCCTCTTCACCAATCAACAAACTGATTTCACGCGGCATAAATTCTTTACCGGTTGGCAGTTTGGCATTAAGGa                 <  1:1646279/100‑1 (MQ=255)
                                                                               aGGATTTCCTGGCTGCTCGGCGGATTTCCCTCTGCACCAATCAACAAACTGATGTCACTCGGCATAAATTCTTTACCGGTTGGCAGTTTGGCATTAAGGa                 >  2:280591/1‑100 (MQ=255)
                                                                                 gATTTCCTGGCTGCTCAGCGGATTTCCCTCTTCACCAATCAACAAACTGATTTCACGCGGCATAAATTCTTTACCGGTTGGCAGTTTGGCATTAAGGAag               <  1:599846/100‑1 (MQ=255)
                                                                                    ttCCTGGCTGCTCGGCGGATTTCCCTCTTCACCAATCAACAAACTGATTTCACGCGGCATAAATTCTTTACCGGTTGGCAGTTTGGCATTAAGGAAGAGc            <  2:221838/100‑1 (MQ=255)
                                                                                      ccTGGCTGCTCGGCGGATTTCCCTCTTCACCAATCAACAAACTGATTTCACGCGGCATAAATTCTTTACCGGTTGGCAGTTTGGCATTAAGGAAGAGCGt          >  2:650243/1‑100 (MQ=255)
                                                                                       cTGGCTGCTCGGCGGATTTCCCTCTTCACCAATCAACAAACTGATTTCACGCGGCATAAATTCTTTACCGGTTGGCAGTTTGGCATTAAGGAAGAGCGtt         <  1:1991903/100‑1 (MQ=255)
                                                                                        tGGCTGCTCGGCGGATTTCCCTCTTCACCAATCAACAAACTGATTTCACGCGGCATAAATTCTTTACCGGTTGGCAGTTTGGCATTAAGGAAGAGCGttt        <  2:1411779/100‑1 (MQ=255)
                                                                                          gctgctCGGCGGATTTCCCTCTTCACCAATCAACAAACTGATTTCACGCGGCATAAATTCTTTACCGGTTGGCAGTTTGGCATTAAGGAAGAGCGTTTCg      <  1:210966/100‑1 (MQ=255)
                                                                                           ctgctCGGCGGATTTCCCTCTTCACCAATCAACAAACTGATTTCACGCGGCATAAATTCTTTACCGGTTGGCAGTTTGGCATTAAGGAAGAGCGTTTCga     >  2:1140817/1‑99 (MQ=255)
                                                                                           ctgctCGGCGGATTTCCCTCTTCACCAATCAACAAACTGATTTCACGCGGCATAAATTCTTTACCGGTTGGCAGTTTGGCATTAAGGAAGAGCGTTTCGc     >  1:3121157/1‑100 (MQ=255)
                                                                                            tgctCGGCGGATTTCCCTCTTCACCAATCAACAAACTGATTTCACGCGGCATAAATTCTTTACCGGTTGGCAGTTTGGCATTAAGGAAGAGCGTTTCGcc    <  2:1100556/100‑1 (MQ=255)
                                                                                             gctCGGCGGATTTCCCTCTTCACCAATCAACAAACTGATTTCACGCGGCATAAATTCTTTACCGGTTGGCAGTTTGGCATTAAGGAAGAGCGTTTCGCCa   >  1:2746835/1‑100 (MQ=255)
                                                                                                 |                                                                                                
AAATAAGGTGGTGAGTGAATCAATCATTTGTGCTGGCGGCTCTGCGACTTCCGATAATATCAGCGATTCACCGCCTTCCAGGATTTCCTGGCTGCTCAGCGGATTTCCCTCTTCACCAATCAACAAACTGATTTCACGCGGCATAAATTCTTTACCGGTTGGCAGTTTGGCATTAAGGAAGAGCGTTTCGCCAA  >  NC_000913/2654448‑2654641

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: