Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,574,749 A→C F85L (TTT→TTG eutQ ← RmlC‑like cupin domain protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,574,7490AC99.0% 687.2 / ‑2.5 194F85L (TTT→TTGeutQRmlC‑like cupin domain protein
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base C (97/94);  total (99/95)
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.03e-01

Only 100 of 194 total aligned reads displayed.

GGTCACCGATTTGAAGCTCGGCTGCATTGCGCCCTGCTCCAGCGACTGTTTTTCCTTCATCACTTTTTCCATCAGCTGTGCGACCAGGCTTTCGGTAAACTGGCCTTCCGGCAGCTGGGCAATGATGGTTTCGCGGATGCGCTGGCTTTCGGTTTTGTCGGCTGGCACGCTGGCGGGCACGGACGCCGTTACC  >  NC_000913/2574653‑2574845
                                                                                                |                                                                                                
ggTCACCGATTTGAAGCTCGGCTGCATTGCGCCCTGCTCCAGCGACTGTTTTTCCTTCATCACTTTTTCCATCAGCTGTGCGACCAGGCTTTCGGTCAAc                                                                                               <  2:3281155/100‑1 (MQ=255)
ggTCACCGATTTGAAGCTCGGCTGCATTGCGCCCTGCTCCAGCGACTGTTTTTCCTTCATCACTTTTTCCATCAGCTGTGCGACCAGGCTTTCGGTCAAc                                                                                               <  1:3092443/100‑1 (MQ=255)
  tCACCGATTTGAAGCTCGGCTGCATTGCGCCCTGCTCCAGCGACTGTTTTTCCTTCATCACTTTTTCCATCAGCTGTGCGACCAGGCTTTCGGTCAACTg                                                                                             <  1:1121260/100‑1 (MQ=255)
   naCCGATTTGAAGCTCGGCTGCATTGCGCCCTGCTCCAGCGACTGTTTTTCCTTCATCACTTTTTCCATCAGCTGTGCGACCAGGCTTTCGGTCAACTgg                                                                                            >  1:4420674/2‑100 (MQ=255)
   cACCGATTTGAAGCTCGGCTGCATTGCGCCCTGCTCCAGCGACTGTTTTTCCTTCATCACTTTTTCCATCAGCTGTGCGACCAGGCTTTCGGTCAACTgg                                                                                            <  2:3924265/100‑1 (MQ=255)
   cACCGATTTGAAGCTCGGCTGCATTGCGCCCTGCTCCAGCGACTGTTTTTCCTTCATCACTTTTTCCATCAGCTGTGCGACCAGGCTTTCGGTCAACTgg                                                                                            <  2:3350532/100‑1 (MQ=255)
   cACCGATTTGAAGCTCGGCTGCATTGCGCCCTGCTCCAGCGACTGTTTTTCCTTCATCACTTTTTCCATCAGCTGTGCGACCAGGCTTTCGGTCAACTgg                                                                                            <  1:2332174/100‑1 (MQ=255)
     ccGATTTGAAGCTCGGCTGCATTGCGCCCTGCTCCAGCGACTGTTTTTCCTTCATCACTTTTTCCATCAGCTGTGCGACCAGGCTTTCGGTCAACTGGcc                                                                                          >  1:5497558/1‑100 (MQ=255)
     ccGATTTGAAGCTCGGCTGCATTGCGCCCTGCTCCAGCGACTGTTTTTCCTTCATCACTTTTTCCATCAGCTGTGCGACCAGGCTTTCGGTCAACTGGcc                                                                                          >  2:4852310/1‑100 (MQ=255)
         tttGAAGCTCGGCTGCATTGCGCCCTGCTCCAGCGACTGTTTTTCCTTCATCACTTTTTCCATCAGCTGTGCGACCAGGCTTTCGGTCAACTGGCCTTcc                                                                                      >  1:4383286/1‑100 (MQ=255)
          ttGAAGCTCGGCTGCATTGCGCCCTGCTCCAGCGACTGTTTTTCCTTCATCACTTTTTCCATCAGCTGTGCGACCAGGCTTTCGGTCAACTGGCCTTCCg                                                                                     <  1:2819353/100‑1 (MQ=255)
           tGAAGCTCGGCTGCATTGCGCCCTGCTCCAGCGACTGTTATTCCTTCATCACTTTTTCCATCAGCTGTGCGACCAGGCTTTCGGTCAACTGGCCTTCCgg                                                                                    >  1:4862380/1‑100 (MQ=255)
            gAAGCTCGGCTGCATTGCGCCCTGCTCCAGCGACTGTTTTTCCTTCATCACTTTTTCCATCAGCTGTGCGACCAGGCTTTCGGTCAACTGGCCTTCCGGc                                                                                   >  1:2935474/1‑100 (MQ=255)
             aaGCTCGGCTGCATTGCGCCCTGCTCCAGCGACTGTTTTTCCTTCATCACTTTTTCCATCAGCTGTGCGACCAGGCTTTCGGTCAACTGGCCTTCCGGCa                                                                                  >  2:5025942/1‑100 (MQ=255)
              aGCTCGGCTGCATTGCGCCCTGCTCCAGCGACTGTTTTTCCTTCATCACTTTTTCCATCAGCTGTGCGACCAGGCTTTCGGTCAACTGGCCTTCCGGCAg                                                                                 >  2:1999498/1‑100 (MQ=255)
                cTCGGCTGCATTGCGCCCTGCTCCAGCGACTGTTTTTCCTTCATCACTTTTTCCATCAGCTGTGCGACCAGGCTTTCGGTCAACTGGCCTTCCGGCAGCt                                                                               <  2:2973951/100‑1 (MQ=255)
                cTCGGCTGCATTGCGCCCTGCTCCAGCGACTGTTTTTCCTTCATCACTTTTTCCATCAGCTGTGCGACCAGGCTTTCGGTCAACTGGCCTTCCGGCAGCt                                                                               >  2:3666494/1‑100 (MQ=255)
                cTCGGCTGCATTGCGCCCTGCTCCAGCGACTGTTTTCCCTACATCACTTTTTCCATCAGCTGTGCGACCAGGCTTTCGGTCAACTGGCCTTCCGGCAGCt                                                                               >  2:5153553/1‑100 (MQ=255)
                   ggCTGCATTGCGCCCTGCTCCAGCGACTGTTTTTCCTTCATCACTTTTTCCATCAGCTGTGCGACCAGGCTTTCGGTCAACTGGCCTTCCGGCAGCTggg                                                                            <  1:54764/100‑1 (MQ=255)
                   ggCTGCATTGCGCCCTGCTCCAGCGACTGTTTTTCCTTCATCACTTTTTCCATCAGCTGTGCGACCAGGCTTTCGGTCAACTGGCCTTCCGGCAGCTggg                                                                            >  1:2666189/1‑100 (MQ=255)
                   ggCTGCATTGCGCCCTGCTCCAGCGACTGTTTTTCCTTCATCACTTTTTCCATCAGCTGTGCGACCAGGCTTTCGGTCAACTGGCCTTCCGGCAGCTggg                                                                            <  2:4257377/100‑1 (MQ=255)
                     cTGCATTGCGCCCTGCTCCAGCGACTGTTTTTCCTTCATCACTTTTTCCATCAGCTGTGCGACCAGGCTTTCGGTCAACTGGCCTTCCGGCAGCTGGGCa                                                                          <  1:176520/100‑1 (MQ=255)
                     cTGCATTGCGCCCTGCTCCAGCGACTGTTTTTCCTTCATCACTTTTTCCATCAGCTGTGCGACCAGGCTTTCGGTCAACTGGCCTTCCGGCAGCTGGGCa                                                                          >  1:1185658/1‑100 (MQ=255)
                        cATTGCGCCCTGCTCCAGCGACTGTTTTTCCTTCATCACTTTTTCCATCAGCTGTGCGACCAGGCTTTCGGTCAACTGGCCTTCCGGCAGCTGGGCAatg                                                                       <  1:1450855/100‑1 (MQ=255)
                         aTTGCGCCCTGCTCCAGCGACTGTTTTTCCTTCATCACTTTTTCCATCAGCTGTGCGACCAGGCTTTCGGTCAACTGGCCTTCCGGCAGCTGGGCAatga                                                                      >  2:1587776/1‑100 (MQ=255)
                         aTTGCGCCCTGCTCCAGCGACTGTTTTTCCTTCATCACTTTTTCCATCAGCTGTGCGACCAGGCTTTCGGTCAACTGGCCTTCCGGCAGCTGGGCAatga                                                                      >  1:4751599/1‑100 (MQ=255)
                          ttGCGCCCTGCTCCAGCGACTGTTTTTCCTTCATCACTTTTTCCATCAGCTGTGCGACCAGGCTTTCGGTCAACTGGCCTTCCGGCAGCTGGGCAatgat                                                                     <  1:2378093/100‑1 (MQ=255)
                           tGCGCCCTGCTCCAGCGACTGTTTTTCCTTCATCACTTTTTCCATCAGCTGTGCGACCAGGCTTTCGGTCAACTGGCCTTCCGGCAGCTGGGCAatgatg                                                                    <  2:581043/100‑1 (MQ=255)
                             cgcCCTGCTCCAGCGACTGTTTTTCCTTCATCACTTTTTCCATCAGCTGTGCGACCAGGCTTTCGGTCAACTGGCCTTCCGGCAGCTGGGCAATGATGGt                                                                  <  2:2312109/100‑1 (MQ=255)
                             cgcCCTGCTCCAGCGACTGTTTTTCCTTCATCACTTTTTCCATCAGCTGTGCGACCAGGCTTTCGGTCAACTGGCCTTCCGGCAGCTGGGCAATGATGGt                                                                  <  2:3121528/100‑1 (MQ=255)
                              gcCCTGCTCCAGCGACTGCTTTTCCTTCATCACTTTTTCCATCAGCTGTGCGACCAGGCTTTCGGTCAACTGGCCTTCCGGCAGCTGGGCAATGATGGtt                                                                 <  1:5224961/100‑1 (MQ=255)
                               cccTGCTCCAGCGACTGTTTTTCCTTCATCACTTTTTCCATCAGCTGTGCGACCAGGCTTTCGGTCAACTGGCCTTCCGGCAGCTGGGCAATGATGGttt                                                                >  1:105427/1‑100 (MQ=255)
                               cccTGCTCCAGCGACTGTTTTTCCTTCATCACTTTTTCCATCAGCTGTGCGACCAGGCTTTCGGTCAACTGGCCTTCCGGCAGCTGGGCAATGATGGttt                                                                <  1:4303288/100‑1 (MQ=255)
                                ccTGCTCCAGCGACTGTTTTTTTTTTTTTTTTTTTTCCATCAGCTGTGCGACCAGGCTTTCGGTTAACTGGCCTTCCGGCGGCTGGGCACTGATGGTTTc                                                               >  2:4755641/1‑100 (MQ=255)
                                  tGCTCCAGCGACTGTTTTTCCTTCATCACTTTTTCCATCAGCTGTGCGACCAGGCTTTCGGTCAACTGGCCTTCCGGCAGCTGGGCAATGATGGTTTcgc                                                             <  1:2833004/100‑1 (MQ=255)
                                   gCTCCAGCGACTGTTTTTCCTTCATCACTTTTTCCATCAGCTGTGCGACCAGGCTTTCGGTCAACTGGCCTTCCGGCAGCTGGGCAATGATGGTTTcgcg                                                            <  1:2966853/100‑1 (MQ=255)
                                    cTCCAGCGACTGTTTTTCCTTCATCACTTTTTCCATCAGCTGTGCGACCAGGCTTTCGGTCAACTGGCCTTCCGGCAGCTGGGCAATGATGGTTTCGCgg                                                           >  1:2999477/1‑100 (MQ=255)
                                    cTCCAGCGACTGTTTTTCCTTCATCACTTTTTCCATCAGCTGTGCGACCAGGCTTTCGGTCAACTGGCCTTCCGGCAGCTGGGCAATGACGGTTTCGCgg                                                           <  2:3871656/100‑1 (MQ=255)
                                     tCCGGCGACTGTTTTTCCTTCATCACTTTTTCCATCAGCTGTGCGACCAGGCTTTCGGTCAACTGGCCTTCCGGCAGCTGGGCAATGATGGTTTCGCGGa                                                          <  2:2246669/100‑1 (MQ=255)
                                     tCCAGCGACTGTTTTTCCTTCATCACTTTTTCCATCAGCTGTGCGACCAGGCTTTCGGTCAACTGGCCTTCCGGCAGCTGGGCAATGATGGTTTCGCGGa                                                          >  1:2672188/1‑100 (MQ=255)
                                     tCCAGCGACTGTTTTTCCTTCATCACTTTTTCCATCAGCTGTGCGACCAGGCTTTCGGTCAACTGGCCTTCCGGCAGCTGGGCAATGATGGTTTCGCGGa                                                          >  1:922224/1‑100 (MQ=255)
                                     tCCAGCGACTGTTTTTCCTTCATCACTTTTTCCATCAGCTGTGCGACCAGGCTTTCGGTCAACTGGCCTTCCGGCAGCTGGGCAATGATGGTTTCGCGGa                                                          >  2:524080/1‑100 (MQ=255)
                                      ccAGCGACTGTTTTTCCTTCATCACTTTTTCCATCAGCTGTGCGACCAGGCTTTCGGTCAACTGGCCTTCCGGCAGCTGCGCAACGATGGTTTCGCGCAt                                                         >  2:5355197/1‑100 (MQ=255)
                                      ccAGCGACTGTTTTTCCTTCATCACCTTTTCCATCAGCTGTGCGACCAGGCTTTCGGTCAACTGGCCTTCCGGCAGCTGGGCGATGATGGTTTCGCGGAt                                                         <  2:2254171/100‑1 (MQ=255)
                                       cAGCGACTGTTTTTCCTTCATCACTTTTTCCATCAGCTGTGCGACCAGGCTTTCGGTCAACTGGCCTTCCGGCAGCTGGGCAATGATGGTTTCGCGGATg                                                        >  1:1920321/1‑100 (MQ=255)
                                       cAGCGACTGTTTTTCCTTCATCACTTTTTCCATCAGCTGTGCGACCAGGCTTTCGGTCAACTGGCCTTCCGGCAGCTGGGCAATGATGGTTTCGCGGATg                                                        <  2:2743248/100‑1 (MQ=255)
                                        aGCGACTGTTTTTCCTTCATCACTTTTTCCATCAGCTGTGCGACCAGGCTTTCGGTCAACTGGCCTTCCGGCAGCTGGGCAATGATGGTTTCGCGGATgc                                                       >  2:1764358/1‑100 (MQ=255)
                                           gACTGTTTTTCCTTCATCACTTTTTCCATCAGCTGTGCGACCAGGCTTTCGGTCAACTGGCCTTCCGGCAGCTGGGCAATGATGGTTTCGCGGATGCGCt                                                    <  1:4384569/100‑1 (MQ=255)
                                           gACTGTTTTTCCTTCATCACTTTTTCCATCAGCTGTGCGACCAGGCTTTCGGTCAACTGGCCTTCCGGCAGCTGGGCAATGATGGTTTCGCGGATGCGCt                                                    <  1:3798582/100‑1 (MQ=255)
                                           gACTGTTTTTCCTTCATCACTTTTTCCATCAGCTGTGCGACCAGGCTTTCGGTCAACTGGCCTTCCGGCAGCTGGGCAATGATGGTTTCGCGGATGCGCt                                                    >  2:424026/1‑100 (MQ=255)
                                              tGTTTTTCCTTCATCACTTTTTCCATCAGCTGTGCGACCAGGCTTTCGGTCAACTGGCCTTCCGGCAGCTGGGCAATGATGGTTTCGCGGATGCGCTGGc                                                 <  2:4209531/100‑1 (MQ=255)
                                                  tttCCTTCATCACTTTTTCCATCAGCTGTGCGACCAGGCTTTCGGTCAACTGGCCTTCCGGCAGCTGGGCAATGATGGTTTCGCGGATGCGCTGGCTTTc                                             <  2:1068282/100‑1 (MQ=255)
                                                   ttCCTTCATCACTTTTTCCATCAGCTGTGCGACCAGGCTTTCGGTCAACTGGCCTTCCGGCAGCTGGGCAATGATGGTTTCGCGGATGCGCTGGCTTTCg                                            >  2:2842779/1‑100 (MQ=255)
                                                    tCCTTCATCACTTTTTCCATCAGCTGTGCGACCAGGCTTTCGGTCAACTGGCCTTCCGGCAGCTGGGCAATGATGGTTTCGCGGATGCGCTGGCTTTCgg                                           <  2:2452658/100‑1 (MQ=255)
                                                      cTTCATCACTTTTTCCATCAGCTGTGCGACCAGGCTTTCGGTCAACTGGCCTTCCGGCAGCTGGGCAATGATGGTTTCGCGGATGCGCTGGCTTTCGGtt                                         >  2:3691981/1‑100 (MQ=255)
                                                       ttCATCACTTTTTCCATCAGCTGTGCGACCAGGCTTTCGGTCAACTGGCCTTCCGGCAGCTGGGCAATGATGGTTTCGCGGATGCGCTGGCTTTCGGttt                                        <  2:4997751/100‑1 (MQ=255)
                                                       ttCATCACTTTTTCCATCAGCTGTGCGACCAGGCTTTCGGTAAACTGGCCTTCCGGCAGCTGGGCAATGATGGTTTCGCGGATGCGCTGGCTTTCGGttt                                        >  1:3221267/1‑100 (MQ=255)
                                                        tcatcaCTTTTTCCATCAGCTGTGCGACCAGGCTTTCGGTCAACTGGCCTTCCGGCAGCTGGGCAATGATGGTTTCGCGGATGCGCTGGCTTTCGGtttt                                       <  1:4358314/100‑1 (MQ=255)
                                                          atcaCTTTTTCCATCAGCTGTGCGACCAGGCTTTCGGTCAACTGGCCTTCCGGCAGCTGGGCAATGATGGTTTCGCGGATGCGCTGGCTTTCGGTTTTGt                                     <  2:4443620/100‑1 (MQ=255)
                                                          atcaCTTTTTCCATCAGCTGTGCGACCAGGCTTTCGGTCAACTGGCCTTCCGGCAGCTGGGCAATGATGGTTTCGCGGATGCGCTGGCTTTCGGTTTTGt                                     >  2:2628952/1‑100 (MQ=255)
                                                           tcaCTTTTTCCATCAGCTGTGCGACCAGGCTTTCGGTCAACTGGCCTTCCGGCAGCTGGGCAATGATGGTTTCGCGGATGCGCTGGCTTTCGGTTTTGTc                                    >  1:3573292/1‑100 (MQ=255)
                                                           tcaCTTTTTCCATCAGCTGTGCGACCAGGCTTTCGGTCAACTGGCCTTCCGGCAGCTGGGCAATGATGGTTTCGCGGATGCGCTGGCTTTCGGTTTTGTc                                    >  1:2046659/1‑100 (MQ=255)
                                                           tcaCTTTTTCCATCAGCTGTGCGACCAGGCTTTCGGTCAACTGGCCTTCCGGCAGCTGGGCAATGATGGTTTCGCGGATGCGCTGGCTTTCGGTTTTGTc                                    >  2:271719/1‑100 (MQ=255)
                                                           tcaCTTTTTCCATCAGCTGTGCGACCAGGCTTTCGGTCAACTGGCCTTCCGGCAGCTGGGCAATGATGGTTTCGCGGATGCGCTGGCTTTCGGTTTTGTc                                    <  1:1368997/100‑1 (MQ=255)
                                                           tcaCTTTTTCCATCAGCTGTGCGACCAGGCTTTCGGTCAACTGGCCTTCCGGCAGCTGGGCAATGATGGTTTCGCGGATGCGCTGGCTTTCGGTTTTGTc                                    >  2:1392263/1‑100 (MQ=255)
                                                            caCTTTTTCCATCAGCTGTGCGACCAGGCTTTCGGTCAACTGGCCTTCCGGCAGCTGGGCAATGATGGTTTCGCGGATGCGCTGGCTTTCGGTTTTGTCg                                   <  2:3039234/100‑1 (MQ=255)
                                                            caCTTTTTCCATCAGCTGTGCGACCAGGCTTTCGGTCAACTGGCCTTCCGGCAGCTGGGCAATGATGGCTTCGCGGATGCGCTGGCTTTCGGTTTTGTCg                                   <  1:4465203/100‑1 (MQ=255)
                                                              cTTTTTCCATCAGCTGTGCGACCAGGCTTTCGGTCAACTGGCCTTCCGGCAGCTGGGCAATGATGGTTTCGCGGATGCGCTGGCTTTCGGTTTTGTCGGc                                 >  2:209397/1‑100 (MQ=255)
                                                                ttttCCATCAGCTGTGCGACCAGGCTTTCGGTCAACTGTGCTTCCGGCAGCTGGGCAATGATGGTTTCGCGGATGCGCTGGCTTTCGCTTTTGTCGGCTg                               <  2:5284528/100‑1 (MQ=255)
                                                                ttttCCATCAGCTGTGCGACCAGGCTTTCGGTCAACTGGCCTTCCGGCAGCTGGGCAATGATGGTTTCGCGGATGCGCTGGCTTTCGGTTTTGTCGGCTg                               <  2:4494459/100‑1 (MQ=255)
                                                                ttttCCATCAGCTGTGCGACCAGGCTTTCGGTCAACTGGCCTTCCGGCAGCTGGGCAATGATGGTTTCGCGGATGCGCTGGCTTTCGGTTTTGTCGGCTg                               <  2:1868194/100‑1 (MQ=255)
                                                                   tCCATCAGCTGTGCGACCAGGCTTTCGGTCAACTGGCCTTCCGGCAGCTGGGCAATGATGGTTTCGCGGATGCGCTGGCTTTCGGTTTTGTCGGCTGGCa                            >  2:2559511/1‑100 (MQ=255)
                                                                   tCCATCAGCTGTGCGACCAGGCTTTCGGTCAACTGGCCTTCCGGCAGCTGGGCAATGATGGTTTCGCGGATGCGCTGGCTTTCGGTTTTGTCGGCTGGCa                            <  1:7030/100‑1 (MQ=255)
                                                                    ccATCAGCTGTGCGACCAGGCTTTCGGTCAACTGGCCTTCCGGCAGCTGGGCAATGATGGTTTCGCGGATGCGCTGGCTTTCGGTTTTGTCGGCTGGCAc                           >  2:3837829/1‑100 (MQ=255)
                                                                     cATCAGCTGTGCGACCAGGCTTTCGGTCAACTGGCCTTCCGGCAGCTGGGCAATGATGGTTTCGCGGATGCGCTGGCTTTCGGTTTTGTCGGCTGGCACg                          >  2:937756/1‑100 (MQ=255)
                                                                       tCAGCTGTGCGACCAGGCTTTCGGTCAACTGGCCTTCCGGCAGCTGGGCAATGATGGTTTCGCGGATGCGCTGGCTTTCGGTTTTGTCGGCTGGCACGCt                        <  2:4109403/100‑1 (MQ=255)
                                                                        cAGCTGTGCGACCAGGCTTTCGGTCAACTGGCCTTCCGGCAGCTGGGCAATGATGGTTTCGCGGATGCGCTGGCTTTCGGTTTTGTCGGCTGGCACGCTg                       <  2:5062063/100‑1 (MQ=255)
                                                                         aGCTGTGCGACCAGGCTTTCGGTCAACTGGCCTTCCGGCAGCTGGGCAATGATGGTTTCGCGGATGCGCTGGCTTTCGGTTTTGTCGGCTGGCACGCTgg                      >  1:2855772/1‑100 (MQ=255)
                                                                         aGCTGTGCGACCAGGCTTTCGGTCAACTGGCCTTCCGGCAGCTGGGCAATGATGGTTTCGCGGATGCGCTGGCTTTCGGTTTTGTCGGCTGGCACGCTgg                      <  1:1571735/100‑1 (MQ=255)
                                                                           cTGTGCGACCAGGCTTTCGGTCAACTGGCCTTCCGGCAGCTGGGCAATGATGGTTTCGCGGATGCGCTGGCTTTCGGTTTTGTCGGCTGGCACGCTGGCg                    <  1:5145358/100‑1 (MQ=255)
                                                                           cTGTGCGACCAGGCTTTCGGTCAACTGGCCTTCCGGCAGCTGGGCAATGATGGTTTCGCGGATGCGCTGGCTTTCGGTTTTGTCGGCTGGCACGCTGGCg                    <  1:242471/100‑1 (MQ=255)
                                                                           cTGTGCGACCAGGCTTTCGGTCAACTGGCCTTCCCGCAGCTGGGCAATGATGGTTTCGCGGATGCGCTGGCTTTCGGTTTTGTCGGCTGGCACGCTGGCg                    <  1:2163129/100‑1 (MQ=255)
                                                                            tgtgCGACCAGGCTTTCGGTCAACTGGCCTTCCGGCAGCTGGGCAATGATGGTTTCGCGGATGCGCTGGCTTTCGGTTTTGTCGGCTGGCACGCTGGCgg                   <  1:4526768/100‑1 (MQ=255)
                                                                            tgtgCGACCAGGCTTTCGGTCAACTGGCCTTCCGGCAGCTGGGCAATGATGGTTTCGCGGATGCGCTGGCTTTCGGTTTTGTCGGCTGGCACGCTGGCgg                   <  1:5542959/100‑1 (MQ=255)
                                                                             gtgCGACCAGGCTTTCGGTCAACTGGCCTTCCGGCAGCTGGGCAATGATGGTTTCGCGGATGCGCTGGCTTTCGGTTTTGTCGGCTGGCACGCTGGCggg                  >  1:359987/1‑100 (MQ=255)
                                                                             gtgCGACCAGGCTTTCGGTCAACTGGCCTTCCGGCAGCTGGGCAATGATGGTTTCGCGGATGCGCTGGCTTTCGGTTTTGTCGGCTGGCACGCTGGCggg                  <  1:3121804/100‑1 (MQ=255)
                                                                             gtgCGACCAGGCTTTCGGTCAACTGGCCTTCCGGCAGCTGGGCAATGATGGTTTCGCGGATGCGCTGGCTTTCGGTTTTGTCGGCTGGCACGCTGGCggg                  <  1:3525584/100‑1 (MQ=255)
                                                                             gtgCGACCAGGCTTTCGGTCAACTGGCCTTCCGGCAGCTGGGCAATGATGGTTTCGCGGATGCGCTGGCTTTCGGTTTTGTCGGCTGGCACGCTGGCggg                  >  1:4596583/1‑100 (MQ=255)
                                                                              tgCGACCAGGCCTTCGGTCAACTGGCCTTCCGGCAGCTGGGCAATGATGGTTTCGCGGATGCGCTGGCTTTCGGCTTTGTCGGCTGGCACGCTGGCGGGc                 <  2:4379115/100‑1 (MQ=255)
                                                                               gCGACCAGGCTTTCGGTCAACTGGCCTTCCGGCAGCTGGGCAATGATGGTTTCGCGGATGCGCTGGCTTTCGGTTTTGTCGGCTGGCACGCTGGCGGGCa                <  2:757834/100‑1 (MQ=255)
                                                                                   ccAGGCTTTCGGTCAACTGGCCTTCCGGCAGCTGGGCAATGATGGTTTCGCGGATGCGCTGGCTTTCGGTTTTGTCGGCTGGCACGCTGGCGGGCACGGa            >  1:324422/1‑100 (MQ=255)
                                                                                   ccAGGCTTTCGGTCAACTGGCCTTCCGGCAGCTGGGCAATGATGGTTTCGCGGATGCGCTGGCTTTCGGTTTTGTCGGCTGGCACGCTGGCGGGCACGGa            >  2:1499485/1‑100 (MQ=255)
                                                                                     aGGCTTTCGGTCAACTGGCCTTCCGGCAGCTGGGCAATGATGGTTTCGCGGATGCGCTGGCTTTCGGTTTTGTCGGCTGGCACGCTGGCGGGCACGGACg          <  1:858071/100‑1 (MQ=255)
                                                                                     aGGCTTTCGGTCAACTGGCCTTCCGGCAGCTGGGCAATGATGGTTTCGCGGATGCGCTGGCTTTCGGTTTTGTCGGCTGGCACGCTGGCGGGCACGGACg          <  2:1147669/100‑1 (MQ=255)
                                                                                      ggCTTTCGGTCAACTGGCCTTCCGGCAGCTGGGCAATGATGGTTTCGCGGATGCGCTGGCTTTCGGTTTTGTCGGCTGGCACGCTGGCGGGCACGGACGc         >  1:58258/1‑100 (MQ=255)
                                                                                        cTTTCGGTCAACTGGCCTTCCGGCAGCTGGGCAATGATGGTTTCGCGGATGCGCTGGCTTTCGGTTTTGTCGGCTGGCACGCTGGCGGGCACGGACGCCg       >  2:291569/1‑100 (MQ=255)
                                                                                         tttCGGTCAACTGGCCTTCCGGCAGCTGGGCAATGATGGTTTCGCGGATGCGCTGGCTTTCGGTTTTGTCGGCTGGCACGCTGGCGGGCACGGACGCCGt      >  1:4139933/1‑100 (MQ=255)
                                                                                          ttCGGTCAACTGGCCTTCCGGCAGCTGGGCAATGATGGTTTCGCGGATGCGCTGGCTTTCGGTTTTGTCGGCTGGCACGCTGGCGGGCACGGACGccagc     <  2:2399962/100‑4 (MQ=255)
                                                                                            cGGTCAACTGGCCTTCCGGCAGCTGGGCAATGATGGTTTCGCGGATGCGCTGGCTTTCGGTTTTGTCGGCTGGCACGCTGGCGGGCACGGACGCCGTTAc   >  2:318319/1‑100 (MQ=255)
                                                                                             ggTCAACTGGCCTTCCGGCAGCTGGGCAATGATGGTTTCGCGGATGCGCTGGCTTTCGGTTTTGTCGGCTGGCACGCTGGCGGGCACGGACGCCGTTAcc  >  1:1290818/1‑100 (MQ=255)
                                                                                                |                                                                                                
GGTCACCGATTTGAAGCTCGGCTGCATTGCGCCCTGCTCCAGCGACTGTTTTTCCTTCATCACTTTTTCCATCAGCTGTGCGACCAGGCTTTCGGTAAACTGGCCTTCCGGCAGCTGGGCAATGATGGTTTCGCGGATGCGCTGGCTTTCGGTTTTGTCGGCTGGCACGCTGGCGGGCACGGACGCCGTTACC  >  NC_000913/2574653‑2574845

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