Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 265,418 A→G 100% R42R (CGT→CGC ykfF ← CP4‑6 prophage; uncharacterized protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913265,4180AG98.2% 346.8 / ‑3.0 108R42R (CGT→CGCykfFCP4‑6 prophage; uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base G (57/49);  total (58/50)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.48e-01

Only 100 of 109 total aligned reads displayed.

CTGTACGGATGCCTGAGGTGCGGATATAGCGGTTAAGACCTTCACCGGCATCCGGCTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCACGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCTTCGAGTGTGATATTGCTGTACGTGGTAGTGACCGCTTGCGCTTGTCTCCGGGTGAAAGG  >  NC_000913/265321‑265513
                                                                                                 |                                                                                               
cTGTACGGATGCCTGAGGTGCGGATATAGCGGTTAAGACCTTCACCGGCATCCGGCTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATcgcg                                                                                               >  2:1417844/1‑100 (MQ=255)
 tGTACGGATGCCTGAGGTGCGGATATAGCGGTTAAGACCTTCACCGGCATCCGGCTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGa                                                                                              >  1:2536151/1‑100 (MQ=255)
 tGTACGGATGCCTGAGGTGCGGATATAGCGGTTAAGACCTTCACCGGCATCCGGCTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGa                                                                                              >  2:2548760/1‑100 (MQ=255)
 tGTACGGATGCCTGAGGTGCGGATATAGCGGTTAAGACCTTCACCGGCATCCGGCTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGa                                                                                              >  1:505788/1‑100 (MQ=255)
  gTACGGATGCCTGAGGTGCGGATATAGCGGTTAAGACCATCACCGGCATCCGGCTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGaa                                                                                             >  2:1864349/1‑100 (MQ=255)
   tACGGATGCCTGAGGTGCGGATATAGCGGTTAAGACCTTCACCGGCATCCGGCTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAac                                                                                            <  2:2258614/100‑1 (MQ=255)
   tACGGATGCCTGAGGTGCGGAAATAGCGGTTAAGACCTTCACCGGCATCCGGCTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAac                                                                                            >  2:1493388/1‑100 (MQ=255)
   tACGGATGCCTGAAGTGCGGATATAGCGGTTAAGACCGACAACGGCATCCGGCTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGaa                                                                                             >  2:1877551/1‑99 (MQ=21)
       gATGCCTGAGGTGCGGATATAGCGGTTAAGACCTTCACCGGCATCCGGCTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAaccacc                                                                                        <  2:1261096/100‑1 (MQ=255)
       gATGCCTGAGGTGCGGATATAGCGGTTAAGACCTTCACCGGCATCCGGCTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAaccacc                                                                                        <  2:2609411/100‑1 (MQ=255)
          gCCTGAGGTGCGGATATAGCGGTTAAGACCTTCACCGGCATCCGGCTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGa                                                                                     >  2:1551632/1‑100 (MQ=255)
           ccTGAGGTGCGGATATAGCGGTTAAGACCTTCACCGGCATCCGGCTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGAc                                                                                    >  1:2723972/1‑100 (MQ=255)
             tGAGGTGCGGATATAGCGGTTAAGACCTTCACCGGCATCCGGCGCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGACgg                                                                                  <  2:1934883/100‑1 (MQ=255)
                ggTGCGGATATAGCGGTTAAGACCTTCACCGGCATCCGGCTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAg                                                                               >  2:1035434/1‑100 (MQ=255)
                ggTGCGGATATAGCGGTTAAGACCTTCACCGGCATCCGGCTCAAAGTTCCATGCCCGCCAAACCATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAg                                                                               >  1:1579591/1‑100 (MQ=255)
                 gTGCGGATATAGCGGTTAAGACCTTCACCGGCATCCGGCTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGt                                                                              <  2:2141777/100‑1 (MQ=255)
                 gTGCGGATATAGCGGTTAAGACCTTCACCGGCATCCGGCTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGt                                                                              >  1:1669014/1‑100 (MQ=255)
                  tGCGTATATAGCGGTTAAGACCTTCACCGTCATCCGGCTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGTg                                                                             <  1:2881485/100‑1 (MQ=255)
                  tGCGGATATAGCGGTTAAGACCTTCACCGGCATCCGGCTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAAACACCAGACGGAAGTg                                                                             >  1:1269077/1‑100 (MQ=255)
                    cGGATATAGCGGTTAAGACCTTCACCGGCATCCGGCTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGTGAc                                                                           <  2:1994602/100‑1 (MQ=255)
                    cGGATATAGCGGTTAAGACCTTCACCGGCATCCGGCTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGTGAc                                                                           >  1:1994625/1‑100 (MQ=255)
                      gATATAGCGGTTAAGACCTTCACCGGCATCCGGCTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGTGACTg                                                                         <  2:1561870/100‑1 (MQ=255)
                        tataGCGTTTAAGAGCCTCACCGGCATCCCGCTCAAAGTTCCAAGCCCGCCCGACCATCCGGCCTTCAGTATCGCCATCCACCAGACGGAAGAGACTga                                                                        <  2:496495/99‑2 (MQ=255)
                         ataGCGGTTAAGACCTTCACCGGCATCCGGCTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGTGACTGccc                                                                      >  1:281533/1‑100 (MQ=255)
                           aGCGGTTAAGACCTTCACCGGCATCCGGCTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGTGACTGCCCTg                                                                    >  2:1085288/1‑100 (MQ=255)
                           aGCGGTTAAGACCTTCACCGGCATCCGGCTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGTGACTGCCCTg                                                                    >  2:2014261/1‑100 (MQ=255)
                            gCGGTTAAGACCTTCACCGGCATCCGGCTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGTGACTGCCCTgg                                                                   <  1:2096795/100‑1 (MQ=255)
                              ggTTAAGACCTTCACCGGCATCCGGCTAGAAGTTGCATGCGCGCCAGACGATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGGGACTGcc                                                                       >  2:2144228/1‑94 (MQ=255)
                               gTTAAGACCTTCACCGGCATCCGGCTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGTGACTGCCCTGgtcg                                                                <  1:2178337/100‑1 (MQ=255)
                               gTTAAGACCTTCACCGGCATCCGGCTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGTGACTGCCCTGgtcg                                                                <  1:2547494/100‑1 (MQ=255)
                                 tAAGACCTTCACCGGCATCCGGCTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGTGACTGCCCTGgtcgtc                                                              >  1:2783269/1‑100 (MQ=255)
                                  aaGACCTTCACCGGCATCCGGCTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCt                                                             <  2:2596465/100‑1 (MQ=255)
                                   aGACCTTCACCGGCATCCGGCTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCtt                                                            >  1:2321929/1‑100 (MQ=255)
                                      ccTTCACCGGCATCCGGCTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCTTCGa                                                         >  2:748222/1‑100 (MQ=255)
                                       cTTCACCGGCATCCGGCTTAAAGTTCCATGTCCGCGAGATCATCCGGCCTTCAGTATCGCGAACCACCAGACGGTAGTGACTGCCCTGGTCGTCTTCGa                                                         <  2:1917876/99‑1 (MQ=25)
                                       cTTCACCGGCATCCGGCTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCTTCGAg                                                        <  1:2746364/100‑1 (MQ=255)
                                       cTTCACCGGCATCCGGCTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCTTCGAg                                                        >  1:2432072/1‑100 (MQ=255)
                                       cTTCACCGGCATCCGGCTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCTTCGAg                                                        <  1:132864/100‑1 (MQ=255)
                                       cTTCACCGGCATCCGGCTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCTTCGAg                                                        <  1:2570558/100‑1 (MQ=255)
                                        ttCACCGGCATCCGGCTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGTGACTGCCGTGGTCGTCTTCGa                                                         >  1:3005214/1‑98 (MQ=255)
                                        ttCACCGGCATCCGGCTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCTTCGAgt                                                       >  1:100787/1‑100 (MQ=255)
                                        ttCACCGGCATCCGGCTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCTTCGAgt                                                       >  2:204686/1‑100 (MQ=255)
                                         tCACCGGCATCCGGCTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCTTCGAgtg                                                      <  1:999364/100‑1 (MQ=255)
                                          cACCGGCATCCGGCTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAACCTCCAGACGGAAGTGACTGCCCTGGTCGTCTTCGAgtgt                                                     >  2:2041784/1‑100 (MQ=255)
                                          cACCGGCATCCGGCTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCTTCGAgtgt                                                     <  1:2270189/100‑1 (MQ=255)
                                          cACCGGCATCCGGCTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCTTCGAgtct                                                     <  1:2041807/100‑3 (MQ=255)
                                            ccGGCATCCGGCTCAAAGTTCCATGTCCGCCAGATTATCCGGTTTTCAGTATCGCGAACCAGCAGACGGAAGTGACTGCCCTGGTCGTCTTCGAGTGTGa                                                   <  1:144077/100‑1 (MQ=25)
                                            ccGGCATCCGGCTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCTTCGAgtagga                                                   >  2:2398992/1‑96 (MQ=255)
                                            ccGGCATCCGGCTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCGCAAACCACCAAACGGAAGTAACTGCCCTGGTCGTCTTCGAGTATGa                                                   >  1:2807727/1‑100 (MQ=25)
                                              ggCATCCGGCTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCTTCGAGTGTGata                                                 >  1:697636/1‑100 (MQ=255)
                                                cATCCGGCTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCTTCGAGTGTGATAtt                                               <  1:2776645/100‑1 (MQ=255)
                                                   ccGGCTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCTTCGAGTGTGATATTg                                              >  2:2819249/1‑98 (MQ=255)
                                                   ccGGCTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCTTCGAGTGTGATATTGCt                                            <  1:1433782/100‑1 (MQ=255)
                                                    cGGCTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCTTCGAGTGTGATATTGCTg                                           <  2:2155251/100‑1 (MQ=255)
                                                     ggCTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGTAACTGCCCTGGTCGTCTTCGAGTGTGATCTTGCTGt                                          >  2:2970788/1‑100 (MQ=255)
                                                      gCTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCTTCGAGTGTGATATTGCTGTa                                         >  1:1580963/1‑100 (MQ=255)
                                                       cTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGTGACCGCCCTGGTCGTCTTCGAGTGTGATATTGCTGTAc                                        >  2:1271238/1‑100 (MQ=255)
                                                         cAAAGTTCCATGCCCGCCAGACCATCTGGTCTTCAGTATCGCGAACCACCAGATGGAAGTGACTGCCCTGGTCGTCTTCGAGTGTGATATTGCTGTACg                                       <  1:2199560/99‑1 (MQ=37)
                                                         cAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCTTCGAGTGTGATATTGCTGTACGt                                      >  2:2169430/1‑100 (MQ=255)
                                                         cAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCTTCGAGTGTGATATTGCTGTACGt                                      >  1:1928171/1‑100 (MQ=255)
                                                          aaaGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCTTCGAGTGTGATATTGATGTACGt                                      >  2:1709534/1‑99 (MQ=255)
                                                           aaGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCTTCGAGTGTGATATTGCTGTACGTgg                                    <  2:505785/100‑1 (MQ=255)
                                                           aaGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCACGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCTTCGAGTGTGATATTGCTGTACGTgg                                    <  2:274073/100‑1 (MQ=255)
                                                            aGTTCCATGCCCGCCAGACTATCCTGCCTTTAGTATCGCGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCTTCGAGTGTGATATTGCTGTACGTGGt                                   <  1:1909338/100‑1 (MQ=38)
                                                            aGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCTTCGAGTGTGATATTGCTGTACGTgat                                   <  1:2041335/100‑3 (MQ=255)
                                                            aGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCTTCGAGTGTGATATTGCTGTACGTGGt                                   <  1:1483357/100‑1 (MQ=255)
                                                            aGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCTTCGAGTGTGATATTGCTGTACGTGGt                                   <  1:145796/100‑1 (MQ=255)
                                                             gTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCTTCGAGTGTGATATTGCTGTACGTGGTa                                  <  1:2113192/100‑1 (MQ=255)
                                                              tcccATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCTTCGAGTGTGATATTGCTGTACGTGGTAg                                 <  1:222452/98‑1 (MQ=255)
                                                               tCCATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCTTCGAGTGTGATATTGCTGTACGTGGTAGt                                >  1:2578629/1‑100 (MQ=255)
                                                                ccATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCTTCGAGTGTGATATTGCTGTACGTGGTAGTg                               >  1:2364447/1‑100 (MQ=255)
                                                                 cATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCTTCGAGTGTGATATTGCTGTACGTGGTAGTGa                              >  1:115361/1‑100 (MQ=255)
                                                                  aTGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAACTACCAGACGGAAGTGACTGCCCTGGTCGTCTTCGAGTGTGATATTGCTGTACGTGGTAGTGAc                             <  1:2545496/100‑1 (MQ=255)
                                                                  aTGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGTGACTGCTCTGGTCGTCTTCGAGTGTGATATTGCTGTACGTGGTAGTggc                             <  1:2046263/100‑3 (MQ=255)
                                                                  aTGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCTTCGAGTGTGATATTGCTGTACGTGGTAGTGAc                             >  1:2017241/1‑100 (MQ=255)
                                                                  aTGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGTGACTGCCCTGGTAGTCTTCGAGTGTGATATTGCTGTACGTGGTAGTGAc                             >  1:108058/1‑100 (MQ=255)
                                                                   tGCCCGCCAGACTATTCTGCCTTCAGTATCGCGTACTACCAGACGGAAGTGACTGCCCTGGTCGTCTTCGAGTGTGATATTGCTGTACGTGGTAGTGAcc                            <  1:1411285/100‑1 (MQ=38)
                                                                   tGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCTTCGAGTGTGATATTGCTGTACGTGGTAGTGAcc                            <  2:820421/100‑1 (MQ=255)
                                                                   tGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGTGACCGCCCTGGTCGTCTTCGAGTGTGATATTGCTGTACGTGGTAGTGAcc                            >  2:2634075/1‑100 (MQ=255)
                                                                    gCCCGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCTTCGAGTGTGATATTGCTGTACGTGGTAGTGACcg                           >  2:1892694/1‑100 (MQ=255)
                                                                    gCCCGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCTTCGAGTGTGATATTGCTGTACGTGGTAGTGACcg                           <  2:2244874/100‑1 (MQ=255)
                                                                     cccGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCTTCGAGTGTGATATTGCTGTACGTGGTAGTGACcgc                          >  1:2290441/1‑100 (MQ=255)
                                                                     cccGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCTTCGAGTGTGATATTGCTGTACGTGGTAGTGACcgc                          >  1:2661394/1‑100 (MQ=255)
                                                                      ccGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCTTCGAGTGTGATATTGCTGTACGTGGTAATGACcgct                         >  1:3045517/1‑100 (MQ=255)
                                                                         ccAGACCATCCGGTCTTCAGTATCGCTAACCATCAGACGGAAGTGACTGCCCTGGTCGTCTGCGAATGTGATATTGCTGTACGTGGTAGTGACcgcttgc                      <  1:236064/100‑1 (MQ=38)
                                                                         ccAGACCATCCGGCCTTCACTATAGCAAAGCACCAAACGGAAGTGACGGCCCTGGTCGTCTTTGAGAGTGATAATGCTCTACGTGGTAGTGACcgctagc                      >  2:1785452/1‑100 (MQ=255)
                                                                            gACCATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCTTCGAGTGTGATATTGCTGTACGTGGTAGTGACcgcttgcgct                   >  2:1908363/1‑100 (MQ=255)
                                                                             aCCATCCGGCCTTCAGTATCGCGAACCACCAGACGGTAGTGACTGCCCTGGTCGTCTTCGAGTGTGATATTGCTGTACGTGGTAGTGACcgcttgcgctt                  <  1:2131704/100‑1 (MQ=255)
                                                                               cATCCGGCCTTCAGTATCGCGAATCACCAGACGGAAGTGACTGCCCTGGTCGTCTTTGAGTGTCATATTGCTGTACGTGGTAGTGACagcttgcgcttg                 <  2:1600834/99‑1 (MQ=255)
                                                                               cATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCTTCGAGTGTGATATTGCTGTACGTGGTAGTGACCGCTTGCGCTTGt                >  2:291214/1‑100 (MQ=255)
                                                                                aTCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGCGACTGCCCTGGTCGTCTTCGAGTGTGATATTGCTGTACGTGGTAGTGACCGCTTGCGCTTGtc               >  1:2141777/1‑100 (MQ=255)
                                                                                  ccGGCCTTCAGTATCGCGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCTTCGAGTGTGATATTGCTGTACGTGGTAGTGACCGCTTGCGCTTGtctc             <  1:1093787/100‑1 (MQ=255)
                                                                                     gCCTTCAGTATCGCGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCTTCGAGTGTGATATTGCTGTACGTGGTAGTGACCGCTTGCGCTTGTCTCCgg          <  1:2112026/100‑1 (MQ=255)
                                                                                      ccTTCAGTATCGCGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCTTCGAGTGTGATATTGCTGTACGTGGTAGTGACCGCTTGCGCTTGTCTCCggg         <  1:1908386/100‑1 (MQ=255)
                                                                                       cTTCAGTATCGCGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCTTCGAGCGTGATATTGCTGTACGTGGTAGTGACCGCTTGCGCTTGTCTCCGGGt        >  1:2618381/1‑100 (MQ=255)
                                                                                        ttCAGTATCGCGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCTTCGAGTGTGATATTGCTGTACGTGGTAGTGACCGCTTGCGCTTGTCTCCGGGTg       >  1:1432621/1‑100 (MQ=255)
                                                                                          cAGTATCGCGAACCCCCAGTCGGAAGTGTCTACCCTGGTCGGCTTCGAGTGTGACATTGCTGTATGTGTTAGTCACCTCTTGCTCTTGTCTCCGGGTGaa     <  2:2864012/100‑1 (MQ=255)
                                                                                          cAGTATCGCGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCTTCGAGTGTGATATTGCTGTACGTGGTAGTGACCGCTTGCGCTTGTCTCCGGGTGaa     <  2:2858691/100‑1 (MQ=255)
                                                                                            gTATCGCGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCTTCGAGTGTGATATTGCTGTACGTGGTAGTGACCGCTTGCGCTTGTCTCCGGGTGAAAg   >  1:2776672/1‑100 (MQ=255)
                                                                                             tATCGCGAACCACCAGACGGAAGTGACTGTCCTGGTCGTCTTCGAGTGTGATATTGCTGTACGTTGTAGTGACCGCTTGCGCTTGTCTCCGGGTGAAAgg  <  2:2783237/100‑1 (MQ=255)
                                                                                                 |                                                                                               
CTGTACGGATGCCTGAGGTGCGGATATAGCGGTTAAGACCTTCACCGGCATCCGGCTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCACGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCTTCGAGTGTGATATTGCTGTACGTGGTAGTGACCGCTTGCGCTTGTCTCCGGGTGAAAGG  >  NC_000913/265321‑265513

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