Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,666,733 T→C H29R (CAC→CGC)  yhjA ← putative cytochrome C peroxidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,666,7330TC85.6% 698.0 / 98.3 271H29R (CAC→CGC) yhjAputative cytochrome C peroxidase
Reads supporting (aligned to +/- strand):  ref base T (22/17);  new base C (136/96);  total (158/113)
Fisher's exact test for biased strand distribution p-value = 8.61e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.04e-01

Only 100 of 281 total aligned reads displayed.

CGCATCCTTTTTCGCGGAGAAAGCCTAAAACCTTATTATTTTCACTGACAGCAGATGCCTGAACATCGGCCTGTTTACTGCGTTTATTATCGTGGTACCACACATAACCAGATAACCCTAAATAGCAAATCGCGACGCCAGCCAGGCCGATCGCGGTAATACGTGAGACCATTTTCATTATTCT  >  NC_000913/3666641‑3666824
                                                                                            |                                                                                           
 gCATCCTTTTTCGCGGAGAAAGCCTAAAACCTTATTATTTTCGCTGACAGCAGATGCCTGAACATCGGCCTGTTTACTGCGTTTNTTATCGTg                                                                                            >  1:5223191/1‑93 (MQ=255)
  cATCCTTTTTCGCGGAGAAAGCCTAAAACCTTATTATTTTCACTGACAGCAGATGCCTGAACATCGGCCTGTTTACTGCGTTTATTATCGCgg                                                                                           >  2:763216/1‑93 (MQ=255)
  cATCCTTTTTCGCGGAGAAAGCCTAAAACCTTATTATTTTCACTGACAGCAGATGCCTGAACATCGGCCTGTTTACTGCGTTTATTATCGCgg                                                                                           >  1:579402/1‑93 (MQ=255)
   aTCCTTTTTCGCGGAGAAAGCCTAAAACCTTATTATTTTCACTGACAGCAGATGCCTGAACATCGGCCTGTTTACTGCGTTTATTATCGTGGt                                                                                          >  2:2884373/1‑93 (MQ=255)
    tCCTTTTTCGCGGAGAAAGCCTAAAACCTTATTATTTTCACTGACAGCAGATGCCTGAACATCGGCCTGTTTACTGCGTTTATTATCGCGGTa                                                                                         >  2:4053298/1‑93 (MQ=255)
    tCCTTTTTCGCGGAGAAAGCCTAAAACCTTATTATTTTCACTGACAGCAGATGCCTGAACATCGGCCTGTTTACTGCGTTTATTATCGAGGTa                                                                                         >  2:4501597/1‑93 (MQ=255)
      cTTTTTCGCGGAGAAAGCCTAAAACCTTATTATTTTCACTGACAGCAGATGCCTGAACATCGGCCTGTTTACTGCGTTTATTATCGCGGTAcc                                                                                       >  2:5482203/1‑93 (MQ=255)
         tttCGCGGAGAAAGCCTAAAACCTTATTATTTTCACTGACAGCAGATGCCTGAACATCGGCCTGTTTACTGCGTTTATTATCGCGGTACcaca                                                                                    <  1:2892651/93‑1 (MQ=255)
           tCGCGGAGAAAGCCTAAAACCTTATTATTTTCACTGACAGCAGATGCCTGAACATCGGCCTGTTTACTGCGTTTATTATCGTGGTACcacaca                                                                                  >  2:3715930/1‑93 (MQ=255)
           tCGCGGAGAAAGCCTAAAACCTTATTATTTTCACTGACAGCAGATGCCTGAACATCGGCCTGTTTACTGCGTTTATTATCGCGGTACcacaca                                                                                  >  1:1902906/1‑93 (MQ=255)
            cgcgGAGAAAGCCTAAAACCTTATTATTTTCACTGACAGCAGATGCCTGAACATCGGCCTGTTTACTGCGTTTATTATCGTGGTACCACACAt                                                                                 <  2:2148360/93‑1 (MQ=255)
            cgcgGAGAAAGCCTAAAACCTTATTATTTTCACTGACAGCAGATGCCTGAACATCGGCCTGTTTACTGCGTTTATTATCGTGGTACCACACAt                                                                                 <  2:1276906/93‑1 (MQ=255)
              cgGAGAAAGCCTAAAACCTTATTATTTTCACTGACAGCAGATGCCTGAACATCGGCCTGTTTACTGCGTTTATTATCGTGGTACCACACATaa                                                                               >  1:1521722/1‑93 (MQ=255)
               ggAGAAAGCCTAAAACCTTATTATTTTTACTGACAGCAGATGCCTGAACATCGGCCTGTTTACTGCGTTTATTATCGCGGTACCACACATAAc                                                                              >  1:949392/1‑93 (MQ=255)
               ggAGAAAGCCTAAAACCTTATTATTTTCACTGACAGCAGATGCCTGAACATCGGCCTGTTTACTGCGTTTATTATCGCGGTACCACACATAAc                                                                              >  1:2928490/1‑93 (MQ=255)
                gagaAAGCCTAAAACCTTATTATTTTCACTGACAGCAGATGCCTGAACATCGGCCTGTTTACTGCGTTTATTATCGTGGTACCACACATAAcc                                                                             >  1:5097209/1‑93 (MQ=255)
                gagaAAGCCTAAAACCTTATTATTTTCACTGACAGCAGATGCCTGAACATCGGCCTGTTTACTGCGTTTATTATCGCGGTACCACACATAAcc                                                                             >  1:5299152/1‑93 (MQ=255)
                gagaAAGCCTAAAACCTTATTATTTTCACTGACAGCAGATGCCTGAACATCGGCCTGTTTACTGCGTTTATTATCGCGGTACCACACATAAcc                                                                             >  2:302423/1‑93 (MQ=255)
                 agaAAGCCTAAAACCTTATTATTTTCACTGACAGCAGATGCCTGAACATCGGCCTGTTTACTGCGTTTATTATCGCGGTACCACACATAACCa                                                                            <  1:3240901/93‑1 (MQ=255)
                 agaAAGCCTAAAACCTTATTATTTTCACTGACAGCAGATGCCTGAACATCGGCCTGTTTACTGCGTTTATTATCGCGGTACCACACATAACCa                                                                            >  1:2512630/1‑93 (MQ=255)
                 agaAAGCCTAAAACCTTATTATTTTCACTGACAGCAGATGCCTGAACATCGGCCTGTTTACTGCGTTTATTATCGCGGTACCACACATAACCa                                                                            >  1:2217397/1‑93 (MQ=255)
                  gaAAGCCTAAAACCTTATTATTTTCACTGACAGCAGATGCCTGAACATCGGCCTGTTTACTGCGTTTATTATCGCGGTACCACACATAACCAg                                                                           <  2:6180010/93‑1 (MQ=255)
                  gaAAGCCTAAAACCTTATTATTTTCACTGACAGCAGATGCCTGAACATCGGCCTGTTTACTGCGTTTATTATCGCGGTACCACACATAACCAg                                                                           >  2:3730876/1‑93 (MQ=255)
                   aaaGCCTAAAACCTTATTATTTTCACTGACAGCAGATGCCTGAACATCGGCCTGTTTACTGCGTTTATTATCGCGGTACCACACATAACCAGa                                                                          <  1:3137381/93‑1 (MQ=255)
                   aaaGCCTAAAACCTTATTATTTTCACTGACAGCAGATGCCTGAACATCGGCCTGTTTACTGCGTTTATTATCGCGGTACCACACATAACCAGa                                                                          >  1:3758560/1‑93 (MQ=255)
                    aaGCCTAAAACCTTATTATTTTCACTGACAGCAGATGCCTGAACATCGGCCTGTTTACTGCGTTTATTATCGTGGTACCACACATAACCAGAt                                                                         >  1:812043/1‑93 (MQ=255)
                    aaGCCTAAAACCTTATTATTTTCACTGACAGCAGATGCCTGAACATCGGCCTGTTTACTGCGTTTATTATCGCGGTACCACACATAACCAGAt                                                                         >  1:5687770/1‑93 (MQ=255)
                    aaGCCTAAAACCTTATTATTTTCACTGACAGCAGATGCCTGAACATCGGCCTGTTTACTGCGTTTATTATCGCGGTACCACACATAACCAGAt                                                                         >  2:3309313/1‑93 (MQ=255)
                     aGCCTAAAACCTTATTATTTTCACTGACAGCAGATGCCTGAACATCGGCCTGTTTACTGCGTTTATTATCGTGGTACCACACATAACCAGATa                                                                        >  2:3235296/1‑93 (MQ=255)
                       ccTAAAACCTTATTATTTTCACTGACAGCAGATGCCTGAACATCGGCCTGTTTACTGCGTTTATTATCGCGGTACCACACATAACCAGATAAc                                                                      >  1:20452/1‑93 (MQ=255)
                         tAAAACCTTATTATTTTCACTGACAGCAGATGCCTGAACATCGGCCTGTTTACTGCGTTTATTATCGTGGTACCACACATAACCAGATAAccc                                                                    <  2:512317/93‑1 (MQ=255)
                         tAAAACCTTATTATTTTCACTGACAGCAGATGCCTGAACATCGGCCTGTTTACTGCGTTTATTATCGCGGTACCACACATAACCAGATAAccc                                                                    >  2:3485599/1‑93 (MQ=255)
                         tAAAACCTTATTATTTTCACTGACAGCAGATGCCTGAACATCGGCCTGTTTACTGCGTTTATTATCGCGGTACCACACATAACCAGATAAccc                                                                    <  1:3587891/93‑1 (MQ=255)
                           aaaCCTTATTATTTTCACTGACAGCAGATGCCTGAACATCGGCCTGTTTACTGCGTTTATTATCGCGGTACCACACATAACCAGATAACCCTa                                                                  >  2:1564427/1‑93 (MQ=255)
                            aaCCTTATTATTTTCACTGACAGCAGATGCCTGAACATCGGCCTGTTTACTGCGTTTATTATCGTGGTACCACACATAACCAGATAACCCTaa                                                                 >  1:1956215/1‑93 (MQ=255)
                             aCCTTATTATTTTCACTGACAGCAGATGCCTGAACATCGGCCTGTTTACTGCGTTTATTATCGCGGTACCACACATAACCAGATAACCCTaaa                                                                >  1:1191629/1‑93 (MQ=255)
                              ccTTATTATTTTCACTGACAGCAGATGCCTGAACATCGGCCTGTTTACTGCGTTTATTATCGCGGTACCACACATAACCAGATAACCCTAAAt                                                               >  2:5209551/1‑93 (MQ=255)
                               cTTATTATTTTCACTGACAGCAGATGCCTGAACATCGGCCTGTTTACTGCGTTTATTATCGCGGTACCACACATAACCAGATAACCCTAAATa                                                              >  1:495491/1‑93 (MQ=255)
                                ttattaTTTTCACTGACAGCAGATGCCTGAACATCGGCCTGTTTACTGCGTTTATTATCGCGGTACCACACATAACCAGATAACCCTAAATAg                                                             >  2:1247383/1‑93 (MQ=255)
                                 tattatTTTCACTGACAGCAGATGCCTGAACATCGGCCTGTTTACTGCGTTTATTATCGCGGTACCACACATAACCAGATAACCCTAAATAGc                                                            >  2:6412193/1‑93 (MQ=255)
                                  attattTTCACTGACAGCAGATGCCTGAACATCGGCCTGTTTACTGCGTTTATTATCGTGGTACCACACATAACCAGATAACCCTAAATAGCa                                                           >  1:3995437/1‑93 (MQ=255)
                                  attattTTCACTGACAGCAGATGCCTGAACATCGGCCTGTTTACTGCGTTTATTATCGCGGTACCACACATAACCAGATAACCCTAAATAGCa                                                           <  1:6389062/93‑1 (MQ=255)
                                      ttttCACTGACAGCAGATGCCTGAACATCGGCCTGTTTACTGCGTTTATTATCGTGGTACCACACATAACCAGATAACCCTAAATAGCAAATc                                                       >  1:2747001/1‑93 (MQ=255)
                                          cACTGACAGCAGATGCCTGAACATCGGCCTGTTTACTGCGTTTATTATCGTGGTACCACACATAACCAGATAACCCTAAATAGCAAATCGCGa                                                   <  2:2402413/93‑1 (MQ=255)
                                          cACTGACAGCAGATGCCTGAACATCGGCCTGTTTACTGCGTTTATTATCGCGGTACCACACATAACCAGATAACCCTAAATAGCAAATCGCGa                                                   <  2:218862/93‑1 (MQ=255)
                                             tGACAGCAGATGCCTGAACATCGGCCTGTTTACTGCGTTTATTATCGCGGTACCACACATAACCAGATAACCCTAAATAGCAAATCGCGACgc                                                <  2:3081550/93‑1 (MQ=255)
                                             tGACAGCAGATGCCTGAACATCGGCCTGTTTACTGCGTTTATTATCGCGGTACCACACATAACCAGATAACCCTAAATAGCAAATCGCGACgc                                                >  2:1893654/1‑93 (MQ=255)
                                                 agcagATGCCTGAACATCGGCCTGTTTACTGCGTTTATTATCGCGGTACCACACATAACCAGATAACCCTAAATAGCAAATCGCGACgccagc                                            <  2:920245/93‑1 (MQ=255)
                                                  gcagATGCCTGAACATCGGCCTGTTTACTGCGTTTATTATCGCGGTACCACACATAACCAGATAACCCTAAATAGCAAATCGCGACgccagcc                                           <  1:1835875/93‑1 (MQ=255)
                                                  gcagATGCCTGAACATCGGCCTGTTTACTGCGTTTATTATCGCGGTACCACACATAACCAGATAACCCTAAATAGCAAATCGCGACgccagcc                                           >  2:849879/1‑93 (MQ=255)
                                                  gcagATGCCTGAACATCGGCCTGTTTACTGCGTTTATTATCGCGGTACCACACATAACCAGATAACCCTAAATAGCAAATCGCGACgccagcc                                           <  1:4730523/93‑1 (MQ=255)
                                                   cagATGCCTGAACATCGGCCTGTTTACTGCGTTTATTATCGTGGTACCACACATAACCAGATAACCCTAAATAGCAAATCGCGACgccagcca                                          <  1:5058321/93‑1 (MQ=255)
                                                   cagATGCCTGAACATCGGCCTGTTTACTGCGTTTATTATCGCGGTACCACACATAACCAGATAACCCTAAATAGCAAATCGCGACgccagcca                                          <  1:6233631/93‑1 (MQ=255)
                                                    agATGCCTGAACATCGGCCTGTTTACTGCGTTTATTATCGCGGTACCACACATAACCAGATAACCCTAAATAGCAAATCGCGACGccagccag                                         <  2:5348870/93‑1 (MQ=255)
                                                    agATGCCTGAACATCGGCCTGTTTACTGCGTTTATTATCGCGGTACCACACATAACCAGATAACCCTAAATAGCAAATCGCGACGccagccag                                         >  1:4227018/1‑93 (MQ=255)
                                                     gATGCCTGAACATCGGCCTGTTTACTGCGTTTATTATCGCGGTACCACACATAACCAGATAACCCTAAATAGCAAATCGCGACGCCAGCCAgg                                        <  2:744999/93‑1 (MQ=255)
                                                     gATGCCTGAACATCGGCCTGTTTACTGCGTTTATTATCGCGGTACCACACATAACCAGATAACCCTAAATAGCAAATCGCGACGCCAGCCAgg                                        >  2:1012791/1‑93 (MQ=255)
                                                     gATGCCTGAACATCGGCCTGTTTACTGCGTTTATTATCGCGGTACCACACATAACCAGATAACCCTAAATAGCAAATCGCGACGCCAGCCAgg                                        >  2:1436524/1‑93 (MQ=255)
                                                     gATGCCTGAACATCGGCCTGTTTACTGCGTTTATTATCGCGGTACCACACATAACCAGATAACCCTAAATAGCAAATCGCGACGCCAGCCAgg                                        >  1:669142/1‑93 (MQ=255)
                                                      aTGCCTGAACATCGGCCTGTTTACTGCGTTTATTATCGCGGTACCACACATAACCAGATAACCCTAAATAGCAAATCGCGACGCCAGCCAGGc                                       <  1:3666041/93‑1 (MQ=255)
                                                       tGCCTGAACATCGGCCTGTTTACTGCGTTTATTATCGCGGTACCACACATAACCAGATAACCCTAAATAGCAAATCGCGACGCCAGCCAGGcc                                      <  1:59914/93‑1 (MQ=255)
                                                         ccTGAACATCGGCCTGTTTACTGCGTTTATTATCGCGGTACCACACATAACCAGATAACCCTAAATAGCAAATCGCGACGCCAGCCAGGCCGa                                    >  1:107218/1‑93 (MQ=255)
                                                         ccTGAACATCGGCCTGTTTACTGCGTTTATTATCGCGGTACCACACATAACCAGATAACCCTAAATAGCAAATCGCGACGCCAGCCAGGCCGa                                    <  2:3616162/93‑1 (MQ=255)
                                                          cTGAACATCGGCCTGTTTACTGCGTTTATTATCGCGGTACCACACATAACCAGATATCCCTAAATAGCAAATCGCGACGCCAGCCAGGCCGAt                                   <  2:4194092/93‑1 (MQ=255)
                                                           tGAACATCGGCCTGTTTACTGCGTTTATTATCGCGGTACCACACATAACCAGATAACCCTAAATAGCAAATCGCGACGCCAGCCAGGCCGATc                                  >  1:1414742/1‑93 (MQ=255)
                                                           tGAACATCGGCCTGTTTACTGCGTTTATTATCGCGGTACCACACATAACCAGATAACCCTAAATAGCAAATCGCGACGCCAGCCAGGCCGATc                                  <  1:2352691/93‑1 (MQ=255)
                                                               cATCGGCCTGTTTACTGCGTTTATTATCGCGGTACCACACATAACCAGATAACCCTAAATAGCAAATCGCGACGCCAGCCAGGCCGATCGCgg                              <  2:1104970/93‑1 (MQ=255)
                                                                aTCGGCCTGTTTACTGCGTTTATTATCGCGGTACCACACATAACCAGATAACCCTAAATAGCAAATCGCGACGCCAGCCAGGCCGATCGCGGt                             >  1:659993/1‑93 (MQ=255)
                                                                aTCGGCCTGTTTACTGCGTTTATTATCGCGGTACCACACATAACCAGATAACCCTAAATAGCAAATCGCGACGCCAGCCAGGCCGATCGCGGt                             <  1:3542200/93‑1 (MQ=255)
                                                                aTCGGCCTGTTTACGGCGTTTATTATCGCGGTACCACACATAACCAGATAACCCTAAATAGCAAATCGCGACGCCAGCCAGGCCGATCGCGGt                             <  1:730289/93‑1 (MQ=255)
                                                                 tCGGCCTGTTTACTGCGTTTATTATCGCGGTACCACACATAACCAGATAACCCTAAATAGCAAATCGCGACGCCAGCCAGGCCGATCGCGGTa                            >  1:1044698/1‑93 (MQ=255)
                                                                 tCGGCCTGTTTACTGCGTTTATTATCGCGGTACCACACATAACCAGATAACCCTAAATAGCAAATCGCGACGCCAGCCAGGCCGATCGCGGTa                            >  2:4113070/1‑93 (MQ=255)
                                                                  cGGCCTGTTTACTGCGTTTATTATCGCGGTACCACACATAACCAGATAACCCTAAATAGCAAATCGCGACGCCAGCCAGGCCGATCGCGGTaa                           <  2:2649174/93‑1 (MQ=255)
                                                                   ggCCTGTTTACTGCGTTTATTATCGCGGTACCACACATAACCAGATAACCCTAAATAGCAAATCGCGACGCCAGCCAGGCCGATCGCGGTAAt                          >  2:2014450/1‑93 (MQ=255)
                                                                    gCCTGTTTACTGCGTTTATTATCGCGGTACCACACATAACCAGATAACCCTAAATAGCAAATCGCGACGCCAGCCAGGCCGATCGCGGTAATa                         <  2:4302733/93‑1 (MQ=255)
                                                                    gCCTGTTTACTGCGTTTATTATCGCGGTACCACACATAACCAGATAACCCTAAATAGCAAATCGCGACGCCAGCCAGGCCGATCGCGGTAATa                         <  2:4254317/93‑1 (MQ=255)
                                                                    gCCTGTTTACTGCGTTTATTATCGCGGTACCACACATAACCAGATAACCCTAAATAGCAAATCGCGACGCCAGCCAGGCCGATCGCGGTAATa                         >  1:4808519/1‑93 (MQ=255)
                                                                     ccTGTTTACTGCGTTTATTATCGCGGTACCACACATAACCAGATAACCCTAAATAGCAAATCGCGACGCCAGCCAGGCCGATCGCGGTAATAc                        >  2:2572516/1‑93 (MQ=255)
                                                                     ccTGTTTACTGCGTTTATTATCGCGGTACCACACATAACCAGATAACCCTAAATAGCAAATCGCGACGCCAGCCAGGCCGATCGCGGCAATAc                        >  2:4540848/1‑93 (MQ=255)
                                                                     ccTGTTTACTGCGTTTATTATCGCGGGACCACACATAACCAGATAACCCTAAATAGCAAATAGCAACGCCAGCCAGGCAGATCGCGGTAATag                        >  2:4355812/1‑92 (MQ=255)
                                                                      cTGTTTACTGCGTTTATTATCGCGGTACCACACATAACCAGATACCCCTAAATAGCAAATCGCGACGCCAGCCAGGCCGATCGCGGTAATACg                       <  2:4945634/93‑1 (MQ=255)
                                                                        gTTTACTGCGTTTATTATCGCGGTACCACACATAACCAGATAACCCTAAATAGCAAATCGCGACGCCAGCCAGGCCGATCGCGGTAATACGTg                     <  1:5567671/93‑1 (MQ=255)
                                                                         tttACTGCGTTTATTATCGCGGTACCACACATAACCAGATAACCCTAAATAGCAAATCGCGACGCCAGCCAGGCCGATCGCGGTAATACGTga                    <  1:6438366/93‑1 (MQ=255)
                                                                          ttACTGCGTTTATTATCGCGGTACCACACATAACCAGATAACCCAAAATAGCAAATCGCGACGCCAGCCAGGCCGATCGCGGTAATACGTgag                   <  1:4513788/93‑1 (MQ=255)
                                                                           tACTGCGTTTATTATCGCGGTACCACACATAACCAGATAACCCTAAATAGCAAATCGCGACGCCAGCCAGGCCGATCGCGGTAATACGTgaga                  >  2:6360214/1‑93 (MQ=255)
                                                                            aCTGCGTTTATTATCGCGGTACCACACATAACCAGATAACCCTAAATAGCAAATCGCGACGCCAGCCAGGCCGATCGCGGTAATACGTGAGAc                 <  1:130256/93‑1 (MQ=255)
                                                                             cTGCGTTTATTATCGCGGTACCACACATAACCAGATAACCCTAAATAGCAAATCGCGACGCCAGCCAGGCCGATCGCGGTAATACGTGAGAcc                <  2:5961010/93‑1 (MQ=255)
                                                                             cTGCGTTTATTATCGCGGTACCACACATAACCAGATAACCCTAAATAGCAAATCGCGACGCCAGCCAGGCCGATCGCGGTAATACGTGAGAcc                >  1:6546687/1‑93 (MQ=255)
                                                                               gCGTTTATTATCGCGGTACCACACATAACCAGATAACCCTAAATAGCAAATCGCGACGCCAGCCAGGCCGATCGCGGTAATACGTGAGACCAt              >  1:6174071/1‑93 (MQ=255)
                                                                                cGTTTATTATCGCGGTACCACACATAACCAGATAACCCTAAATAGCAAATCGCGACGCCAGCCAGGCCGATCGCGGTAATACGTGAGACCAtt             <  2:3839266/93‑1 (MQ=255)
                                                                                cGTTTATTATCGCGGTACCACACATAACCAGATAACCCTAAATAGCAAATCGCGACGCCAGCCAGGCCGATCGCGGTAATACGTGAGACCAtt             <  2:5448567/93‑1 (MQ=255)
                                                                                 gTTTATTATCGCGGTACCACACATAACCAGATAACCCTAAATAGCAAATCGCGACGCCAGCCAGGCCGATCGCGGTAATACGTGAGACCAttt            <  2:1357126/93‑1 (MQ=255)
                                                                                  tttATTATCGTGGTACCACACATAACCAGATAACCCTAAATAGCAAATCGCGACGCCAGCCAGGCCGATCGCGGTAATACGTGAGACCAtttt           >  1:4052070/1‑93 (MQ=255)
                                                                                  tttATTATCGCGGTACCACACATAACCAGATAACCCTAAATAGCAAATCGCGACGCCAGCCAGGCCGATCGCGGTAATACGTGAGACCAtttt           >  1:4371213/1‑93 (MQ=255)
                                                                                     attatCGCGGTACCACACATAACCAGATAACCCTAAATAGCAAATCGCGACGCCAGCCAGGCCGATCGCGGTAATACGTGAGACCATTTTCat        <  2:5830623/93‑1 (MQ=255)
                                                                                      ttatCGCGGTACCACACATAACCAGATAACCCTAAATAGCAAATCGCGACGCCAGCCAGGCCGATCGCGGTAATACGTGAGACCATTTTCatt       <  2:2487754/93‑1 (MQ=255)
                                                                                        atCGCGGTACCACACATAACCAGATAACCCTAAATAGCAAATCGCGACGCCAGCCAGGCCGATCGCGGTAATACGTGAGACCATTTTCattat     >  2:2696551/1‑93 (MQ=255)
                                                                                         tCGCGGTACCACACATAACCAGATAACCCTAAATAGCAAATCGCGACGCCTGCCAGGCCGATCGCGGTAATACGTGAGACCATTTTCattatt    <  2:4668120/93‑1 (MQ=255)
                                                                                         tCGCGGTACCACACATAACCAGATAACCCTAAATAGCAAATCGCGACGCCAGCCAGGCCGATCGCGGTAATACGTGAGACCATTTTCattatt    <  1:1547461/93‑1 (MQ=255)
                                                                                         tCGCGGTACAACACATAACCAGATAACCCTAAATAGCAACTCGCGACGCCAGCCAGGCCGATCGCGGTAATACGTGAGACCATTTTCattatt    <  1:4083908/93‑1 (MQ=255)
                                                                                            |                                                                                           
CGCATCCTTTTTCGCGGAGAAAGCCTAAAACCTTATTATTTTCACTGACAGCAGATGCCTGAACATCGGCCTGTTTACTGCGTTTATTATCGTGGTACCACACATAACCAGATAACCCTAAATAGCAAATCGCGACGCCAGCCAGGCCGATCGCGGTAATACGTGAGACCATTTTCATTATTCT  >  NC_000913/3666641‑3666824

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: