Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,755,517 T→C intergenic (‑63/+127) yiaY ← / ← selB putative oxidoreductase/selenocysteinyl‑tRNA‑specific translation factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,755,5170TC99.3% 461.1 / ‑4.8 142intergenic (‑63/+127)yiaY/selBputative oxidoreductase/selenocysteinyl‑tRNA‑specific translation factor
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (5/136);  total (5/137)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.55e-01

Only 100 of 142 total aligned reads displayed.

AAGGAATAAAGAACGTTGAAGCTGCCATAACACACTCCTGAAAGTGTAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTAAATTCATAAGAAT  >  NC_000913/3755427‑3755604
                                                                                          |                                                                                       
aaGGAATAAAGAACGTTGAAGCTGCCATAACACACTCCTGAAAGTGTAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATACGa                                                                                       <  2:3255817/93‑1 (MQ=255)
 aGGAATAAAGAACGTTGAAGCTGCCATAACACACTCCTGAAAGTGTAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATACGAt                                                                                      <  2:1150055/93‑1 (MQ=255)
 aGGAATAAAGAACGTTGAAGCTGCCATAACACACTCCTGAAAGTGTAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATACGAt                                                                                      <  2:3591832/93‑1 (MQ=255)
   gAATAAAGAACGTTGAAGCTGCCATAACACACTCCTGAAAGTGTAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATACGATCa                                                                                    <  1:236668/93‑1 (MQ=255)
      tAAAGAACGTTGAAGCTGCCATAACACACTCCTGAAAGTGTAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATACGATCAGta                                                                                 <  1:1006412/93‑1 (MQ=255)
      tAAAGAACGTTGAAGCTGCCATAACACACTCCTGAAAGTGTAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATACGATCAGta                                                                                 <  2:743388/93‑1 (MQ=255)
        aaGAACGTTGAAGCTGCCATAACACACTCCTGAAAGTGTAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATACGATCAGtatt                                                                               <  2:3956151/93‑1 (MQ=255)
         aGAACTTTGAACCTGCCATAACACACTCCTGAAAGTGTAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATACGATCAGtatta                                                                              <  1:2972644/93‑1 (MQ=255)
         aGAACGTTGAAGCTGCCATAACACACTCCTGAAAGTGTAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATACGATCAGtatta                                                                              <  1:3663086/93‑1 (MQ=255)
         aGAACGTTGAAGCTGCCATAACACACTCCTGAAAGTGTAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATACGATCAGtatta                                                                              <  2:2083159/93‑1 (MQ=255)
         aGAACGTTGAAGCTGCCATAACACACTCCTGAAAGTGTAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATACGATCAGtatta                                                                              <  1:342425/93‑1 (MQ=255)
         aGAACGTGGAAGCTGCCATAACACACTCCTGAAAGTGTAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATACGATCAGtatta                                                                              <  1:4709208/93‑1 (MQ=255)
          agaCGTGGAGGCTGCATTAACACACTCCTGAAAGTGTAAAAAAGGGCTGGCCGCTTACTGAGGCTTTTCATTCGAAAATACGATCAGtattat                                                                             <  1:3976885/91‑1 (MQ=255)
          aaaCGTTGAAGCTGCCATAACACACTCCTGAAAGTGTAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATACGATCAGtattat                                                                             <  1:3976884/92‑1 (MQ=255)
           accGTTGAAGCTGCCATAACACACTCCTGAAAGTGTAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATACGATCAGTattatt                                                                            <  1:186122/91‑1 (MQ=255)
           aaCGTTGAAGCTGCCATAACACACTCCTGAAAGTGTAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATACGATCAGTattatt                                                                            <  2:619390/93‑1 (MQ=255)
           aaCGTTGAAGCTGCCATAACACACTCCTGAAAGTGTAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATACGATCAGTattatt                                                                            <  2:353585/93‑1 (MQ=255)
            aCGTTGAAGCTGCCATAACACACTCCTGAAAGTGTAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATACGATCAGTATTATTg                                                                           >  2:4630521/1‑93 (MQ=255)
             cGTTGAAGCTGCCATAACACACTCCTGAAAGTGTAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATACGATCAGTATTATTGa                                                                          <  2:1650693/93‑1 (MQ=255)
              gTTGAAGCTGCCATAACACACTCCTGAAAGTGTAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATACGATCAGTATTATTGAt                                                                         <  2:4572420/93‑1 (MQ=255)
               ttGAAGCTGCCATAACACACTCCTGAAAGTGTAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATACGATCAGTATTATTGAtt                                                                        <  2:1780071/93‑1 (MQ=255)
                tGAAGCTGCCATAACACACTCCTGAAAGTGTAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATACGATCAGTATTATTGATTa                                                                       <  2:2987012/93‑1 (MQ=255)
                tGAAGCTGCCATAACACACTCCTGAAAGTGTAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATACGATCAGTATTATTGATTa                                                                       <  2:4232340/93‑1 (MQ=255)
                  aaGCTGCCATAACACACTCCTGAAAGTGTAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATACGATCAGTATTATTGATTATc                                                                     <  2:2542547/93‑1 (MQ=255)
                  aaGCTGCCATAACACACTCCTGAAAGTGTAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATACGATCAGTATTATTGATTATc                                                                     <  2:3516336/93‑1 (MQ=255)
                   aGCTGCCATAACACACTCCTGAAAGTGTAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATACGATCAGTATTATTGATTATCa                                                                    <  2:2476007/93‑1 (MQ=255)
                    gCTGCCATAACACACTCCTGAAAGTGTAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATACGATCAGTATTATTGATTATCaa                                                                   <  1:3096237/93‑1 (MQ=255)
                    gCTGCCATAACACACTCCTGAAAGTGTAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATACGATCAGTATTATTGATTATCaa                                                                   <  1:3096238/93‑1 (MQ=255)
                    gCTGCCATAACACACTCCTGAAAGTGTAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATACGATCAGTATTATTGATTATCaa                                                                   <  2:1040857/93‑1 (MQ=255)
                      tGCCATAACACACGCCTGAAAGTGTAAAAAAGGGCGTGCAGCTTACTGAGGATTTTCATTCGAAAATACGTACAGTATTATTGATTATCaaaa                                                                 <  2:3987621/93‑1 (MQ=255)
                       gCCATAACACACTCCTGAAAGTGTAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATACGATCAGTATTATTGATTATCAAAAt                                                                <  2:1374120/93‑1 (MQ=255)
                       gCCATAACACACTCCTGAAAGTGTAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATACGATCAGTATTATTGATTATCAAAAt                                                                <  1:1379936/93‑1 (MQ=255)
                       gCCATAACACACTCCTGAAAGTGTAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATACGATCAGTATTATTGATTATCAAAAt                                                                <  2:2315225/93‑1 (MQ=255)
                        ccATAACACACTCCTGAAAGTGTAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATACGATCAGTATTATTGATTATCAAAAtt                                                               <  1:3038620/93‑1 (MQ=255)
                         cATAACACACTCCTGAAAGTGTAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATACGATCAGTATTATTGATTATCAAAATTa                                                              <  2:1373106/93‑1 (MQ=255)
                          aTAACACACTCCTGAAAGTGTAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATACGATCAGTATTATTGATTATCAAAATTaa                                                             <  2:675909/93‑1 (MQ=255)
                          aTAACACACTCCTGAAAGTGTAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATACGATCAGTATTATTGATTATCAAAATTaa                                                             <  2:4670230/93‑1 (MQ=255)
                           tAACACACTCCTGAAAGTGTAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATACGATCAGTATTATTGATTATCAAAATTAAt                                                            <  2:4647650/93‑1 (MQ=255)
                            aaCACACTCCTGAAAGTGTAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATACGATCAGTATTATTGATTATCAAAATTAATc                                                           <  2:1792333/93‑1 (MQ=255)
                              cactcTACTGCAAGTGTAAGAACGGGCGGGCAGCATACAGAAAATTTTCAATCGAAAATACGATCAGTATTATTGATTATCAGAATTAATCta                                                         <  1:4223555/93‑1 (MQ=255)
                              cacacTCCTGAAAGTGTAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATACGATCAGTATTATTGATTATCAAAATTAATCta                                                         <  2:4685908/93‑1 (MQ=255)
                              cacacTCCTGAAAGTGTAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATACGATCAGTATTATTGATTATCAAAATTAATCta                                                         <  1:1316088/93‑1 (MQ=255)
                               acacTCCTGAAAGTGGAAAAAGGGGCGGGCATCTTACTGAGGATTTTCATTCGAAAATACGATCAGTATTATTGATTATCAAAATTAATCtaa                                                        <  1:2978114/93‑1 (MQ=255)
                                cacTCCTGAAAGTGTAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATACGATCAGTATTATTGATTATCAAAATTAATCtaat                                                       <  2:4593394/93‑1 (MQ=255)
                                     cTGAAAGTGTAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATACGATCAGTATTATTGATTATCAAAATTAATCTAATaaaaa                                                  <  2:2088853/93‑1 (MQ=255)
                                       gAAAGTGTAAAAAAGGGCGGGCAGCTTGCTGAGGATTTTCATTCGAAAATACGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGc                                                <  1:23602/93‑1 (MQ=255)
                                        aaaGTGTAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATACGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGca                                               <  2:3374426/93‑1 (MQ=255)
                                        aaaGTGTAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATACGATCAGTATTATTGATATCAAAATTAATCTAATAAAAAGcac                                              <  1:2286987/93‑1 (MQ=255)
                                         aaGTGTAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATACGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGcac                                              <  2:1696130/93‑1 (MQ=255)
                                         aaGTGTAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATACGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGcac                                              <  2:2861266/93‑1 (MQ=255)
                                         aaGTGTAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATACGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGcac                                              <  1:550559/93‑1 (MQ=255)
                                          aGTGTAAAAAAGGGCGAGCAGCTTACTGAGGATTTTCATTCGAAAATACGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGcaca                                             <  2:3173475/93‑1 (MQ=255)
                                           gtgtAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATACGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACAt                                            >  1:2294920/1‑93 (MQ=255)
                                           gtgtAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATACGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACAt                                            <  2:1412878/93‑1 (MQ=255)
                                           gtgtAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATACGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACAt                                            <  2:2306968/93‑1 (MQ=255)
                                            tgtAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATACGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACAtt                                           <  1:3498971/93‑1 (MQ=255)
                                               aaaaaaGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATACGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTa                                        <  1:3372636/93‑1 (MQ=255)
                                               aaaaaaGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATACGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTa                                        <  1:160044/93‑1 (MQ=255)
                                               aaaaaaGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATACGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTa                                        <  1:77122/93‑1 (MQ=255)
                                                aaaaaGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATACGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTAt                                       <  2:578761/93‑1 (MQ=255)
                                                 aaaaGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATACGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTAtt                                      <  1:4370570/93‑1 (MQ=255)
                                                 aaaaGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATACGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTAtt                                      <  2:2853926/93‑1 (MQ=255)
                                                 aaaaGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATACGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTAtt                                      <  1:2041542/93‑1 (MQ=255)
                                                   aaGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATACGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCa                                    <  1:2120790/93‑1 (MQ=255)
                                                    aGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATACGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCAt                                   >  1:2589783/1‑93 (MQ=255)
                                                         gggcAGCTTACTGAGGATTTTCATTCGAAAATACGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAAc                              <  1:1082192/93‑1 (MQ=255)
                                                           gcAGCTTACTGAGGATTTTCATTCGAAAATACGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACaa                            <  1:3529133/93‑1 (MQ=255)
                                                            cAGCTTACTGAGGATTTTCATTCGAAAATACGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAAt                           <  1:2654507/93‑1 (MQ=255)
                                                            cAGCTTACTGAGGATTTTCATTCGAAAATACGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAAt                           <  1:2453917/93‑1 (MQ=255)
                                                             aGCTTACTGAGGATTTTCATTCGAAAATACGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAAtg                          <  1:1473451/93‑1 (MQ=255)
                                                             aGCTTACTGAGGATTTTCATTCGAAAATACGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAAtg                          <  2:3049690/93‑1 (MQ=255)
                                                              gcgtACTGAGGATTTTCATTCGAAAATACGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAAtgt                         <  1:1510137/90‑1 (MQ=255)
                                                                ttACTGAGGATTTTCATTCGAAAATACGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGc                       <  2:3439106/93‑1 (MQ=255)
                                                                 tACTGAGGATTTTCATTCGAAAATACGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCt                      <  2:2878769/93‑1 (MQ=255)
                                                                   cTGAGGATTTTCATTCGAAAATACGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCttt                    <  1:619381/93‑1 (MQ=255)
                                                                   cTGAGGATTTTCATTCGAAAATACGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCttt                    <  2:4265906/93‑1 (MQ=255)
                                                                   cTGAGGATTTTCATTCGAAAATACGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCttt                    <  2:1758394/93‑1 (MQ=255)
                                                                   cTGAGGATTTTCATTCGAAAATACGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCttt                    <  1:220004/93‑1 (MQ=255)
                                                                    tGAGGATTTTCATTCGAAAATACGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCtttt                   <  1:4590454/93‑1 (MQ=255)
                                                                     gAGGATTTTCATTCGAAAATACGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCttttt                  <  1:536311/93‑1 (MQ=255)
                                                                     gAGGATTTTCATTCGAAAATACGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCttttt                  <  2:4627141/93‑1 (MQ=255)
                                                                      aGGATTTTCATTCGAAAATACGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCtttttt                 <  1:1286643/93‑1 (MQ=255)
                                                                       ggATTTTCATTCGAAAATACGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCttttttt                <  2:141360/93‑1 (MQ=255)
                                                                           tttCATTCGAAAATACGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTAAAt            <  2:4392849/93‑1 (MQ=255)
                                                                            ttCATTCGAAAATACGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTAAAtt           <  2:410860/93‑1 (MQ=255)
                                                                             tCATTCGAAAATACGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTAAATTc          <  2:3762319/93‑1 (MQ=255)
                                                                               aTTCGAAAATACGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTAAATTCAt        <  2:1218647/93‑1 (MQ=255)
                                                                                ttCGAAAATACGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTAAATTCATa       <  1:3235829/93‑1 (MQ=255)
                                                                                ttCGAAAATACGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTAAATTCATa       <  1:3235828/93‑1 (MQ=255)
                                                                                 tCGAAAATACGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTAAATTCATaa      <  1:4736963/93‑1 (MQ=255)
                                                                                  cGAAAATACGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTAAATTCATAAg     <  1:4302946/93‑1 (MQ=255)
                                                                                  cGAAAATACGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTAAATTCATAAg     <  2:2132070/93‑1 (MQ=255)
                                                                                   gAAAATACGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTTCTTTTTTTAAATTCATAAGa    <  1:1864566/93‑1 (MQ=255)
                                                                                   gAAAATACGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTAAATTCATAAGa    <  2:1941722/93‑1 (MQ=255)
                                                                                    aaaaTACGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTAAATTCATAAGaa   <  2:1212228/93‑1 (MQ=255)
                                                                                    aaaaTACGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTAAATTCATAAGaa   >  1:1655875/1‑93 (MQ=255)
                                                                                    aaaaTACGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTAAATTCATAAGaa   <  2:2131345/93‑1 (MQ=255)
                                                                                    aaaaTACGAGCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTAAATTCATAAGaa   <  1:1160699/93‑1 (MQ=255)
                                                                                     aaaTACGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTAAATTCATAAGAAt  <  1:4630521/93‑1 (MQ=255)
                                                                                     aaaTACGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTAAATTCATAAGAAt  <  2:2444406/93‑1 (MQ=255)
                                                                                          |                                                                                       
AAGGAATAAAGAACGTTGAAGCTGCCATAACACACTCCTGAAAGTGTAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTAAATTCATAAGAAT  >  NC_000913/3755427‑3755604

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