Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,837,763 T→C H91H (CAT→CAC b1757 → putative thiosulfate sulfur transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,837,7630TC86.0% 608.5 / 84.3 229H91H (CAT→CACb1757putative thiosulfate sulfur transferase
Reads supporting (aligned to +/- strand):  ref base T (15/17);  new base C (99/98);  total (114/115)
Fisher's exact test for biased strand distribution p-value = 8.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.35e-01

Only 100 of 234 total aligned reads displayed.

CCTTCTGGTCATGAACTTGCCGCCTTAAACCTCTCTGCCAGCTGGCTTGACAAAATGAGCACCGAACAGCTCAACGCGTGGATCAAGCAACATAACCTGAAAACCGATGCTCCGGTGGCGCTGTACGGTAATGACAAAGATGTCGACGCCGTCAAAACGCGACTGCAAAAAGCAGGTTTAAC  >  NC_000913/1837671‑1837852
                                                                                            |                                                                                         
ccTTCTGGTCATGAACTTGCCGCCTTAAACCTCTCTGCCAGCTGGCTTGACAAAATGAGCACCGAACAGCTCAACGCGTGGATCAAGCAACAt                                                                                           >  2:3174084/1‑93 (MQ=255)
ccTTCTGGTCATGAACTTGCCGCCTTAAACCTCTCTGCCAGCTGGCTTGACAAAATGAGCACCGAACAGCTCAACGCGTGGATCAAGCAACAt                                                                                           >  2:387039/1‑93 (MQ=255)
 cTTCTGGTCATGAACTTGCCGCCTTAAACCTCTCTGCCAGCTGGCTTGACAAAATGAGCACCGAACAGCTCAACGCGTGGATCAAGCAACATa                                                                                          >  1:4267860/1‑93 (MQ=255)
  ttCTGGTCATGAACTTGCCGCCTTAAACCTCTCTGCCAGCTGGCTTGACAAAATGAGCACCGAACAGCTCAACGCGTGGATCAAGCAACACaa                                                                                         >  1:3777049/1‑93 (MQ=255)
  ttCTGGTCATGAACTTGCCGCCTTAAACCTCTCTGCCAGCTGGCTTGACAAAATGAGCACCGAACAGCTCAACGCGTGGATCAAGCAACACaa                                                                                         <  1:3965463/93‑1 (MQ=255)
   tCTGGTCATGAACTTGCCGCCTTAAACCTCTCTGCCAGCTGGCTTGACAAAATGAGCACCGAACAGCTCAACGCGTGGATCAAGCAACATAAc                                                                                        >  2:2135451/1‑93 (MQ=255)
   tCTGGTCATGAACTTGCCGCCTTAAACCTCTCTGCCAGCTGGCTTGACAAAATGAGCACCGAACAGCTCAACGCGTGGATCAAGCAACACAAc                                                                                        <  2:616275/93‑1 (MQ=255)
     tGGTCATGAACTTGCCGCCTTAAACCTCTCTGCCAGCTGGCTTGACAAAATGAGCACCGAACAGCTCAACGCGTGGATCAAGCAACATAACCt                                                                                      >  2:17606/1‑93 (MQ=255)
     tGGTCATGAACTTGCCGCCTTAAACCTCTCTGCCAGCTGGCTTGACAAAATGAGCACCGAACAGCTCAACGCGTGGATCAAGCAACACAACCt                                                                                      >  1:3496456/1‑93 (MQ=255)
      ggTCATGAACTTGCCGCCTTAAACCTCTCTGCCAGCTGGCTTGACAAAATGAGCACCGAACAGCTCAACGCGTGGATCAAGCAACACAACCTg                                                                                     >  2:1731999/1‑93 (MQ=255)
        tCATGAACTTGCCGCCTTAACCCTCTCTGCCAGCTGGCTTGACAAAATGAGCACCGAACAGCTCAACGCGTGGATCAAGCAACACAACCTGaa                                                                                   <  1:2050921/93‑1 (MQ=255)
         cATGAACTTGCCGCCTTAAACCTCTCTGCCAGCTGGCTTGACAAAATGAGCACCGAACAGCTCAACGCGTGGATCAAGCAACACAACCTGaaa                                                                                  <  1:2399743/93‑1 (MQ=255)
           tGAACTTGCCGCCTTAAACCTCTCTGCCAGCTGGCTTGACAAAATGAGCACCGAACAGCTCAACGCGTGGATCAAGCAACACAACCTGAAAAc                                                                                <  2:5073816/93‑1 (MQ=255)
           tGAACTTGCCGCCTTAAACCTCTCTGCCAGCTGGCTTGACAAAATGAGCACCGAACAGCTCAACGCGTGGATCAAGCAACACAACCTGAAAAc                                                                                >  1:1897841/1‑93 (MQ=255)
           tGAACTTGCCGCCTTAAACCTCTCTGCCAGCTGGCTTGACAAAATGAGCACCGAACAGCTCAACGCGTGGATCAAGCAACACAACCTGAAAAc                                                                                >  2:6603848/1‑93 (MQ=255)
             aaCTTGCCGCCTTAAACCTCTCTGCCAGCTGGCTTGACAAAATGAGCACCGAACAGCTCAACGCGTGGATCAAGCAACACAACCTGAAAACCg                                                                              >  2:1909636/1‑93 (MQ=255)
             aaCTTGCCGCCTTAAACCTCTCTGCCAGCTGGCTTGACAAAATGAGCACCGAACAGCTCAACGCGTGGATCAAGCAACACAACCTGAAAACCg                                                                              >  2:2544097/1‑93 (MQ=255)
               cTTGCCGCCTTAAACCTCTCTGCCAGCTGGCTTGACAAAATGAGCACCGAACAGCTCAACGCGTGGATCAAGCAACACAACCTGAAAACCGAt                                                                            >  2:1227354/1‑93 (MQ=255)
                 tGCCGCCTTAAACCTCTCTGCCAGCTGGCTTGACAAAATGAGCACCGAACAGCTCAACGCGTGGATCAAGCAACACAACCTGAAAACCGATGc                                                                          >  1:3453937/1‑93 (MQ=255)
                 tGCCGCCTTAAACCTCTCTGCCAGCTGGCTTGACAAAATGAGCACCGAACAGCTCAACGCGTGGATCAAGCAACACAACCTGAAAACCGATGc                                                                          <  1:6503146/93‑1 (MQ=255)
                  gccgccTTAAACCTCTCTGCCAGCTGGCTTGACAAAATGAGCACCGAACAGCTCAACGCGTGGATCAAGCAACATAACCTGAAAACCGATGCt                                                                         <  1:3113558/93‑1 (MQ=255)
                  gccgccTTAAACCTCTCTGCCAGCTGGCTTGACAAAATGAGCACCGAACAGCTCAACGCGTGGATCAAGCAACATAACCTGAAAACCGATGCt                                                                         <  2:4307149/93‑1 (MQ=255)
                   ccgccTTAAACCTCTCTGCCAGCTGGCTTGACAAAATGAGCACCGAACAGCTCAACGCGTGGATCAAGCAACATAACCTGAAAACCGATGCTc                                                                        <  2:6241367/93‑1 (MQ=255)
                   ccgccTTAAACCTCTCTGCCAGCTGGCTTGACAAAATGAGCACCGAACAGCTCAACGCGTGGATCAAGCAACACAACCTGAAAACCGATGCTc                                                                        <  1:3039892/93‑1 (MQ=255)
                    cgccTTAAACCTCTCTGCCAGCTGGCTTGACAAAATGAGCACCGAACAGCTCAACGCGTGGATCAAGCAACATAACCTGAAAACCGATGCTcc                                                                       <  2:3945149/93‑1 (MQ=255)
                     gccTTAAACCTCTCTGCCAGCTGGCTTGACAAAATGAGCACCGAACAGCTCAACGCGTGGATCAAGCAACACAACCTGAAAACCGATGCTCCg                                                                      >  1:4799876/1‑93 (MQ=255)
                      ccTTAAACCTCTCTGCCAGCTGGCTTGACAAAATGAGCACCGAACAGCTCAACGCGTGGATCAAGCAACATAACCTGAAAACCGATGCTCCgg                                                                     <  1:211689/93‑1 (MQ=255)
                      ccTTAAACCTCTCTGCCAGCTGGCTTGACAAAATGAGCACCGAACAGCTCAACGCGTGGATCAAGCAACACAACCTGAAAACCGATGCTCCgg                                                                     >  2:203140/1‑93 (MQ=255)
                      ccTTAAACCTCTCTGCCAGCTGGCTTGACAAAATGAGCACCGAACAGCTCAACGCGTGGATCAAGCAACACAACCTGAAAACCGATGCTCCgg                                                                     <  1:4641517/93‑1 (MQ=255)
                       cTTAAACCTCTCTGCCAGCTGGCTTGACAAAATGAGCACCGAACAGCTCAACGCGTGGATCAAGCAACATAACCTGAAAACCGATGCTCCGGt                                                                    >  1:5092702/1‑93 (MQ=255)
                        ttAAACCTCTCTGCCAGCTGGCTTGACAAAATGAGCACCGAACAGCTCAACGCGTGGATCAAGCAACACAACCTGAAAACCGATGCTCCGGTg                                                                   <  2:1046456/93‑1 (MQ=255)
                         tAAACCTCTCTGCCAGCTGGCTTGACAAAATGAGCACCGAACAGCTCAACGCGTGGATCAAGCAACACAACCTGAAAACCGATGCTCCGGTgg                                                                  >  1:5829600/1‑93 (MQ=255)
                          aaaCCTCTCTGCCAGCTGGCTTGACAAAATGAGCACCGAACAGCTCAACGCGTGGATCAAGCAACACAACCTGAAAACCGATGCTCCGGTGgc                                                                 >  1:1341448/1‑93 (MQ=255)
                           aaCCTCTCTGCCAGCTGGCTTGACAAAATGAGCACCGAACAGCTCAACGCGTGGATCAAGCAACACAACCTGAAAACCGATGCTCCGGTGgcg                                                                >  1:4024885/1‑93 (MQ=255)
                             ccTCTCTGCCAGCTGGCTTGACAAAATGAGCACCGAACAGCTCAACGCGTGGATCAAGCAACACAACCTGAAAACCGATGCTCCGGTGGCGCt                                                              <  1:4601911/93‑1 (MQ=255)
                              ctctctGCCAGCTGGCTTGACAAAATGAGCACCGAACAGCTCAACGCGTGGATCAAGCAACACAACCTGAAAACCGATGCTCCGGTGGCGCTg                                                             >  2:2840708/1‑93 (MQ=255)
                                ctctGCCAGCTGGCTTGACAAAATGAGCACCGAACAGCTCAACGCGTGGATCAAGCAACACAACCTGAAAACCGATGCTCCGGTGGCGCTGTa                                                           >  1:1749558/1‑93 (MQ=255)
                                ctctGCCAGCTGGCTTGACAAAATGAGCACCGAACAGCTCAACGCGTGGATCAAGCAACACAACCTGAAAACCGATGCTCCGGTGGCGCTGTa                                                           >  1:1043315/1‑93 (MQ=255)
                                 tctGCCAGCTGGCTTGACAAAATGAGCACCGAACAGCTCAACGCGTGGATCAAGCAACACAACCTGAAAACCGATGCTCCGGTGGCGCTGTAc                                                          <  2:1841544/93‑1 (MQ=255)
                                 tctGCCAGCTGGCTTGACAAAATGAGCACCGAACAGCTCAACGCGTGGATCAAGCAACACAACCTGAAAACCGATGCTCCGGTGGCGCTGTAc                                                          <  2:1071520/93‑1 (MQ=255)
                                   tGCCAGCTGGCTTGACAAAATGAGCACCGAACAGCTCAACGCGTGGATCAAGCAACACAACCTGAAAACCGATGCTCCGGTGGCGCTGTACgg                                                        <  2:1949110/93‑1 (MQ=255)
                                    gCCAGCTGGCTTGACAAAATGAGCACCGAACAGCTCAACGCGTGGATCAAGCAACATAACCTGAAAACCGATGCTCCGGTGGCGCTGTACGGt                                                       <  2:1445776/93‑1 (MQ=255)
                                    gCCAGCTGGCTTGACAAAATGAGCACCGAACAGCTCAACGCGTGGATCAAGCAACACAACCTGAAAACCGATGCTCCGGTGGCGCTGTACgga                                                       >  1:2265305/1‑92 (MQ=255)
                                    gCCAGCTGGCTTGACAAAATGAGCACCGAACAGCTCAACGCGTGGATCAAGCAACACAACCTGAAAACCGATGCTCCGGTGGCGCTGTACGGt                                                       <  1:1056797/93‑1 (MQ=255)
                                    gCCAGCTGGCTTGACAAAATGAGCACCGAACAGCTCAACGCGTGGATCAAGCAACACAACCTGAAAACCGATGCTCCGGTGGCGCTGTACGGt                                                       >  1:1626487/1‑93 (MQ=255)
                                      cAGCTGGCTTGACAAAATGAGCACCGAACAGCTCAACGCGTGGATCAAGCAACATAACCTGAAAACCGATGCTCCGGTGGCGCTGTACGGTaa                                                     <  2:6045536/93‑1 (MQ=255)
                                      cAGCTGGCTTGACAAAATGAGCACCGAACAGCTCAACGCGTGGATCAAGCAACACAACCTGAAAACCGATGCTCCGGTGGCGCTGTACGGTaa                                                     <  1:776061/93‑1 (MQ=255)
                                       aGCTGGCTTGACAAAATGAGCACCGAACAGCTCAACGCGTGGATCAAGCAACATAACCTGAAAACCGATGCTCCGGTGGCGCTGTACGGTAAt                                                    <  1:2033943/93‑1 (MQ=255)
                                       aGCTGGCTTGACAAAATGAGCACCGAACAGCTCAACGCGTGGATCAAGCAACACAACCTGAAAACCGATGCTCCGGTGGCGCTGTACGGTAAt                                                    >  2:1914821/1‑93 (MQ=255)
                                       aGCTGGCTTGACAAAATGAGCACCGAACAGCTCAACGCGTGGATCAAGCAACACAACCTGAAAACCGATGCTCCGGTGGCGCTGTACGGTAAt                                                    >  2:3591290/1‑93 (MQ=255)
                                        gCTGGCTTGACAAAATGAGCACCGAACAGCTCAACGCGTGGATCAAGCAACACAACCTGAAAACCGATGCTCCGGTGGCGCTGTACGGTAATg                                                   <  2:6313246/93‑1 (MQ=255)
                                         cTGGCTTGACAAAATGAGCACCGAACAGCTCAACGCGTGGATCAAGCAACACAACCTGAAAACCGATGCTCCGGTGGCGCTGTACGGTAATGa                                                  <  1:4713145/93‑1 (MQ=255)
                                           ggCTTGACAAAATGAGCACCGAACAGCTCAACGCGTGGATCAAGCAACACAACCTGAAAACCGATGCTCCGGTGGCGCTGTACGGTAATGACa                                                >  2:5196753/1‑93 (MQ=255)
                                           ggCTTGACAAAATGAGCACCGAACAGCTCAACGCGTGGATCAAGCAACACAACCTGAAAACCGATGCTCCGGTGGCGCTGTACGGTAATGACa                                                <  2:3681729/93‑1 (MQ=255)
                                            gCTTGACAAAATGAGCACCGAACAGCTCAACGCGTGGATCAAGCAACACAACCTGAAAACCGATGCTCCGGTGGCGCTGTACGGTAATGACaa                                               <  1:5595836/93‑1 (MQ=255)
                                            gCTTGACAAAATGAGCACCGAACAGCTCAACGCGTGGATCAAGCAACACAACCTGAAAACCGATGCTCCGGTGGCGCTGTACGGTAATGACaa                                               <  1:4506324/93‑1 (MQ=255)
                                            gCTTGACAAAATGAGCACCGAACAGCTCAACGCGTGGATCAAGCAACACAACCTGAAAACCGATGCTCCGGTGGCGCTGTACGGTAATGACaa                                               <  2:6121289/93‑1 (MQ=255)
                                             cTTGACAAAATGAGCACCGAACAGCTCAACGCGTGGATCAAGCAACACAACCTGAAAACCGATGCTCCGGTGGCGCTGTACGGTAATGACaaa                                              >  1:5016405/1‑93 (MQ=255)
                                             cTTGACAAAATGAGCACCGAACAGCTCAACGCGTGGATCAAGCAACACAACCTGAAAACCGATGCTCCGGTGGCGCTGTACGGTAATGACaaa                                              <  2:4733569/93‑1 (MQ=255)
                                              ttGACAAAATGAGCACCGAACAGCTCAACGCGTGGATCAAGCAACACAACCTGAAAACCGATGCTCCGGTGGCGCTGTACGGTAATGACAAAg                                             <  1:423221/93‑1 (MQ=255)
                                               tGACAAAATGAGCACCGAACAGCTCAACGCGTGGATCAAGCAACACAACCTGAAAACCGATGCTCCGGTGGCGCTGTACGGTAATGACAAAGa                                            <  1:549810/93‑1 (MQ=255)
                                               tGACAAAATGAGCACCGAACAGCTCAACGCGTGGATCAAGCAACACAACCTGAAAACCGATGCTCCGGTGGCGCTGTACGGTAATGACAAAGa                                            <  2:5682640/93‑1 (MQ=255)
                                               tGACAAAATGAGCACCGAACAGCTCAACGCGTGGATCAAGCAACACAACCTGAAAACCGATGCTCCGGTGGCGCTGTACGGTAATGACAAAGa                                            <  1:1930674/93‑1 (MQ=255)
                                                gACAAAATGAGCACCGAACAGCTCAACGCGTGGATCAAGCAACACAACCTGAAAACCGATGCTCCGGTGGCGCTGTACGGTAATGACAAAGAt                                           <  2:5728759/93‑1 (MQ=255)
                                                gACAAAATGAGCACCGAACAGCTCAACGCGTGGATCAAGCAACACAACCTGAAAACCGATGCTCCGGTGGCGCTGTACGGTAATGACAAAGAt                                           >  2:1330025/1‑93 (MQ=255)
                                                  cAAAATGAGCACCGAACAGCTCAACGCGTGGATCAAGCAACATAACCTGAAAACCGATGCTCCGGTGGCGCTGTACGGTAATGACAAAGATGt                                         <  1:1113405/93‑1 (MQ=255)
                                                   aaaaTGAGCACCGAACAGCTCAACGCGTGGATCAAGCAACACAACCTGAAAACCGATGCTCCGGTGGCGCTGTACGGTAATGACAAAGATGTc                                        <  1:6037729/93‑1 (MQ=255)
                                                    aaaTGAGCACCGAACAGCTCAACGCGTGGATCAAGCAACACAACCTGAAAACCGATGCTCCGGTGGCGCTGTACGGTAATGACAAAGATGTCg                                       >  1:5600260/1‑93 (MQ=255)
                                                    aaaTGAGCACCGAACAGCTCAACGCGTGGATCAAGCAACACAACCTGAAAACCGATGCTCCGGTGGCGCTGTACGGTAATGACAAAGATGTCg                                       <  2:6580848/93‑1 (MQ=255)
                                                     aaTGAGCACCGAACAGCTCAACGCGTGGATCAAGCAACACAACCTGAAAACCGATGCTCCGGTGGCGCTGTACGGTAATGACAAAGATGTCGa                                      <  2:1174996/93‑1 (MQ=255)
                                                            aCCGAACAGCTCAACGCGTGGATCAAGCAACACAACCTGAAAACCGATGCTCCGGTGGCGCTGTACGGTAATGACAAAGATGTCGACGCCGTc                               >  1:92340/1‑93 (MQ=255)
                                                               gAACAGCTCAACGCGTGGATCAAGCAACATAACCTGAAAACCGATGCTCCGGTGGCGCTGTACGGTAATGACAAAGATGTCGACGCCGTCaaa                            >  1:2768581/1‑93 (MQ=255)
                                                               gAACAGCTCAACGCGTGGATCAAGCAACACAACCTGAAAACCGATGCTCCGGTGGCGCTGTACGGTAATGACAAAGATGTCGACGCCGTCaaa                            >  1:5204417/1‑93 (MQ=255)
                                                                aaCAGCTCAACGCGTGGATCAAGCAACACAACCTGAAAACCGATGCTCCGGTGGCGCTGTACGGTAATGACAAAGATGTCGACGCCGTCaaaa                           >  1:444648/1‑93 (MQ=255)
                                                                aaCAGCTCAACGCGTGGATCAAGCAACACAACCTGAAAACCGATGCTCCGGTGGCGCTGTACGGTAATGACAAAGATGTCGACGCCGTCaaaa                           >  2:6278657/1‑93 (MQ=255)
                                                                aaCAGCTCAACGCGTGGATCAAGCAACACAACCTGAAAACCGATGCTCCGGTGGCGCTGTACGGTAATGACAAAGATGTCGACGCCGTCaaaa                           >  2:4509183/1‑93 (MQ=255)
                                                                  cAGCTCAACCCGTGGATCAAGCAACACAACCTGAAAACCGATGCTCCGGTGGCGCTGTACGGTAATGACAAAGATGTCGACGCCGTCAAAAcg                         >  1:5302358/1‑93 (MQ=255)
                                                                      tCAACGCTTGGATCAAGCAACACAACCTGAAAACCGATGCTCCGGTGGCGCTGTACGGTAATGACAAAGATGTCGACGCCGTCAAAACGCGAc                     <  1:2572587/93‑1 (MQ=255)
                                                                        aaCGCGTGGATCAAGCAACACAACCTGAAAACCGATGCTCCGGTGGCGCTGTACGGTAATGACAAAGATGTCGACGCCGTCAAAACGCGACTg                   >  1:2502822/1‑93 (MQ=255)
                                                                        aaCGCGTGGATCAAGCAACACAACCTGAAAACCGATGCTCCGGTGGCGCTGTACGGTAATGACAAAGATGTCGACGCCGTCAAAACGCGACTg                   >  2:590459/1‑93 (MQ=255)
                                                                         aCGCGTGGATCAAGCAACACAACCTGAAAACCGATGCTCCGGTGGCGCTGTACGGTAATGACAAAGATGTCGACGCCGTCAAAACGCGACTGc                  >  1:617459/1‑93 (MQ=255)
                                                                          cgcgTGGATCAAGCAACATAACCTGAAAACCGATGCTCCGGTGGCGCTGTACGGTAATGACAAAGATGTCGACGCCGTCAAAACGCGACTGCa                 >  1:1565560/1‑93 (MQ=255)
                                                                          cgcgTGGATCAAGCAACACAACCTGAAAACCGATGCTCCGGTGGCGCTGTACGGTAATGACAAAGATGTCGACGCCGTCAAAACGCGACTGCa                 >  2:968176/1‑93 (MQ=255)
                                                                          cgcgTGGATCAAGCAACACAACCTGAAAACCGATGCTCCGGTGGCGCTGTACGGTAATGACAAAGATGTCGACGCCGTCAAAACGCGACTGCa                 <  1:1888286/93‑1 (MQ=255)
                                                                             gTGGATCAAGCAACACAACCTGAAAACCGATGCTCCGGTGGCGCTGTACGGTAATGACAAAGATGTCGACGCCGTCAAAACGCGACTGCaaaa              >  1:4158948/1‑93 (MQ=255)
                                                                              tGGATCAAGCAACACAACCTGAAAACCGATGCTCCGGTGGCGCTGTACGGTAATGACAAAGATGTCGACGCCGTCAAAACGCGACTGCaaaaa             >  2:5286731/1‑93 (MQ=255)
                                                                              tGGATCAAGCAACACAACCTGAAAACCGATGCTCCGGTGGCGCTGTACGGTAATGACAAAGATGTCGACGCCGTCAAAACGCGACTGCaaaaa             <  2:3207482/93‑1 (MQ=255)
                                                                              tGGATCAAGCAACACAACCTGAAAACCGATGCTCCGGTGGCGCTGTACGGTAATGACAAAGATGTCGACGCCGTCAAAACGCGACTGCaaaaa             <  2:6248115/93‑1 (MQ=255)
                                                                               ggATCAAGCAACATAACCTGAAAACCGATGCTCCGGTGGCGCTGTACGGTAATGACAAAGATGTCGACGCCGTCAAAACGCGACTGCAAAAAg            >  1:2301443/1‑93 (MQ=255)
                                                                               ggATCAAGCAACACAACCTGAAAACCGATGCTCCGGTGGCGCTGTACGGTAATGACAAAGATGTCGACGCCGTCAAAACGCGACTGCAAAAAg            <  1:5686335/93‑1 (MQ=255)
                                                                               ggATCAAGCAACACAACCTGAAAACCGATGCTCCGGTGGCGCTGTACGGTAATGACAAAGATGTCGACGCCGTCAAAACGCGACTGCAAAAAg            <  2:3292796/93‑1 (MQ=255)
                                                                               ggATCAAGCAACACAACCTGAAAACCGATGCTCCGGTGGCGCTGTACGGTAATGACAAAGATGTCGACGCCGTCAAAACGCGACTGCAAAAAg            <  1:6244765/93‑1 (MQ=255)
                                                                                 aTCAAGCAACACAACCTGAAAACCGATGCTCCGGTGGCGCTGTACGGTAATGACAAAGATGTCGACGCCGTCAAAACGCGACTGCAAAAAGCa          >  2:715820/1‑93 (MQ=255)
                                                                                 aTCAAGCAACACAACCTGAAAACCGATGCTCCGGTGGCGCTGTACGGTAATGACAAAGATGTCGACGCCGTCAAAACGCGACTGCAAAAAGCa          <  2:2497675/93‑1 (MQ=255)
                                                                                   cAAGCAACACAACCTGAAAACCGATGCTCCGGTGGCGCTGTACGGTAATGACAAAGATGTCGACGCCGTCAAAACGCGACTGCAAAAAGCAgg        >  1:3210510/1‑93 (MQ=255)
                                                                                   cAAGCAACACAACCTGAAAACCGATGCTCCGGTGGCGCTGTACGGTAATGACAAAGATGTCGACGCCGTCAAAACGCGACTGCAAAAAGCAgg        <  2:5265583/93‑1 (MQ=255)
                                                                                   cAAGCAACACAACCTGAAAACCGATGCTCCGGTGGCGCTGTACGGTAATGACAAAGATGTCGACGCCGTAAAAACGCGACTGCAAAAAGCAgg        <  1:5350865/93‑1 (MQ=255)
                                                                                     aGCAACACAACCTGAAAACCGATGCTCCGGTGGCGCTGTACGGTAATGACAAAGATGTCGACGCCGTCAAAACGCGACTGCAAAAAGCAGGtt      >  1:8018/1‑93 (MQ=255)
                                                                                      gCAACATAACCTGAAAACCGATGCTCCGGTGGCGCTGTACGGTAATGACAAAGATGTCGACGCCGTCAAAACGCGACTGCAAAAAGCAGGttt     <  1:5995315/93‑1 (MQ=255)
                                                                                         aCATAACCTGAAAACCGATGCTCCGGTGGCGCTGTACGGTAATGACAAAGATGTCGACGCCGTCAAAACGCGACTGCAAAAAGCAGGTTTAAc  <  2:4203816/93‑1 (MQ=255)
                                                                                            |                                                                                         
CCTTCTGGTCATGAACTTGCCGCCTTAAACCTCTCTGCCAGCTGGCTTGACAAAATGAGCACCGAACAGCTCAACGCGTGGATCAAGCAACATAACCTGAAAACCGATGCTCCGGTGGCGCTGTACGGTAATGACAAAGATGTCGACGCCGTCAAAACGCGACTGCAAAAAGCAGGTTTAAC  >  NC_000913/1837671‑1837852

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