Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,004,812 T→C R448R (CGT→CGC recQ → ATP‑dependent DNA helicase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,004,8120TC99.6% 752.8 / ‑4.6 234R448R (CGT→CGCrecQATP‑dependent DNA helicase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (116/117);  minor base A (0/1);  total (116/118)
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.52e-01

Only 100 of 242 total aligned reads displayed.

TCAGATTGCCCTTTCCACCATTGGTCGTGTGAATCAGCGGTTTGGGATGGGTTATGTGGTGGAAGTGATTCGTGGTGCTAATAACCAGCGTATCCGCGACTATGGTCATGACAAACTGAAAGTCTATGGCATGGGCCGTGATAAAAGCCATGAACATTGGGTGAGCGTGATCCGCCAGCT  >  NC_000913/4004722‑4004901
                                                                                          |                                                                                         
tCAGATTGCCCTTTCCACCATTGGTCGTGTGAATCAGCGGTTTGGGATGGGTTATGTGGTGGAAGTGATTCGTGGTGCTAATAACCAGCGCAt                                                                                         >  2:2528971/1‑93 (MQ=255)
tCAGATTGCCCTTTCCACCATTGGTCGTGTGAATCAGCGGTTTGGGATGGGTTATGTGGTGGAAGTGATTCGCGGTGCTAATAACCAGCGCAt                                                                                         >  2:859023/1‑93 (MQ=255)
tCAGATTGCCCTTTCCACCATTGGTCGTGTGAATCAGCGGGTTGGGATGGGTTATGTGGTGGAAGTGATTCGTGGTGCTAATAACCAGCGCAt                                                                                         >  1:1825213/1‑93 (MQ=255)
 cAGATTGCCCTTTCCACCATTGGTCGTGTGAATCAGCGGTTTGGGATGGGTTATGTGGTGGAAGTGATTCGTGGTGCTAATAACCAGCGCATc                                                                                        <  1:74630/93‑1 (MQ=255)
 cAGATTGCCCTTTCCACCATTGGTCGTGTGAATCAGCGGTTTGGGATGGGTTATGTGGTGGAAGTGATTCGTGGTGCTAATAACCAGCGCATc                                                                                        <  2:4557533/93‑1 (MQ=255)
   gATTGCCCTTTCCACCATTGGTCGTGTGAATCAGCGGTTTGGGATGGGTTATGTGGTGGAAGTGATTCGTGGTGCTAATAACCAGCGCATCcg                                                                                      >  2:2259380/1‑93 (MQ=255)
     ttGCCCTTTCCACCATTGGTCGTGTGAATCAGCGGTTTGGGATGGGTTATGTGGTGGAAGTGATTCGTGGTGCTAATAACCAGCGCATCcgcg                                                                                    >  1:1134889/1‑93 (MQ=255)
     ttGCCCTTTCCACCATTGGTCGTGTGAATCAGCGGTTTGGGATGGGTTATGTGGTGGAAGTGATTCGTGGTGCTAATAACCAGCGCATCcgcg                                                                                    >  1:5678084/1‑93 (MQ=255)
     ttGCCCTTTCCACCATTGGTCGTGTGAATCAGCGGTTTGGGATGGGTTATGTGGTGGAAGTGATTCGTGGTGCTAATAACCAGCACATCcgcg                                                                                    <  2:6155956/93‑1 (MQ=255)
      tGCCCTTTCCACCATTGGTCGTGTGAATCAGCGGTTTGGGATGGGTTATGTGGTGGAAGTGATTCGTGGTGCTAATAACCAGCGCATCCGCGa                                                                                   <  2:584043/93‑1 (MQ=255)
      tGCCCTTTCCACCATTGGCCGTGTGAATCAGCGGTTTGGGATGGGTTATGTGGTGGAAGTGATTCGTGGTGCTAATAACCAGCGCATCCGCGa                                                                                   <  2:1827929/93‑1 (MQ=255)
       gCCCTTTCCACCATTGGTCGTGTGAATCAGCGGTTTGGGATGGGTTATGTGGTGGAAGTGATTCGTGGTGCTAATAACCAGCGCATCCGCGAc                                                                                  <  1:61342/93‑1 (MQ=255)
        cccTTTCCACCATTGGTCGTGTGAATCAGCGGTTTGGGATGGGTTATGTGGTGGAAGTGATTCGTGGTGCTAATAACCAGCGCATCCGCGACt                                                                                 <  1:3222706/93‑1 (MQ=255)
         ccTTTCCACCATTGGTCGTGTGAATCAGCGGTTTGGGATGGGTTATGTGGTGGAAGTGATTCGTGGTGCTAATAACCAGCGCATCCGCGACTa                                                                                >  1:519018/1‑93 (MQ=255)
         ccTTTCCACCATTGGTCGTGTGAATCAGCGGTTTGGGATGGGTTATGTGGTGGAAGTGATTCGTGGTGCTAATAACCAGCGCATCCGCGACTa                                                                                <  2:5291436/93‑1 (MQ=255)
          cTTTCCACCATTGGTCGTGTGAATCAGCGGTTTGGGATGGGTTATGTGGTGGAAGTGATTCGTGGTGCTAATAACCAGCGCATCCGCGACTAt                                                                               <  2:5644715/93‑1 (MQ=255)
          cTTTCCACCATTGGTCGTGTGAATCAGCGGCTTGGGATGGGTTATGTGGTGGAAGTGATTCGTGGTGCTAATAACCAGCGCATCCGCGACTAt                                                                               >  1:5782918/1‑93 (MQ=255)
           tttCCACCATTGGTCGTGTGAATCAGCGGTTTGGGATGGGTTATGTGGTGGAAGTGATTCGTGGTGCTAATAACCAGCGCATACGCGACTATg                                                                              >  2:3073812/1‑93 (MQ=255)
            ttCCACCATTGGTCGTGTGAATCAGCGGTTTGGGATGGGTTATGTGGTGGAAGTGATTCGTGGTGCTAATAACCAGCGCATCCGCGACTATgg                                                                             >  1:6215690/1‑93 (MQ=255)
             tCCACCATTGGTCGTGTGAATCAGCGGTTTGGGATGGGTTATGTGGTGGAAGTGATTCGTGGTGCTAATAACCAGCGCATCCGCGACTATGGt                                                                            <  2:520977/93‑1 (MQ=255)
              ccaccaTTGGTCGTGTGAATCAGCGGTTTGGGATGGGTTATGTGGTGGAAGTGATTCGTGGTGCTAATAACCAGCGCATCCGCGACTATGGTc                                                                           >  1:3673566/1‑93 (MQ=255)
              ccaccaTTGGTCGTGTGAATCAGCGGTTTGGGATGGGTTATGTGGTGGAAGTGATTCGTGGTGCTAATAACCAGCGCATCCGCGACTATGGTc                                                                           >  1:2352357/1‑93 (MQ=255)
              ccaccaTTGGTCGTGTGAATCAGCGGTTTGGGATGGGTTATGTGGTGGAAGTGATTCGTGGTGCTAATAACCAGCGCATCCGCGACTATGGTc                                                                           <  1:1492850/93‑1 (MQ=255)
              ccaccaTTGGTCGTGTGAATCAGCGGTTTGGGATGGGTTATGTGGTGGAAGTGATTCGTGGTGCTAATAACCAGCGCATCCGCGACTATGGTc                                                                           <  2:3384712/93‑1 (MQ=255)
               caccaTTGGTCGTGTGAATCAGCGGTTTGGGATGGGTTATGTGGTGGAAGTGATTCGTGGTGCTAATAACCAGCGCATCCGCGACTATGGTCa                                                                          >  2:1572688/1‑93 (MQ=255)
                accaTTGGTCGTGTGAATCAGCGGTTTGGGATGGGTTATGTGGTGGAAGTGATTCGTGGTGCTAATAACCAGCGCATCCGCGACTATGGTCAt                                                                         <  1:726790/93‑1 (MQ=255)
                accaTTGGTCGTGTGAATCAGCGGTTTGGGATGGGTTATGTGGTGGAAGTGATTCGTGGTGCTAATAACCAGCGCATCCGCGACTATGGTCAt                                                                         >  2:1781216/1‑93 (MQ=255)
                  caTTGGTCGTGTGAATCAGCGGTTTGGGATGGGTTATGTGGTGGAAGTGATTCGTGGTGCTAATAACCAGCGCATCCGCGACTATGGTCATGa                                                                       <  1:1975440/93‑1 (MQ=255)
                   aTTGGTCGTGTGAATCAGCGGTTTGGGATGGGTTATGTGGTGGAAGTGATTCGTGGTGCTAATAACCAGCGCATCCGCGACTATGGTCATGAc                                                                      >  1:5495568/1‑93 (MQ=255)
                   aTTGGTCGTGTGAATCAGCGGTTTGGGATGGGTTATGTGGTGGAAGTGATTCGTGGTGCTAATAACCAGCGCATCCGCGACTATGGTCATGAc                                                                      <  1:5069442/93‑1 (MQ=255)
                     tGGTCGTGTGAATCAGCGGTTTGGGATGGGTTATGTGGTGGAAGTGATTCGTGGTGCTAATAACCAGCGCATCCGCGACTATGGTCATGACaa                                                                    <  1:3998508/93‑1 (MQ=255)
                      ggTCGTGTGAATCAGCGGTTTGGGATGGGTTATGTGGTGGAAGTGATTCGTGGTGCTAATAACCAGCGCATCCGCGACTATGGTCATGACaaa                                                                   <  2:4022612/93‑1 (MQ=255)
                        tCGTGTGAATCAGCGGTTTGGGATGGGTTATGTGGTGGAAGTGATTCGTGGTGCTAATAACCAGCGCATCCGCGACTATGGTCATGACAAACt                                                                 <  2:1732574/93‑1 (MQ=255)
                          gtgtGAATCAGCGGTTTGGGATGGGTTATGTGGTGGAAGTGATTCGTGGTGCTAATAACCAGCGCATCCGCGACTATGGTCATGACAAACTGa                                                               <  2:2623784/93‑1 (MQ=255)
                          gtgtGAATCAGCGGTTTGGGATGGGTTATGTGGTGGAAGTGATTCGTGGTGCTAATAACCAGCGCATCCGCGACTATGGTCATGACAAACTGa                                                               <  1:853790/93‑1 (MQ=255)
                            gtgAATCAGCGGTTTGGGATGGGTTATGTGGTGGAAGTGATTCGTGGTGCTAATAACCAGCGCATCCGCGACTATGGTCATGACAAACTGaaa                                                             <  1:2550057/93‑1 (MQ=255)
                             tgAATCAGCGGTTTGGGATGGGTTATGTGGTGGAAGTGATTCGTGGTGCTAATAACCAGCGCATCCGCGACTATGGTCATGACAAACTGAAAg                                                            <  2:3719841/93‑1 (MQ=255)
                              gAATCAGCGGTTTGGGATGGGTTATGTGGTGGAAGTGATTCGTGGTGCTAATAACCAGCGCATCCGCGACTATGGTCATGACAAACTGAAAGt                                                           >  2:3540808/1‑93 (MQ=255)
                              gAATCAGCGGTTTGGGATGGGTTATGTGGTGGAAGTGATTCGTGGTGCTAATAACCAGCGCATCCGCGACTATGGTCATGACAAACTGAAAGt                                                           >  2:2114309/1‑93 (MQ=255)
                                aTCAGCGGTTTGGGATGGGTTATGTGGTGGAAGTGATTCGTGGTGCTAATAACCAGCGCATCCGCGACTATGGTCATGACAAACTGAAAGTCt                                                         <  2:2373372/93‑1 (MQ=255)
                                 tCAGCGGTTTGGGATGGGTTATGTGGTGGAAGTGATTCGTGGTGCTAATAACCAGCGCATCCGCGACTATGGTCATGACAAACTGAAAGTCTa                                                        <  1:3536929/93‑1 (MQ=255)
                                 tCAGCGGTTTGGGATGGGTTATGTGGTGGAAGTGATTCGTGGTGCTAATAACCAGCGCATCCGCGACTATGGTCATGACAAACTGAAAGTCTa                                                        >  1:1415964/1‑93 (MQ=255)
                                 tCAGCGGTTTGGGATGGGTTATGTGGTGGAAGTGATTCGTGGTGCTAATAACCAGCGCATCCGCGACTATGGTCATGACAAACTGAAAGTCTa                                                        <  1:5556812/93‑1 (MQ=255)
                                 tCAGCGGTTTGGGATGGGTTATGTGGTGGAAGTGATTCGTGGTGCTAATAACCAGCGCATCCGCGACTATGGTCATGACAAACTGAAAGTCTa                                                        <  2:1676305/93‑1 (MQ=255)
                                    gCGGTTTGGGATGGGTTATGTGGTGGAAGTGATTCGTGGTGCTAATAACCAGCGCATCCGCGACTATGGTCATGACAAACTGAAAGTCTATgg                                                     <  1:3464718/93‑1 (MQ=255)
                                        tttGGGATGGGTTATGTGGTGGAAGTGATTCGTGGTGCTAATAACCAGCGCATCCGCGACTATGGTCATGACAAACTGAAAGTCTATGGCATg                                                 >  1:1190907/1‑93 (MQ=255)
                                           gggATGGGTTATGTGGTGGAAGTGATTCGTGGTGCTAATAACCAGCGCATCCGCGACTATGGTCATGACAAACTGAAAGTCTATGGCATGGGc                                              >  2:5788524/1‑93 (MQ=255)
                                           gggATGGGTTATGTGGTGGAAGTGATTCGTGGTGCTAATAACCAGCGCATCCGCGACTATGGTCATGACAAACTGAAAGTCTATGGCATGGGc                                              <  2:3361650/93‑1 (MQ=255)
                                            ggATGGGTTATGTGGTGGAAGTGATTCGTGGTGCTAATAACCAGCGCATCCGCGACTATGGTCATGACAAACTGAAAGTCTATGGCATGGGcc                                             >  1:2880954/1‑93 (MQ=255)
                                            ggATGGGTTATGTGGTGGAAGTGATTCGTGGTGCTAATAACCAGCGCATCCGCGACTATGGTCATGACAAACTGAAAGTCTATGGCATGGGcc                                             <  2:4883193/93‑1 (MQ=255)
                                                  gTTATGTGGTGGAAGTGATTCGTGGTGCTAATAACCAGCGCATCCGCGACTATGGTCATGACAAACTGAAAGTCTATGGCATGGGCCGTGATa                                       >  1:3395212/1‑93 (MQ=255)
                                                  gTTATGTGGTGGAAGTGATTCGTGGTGCTAATAACCAGCGCATCCGCGACTATGGTCATGACAAACTGAAAGTCTATGGCATGGGCCGTGATa                                       >  2:4450130/1‑93 (MQ=255)
                                                   ttATGTGGTGGAAGTGATTCGTGGTGCTAATAACCAGCGCATCCGCGACTATGGTCATGACAAACTGAAAGTCTATGGCATGGGCCGTGATaa                                      >  2:5469532/1‑93 (MQ=255)
                                                    tATGTGGTGGAAGTGATTCGTGGTGCTAATAACCAGCGCATCCGCGACTATGGTCATGACAAACTGAAAGTCTATGGCATGGGCCGTGATaaa                                     >  1:1667334/1‑93 (MQ=255)
                                                      tgtgGTGGAAGTGATTCGTGGTGCTAATAACCAGCGCATCCGCGACTATGGTCATGACAAACTGAAAGTCTATGGCATGGGCCGTGATAAAAg                                   <  1:1767291/93‑1 (MQ=255)
                                                       gtggtgGAAGTGATTCGTGGTGCTAATAACCAGCGCATCCGCGACTATGGTCATGACAAACTGAAAGTCTATGGCATGGGCCGTGATAAAAGc                                  <  2:2633775/93‑1 (MQ=255)
                                                       gtggtgGAAGTGATTCGTGGTGCTAATAACCAGCGCATCCGCGACTATGGTCATGACAAACTGAAAGTCTATGGCATGGGCCGTGATAAAAGc                                  <  1:4868794/93‑1 (MQ=255)
                                                          gtggAAGTGATTCGTGGTGCTAATAACCAGCGCATCCGCGACTATGGTCATGACAAACTGAAAGTCTATGGCATGGGCCGTGATAAAAGCCAt                               <  1:1040223/93‑1 (MQ=255)
                                                           tggAAGTGATTCGTGGTGCTAATAACCAGCGCATCCGCGACTATGGTCATGACAAACTGAAAGTCTATGGCATGGGCCGTGATAAAAGCCATg                              <  2:3247104/93‑1 (MQ=255)
                                                           tggAAGTGATTCGTGGTGCTAATAACCAGCGCATCCGCGACTATGGTCATGACAAACTGAAAGTCTATGGCATGGGCCGTGATAAAAGCCATg                              >  2:2327904/1‑93 (MQ=255)
                                                             gAAGTGATTCGTGGTGCTAATAACCAGCGCATCCGCGACTATGGTCATGACAAACTGAAAGTCTATGGCATGGGCCGTGATAAAAGCCATGaa                            >  1:3631226/1‑93 (MQ=255)
                                                             gAAGTGATTCGTGGTGCTAATAACCAGCGCATCCGCGACTATGGTCATGACAAACTGAAAGTCTATGGCATGGGCCGTGATAAAAGCCATGaa                            <  2:2481855/93‑1 (MQ=255)
                                                              aaGTGATTCGTGGTGCTAATAACCAGCGCATCCGCGACTATGGTCATGACAAACTGAAAGTCTATGGCATGGGCCGTGATAAAAGCCATGAAc                           >  1:1365114/1‑93 (MQ=255)
                                                                gTGATTCGTGGTGCTAATAACCAGCGCATCCGCGACTATGGTCATGACAAACTGAAAGTCTATGGCATGGGCCGTGATAAAAGCCATGAACAt                         >  2:3833821/1‑93 (MQ=255)
                                                                gTGATTCGTGGTGCTAATAACCAGCGCATCCGCGACTATGGTCATGACAAACTGAAAGTCTATGGCATGGGCCGTGATAAAAGCCATGAACAt                         >  1:121320/1‑93 (MQ=255)
                                                                  gATTCGTGGTGCTAATAACCAGCGCATCCGCGACTATGGTCATGACAAACTGAAAGTCTATGGCATGGGCCGTGATAAAAGCCATGAACATTg                       >  2:4622974/1‑93 (MQ=255)
                                                                  gATTCGTGGTGCTAATAACCAGCGCATCCGCGACTATGGTCATGACAAACTGAAAGTCTATGGCATGGGCCGTGATAAAAGCCATGAACATTg                       <  2:2076310/93‑1 (MQ=255)
                                                                  gATTCGTGGTGCTAATAACCAGCGCATCCGCGACTATGGTCATGACAAACTGAAAGTCTATGGCATGGGCCGTGATAAAAGCCATGAACATTg                       >  1:5133370/1‑93 (MQ=255)
                                                                  gATTCGTGGTGCTAATAACCAGCGCATCCGCGACTATGGTCATGACAAACTGAAAGTCTATGGCATGGGCCGTGATAAAAGCCATGAACATTg                       >  2:4835530/1‑93 (MQ=255)
                                                                  gATTCGTGGTGCTAATAACCAGCGCATCCGCGACTATGGTCATGACAAACTGAAAGTCTATGGCATGGGCCGTGATAAAAGCCATGAACATTg                       <  2:4810074/93‑1 (MQ=255)
                                                                     tCGTGGTGCTAATAACCAGCGCATCCGCGACTATGGTCATGACAAACTGAAAGTCTATGGCATGGGCCGTGATAAAAGCCATGAACATTGGGt                    >  1:2910997/1‑93 (MQ=255)
                                                                      cGTGGTGCTAATAACCAGCGCATCCGCGACTATGGTCATGACAAACTGAAAGTCTATGGCATGGGCCGTGATAAAAGCCATGAACATTGGGTg                   <  2:1450895/93‑1 (MQ=255)
                                                                      cGTGGTGCTAATAACCAGCGCATCCGCGACTATGGTCATGACAAACTGAAAGTCTATGGCATGGGCCGTGATAAAAGCCATGAACATTGGGTg                   >  1:3750858/1‑93 (MQ=255)
                                                                       gtggtgCTAATAACCAGCGCATCCGCGACTATGGTCATGACAAACTGAAAGTCTATGGCATGGGCCGTGATAAAAGCCATGAACATTGGGTGa                  <  2:4215022/93‑1 (MQ=255)
                                                                       gtggtgCTAATAACCAGCGCATCCGCGACTATGGTCATGACAAACTGAAAGTCTATGGCATGGGCCGTGATAAAAGCCATGAACATTGGGTGa                  >  2:4258319/1‑93 (MQ=255)
                                                                       gtggtgCTAATAACCAGCGCATCCGCGACTATGGTCATGACAAACTGAAAGTCTATGGCATGGGCCGTGATAAAAGCCATGAACATTGGGTGa                  >  1:2654682/1‑93 (MQ=255)
                                                                       gtggtgCTAATAACCAGCGCATCCGCGACTATGGTCATGACAAACTGAAAGTCTATGGCATGGGCCGTGATAAAAGCCATGAACATTGGGTGa                  <  1:426018/93‑1 (MQ=255)
                                                                        tggtgCTAATAACCAGCGCATCCGCGACTATGGTCATGACAAACTGAAAGTCTATGGCATGGGCCGTGATAAAAGCCATGAACATTGGGTGAg                 >  1:4315170/1‑93 (MQ=255)
                                                                        tggtgCTAATAACCAGCGCATCCGCGACTATGGTCATGACAAACTGAAAGTCTATGGCATGGGCCGTGATAAAAGCCATGAACATTGGGTGAg                 >  2:2416106/1‑93 (MQ=255)
                                                                          gtgCTAATAACCAGCGCATCCGCGACTATGGTCATGACAAACTGAAAGTCTATGGCATGGGCCGTGATAAAAGCCATGAACATTGGGTGAGCg               >  1:1058371/1‑93 (MQ=255)
                                                                           tgCTAATAACCAGCGCATCCGCGACTATGGTCATGACAAACTGAAAGTCTATGGCATGGGCCGTGATAAAAGCCATGAACATTGGGTGAGCGt              >  1:4593378/1‑93 (MQ=255)
                                                                           tgCTAATAACCAGCGCATCCGCGACTATGGTCATGACAAACTGAAAGTCTATGGCATGGGCCGTGATAAAAGCCATGAACATTGGGTGAGCGt              <  2:8888/93‑1 (MQ=255)
                                                                           tgCTAATAACCAGCGCATCCGCGACTATGGTCATGACAAACTGAAAGTCTATGGCATGGGCCGTGATAAAAGCCATGAACATTGGGTGAGCGt              <  2:1590839/93‑1 (MQ=255)
                                                                           tgCTAATAACCAGCGCATCCGCGACTATGGTCATGACAAACTGAAAGTCTATGGCATGGGCCGTGATAAAAGCCATGAACATTGGGTGAGCGt              >  1:5323921/1‑93 (MQ=255)
                                                                            gCTAATAACCAGCGCATCCGCGACTATGGTCATGACAAACTGAAAGTCTATGGCATGGGCCGTGATAAAAGCCATGAACATTGGGTGAGCGTg             >  1:2702180/1‑93 (MQ=255)
                                                                             cTAATAACCAGCGCATCCGCGACTATGGTCATGACAAACTGAAAGTCTATGGCATGGGCCGTGATAAAAGCCATGAACATTGGGTGAGCGTGa            <  1:2259380/93‑1 (MQ=255)
                                                                             cTAATAACCAGCGCATCCGCGACTATGGTCATGACAAACTGAAAGTCTATGGCATGGGCCGTGATAAAAGCCATGAACATTGGGTGAGCGTGa            <  2:649910/93‑1 (MQ=255)
                                                                              taataaCCAGCGCATCCGCGACTATGGTCATGACAAACTGAAAGTCTATGGCATGGGCCGTGATAAAAGCCATGAACATTGGGTGAGCGTGAt           >  1:896443/1‑93 (MQ=255)
                                                                               aataaCCAGCGCATCCGCGACTATGGTCATGACAAACTGAAAGTCTATGGCATGGGCCGTGATAAAAGCCATGAACATTGGGTGAGCGTGATc          >  2:4525320/1‑93 (MQ=255)
                                                                               aataaCCAGCGCATCCGCGACTATGGTCATGACAAACTGAAAGTCTATGGCATGGGCCGTGATAAAAGCCATGAACATTGGGTGAGCGTGATc          >  2:507866/1‑93 (MQ=255)
                                                                               aataaCCAGCGCATCCGCGACTATGGTCATGACAAACTGAAAGTCTATGGCATGGGCCGTGATAAAAGCCATGAACATTGGGTGAGCGTGATc          >  1:190442/1‑93 (MQ=255)
                                                                                ataaCCAGCGCATCCGCGACTATGGTCATGACAAACTGAAAGTCTATGGCATGGGCCGTGATAAAAGCCATGAACATTGGGTGAGCGTGATcc         >  1:5930776/1‑93 (MQ=255)
                                                                                ataaCCAGCGCATCCGCGACTATGGTCATGACAAACTGAAAGTCTATGGCATGGGCCGTGATAAAAGCCATGAACATTGGGTGAGCGTGATcc         >  1:4172560/1‑93 (MQ=255)
                                                                                  aaCCAGCGCATCCGCGACTATGGTCATGACAAACTGAAAGTCTATGGCATGGGCCGTGATAAAAGCCATGAACATTGGGTGAGCGTGATCCGc       <  1:296157/93‑1 (MQ=255)
                                                                                   aCCAGCGCATCCGCGACTATGGTCATGACAAACTGAAAGTCTATGGCATGGGCCGTGATAAAAGCCATGAACATTGGGTGAGCGTGATCCGcc      >  2:1153896/1‑93 (MQ=255)
                                                                                   aCCAGCGCATCCGCGACTATGGTCATGACAAACTGAAAGTCTATGGCATGGGCCGTGATAAAAGCCATGAACATTGGGTGAGCGTGATCCGcc      >  1:6041194/1‑93 (MQ=255)
                                                                                   aCCAGCGCATCCGCGACTATGGTCATGACAAACTGAAAGTCTATGGCATGGGCCGTGATAAAAGCCATGAACATTGGGTGAGCGTGATCCGcc      <  1:1277214/93‑1 (MQ=255)
                                                                                      aGCGCATCCGCGACTATGGTCATGACAAACTGAAAGTCTATGGCATGGGCCGTGATAAAAGCCATGAACATTGGGTGAGCGTGATCCGCCAGc   <  2:1339146/93‑1 (MQ=255)
                                                                                       gcgcATCCGCGACTATGGTCATGACAAACTGAAAGTCTATGGCATGGGCCGTGATAAAAGCCATGAACATTGGGTGAGCGTGATCCGCCAGCt  >  2:120518/1‑93 (MQ=255)
                                                                                       gcgcATCCGCGACTATGGTCATGACAAACTGAAAGTCTATGGCATGGGCCGTGATAAAAGCCATGAACATTGGGTGAGCGTGATCCGCCAGCt  <  1:2066960/93‑1 (MQ=255)
                                                                                          |                                                                                         
TCAGATTGCCCTTTCCACCATTGGTCGTGTGAATCAGCGGTTTGGGATGGGTTATGTGGTGGAAGTGATTCGTGGTGCTAATAACCAGCGTATCCGCGACTATGGTCATGACAAACTGAAAGTCTATGGCATGGGCCGTGATAAAAGCCATGAACATTGGGTGAGCGTGATCCGCCAGCT  >  NC_000913/4004722‑4004901

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