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.7% 1077.0 / ‑5.6 312R80H (CGC→CAC) yadShypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (159/151);  total (160/152)
Fisher's exact test for biased strand distribution p-value = 4.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.95e-01

Only 100 of 340 total aligned reads displayed.

AAGGCTTTATTCACGCCAATGCCGACAAACACCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGCGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCAAAATACCGGGCCGTGATCCAGC  >  NC_000913/177296‑177475
                                                                                          |                                                                                         
aaGGCTTTATTCACGCCAATGCCGACAAACACCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGgtgc                                                                                         >  1:1286999/1‑93 (MQ=255)
aaGGCTTTATTCACGCCAATGCCGACAAACACCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGgtgc                                                                                         <  2:3443647/93‑1 (MQ=255)
aaGGCTTTATTCACGCCAATGCCGACAAACACCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGgtgc                                                                                         >  2:6250506/1‑93 (MQ=255)
 aGGCTTTATTCACGCCAATGCCGACAAACACCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGTGCa                                                                                        >  1:6297241/1‑93 (MQ=255)
  ggCTTTATTCACGCCAATGCCGACAAACACCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGTGCAc                                                                                       <  1:719441/93‑1 (MQ=255)
   gCTTTATTCACGCCAATGCCGACAAACACCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGTGCAcc                                                                                      >  1:1606729/1‑93 (MQ=255)
   gCTTTATTCACGCCAATGCCGACAAACACCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGTGCAcc                                                                                      <  2:6786040/93‑1 (MQ=255)
   gCTTTATTCACGCCAATGCCGACAAACACCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGTGCAcc                                                                                      >  2:2442286/1‑93 (MQ=255)
   gCTTTATTCACGCCAATGCCGACAAACACCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGTGCAcc                                                                                      >  1:3511158/1‑93 (MQ=255)
    cTTTATTCACGCCAATGCCGACAAACACCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGTGCACCa                                                                                     <  2:3004469/93‑1 (MQ=255)
     tttATTCACGCCAATGCCGACAAACACCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGTGcagcag                                                                                    >  2:2066911/1‑93 (MQ=255)
     tttATTCACGCCAATGCCGACAAACACCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGTGCACCAg                                                                                    >  1:7227550/1‑93 (MQ=255)
       tATTCACGCCAATGCCGACAAACACCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGTGCACCAGCa                                                                                  <  1:6415531/93‑1 (MQ=255)
        aTTCACGCCAATGCCGACAAACACCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGTGCACCAGCAc                                                                                 <  2:7153398/93‑1 (MQ=255)
         ttCACGCCAATGCCGACAAACACCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGTGCACCAGCACg                                                                                >  2:1974283/1‑93 (MQ=255)
         ttCACGCCAATGCCGACAAACACCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGTGCACCAGCACg                                                                                >  1:5303758/1‑93 (MQ=255)
            aCGCCAATGCCGACAAACACCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGTGCACCAGCACGATg                                                                             >  2:3664772/1‑93 (MQ=255)
             cGCCAATGCCGACAAACACCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGTGCACCAGCACGATgg                                                                            >  2:5780959/1‑93 (MQ=255)
             cGCCAATGCCGACAAACACCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGTGCACCAGCACGATgg                                                                            <  2:3514822/93‑1 (MQ=255)
              gCCAATGCCGACAAACACCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGTGCACCAGCACGATGGt                                                                           <  2:3287269/93‑1 (MQ=255)
              cccAATGCCGACAAACACCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGTGCACCAGCACGATGGt                                                                           <  1:5613099/92‑1 (MQ=255)
               ccAATGCCGACAAACACCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGTGCACCAGCACGATGGTc                                                                          >  2:1763964/1‑93 (MQ=255)
               ccAATGCCGACAAACACCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGTGAACCAGCACGATGGTc                                                                          >  2:5531988/1‑93 (MQ=255)
               ccAATGCCGACAAACACCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGCCTGGGGCACCAGCACGATGGTc                                                                          >  1:2687281/1‑93 (MQ=255)
                 aaTGCCGACAAACACCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGTGCACCAGCACGATGGTCAg                                                                        <  1:5798999/93‑1 (MQ=255)
                   tGCCGACAAACACCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGTGCACCAGCACGATGGTCAGCa                                                                      <  1:4165221/93‑1 (MQ=255)
                   tGCCGACAAACACCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGTGCACCAGCACGATGGTCAGCa                                                                      >  1:4385496/1‑93 (MQ=255)
                   tGCCGACAAACACCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGTGCACCAGCACGATGGTCAGCa                                                                      <  1:4885758/93‑1 (MQ=255)
                      cGACAAACACCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGTGCACCAGCACGATGGTCAGCATGc                                                                   <  1:3765856/93‑1 (MQ=255)
                      cGACAAACACCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGTGCACCAGCACGATGGTCAGCATGc                                                                   <  2:1367490/93‑1 (MQ=255)
                       gACAAACACCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGTGCACCAGCACGATGGTCAGCATGCt                                                                  <  1:7083031/93‑1 (MQ=255)
                        aCAAACACCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGTGCACCAGCACGATGGTCAGCATGCTg                                                                 >  1:2462876/1‑93 (MQ=255)
                        aCAAACACCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGTGCACCAGCACGATGGTCAGCATGCTg                                                                 <  1:1936584/93‑1 (MQ=255)
                         cAAACACCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGTGCACCAGCACGATGGTCAGCATGCTgg                                                                >  2:4145969/1‑93 (MQ=255)
                          aaaCACCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGTGCACCAGCACGATGGTCAGCATGCTGGt                                                               >  2:1849803/1‑93 (MQ=255)
                           aaCACCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGTGCACCAGCACGATGGTCAGCATGCTGGTg                                                              >  1:359581/1‑93 (MQ=255)
                              acCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGTGCACCAGCACGATGGTCAGCATGCTGGTGAcc                                                           >  2:1907711/1‑93 (MQ=255)
                                cGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGCGCACCAGCACGATGGTCAGCATGCTGGTGACCAt                                                         <  1:5719373/93‑1 (MQ=255)
                                 gCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGTGCACCAGCACGATGGTCAGCATGCTGGTGACCAtt                                                        <  2:6346754/93‑1 (MQ=255)
                                  ccAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGTGCACCAGCACGATGGTCAGCATGCTGGTGACCATTg                                                       >  2:423774/1‑93 (MQ=255)
                                   cAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGTGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGc                                                      <  1:209627/93‑1 (MQ=255)
                                    aGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGTGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCa                                                     >  1:3001420/1‑93 (MQ=255)
                                     gACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGTGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCaa                                                    >  2:1470422/1‑93 (MQ=255)
                                     gACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGTGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCaa                                                    >  1:2941302/1‑93 (MQ=255)
                                     gACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGTGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCaa                                                    >  1:7277629/1‑93 (MQ=255)
                                     gACCAACGGCGTCCAGCACCAGCAACATCCATTTTGGTAAGCGTCTTGGCTGGTGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCaa                                                    >  2:4874142/1‑93 (MQ=255)
                                      aCCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGTGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAAc                                                   <  2:2719344/93‑1 (MQ=255)
                                       ccAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGTGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACg                                                  >  1:3554920/1‑93 (MQ=255)
                                        cAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGTGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGa                                                 <  2:2525103/93‑1 (MQ=255)
                                         aaCGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGTGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGAc                                                >  2:2825945/1‑93 (MQ=255)
                                         aaCGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGTGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGAc                                                >  1:5524127/1‑93 (MQ=255)
                                         aaCGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGTGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGAc                                                >  2:3892640/1‑93 (MQ=255)
                                         aaCGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGTGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGAc                                                >  2:5145350/1‑93 (MQ=255)
                                           cGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGTGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCa                                              <  2:1242566/93‑1 (MQ=255)
                                             gCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGTGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGa                                            <  1:6227370/93‑1 (MQ=255)
                                             gCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGTGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGa                                            <  2:6971521/93‑1 (MQ=255)
                                              cGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGTGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGAt                                           <  1:4648870/93‑1 (MQ=255)
                                               gTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGTGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATc                                          <  2:7932956/93‑1 (MQ=255)
                                               gTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGTGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATc                                          <  1:3159150/93‑1 (MQ=255)
                                               gTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGTGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATc                                          <  1:1716087/93‑1 (MQ=255)
                                                  cAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGTGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGt                                       >  2:5699763/1‑93 (MQ=255)
                                                  cAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGTGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGt                                       <  1:457863/93‑1 (MQ=255)
                                                    gCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGTGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTgg                                     <  2:7808314/93‑1 (MQ=255)
                                                    gCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGTGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTgg                                     >  1:7861460/1‑93 (MQ=255)
                                                      aCCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGTGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGa                                   >  1:3105410/1‑93 (MQ=255)
                                                       ccGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGTGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGAt                                  <  1:2257667/93‑1 (MQ=255)
                                                        cGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGTGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATc                                 <  1:800652/93‑1 (MQ=255)
                                                        cGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGTGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATc                                 >  1:7629695/1‑93 (MQ=255)
                                                        cGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGTGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATc                                 <  2:4579290/93‑1 (MQ=255)
                                                          gCAACATCCATTTTGGTAAGCGTCTTGGCTGGTGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCtt                               <  2:737656/93‑1 (MQ=255)
                                                            aaCATCCATTTTGGTAAGCGTCTTGGCTGGTGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATTTTTc                             >  2:5392254/1‑93 (MQ=255)
                                                             aCATCCATTTTGGTAAGCGTCTTGGCTGGTGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCa                            <  2:306900/93‑1 (MQ=255)
                                                               aTCCATTTTGGTAAGCGTCTTGGCTGGTGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCAcc                          >  2:2276203/1‑93 (MQ=255)
                                                               aTCCATTTTGGTAAGCGTCTTGGCTGGTGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCAcc                          >  1:7658831/1‑93 (MQ=255)
                                                                tCCATTTTGGTAAGCGTCTTGGCTGGTGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCAccc                         >  1:1822633/1‑93 (MQ=255)
                                                                 ccATTTTGGTAAGCGTCTTGGCTGGTGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCa                        >  2:6187367/1‑93 (MQ=255)
                                                                 ccATTTTGGTAAGCGTCTTGGCTGGTGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCa                        <  2:4811269/93‑1 (MQ=255)
                                                                  cATTTTGGTAAGCGTCTTGGCTGGTGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCaa                       <  2:6425490/93‑1 (MQ=255)
                                                                  cATTTTGGTAAGCGTCTTGGCTGGTGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCaa                       <  1:5868392/93‑1 (MQ=255)
                                                                    ttttGGTAAGCGTCTTGGCTGGTGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCaaaa                     >  2:5941857/1‑93 (MQ=255)
                                                                    ttttGGTAAGCGTCTTGGCTGGTGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCaaaa                     <  2:7201309/93‑1 (MQ=255)
                                                                    ttttGGTAAGCGTCTTGGCTGGTGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCaaaa                     <  1:3325523/93‑1 (MQ=255)
                                                                    ttttGGTAAGCGTCTTGGCTGGTGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCaaaa                     >  1:3934395/1‑93 (MQ=255)
                                                                     tttGGTAAGCGTCTTGGCTGGTGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCAAAAt                    >  2:1657695/1‑93 (MQ=255)
                                                                      ttGGTAAGCGTCTTGGCTGGTGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCAAAATa                   >  1:4483837/1‑93 (MQ=255)
                                                                      ttGGTAAGCGTCGTGGCTGGTGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCAAAATa                   <  2:4060764/93‑1 (MQ=255)
                                                                       tGGTAAGCGTCTTGGCTGGTGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCAAAATAc                  >  2:4372032/1‑93 (MQ=255)
                                                                            aGCGTCTTGGCTGGTGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCAAAATACCGGGc             <  1:1007083/93‑1 (MQ=255)
                                                                             gCGTCTTGGCTGGTGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCAAAATACCGGGcc            >  2:2502822/1‑93 (MQ=255)
                                                                             gCGTCTTGGCTGGTGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCAAAATACCGGGcc            >  1:14244/1‑93 (MQ=255)
                                                                             gCGTCTTGGCTGGTGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCAAAATACCGGGcc            >  1:1161221/1‑93 (MQ=255)
                                                                             gCGTCTTGGCTGGTGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCAAAATACCGGGcc            <  1:933400/93‑1 (MQ=255)
                                                                               gTCTTGGCTGGTGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCAAAATACCGGGCCGt          >  2:1052675/1‑93 (MQ=255)
                                                                                tCTTGGCTGGTGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCAAAATACCGGGCCGTg         <  2:7471775/93‑1 (MQ=255)
                                                                                tCTTGGCTGGTGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCAAAATACCGGGCCGTg         <  2:5358831/93‑1 (MQ=255)
                                                                                 cTTGGCTGGTGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCAAAATACCGGGCCGTGa        <  1:5489631/93‑1 (MQ=255)
                                                                                 cTTGGCTGGTGAACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCAAAATACCGGGGCGTGa        >  2:505109/1‑93 (MQ=255)
                                                                                   tggctggtGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGCGCTTTCACCAAAAATAACGGGCCTTGATc      >  1:2346185/1‑93 (MQ=255)
                                                                                     gctggtGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCAAAATACCGGGCCGTGATCCa    >  1:7745551/1‑93 (MQ=255)
                                                                                       tggtGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCAAAATACCGGGCCGTGATCCAgc  <  2:4628318/93‑1 (MQ=255)
                                                                                          |                                                                                         
AAGGCTTTATTCACGCCAATGCCGACAAACACCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGCGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCAAAATACCGGGCCGTGATCCAGC  >  NC_000913/177296‑177475

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