Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,863,429 T→C M158V (ATG→GTG)  yidE ← putative transport protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,863,4290TC97.7% 554.3 / ‑0.6 175M158V (ATG→GTG) yidEputative transport protein
Reads supporting (aligned to +/- strand):  ref base T (3/1);  new base C (78/93);  total (81/94)
Fisher's exact test for biased strand distribution p-value = 3.37e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 181 total aligned reads displayed.

TCTCGACATTGACGCGGAAAATAACCCGCAACATCCACATGGTGAACAAAATCCCGCAAATGCCGAATGGATACGCCATCGCGTAACTCATCCCCATCTGATCGACCATTTCCATTGGTGTACCCAGGTCGCGCAAAATCTGCTGCCCTGCCCCCAGCGCTGGCGTATTGGTAACCGCAC  >  NC_000913/3863339‑3863518
                                                                                          |                                                                                         
tctcGACATTGACGCGGAAAATAACCCGCAACATCCACATGGTGAACAAAATCCCGCAAATGCCGAATGGATACGCCATCGCGTAACTCAccc                                                                                         <  1:3669527/93‑1 (MQ=255)
   cGACATTGACGCGGAAAATAACCCGCAACATCCACATGGTGAACAAAATCCCGCAAATGCCGAATGGATACGCCATCGCGTAACTCACCCCCa                                                                                      <  2:2858975/93‑1 (MQ=255)
   cGACATTGACGCGGAAAATAACCCGCAACATCCACATGGTGAACAAAATCCCGCAAATGCCGAATGGATACGCCATCGCGTAACTCACCCCCa                                                                                      <  2:2545365/93‑1 (MQ=255)
    gACATTGCCGCGGAAAATACCCCGCAACATCCACATGGTGAACAAAATCCCGCAAATGCCGAATGGATACGCATCGCGTAACTCACCCCCATc                                                                                    <  1:3326960/93‑1 (MQ=255)
    gACATTGACGCGGAAAATAACCCGCAACATCCACATGGTGAACAAAATCCCGCAAATGCCGAATGGATACGCCATCGCGTAACTCACCCCCAt                                                                                     >  1:3619380/1‑93 (MQ=255)
     gcATTGACGCGGAAAATAACCCGCAACATCCACATGGTGAACAAAATCCCGCAAATGCCGAATGGATACGCCATCGCGTAACTCACCCCCATc                                                                                    <  1:1534747/92‑1 (MQ=255)
     aCATTGACGCGGAAAATAACCCGCAACATCCACATGGTGAACAAAATCCCGCAAATGCCGAATGGATACGCCATCGCGTAACTCACCCCCATc                                                                                    <  1:3324667/93‑1 (MQ=255)
     aCATTGACGCGGAAAATAACCCGCAACATCCACATGGTGAACAAAATCCCGCAAATGCCGAATGGATACGCCATCGCGTAACTCACCCCCATc                                                                                    >  2:3104860/1‑93 (MQ=255)
     aCATTGACGCGGAAAATAACACGCAACATCCACATGGTGAACAAAATCCCGCAAATGCCGAATGGATACGCCATCGCGTAACTCACCCCCATc                                                                                    <  1:3246372/93‑1 (MQ=255)
      cATTGACGCGGAAAATAACCCGCAACATCCACATGGTGAACAAAATCCCGCAAATGCCGAATGGATACGCCATCGCGTAACTCACCCCCATCt                                                                                   >  1:4121109/1‑93 (MQ=255)
      cATTGACGCGAAAAATAACCCCCAACATCCACATGGTGAACAAAATCCCGCAAATGCCGAATGGATACGCCATCGCGTAACTCACCCCCATCt                                                                                   <  2:4072963/93‑1 (MQ=255)
       aTTGACGCGGAAAATAACCCGCAACATCCACATGGTGAACAAAATCCCGCAAATGCCGAATGGATACGCCATCGCGTAACTCACCCCCATCTg                                                                                  >  2:2559584/1‑93 (MQ=255)
       aTTGACGCGGAAAATAACCCGCAACATCCACATGGTGAACAAAATCCCGCAAATGCCGAATGGATACGCCATCGCGTAACTCACCCCCATCTg                                                                                  >  1:1001747/1‑93 (MQ=255)
       aTTGACGCGGAAAATAACCCGCAACATCCACATGGTGAACAAAATCCCGCAAATGCCGAATGGATACGCCATCGCGTAACTCACCCCCATCTg                                                                                  >  1:1158863/1‑93 (MQ=255)
        ttGACGCGGAAAATAACCCGCAACATCCACATGGTGAACAAAATCCCGCAAATGCCGAATGGATACGCCATCGCGTAACTCACCCCCATCTGa                                                                                 >  2:2933544/1‑93 (MQ=255)
        ttGACGCGGAAAATAACCCGCAACATCCACATGGTGAACAAAATCCCGCAAATGCCGAATGGATACGCCATCGCGTAACTCACCCCCATCTGa                                                                                 <  1:3806908/93‑1 (MQ=255)
         tGACGCGGAAAATAACCCGCAACATCCACATGGTGAACAAAATCCCGCAAATGCCGAATGGATACGCCATCGCGTAACTCACCCCCATCTGAt                                                                                <  2:3550948/93‑1 (MQ=255)
         tGACGCGGAAAATAACCCGCAACATCCACATGGTGAACAAAATCCCGCAAATGCCGAATGGATACGCCATCGCGTAACTCACCCCCATCTGAt                                                                                <  2:1874042/93‑1 (MQ=255)
         tGACGCGGAAAATAACCGCAACATCCACATGGTGAACAAAATCCCGCAAATGCCGAATGGATACGCCATCGCGTAACTCACCCCCATCTGATc                                                                               <  1:3469361/93‑1 (MQ=255)
          gACGCGGAAAATAACCCGCAACATCCACATGGTGAACAAAATCCCGCAAATGCCGAATGGATACGCCATCGCGTAACTCACCCCCATCTGATc                                                                               <  1:2535377/93‑1 (MQ=255)
             gcgGAAAATAACCCGCAACATCCACATGGTGAACAAAATCCCGCAAATGCCGAATGGATACGCCATCGCGTAACTCACCCCCATCTGATCGAc                                                                            >  2:702840/1‑93 (MQ=255)
               ggAAAATAACCCGCAACATCCACATGGTGAACAAAATCCCGCAAATGCCGAATGGATACGCCATCGCGTAACTCACCCCCATCTGATCGACCa                                                                          <  2:3768261/93‑1 (MQ=255)
                gAAAATAACCCGCAACATCCACATGGTGAACAAAATCCCGCAAATGCCGAATGGATACGCCATCGCGTAACTCACCCCCATCTGATCGACCAt                                                                         >  1:1176820/1‑93 (MQ=255)
                gAAAATAACCCGCAACATCCACATGGTCAACAAAATCCCGCAAATGCCGAATGGATACGCCATCGCGTAACTCACCCCCATCTGATCGACCAt                                                                         <  1:2129623/93‑1 (MQ=255)
                 aaaaTAACCCGCAACATCCACATGGTGAACAAAATCCCGCAAATGCCGAATGGATACGCCATCGCGTAACTCACCCCCATCTGATCGACCAtt                                                                        <  2:432944/93‑1 (MQ=255)
                    aTAACCCGCAACATCCACATGGTGAACAAAATCCCGCAAATGCCGAATGGATACGCCATCGCGTAACTCATCCCCATCTGATCGACCATTTcc                                                                     >  1:1996712/1‑93 (MQ=255)
                    aTAACCCGCAACATCCACATGGTGAACAAAATCCCGCAAATGCCGAATGGATACGCCATCGCGTAACTCACCCCCATCTGATCGACCATTTcc                                                                     >  1:2436300/1‑93 (MQ=255)
                     tAACCCGCAACATCCACATGGTGAACAAAATCCCGCAAATGCCGAATGGATACGCCATCGCGTAACTCACCCCCATCTGATCGACCATTTCCa                                                                    >  2:4015535/1‑93 (MQ=255)
                     tAACCCGCAACATCCACATGGTGAACAAAATCCCGCAAATGCCGAATGGATACGCCATCGCGTAACTCACCCCCATCTGATCGACCATTTCCa                                                                    >  2:22704/1‑93 (MQ=255)
                      aaCCCGCAACATCCACATGGTGAACAAAATCCCGCAAATGCCGAATGGATACGCCATCGCGTAACTCACCCCCATCTGATCGACCATTTCCAt                                                                   >  1:1274868/1‑93 (MQ=255)
                       aCCCGCAACATCCACATGGTGAACAAAATCCCGCAAATGCCGAATGGATACGCCATCGCGTAACTCACCCCCATCTGATCGACCATTTCCAtt                                                                  <  1:2896903/93‑1 (MQ=255)
                       aCCCGCAACATCCACATGGTGAACAAAATCCCGCAAATGCCGAATGGATACGCCATCGCGTAACTCACCCCCATCTGATCGACCATTTCCAtt                                                                  <  1:247117/93‑1 (MQ=255)
                       aCCCGCAACATCCACATGGTGAACAAAATCCCGCAAATGCCGAATGGATACGCCATCGCGTAACTCACCCCCATCTGATCGACCATTTCCAtt                                                                  >  2:2106088/1‑93 (MQ=255)
                         ccGCAACATCCACATGGTGAACAAAATCCCGCAAATGCCGAATGGATACGCCATCGCGTAACTCACCCCCATCTGTTCGACCATTTCCATTgg                                                                >  1:4215117/1‑93 (MQ=255)
                         ccGCAACATCCACATGGTGAACAAAATCCCGCAAATGCCGAATGGATACGCCATCGCGTAACTCACCCCCATCTGATCGACCATTTCCATTgg                                                                <  1:1675739/93‑1 (MQ=255)
                         ccGCAACATCCACATGGTGAACAAAATCCCGCAAATGCCGAATGGATACGCCATCGCGTAACTCACCCCCATCTGATCGACCATTTCCATTgg                                                                >  1:3872037/1‑93 (MQ=255)
                          cGCAACATCCACATGGTGAACAAAATCCCGCAAATGCCGAATGGATACGCCATCGCGTAACTCACCCCCATCTGATCGACCATTTCCATTGgt                                                               >  2:417380/1‑93 (MQ=255)
                          cGCAACATCCACATGGTGAACAAAATCCCGCAAATGCCGAATGGATACGCCATCGCGTAACTCACCCCCATCTGATCGACCATTTCCATTGgt                                                               >  1:2080070/1‑93 (MQ=255)
                          cGCAACATCCACATGGTGAACAAAATCCCGCAAATGCCGAATGGATACGCCATCGCGTAACTCACCCCCATCTGATCGACCATTTCCATTGgt                                                               >  1:100302/1‑93 (MQ=255)
                            cAACATCCACATGGTGAACAAAATCCCTCAAATGCCGAATGGATACGCCATCGCGTAACTCACCCCCATCTGATCGACCATTTCCATTGgtgt                                                             <  2:1986093/93‑1 (MQ=255)
                            cAACATCCACATGGTGAACAAAATCCCGCAAATGCCGAATGGATACGCCATCGCGTAACTCACCCCCATCTGATCGACCATTTCCATTGgtgt                                                             <  2:3298651/93‑1 (MQ=255)
                            cAACATCCACATGGTGAACAAAATCCCGCAAATGCCGAATGGATACGCCATCGCGTAACTCACCCCCATCTGATCGACCATTTCCATTGgtgt                                                             >  1:3687305/1‑93 (MQ=255)
                              aCATCCACATGGTGAACAAAATCCCGCAAATGCCGAATGGATACGCCATCGCGTAACTCACCCCCATCTGATCGACCATTTCCATTGGTGTAc                                                           >  2:4666820/1‑93 (MQ=255)
                              aCATCCACATGGTGAACAAAATCCCGCAAATGCCGAATGGATACGCCATCGCGTAACTCACCCCCATCTGATCGACCATTTCCATTGGTGTAc                                                           <  1:3943517/93‑1 (MQ=255)
                                aTCCACATGGTGAACAAAATCCCGCAAATGCCGAATGGATACGCCATCGCGTAACTCACCCCCATCTGATCGACCATTTCCATTGGTGTAccc                                                         <  2:173173/93‑1 (MQ=255)
                                aTCCACATGGTGAACAAAATCCCGCAAATGCCGAATGGATACGCCATCGCGTAACTCACCCCCATCTGATCGACCATTTCCATTGGTGTAccc                                                         <  2:173172/93‑1 (MQ=255)
                                aTCCACATGGTGAACAAAATCCCGCAAATGCCGAATGGATACGCCATCGCGTAACTCACCCCCATCTGATCGACCATTTCCATTGGTGTAccc                                                         <  1:4087080/93‑1 (MQ=255)
                                 tCCACATGGTGAACAAAATCCCGCAAATGCCGAATGGATACGCCATCGCGTAACTCACCCCCATCTGATCGACCATTTCCATTGGTGTACCCa                                                        >  2:1080816/1‑93 (MQ=255)
                                 tCCACATGGTGAACAAAATCCCGCAAATGCCGAATGGATACGCCATCGCGTAACTCACCCCCATCTGATCGACCATTTCCATTGGTGTACCCa                                                        >  2:420149/1‑93 (MQ=255)
                                 tCCACATGGTGAACAAAATCCCGCAAATGCCGAATGGATACGCCATCGCGTAACTCACCCCCATCTGATCGACCATTTCCATTGGTGTACCCa                                                        <  1:2282853/93‑1 (MQ=255)
                                  ccACATGGTGAACAAAATCCCGCAAATGCCGAATGGATACGCCATCGCGTAACTCGCCCCCATCTGATCGACCATTTCCATTGGTGTACCCAg                                                       <  2:373807/93‑1 (MQ=255)
                                  ccACATGGTGAACAAAATCCCGCAAATGCCGAATGGATACGCCATCGCGTAACTCACCCCCATCTGATCGACCATTTCCATTGGTGTACCCAg                                                       <  2:3507077/93‑1 (MQ=255)
                                  ccACATGGTGAACAAAATCCCGCAAATGCCGAATGGATACGCCATCGCGTAACTCACCCCCATCTGATCGACCATTTCCATTGGTGTACCCAg                                                       <  2:2469486/93‑1 (MQ=255)
                                    acaTGGTGAACAAAATCCCGCAAATGCCGAATGGATACGCCATCGCGTAACTCACCCCCATCTGATCGACCATTTCCATTGGTGTACCCAGGt                                                     >  1:542055/1‑93 (MQ=255)
                                    acaTGGTGAACAAAATCCCGCAAATGCCGAATGGATACGCCATCGCGTAACTCACCCCCATCTGATCGACCATTTCCATTGGTGTACCCAGGt                                                     >  1:969212/1‑93 (MQ=255)
                                       tGGTGAACAAAATCCCGCAAATGCCGAATGGATACGCCATCGCGTAACTCACCCCCATCTGATCGACCATTTCCATTGGTGTACCCAGGTcgc                                                  >  1:1732367/1‑93 (MQ=255)
                                       tGGTGAACAAAATCCCGCAAATGCCGAATGGATACGCCATCGCGTAACTCACCCCCATCTGATCGACCATTTCCATTGGTGTACCCAGGTcgc                                                  >  2:3084098/1‑93 (MQ=255)
                                          tGAACAAAATCCCGCAAATGCCGAATGGATACGCCATCGCGTAACTCACCCCCATCTGATTGACCATTTCCATTGGTGTACCCAGGTCGCGCa                                               >  2:1579167/1‑93 (MQ=255)
                                            aaCAAAATCCCGCAAATGCCGAATGGATACGCCATCGCGTAACTCACCCCCATCTGATCGACCATTTCCATTGGTGTACCCAGGTCGCGCaaa                                             >  1:4134246/1‑93 (MQ=255)
                                              cAAAATCCCGCAAATGCCGAATGGATACGCCATCGCGTAACTCACCCCCATCTGATCGACCATTTCCATTGGTGTACCCAGGTCGCGCAAAAt                                           <  1:1382928/93‑1 (MQ=255)
                                                aaaCCCCGCAAATGCCGAATGGATACGCCATCGCGTAACTCACCCCCATCTGATCGACCATTTCCATTGGTGTACCCAGGTCGCGCAAAATct                                         <  2:3684507/93‑1 (MQ=255)
                                                  aTCCCGCAAATGCCGAATGGATACGCCATCGCGTAACTCACCCCCATCTGATCGACCATTTCCATTGGTGTACCCAGGTCGCGCAAAATctgc                                       <  2:1435597/93‑1 (MQ=255)
                                                  aTCCCGCAAATGCCGAATGGATACACCATCGCGTAACTCACCCCCATCTGATCGACCATTTCCATTGGTGTACCCAGGTCGCGCAAAATctgc                                       <  2:2398384/93‑1 (MQ=255)
                                                   tCCCGCAAATGCCGAATGGATACGCCATCGCGTAACTCACCCCCATCTGATCGACCATTTCCATTGGTGTACCCAGGTCGCGCAAAATctgct                                      >  2:4394157/1‑93 (MQ=255)
                                                   tCCCGCAAATGCCGAATGGATACGCCATCGCGTAACTCACCCCCATCTGATCGACCATTTCCATTGGTGTACCCAGGTCGCGCAAAATctgct                                      <  1:2014638/93‑1 (MQ=255)
                                                   tCCCGCAAATGCCGAATGGATACGCCATCGCGTAACTCACCCCCATCTGATCGACCATTTCCATTGGTGTACCCAGGTCGCGCAAAATctgct                                      <  1:246312/93‑1 (MQ=255)
                                                    cccGCAAATGCCGAATGGATACGCCATCGCGTAACTCACCCCCATCTGATCGACCATTTCCATTGGTGTACCCAGGTCGCGCAAAATctgctg                                     >  2:842765/1‑93 (MQ=255)
                                                     ccGCAAATGCCGAATGGATACGCCATCGCGTAACTCACCCCCATCTGATCGACCATTTCCATTGGTGTACCCAGGTCGCGCAAAATCtgctgc                                    <  1:4152676/93‑1 (MQ=255)
                                                      cGCAAATGCCGAATGGATACGCCATCGCGTAACTCACCCCCATCTGATCGACCATTTCCATTGGTGTACCCAGGTCGCGCAAAATCTGctgcc                                   <  2:3170945/93‑1 (MQ=255)
                                                      cGCAAATGCCGAATGGATACGCCATCGCGTAACTCACCCCCATCTGATCGACCATTACCATTGGTGTACCCAGGTCGCGCAAAATCTGctgcc                                   <  2:3170944/93‑1 (MQ=255)
                                                       gCAAATGCCGAATGGATACGCCATCGCGTAACTCACCCCCATCTGATCGACCATTTCCATTGGTGTACCCAGGTCGCGCAAAATCTGctgccc                                  >  1:3071415/1‑93 (MQ=255)
                                                       gCAAATGCCGAATGGATACGCCATCGCGTAACTCACCCCCATCTGATCGACCATTTCCATTGGTGTACCCAGGTCGCGCAAAATCTGctgccc                                  <  1:60341/93‑1 (MQ=255)
                                                       gCAAATGCCGAATGGATACGCCATCGCGTAACTCACCCCCATCTGATCGACCATTTCCATTGGTGTACCCAGGTCGCGCAAAATCTGctgccc                                  >  2:4738884/1‑93 (MQ=255)
                                                       gCAAATGCCGAATGGATACGCCATCGCGTAACTCACCCCCATCTGATCGACCATTTCCATTGGTGTACCCAGGTCGCGCAAAATCTGctgccc                                  >  1:420346/1‑93 (MQ=255)
                                                         aaaTGCCGAATGGATACGCCATCGCGTAACTCACCCCCATCTGATCGACCATTTCCATTGGTGTACCCATGTCGCGCAAAATCTGctgccctg                                >  2:2617231/1‑93 (MQ=255)
                                                          aaTGCCGAATGGATACGCCATCGCGTAACTCACCCCCATCTGATCGACCATTTCCATTGGTGTACCCAGGTCGCGCAAAATCTGctgccctgc                               >  2:2323693/1‑93 (MQ=255)
                                                             gCCGAATGGATACGCCATCGCGTAACTCACCCCCATCTGATCGACCATTTCCATTGGTGTACCCAGGTCGCGCAAAATCTGCTGCCCTGcccc                            <  1:2158404/93‑1 (MQ=255)
                                                              ccGAATGGATACGCCATCGCGTAACTCACCCCCATCTGATCGACCATTTCCATTGGTGTACCCAGGTCGCGCAAAATCTGCTGCCCTGccccc                           <  1:4439843/93‑1 (MQ=255)
                                                              ccGAATGGATACGCCATCGCGTAACTCACCCCCATCTGATCGACCATTTCCATTGGTGTACCCAGGTCGCGCAAAATCTGCTGCCCTGccccc                           <  1:1860281/93‑1 (MQ=255)
                                                               cGAATGGATACGCCATCGCGTAACTCACCCCCATCTGATCGACCATTTCCATTGGTGTACCCAGGTCGCGCAAAATCTGCTGCCCTGCCCCCa                          <  2:3061471/93‑1 (MQ=255)
                                                                gAATGGATACGCCATCGCGTAACTCACCCCCATCTGATCGACCATTTCCATTGGTGTACCCAGGTCGCGCAAAATCTGCTGCCCTGCCCCCAg                         >  1:891540/1‑93 (MQ=255)
                                                                gAATGGATACGCCATCGCGTAACTCACCCCCATCTGATCGACCATTTCCATTGGTGTACCCAGGTCGCGCAAAATCTGCTGCCCTGCCCCCAg                         >  1:1826530/1‑93 (MQ=255)
                                                                gAATGGATACGCCATCGAGTAACTCACCCCCATCTGATCGACCATTTCCATTGGTGTACCCAGGTCGCGCAAAATCTGCTGCCCTGCCCCCAg                         >  2:1329156/1‑93 (MQ=255)
                                                                  aTGGATACGCCATCGCGTAACTCACCCCCCTCTGATCGACCATTTCCATTGGTGTACCCAGGTCGCGCAAAATCTGCTGCCCTGCCCCCAgcg                       >  2:2682436/1‑93 (MQ=255)
                                                                  aTGGATACGCCATCGCGTAACTCACCCCCATCTGATCGACCATTTCCATTGGTGTACCCAGGTCGCGCAAAATCTGCTGCCCTGCCCCCAgcg                       >  1:3130902/1‑93 (MQ=255)
                                                                   tGGATACGCCATCGCGTAACTCACCCCCATCTGATCGACCATTTCCATTGGTGTACCCAGGTCGCGCAAAATCTGCTGCCCTGCCCCCAgcgc                      >  1:2691210/1‑93 (MQ=255)
                                                                   tGGATACGCCATCGCGTAACTCACCCCCATCTGATCGACCATTTCCATTGGTGTACCCAGGTCGCGCAAAATCTGCTGCCCTGCCCCCAgcgc                      <  1:4279367/93‑1 (MQ=255)
                                                                      aTACGCCATCGCGTAACTCACCCCCATCTGATCGACCATTTCCATTGGTGTACCCAGGTCGCGCAAAATCTGCTGCCCTGCCCCCAGCGCTgg                   <  2:1921635/93‑1 (MQ=255)
                                                                       tACGCCATCGCGTAACTCACCCCCATCTGATCGACCATTTCCATTGGTGTACCCAGGTCGCGCAAAATCTGCTGCCCTGCCCCCAGCGCTGGc                  >  2:954437/1‑93 (MQ=255)
                                                                         cGCCATCGCGTAACTCCCCCCCATCTGATCGACCATTTCCATTGGTGTACCCAGGTCGCGCAAAATCTGCTGCCCTGCCCCCAGCGCTGGCGt                <  2:2520131/93‑1 (MQ=255)
                                                                          gCCATCGCGTAACTCACCCCCATCTGATCGACCATGTCCATTGGTGTACCCAGGTCGCGCAAAATCTGCTGCCCTGCCCCCAGCGCTGGCGTa               >  2:4675863/1‑93 (MQ=255)
                                                                           ccATCGCGTAACTCACCCCCATCTGATCGACCATTTCCATTGGTGTACCCAGGTCGCGCAAAATCTGCTGCCCTGCCCCCAGCGCTGGCGTAt              >  1:2993932/1‑93 (MQ=255)
                                                                             aTCGCGTAACTCACCCCCATCTGATCGACCATTTCCATTGGTGTACCCAGGTCGCGCAAAATCTGCTGCCCTGCCCCCAGCGCTGGCGTATTg            >  2:2772652/1‑93 (MQ=255)
                                                                                  gTAACTCACCCCCATCTGATCGACCATTTCCATTGGTGTACCCAGGTCGCGCAAAATCTGCTGCCCTGCCCCCAGCGCTGGCGTATTGGTAAc       <  1:2954659/93‑1 (MQ=255)
                                                                                  gTAACTCACCCCCATCTGATCGACCATTTCCATTGGTGTACCCAGGTCGCGCAAAATCTGCTGCCCTGCCCCCAGCGCTGGCGTATTGGTAAc       <  1:1492007/93‑1 (MQ=255)
                                                                                   tAACTCACCCCCATCTGATCGACCATTTCCATTGGTGTACCCAGGTCGCGCAAAATCTGCTGCCCTGCCCCCAGCGCTGGCGTATTGGTAAcc      >  2:3937320/1‑93 (MQ=255)
                                                                                   tAACTCACCCCCATCTGATCGACCATTTCCATTGGTGTACCCAGGTCGCGCAAAATCTGCTGCCCTGCCCCCAGCGCTGGCGTATTGGTAAcc      <  1:2111877/93‑1 (MQ=255)
                                                                                    aaCTCACCCCCATCTGATCGACCATTTCCATTGGTGTACCCAGGTCGCGCAAAATCTGCTGCCCTGCCCCCAGCGCTGGCGTATTGGTAACCg     <  2:1467015/93‑1 (MQ=255)
                                                                                     aCTCACCCCCATCTGATCGACCATTTCCATTGGTGTACCCAGGTCGCGCAAAATCTGCTGCCCTTCCCCCAGCGCTGGCGTATTGGTAACCGc    >  2:4332058/1‑93 (MQ=255)
                                                                                       tCACCCCCATCTGATCGACCATTTCCATTGGTGTACCCAGGTCGCGCAAAATCTGCTGCCCTGCCCCCAGCGCTGGCGTATTGGTAACCGCAc  <  1:2060675/93‑1 (MQ=255)
                                                                                          |                                                                                         
TCTCGACATTGACGCGGAAAATAACCCGCAACATCCACATGGTGAACAAAATCCCGCAAATGCCGAATGGATACGCCATCGCGTAACTCATCCCCATCTGATCGACCATTTCCATTGGTGTACCCAGGTCGCGCAAAATCTGCTGCCCTGCCCCCAGCGCTGGCGTATTGGTAACCGCAC  >  NC_000913/3863339‑3863518

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: