Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,564,525 G→A 100% M1M (ATG→ATA) † yffS → CPZ‑55 prophage; uncharacterized protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,564,5250GA97.3% 479.1 / 1.7 147M1M (ATG→ATA) †yffSCPZ‑55 prophage; uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base G (1/3);  new base A (80/63);  total (81/66)
Fisher's exact test for biased strand distribution p-value = 3.26e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.59e-01

Only 100 of 147 total aligned reads displayed.

GTGATCTGCGTCAATATTCCCTTCCGTTACGCTTACCCTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATGAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGAT‑GGAAGGACAATATAAGGGCGCTCAACGTGGGGCCAAATGCGAGGAGATCGC  >  NC_000913/2564428‑2564620
                                                                                                 |                                                                                                
gtgATCTGCGTCAATATTCCCTTCCGTTACGCTTACCCTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAg                                                                                                <  2:4032588/100‑1 (MQ=255)
  gATCTGCGTCAATATTCCCTTCCGTTACGCTTACCCTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAAGCATAAGtt                                                                                              >  2:2221145/1‑100 (MQ=255)
   aTCTGCGTCAATATTCCCTTCCGTTACGCTTACCCTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTa                                                                                             >  2:2727785/1‑100 (MQ=255)
    tCTGCGTCAATATTCCCTTCCGTTACGCTTACCCTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTAc                                                                                            >  1:511233/1‑100 (MQ=255)
    tCTGCGTCAATATTCCCTTCCGTTACGCTTACCCTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTAc                                                                                            <  1:1461411/100‑1 (MQ=255)
     cTGCGTCAATATTCCCTTCCGTTACGCTTACCCTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACg                                                                                           <  2:1328139/100‑1 (MQ=255)
      tGCGTCAATATTCCCTTCCGTTACGCTTACCCTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGa                                                                                          <  1:3550751/100‑1 (MQ=255)
      tGCGTCAATATTCCCTTCCGTTACGCTTACCCTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGa                                                                                          <  1:1730426/100‑1 (MQ=255)
       ncGTCAATATTCCCTTCCGTTACGCTTACCCTCCACCAATACGCGTGTTAACCTCCCAATGGATGCTATGAGATGGGAGATAAAAATCATAAGTTACGaa                                                                                         >  1:2054661/2‑100 (MQ=255)
       gCGTCAATATTCCCTTCCGTTACGCTTACCCTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGaa                                                                                         <  1:3532440/100‑1 (MQ=255)
       gCGTCAATATTCCCTTCCGTTACGCTTACCCTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGaa                                                                                         >  1:200945/1‑100 (MQ=255)
       gCGTCAATATTCCCTTCCGTTACGCTTACCCTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGaa                                                                                         >  2:2695464/1‑100 (MQ=255)
        cGTCAATATTCCCTTCCGTTACGCTTACCCTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGaaa                                                                                        >  1:448995/1‑100 (MQ=255)
          tCAATATTCCCTTCCGTTACGCTTACCCTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAAt                                                                                       <  1:3018699/99‑1 (MQ=255)
             atatTCCCTTCCGTTACGCTTACCCTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCaaa                                                                                   >  1:1496401/1‑100 (MQ=255)
             atatTCCCTTCCGTTACGCTTACCCTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCaaa                                                                                   >  1:941763/1‑100 (MQ=255)
               atTCCCTTCCGTTACGCTTACCCTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCaaaa                                                                                  <  2:658607/99‑1 (MQ=255)
                  cgcaaccGTTACGCTTACCCTCCAACAAAACGGGTGTTAACCTCCCAATGTATTCTATGAAAAAAAAAAAAAAAATCATAAGTTACGAAATCAAAATTTg                                                                              <  2:2967576/95‑1 (MQ=255)
                   ccTTCCGTTACGCTTACCCTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATTGGAGATAAAAACCATAAGTTACGAAATCAAAATTTGc                                                                             >  2:60476/1‑100 (MQ=255)
                   ccTTCCGTTACGCTTACCCTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGTAGATAAAAATCATAAGTTACGAAATCAAAATTTGc                                                                             >  1:1570050/1‑100 (MQ=255)
                    cTTCCGTTACGCTTACCCTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCg                                                                            <  1:2181292/100‑1 (MQ=255)
                    cTTCCGTTACGCTTACCCTCCACCAATACGCGCGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCg                                                                            >  2:1897129/1‑100 (MQ=255)
                      tCCGTTACGCTTACCCTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGAc                                                                          <  1:4065453/100‑1 (MQ=255)
                       ccGTTACGCTTACCCTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCg                                                                            <  2:3824630/97‑1 (MQ=255)
                            aCGCTTACCCTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTc                                                                    <  1:2126028/100‑1 (MQ=255)
                            aCGCTTACCCTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTc                                                                    <  2:1123962/100‑1 (MQ=255)
                            aCGCTTACCCTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTc                                                                    <  2:794985/100‑1 (MQ=255)
                             cGCTTACCCTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCa                                                                   <  2:3650778/100‑1 (MQ=255)
                               cTTACCCTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGTGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCaaa                                                                 >  1:897429/1‑100 (MQ=255)
                                 tACCCTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAgg                                                               >  1:110708/1‑100 (MQ=255)
                                   cccTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGgcg                                                             >  1:3139343/1‑100 (MQ=255)
                                   cccTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGgcg                                                             >  1:3234511/1‑100 (MQ=255)
                                   cccTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGgcg                                                             >  2:299581/1‑100 (MQ=255)
                                    ccTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGgcgc                                                            >  2:1034313/1‑100 (MQ=255)
                                    ccTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGgcgc                                                            >  1:2461909/1‑100 (MQ=255)
                                    ccTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGgcgc                                                            >  2:2124357/1‑100 (MQ=255)
                                     cTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCt                                                           >  1:3907311/1‑100 (MQ=255)
                                      tCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGgcgc                                                            <  2:3865469/98‑1 (MQ=255)
                                        caccaATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCg                                                        <  2:3707243/100‑1 (MQ=255)
                                        caccaATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCg                                                        <  2:4078644/100‑1 (MQ=255)
                                        caccaATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCg                                                        >  1:1135046/1‑100 (MQ=255)
                                         accaATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGa                                                       >  1:2194679/1‑100 (MQ=255)
                                          ccaATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGAt                                                      >  1:1873777/1‑100 (MQ=255)
                                           caATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGAT‑g                                                    <  1:3026132/100‑1 (MQ=255)
                                           caATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGAT‑g                                                    <  2:2220397/100‑1 (MQ=255)
                                            aaTACGCGTGTTAACCTCCCAATAGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGAT‑gg                                                   <  1:4031189/100‑1 (MQ=255)
                                               aCGCGTGTTAACCTCCCAATGGATTCTATTAGATGGTAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGAT‑GGAAg                                                <  1:478647/100‑1 (MQ=255)
                                               aCGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCGAAGGCGCTGCGAT‑GGAAg                                                <  2:2809791/100‑1 (MQ=255)
                                               aCGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGAT‑GGaa                                                 <  2:3051011/99‑1 (MQ=255)
                                               aCGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGAT‑GGAAg                                                <  1:2124370/100‑1 (MQ=255)
                                               aCGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGAT‑GGAAg                                                <  2:1771295/100‑1 (MQ=255)
                                                cgcgTGTTAACCTCCCAATGGATTCCATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGAT‑GGaa                                                 >  2:2011433/1‑98 (MQ=255)
                                                  cgTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTAAAAGGCGCTGCGAT‑GGAAGGAc                                             >  2:1951834/1‑100 (MQ=255)
                                                   gtgtTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGAT‑GGAAGGACa                                            <  2:2760023/100‑1 (MQ=255)
                                                    tgtTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGAT‑GGAAGGACaa                                           >  2:470028/1‑100 (MQ=255)
                                                     gtTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCGAAATTTGCGACATCCTCAAAGGCGCTGCGAT‑GGAAGGACAat                                          <  1:2797839/100‑1 (MQ=255)
                                                     gtTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGAT‑GGAAGGACAat                                          <  1:1113275/100‑1 (MQ=255)
                                                        aaCCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGAT‑GGAAGGACAAtata                                       >  1:1205046/1‑100 (MQ=255)
                                                         aCCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATGAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGAGATCGGAAgagcacacgt                                       >  1:4015856/1‑91 (MQ=255)
                                                          ccTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGAT‑GGAAGGACAATATAAg                                     >  1:1047361/1‑100 (MQ=255)
                                                          ccTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGAT‑GGAAGGACAATATAAg                                     >  1:1231469/1‑100 (MQ=255)
                                                           cTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGAT‑GGAAGGACAATATAAgg                                    >  1:2489201/1‑100 (MQ=255)
                                                           cTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGAT‑GGAAGGACAATATAAgg                                    >  2:464101/1‑100 (MQ=255)
                                                           cTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGAT‑GGAAGGACAATATAAg                                     <  2:2582400/99‑1 (MQ=255)
                                                            tCCCAATGGATTCTATGAGATGGGAGATAAAAATTATAAGTTACGAAATCAAAATTTGCGACATCTTCAAAGTCGTTGCGAT‑GGAAGGACAATATAAggg                                   <  2:3358478/100‑1 (MQ=255)
                                                            tCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGAT‑GGAAGGACAATATAAggg                                   <  2:4029299/100‑1 (MQ=255)
                                                             cacAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGAT‑GGAAGGACAATATAAggg                                   <  2:2341868/97‑1 (MQ=255)
                                                              ccAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGAT‑GGAAGGACAATATAAGGgcg                                 >  2:3398331/1‑100 (MQ=255)
                                                              ccAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGAT‑GGAAGGACAATATAAGGgcg                                 >  1:3834017/1‑100 (MQ=255)
                                                              ccAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGAT‑GGAAGGACAATATAAGGgcg                                 >  1:2768001/1‑100 (MQ=255)
                                                              ccAATGGATTCTATGAGATGGGAGATAAAAAGCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGAT‑GGAAGGACAATATAAGGgcg                                 <  2:454964/100‑1 (MQ=255)
                                                               cAATGGATTCTATGAGATGGGAGATAAAAATCATGAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGAT‑GGAAGGACAATATAAGGgcgc                                <  1:860409/100‑1 (MQ=255)
                                                                aaTGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGAT‑GGAAGGACAATATAAGGGCGCt                               <  2:1056965/100‑1 (MQ=255)
                                                                  tGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGAT‑GGAAGGACAATATAAGGGCGCTCa                             >  1:2419576/1‑100 (MQ=255)
                                                                    gATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGAT‑GGAAGGACAATATAAGGGCGCTCAAc                           >  1:207005/1‑100 (MQ=255)
                                                                        cTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGAT‑GGAAGGACAATATAAGGGCGCTCAACGTgg                       >  2:2268332/1‑100 (MQ=255)
                                                                        cTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGAT‑GGAAGGACAATATAAGGGCGCTCAACGTgg                       >  2:1479179/1‑100 (MQ=255)
                                                                        cTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGAT‑GGAAGGACAATATAAGGGCGCTAAACGTgg                       >  1:3793897/1‑100 (MQ=255)
                                                                         tATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGAT‑GGAAGGACAATATAAGGGCGCTCAACGTggg                      >  1:1942027/1‑100 (MQ=255)
                                                                         tATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGAT‑GGAAGGACAATATAAGGGCGCTCAACGTggg                      <  1:4123658/100‑1 (MQ=255)
                                                                         tATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGAT‑GGAAGGACAATATAAGGGCGCTCAACGTggg                      >  1:60525/1‑100 (MQ=255)
                                                                          aTGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGAT‑GGAAGGACAATATAAGGGCGCTCAACGTgggg                     >  2:1536780/1‑100 (MQ=255)
                                                                           tGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGAT‑GGAAGGACAATATAAGGGCGCTCAACGTGGGGc                    >  1:3267743/1‑100 (MQ=255)
                                                                             agaTGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGAT‑GGAAGGACAATATAAGGGCGCTCAACGTGGGGCCa                  >  1:1874567/1‑100 (MQ=255)
                                                                             agaTGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGAT‑GGAAGGACAATATAAGGGCGCTCAACGTGGGGCCa                  >  1:2080761/1‑100 (MQ=255)
                                                                              gaTGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGAT‑GGAAGGACAATATAAGGGCGCTCAACGTGGGGCCaa                 >  1:850205/1‑100 (MQ=255)
                                                                              gaTGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGAT‑GGAAGGACAATATAAGGGCGCTCAACGTGGGGCCaa                 >  1:2314825/1‑100 (MQ=255)
                                                                                tGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGAT‑GGAAGGACAATATAAGGGCGCTCAACGTGGGGCCAAAt               >  2:2135144/1‑100 (MQ=255)
                                                                                 gggAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGAT‑GGAAGGACAATATAAGGGCGCTCAACGTGGGGCCAAATg              <  2:1738572/100‑1 (MQ=255)
                                                                                 gggAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGAT‑GGAAGGACAATATAAGGGCGCTCAACGTGGGGCCAAATg              >  2:1131935/1‑100 (MQ=255)
                                                                                   gagaTAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGAT‑GGAAGGACAATATAAGGGCGCTCAACGTGGGGCCAAATGCg            <  2:1380335/100‑1 (MQ=255)
                                                                                     gaTAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGAT‑GGAAGGACAATATAAGGGCGCTCAACGTGGGGCCAAATGCgag          >  1:3283749/1‑100 (MQ=255)
                                                                                     gaTAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGAT‑GGAAGGACAATATAAGGGCGCTCAACGTGGGGCCAAATGCgag          <  2:4128805/100‑1 (MQ=255)
                                                                                       tAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGAT‑GGAAGGACAATATAAGGGCGCTCAACGTGGGGCCAAATGCgagga        >  1:659297/1‑100 (MQ=255)
                                                                                       tAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGAT‑GGAAGGACAATATAAGGGCGCTCAACGTGGGGCCAAATGCgagga        >  1:749140/1‑100 (MQ=255)
                                                                                        aaaaaTCATAAGTTACGAAATCAAAATTTGCGACATCCTTAAAGGCGCTGCGAT‑GGAAGGACAATATAAGGGCGCTCAACGTGGGGACAAATg              >  2:2500894/1‑93 (MQ=255)
                                                                                            aTCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGAT‑GGAAGGACAATATAAGGGCGCTCAACGTGGGGCCAAATGCGAGGAGATc    <  1:3089056/99‑1 (MQ=255)
                                                                                            aTCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGAT‑GGAAGGACAATATAAGGGCGCTCAACGTGGGGCCAAATGCGAGGAGATCg   <  1:3640453/100‑1 (MQ=255)
                                                                                            aTCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGAT‑GGAAGGACAATATAAGGGCGCTCAACGTGGGGCCAAATGCGAGGAGATCg   >  1:4050772/1‑100 (MQ=255)
                                                                                             tCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGAT‑GGAAGGACAATATAAGGGCGCTCAACGTGGGGCCAAATGCGAGGAGATCGc  <  1:1732390/100‑1 (MQ=255)
                                                                                                 |                                                                                                
GTGATCTGCGTCAATATTCCCTTCCGTTACGCTTACCCTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATGAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGAT‑GGAAGGACAATATAAGGGCGCTCAACGTGGGGCCAAATGCGAGGAGATCGC  >  NC_000913/2564428‑2564620

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