Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 177,386 C→T R80H (CGC→CAC)  yadS ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913177,3860CT99.6% 838.0 / NA 232R80H (CGC→CAC) yadShypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (101/130);  minor base G (1/0);  total (102/130)
Fisher's exact test for biased strand distribution p-value = 4.40e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 248 total aligned reads displayed.

AAGGCTTTATTCACGCCAATGCCGACAAACACCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGCGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCAAAATACCGGGCCGTGATCCAGC  >  NC_000913/177296‑177475
                                                                                          |                                                                                         
    cTTTATTCACGCCAATGCCGACAAACACCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGTGCACCa                                                                                     >  1:2073092/1‑93 (MQ=255)
    cTTTATTCACGCCAATGCCGACAAACACCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGTGCACCa                                                                                     <  1:5343333/93‑1 (MQ=255)
     tttATTCACGCCAATGCCGACAAACACCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGTGCACCAg                                                                                    >  1:1244243/1‑93 (MQ=255)
       tATTCACGCCAATGCCGACAAACACCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGTGCACCAGCa                                                                                  <  1:3417461/93‑1 (MQ=255)
       tATTCACGCCAATGCCGACAAACACCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGTGCACCAGCa                                                                                  <  1:5423813/93‑1 (MQ=255)
         ttCACGCCAATGCCGACAAACACCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGTGCACCAGCACg                                                                                >  1:3054251/1‑93 (MQ=255)
           cACGCCAATGCCGACAAACACCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGTGCACCAGCACGAt                                                                              >  1:6188791/1‑93 (MQ=255)
           cACGCCAATGCCGACAAACACCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGTGCACCAGCACGAt                                                                              <  2:4276164/93‑1 (MQ=255)
           cACGCCAATGCCGACAAACACCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGTGCACCAGCACGAt                                                                              <  2:432832/93‑1 (MQ=255)
           cACGCCAATGCCGACAAACACCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGTGCACCAGCACGAt                                                                              >  1:2966097/1‑93 (MQ=255)
           cACGCCAAGACGGACAAACACCGCCACACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGTGCACCAGCACGAt                                                                              <  2:2207602/93‑1 (MQ=255)
            aCGCCAATGCGGACAAACACCGCCAGACCAACGGCGGCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGTGCACCAGCACGATg                                                                             <  2:2845566/93‑1 (MQ=255)
             cGCCAATGCCGACAAACACCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGTGCACCAGCACGATgg                                                                            <  2:6047925/93‑1 (MQ=255)
              gCCAATGCCGACAAACACCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGTGCACCAGCACGATGGt                                                                           >  2:1487866/1‑93 (MQ=255)
              gCCAATGCCGACAAACACCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGTGCACCAGCACGATGGt                                                                           >  2:443414/1‑93 (MQ=255)
               ccAATGCCGACAAACACCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGTGCACCAGCACGATGGTc                                                                          <  1:575307/93‑1 (MQ=255)
                cAATGCCGACAAACACCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGTGCACCAGCACGATGGTCa                                                                         >  1:2231745/1‑93 (MQ=255)
                cAATGCCGACAAACACCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGTGCACCAGCACGATGGTCa                                                                         <  2:3003806/93‑1 (MQ=255)
                   tGCCGACAAACACCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGTGCACCAGCACGATGGTCAGCa                                                                      <  2:4049647/93‑1 (MQ=255)
                      cGACAAACACCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGTGCACCAGCACGATGGTCAGCATGc                                                                   >  2:5269826/1‑93 (MQ=255)
                      cGACAAACACCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGTGCACCAGCACGATGGTCAGCATGc                                                                   <  1:3475735/93‑1 (MQ=255)
                      cGACAAACACCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGTGCACCAGCACGATGGTCAGCATGc                                                                   <  1:4523647/93‑1 (MQ=255)
                       gACAAACACCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGTGCACCAGCACGATGGTCAGCATGCt                                                                  <  1:4939600/93‑1 (MQ=255)
                        aCAAACACCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGTGCACCAGCACGATGGTCAGCATGCTg                                                                 >  1:234399/1‑93 (MQ=255)
                         cAAACACCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGTGCACCAGCACGATGGTCAGCATGCTgg                                                                <  1:964902/93‑1 (MQ=255)
                         cAAACACCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGTGCACCAGCACGATGGTCAGCATGCTgg                                                                <  2:5679261/93‑1 (MQ=255)
                          aaaCACCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGTGCACCAGCACGATGGTCAGCATGCTGGt                                                               <  1:3780763/93‑1 (MQ=255)
                          aaaCACCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGTGCACCAGCACGATGGTCAGCATGCTGGt                                                               >  1:4732946/1‑93 (MQ=255)
                           aaCACCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGTGCACCAGCACGATGGTCAGCATGCTGGTg                                                              >  1:1997564/1‑93 (MQ=255)
                           aaCACCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGTGCACCAGCACGATGGTCAGCATGCTGGTg                                                              >  1:4649479/1‑93 (MQ=255)
                            acacCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGTGCACCAGCACGATGGTCAGCATGCTGGTGa                                                             >  2:5648873/1‑93 (MQ=255)
                            acacCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGTGCACCAGCACGATGGTCAGCATGCTGGTGa                                                             <  2:5066343/93‑1 (MQ=255)
                             cacCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGTGCACCAGCACGATGGTCAGCATGCTGGTGAc                                                            <  2:5750733/93‑1 (MQ=255)
                              cccGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGTGCACCAGCACGATGGTCAGCATGCTGGTGAcc                                                           <  2:2321089/92‑1 (MQ=255)
                              acCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGTGCACCAGCACGATGGTCAGCATGCTGGTGAcc                                                           >  2:1603405/1‑93 (MQ=255)
                               ccGCCAGACCAACGGCGTCCAGCACCGGCAATATCCATTTTGGTAAGCGTCTTGGCTGGTGCACCAGCACGATGGTCAGCATGCTGGTGACCa                                                          >  1:3670645/1‑93 (MQ=255)
                               ccGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGTGCACCAGCACGATGGTCAGCATGCTGGTGACCa                                                          >  2:4633771/1‑93 (MQ=255)
                               ccGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGTGCACCAGCACGATGGTCAGCATGCTGGTGACCa                                                          <  2:195546/93‑1 (MQ=255)
                               ccGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGTGCACCAGCACGATGGTCAGCATGCTGGTGACCa                                                          <  1:3998875/93‑1 (MQ=255)
                               ccGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGTGCACCAGCACGATGGTCAGCATGCTGGTGACCa                                                          <  1:5170595/93‑1 (MQ=255)
                               ccGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGTGCACCAGCACGATGGTCAGCATGCTGGTGACCa                                                          >  1:3231773/1‑93 (MQ=255)
                                cGCCAGACCAACGGCGTCGAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGTGCACCAGCACGATGGTCAGCATGCTGGTGACCAt                                                         <  2:311273/93‑1 (MQ=255)
                                 gCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGTGCACCAGCACGATGGTCAGCATGCTGGTGACCAtt                                                        >  1:5642055/1‑93 (MQ=255)
                                  ccAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGTGCACCAGCACGATGGTCAGCATGCTGGTGACCATTg                                                       <  1:426223/93‑1 (MQ=255)
                                     gACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGTGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCaa                                                    <  1:2300892/93‑1 (MQ=255)
                                     gACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGTGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCaa                                                    <  2:2963043/93‑1 (MQ=255)
                                      aCCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGTGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAAc                                                   <  2:5905789/93‑1 (MQ=255)
                                        cAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGTGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGa                                                 >  1:4211721/1‑93 (MQ=255)
                                        cAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGTGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGa                                                 >  2:1894188/1‑93 (MQ=255)
                                         aaCGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGTGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGAc                                                >  2:1186226/1‑93 (MQ=255)
                                           cGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGTGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCa                                              <  2:3400015/93‑1 (MQ=255)
                                            ggCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGTGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAg                                             >  1:5950792/1‑93 (MQ=255)
                                             gCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGTGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGa                                            <  2:4002427/93‑1 (MQ=255)
                                             gCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGTGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGa                                            <  1:2555978/93‑1 (MQ=255)
                                             gCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGTGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGa                                            <  2:3627209/93‑1 (MQ=255)
                                              cGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGTGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGAt                                           <  1:1575567/93‑1 (MQ=255)
                                               gTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGTGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATc                                          <  1:4328111/93‑1 (MQ=255)
                                               gTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGTGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATc                                          >  1:1962628/1‑93 (MQ=255)
                                                tCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGTGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCg                                         >  1:2114278/1‑93 (MQ=255)
                                                 ccAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGTGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCgg                                        <  1:3094878/93‑1 (MQ=255)
                                                  cAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGTGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGt                                       <  1:1305237/93‑1 (MQ=255)
                                                  cAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGTGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGt                                       <  2:3889333/93‑1 (MQ=255)
                                                  cAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGTGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGt                                       >  2:5986338/1‑93 (MQ=255)
                                                   aGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGTGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTg                                      <  1:4349457/93‑1 (MQ=255)
                                                    gCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGTGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTgg                                     <  2:275779/93‑1 (MQ=255)
                                                     cACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGTGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTggg                                    <  1:5751039/93‑1 (MQ=255)
                                                      aCCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGTGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGa                                   >  2:3162506/1‑93 (MQ=255)
                                                       ccGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGTGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGAt                                  <  1:1452117/93‑1 (MQ=255)
                                                        cGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGTGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATc                                 <  2:2368224/93‑1 (MQ=255)
                                                         ggCAACATCCATTTTGGTAAGCGTCTTGGCTGGTGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCt                                >  2:5541210/1‑93 (MQ=255)
                                                          gCAACATCCATTTTGGTAAGCGTCTTGGCTGGTGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCtt                               <  2:2471341/93‑1 (MQ=255)
                                                          gCAACATCCATTTTGGTAAGCGTCTTGGCTGGTGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCtt                               <  2:5448929/93‑1 (MQ=255)
                                                           cAACATCCATTTTGGTAAGCGTCTTGGCTGGTGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCttt                              <  1:313767/93‑1 (MQ=255)
                                                               aTCCATTTTGGTAAGCGTCTTGGCTGGTGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCAcc                          <  1:5485283/93‑1 (MQ=255)
                                                                  cATTTTGGTAAGCGTCTTGGCTGGTGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCaa                       >  2:4897274/1‑93 (MQ=255)
                                                                  cATTTTGGTAAGCGTCTTGGCTGGTGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCaa                       >  1:748692/1‑93 (MQ=255)
                                                                  cATTTTGGTAAGCGTCTTGGCTGGTGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCaa                       >  2:3037754/1‑93 (MQ=255)
                                                                    ttttGGTAAGCGTCTTGGCTGGTGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCaaaa                     <  2:4351433/93‑1 (MQ=255)
                                                                       tGGTAAGCGTCTTGGCTGGTGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCAAAATAc                  <  2:4768504/93‑1 (MQ=255)
                                                                        ggTAAGCGTCTTGGCTGGTGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCAAAATAcc                 <  1:824745/93‑1 (MQ=255)
                                                                        ggTAAGCGTCTTGGCTGGTGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCAAAATAcc                 <  1:2472215/93‑1 (MQ=255)
                                                                          tAAGCGTCTTGGCTGGTGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCAAAATACCgg               <  2:1762657/93‑1 (MQ=255)
                                                                           aaGCGTCTTGGCTGGTGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCAAAATACCggg              <  2:852460/93‑1 (MQ=255)
                                                                           aaGCGTCTTGGCTGGTGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCAAAATACCggg              <  1:2670102/93‑1 (MQ=255)
                                                                           aaGCGTCTTGGCTGGTGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCAAAATACCggg              >  2:5827212/1‑93 (MQ=255)
                                                                           aaGCGTCTTGGCTGGTGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCAAAATACCggg              <  1:3625834/93‑1 (MQ=255)
                                                                               gTCTTGGCTGGTGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCAAAATACCGGGCCgg          >  2:1136796/1‑92 (MQ=255)
                                                                               gTCTTGGCTGGTGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCAAAATACCGGGCCGt          <  1:1351148/93‑1 (MQ=255)
                                                                                tCTTGGCTGGTGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCAAAATACCGGGCCGTg         >  2:2792153/1‑93 (MQ=255)
                                                                                tCTTGGCTGGTGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCAAAATACCGGGCCGTg         >  1:3916362/1‑93 (MQ=255)
                                                                                tCTTGGCTGGTGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCAAAATACCGGGCCGTg         >  2:1304684/1‑93 (MQ=255)
                                                                                tCTTGGCTGGGGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCCTTCACCCAAAATACCGGGCCGTg         >  2:76113/1‑93 (MQ=255)
                                                                                  ttGGCTGGTGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGTATCTTTCACCCAAAATACCGGGCCGTGAt       <  2:6087873/93‑1 (MQ=255)
                                                                                   tggctggtGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCAAAATACCGGGCTGTGATc      >  1:2814357/1‑93 (MQ=255)
                                                                                    ggctggtGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCAAAATACCGGGCCGTGATcc     <  1:5266797/93‑1 (MQ=255)
                                                                                    ggctggtGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCAAAATACCGGGCCGTGATcc     <  1:5850556/93‑1 (MQ=255)
                                                                                     gctggtGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCAAAATACCGGGCCGTGATCCa    >  1:5687412/1‑93 (MQ=255)
                                                                                     gctggtGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCAAAATACCGGGCCGTGATCCa    <  1:6140589/93‑1 (MQ=255)
                                                                                       tggtGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCAAAATACCGGGCCGTGATCCAgc  >  1:4845523/1‑93 (MQ=255)
                                                                                       tggtGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCAAAATACCGGGCCGTGATCCAgc  >  2:4254521/1‑93 (MQ=255)
                                                                                          |                                                                                         
AAGGCTTTATTCACGCCAATGCCGACAAACACCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGCGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCAAAATACCGGGCCGTGATCCAGC  >  NC_000913/177296‑177475

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