Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,724,344 T→C P104P (CCA→CCG sodC ← superoxide dismutase, Cu, Zn, periplasmic

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,724,3440TC99.0% 342.2 / ‑4.6 98P104P (CCA→CCGsodCsuperoxide dismutase, Cu, Zn, periplasmic
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (47/50);  total (48/50)
Fisher's exact test for biased strand distribution p-value = 4.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.48e-01

CATATTATCGCCGCCAACGTGGACCATCAGCGCTTTGTCTTTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTTGGTGGCTGGCTGGCAGCTTCCTTTGGCATGAATATGGAAGCCATGTTCACCGGGGGGTAATGCTTTC  >  NC_000913/1724197‑1724488
                                                                                                                                                   |                                                                                                                                                
cATATTATCGCCGCCAACGTGGACCATCAGCGCTTTGTCTTTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCCGGc                                                                                                                                               >  1:1215708/1‑151 (MQ=255)
 atatTATCGCCGCCAACGTGGACCATCAGCGCTTTGTCTTTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCCGGcc                                                                                                                                              <  2:663824/151‑1 (MQ=255)
 atatTATCGCCGCCAACGTGGACCATCAGCGCTTTGTCTTTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCCGGcc                                                                                                                                              >  2:983610/1‑151 (MQ=255)
 atatTATCGCCGCCAACGTGGACCATCAGCGCTTTGTCTTTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCCGGcc                                                                                                                                              <  2:982452/151‑1 (MQ=255)
    ttatCGCCGCCAACGTGGACCATCAGCGCTTTGTCTTTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCCGGCCCtt                                                                                                                                           >  1:1663229/1‑151 (MQ=255)
      atCGCCGCCAACGTGGACCATCAGCGCTTTGTCTTTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCCGGCCCTTc                                                                                                                                          <  1:193255/150‑1 (MQ=255)
        cgccgcCAACGTGGACCATCAGCGCTTTGTCTTTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCCGGCCCTTCATg                                                                                                                                       >  1:1023541/1‑151 (MQ=255)
        cgccgcCAACGTGGACCATCAGCGCTTTGTCTTTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCCGGCCCTTCATg                                                                                                                                       <  2:636741/151‑1 (MQ=255)
         gccgccAACGTGGACCATCAGCGCTTTGTCTTTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCCGGCCCTTCATGt                                                                                                                                      <  1:1160442/151‑1 (MQ=255)
           cgccAACGTGGACCATCAGCGCTTTGTCTTTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCCGGCCCTTCATGttt                                                                                                                                    <  1:1226112/151‑1 (MQ=255)
              cAACGTGGACCATCAGCGCTTTGTCTTTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCCGGCCCTTCATGTTTAcc                                                                                                                                 >  2:774943/1‑151 (MQ=255)
               aaCGTGGACCATCAGCGCTTTGTCTTTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCCGGCCCTTCATGTTTAcct                                                                                                                                <  1:1358046/151‑2 (MQ=255)
               aaCGTGGACCATCAGCGCTTTGTCTTTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCCGGCCCTTCATGTTTACCg                                                                                                                                >  2:1204166/1‑151 (MQ=255)
               aaCGTGGACCATCAGCGCTTTGTCTTTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCCGGCCCTTCATGTTTACCg                                                                                                                                <  1:172648/151‑1 (MQ=255)
                aCGTGGACCATCAGCGCTTTGTCTTTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCCGGCCCTTCATGTTTACCgg                                                                                                                               <  2:1149909/151‑1 (MQ=255)
                   tGGACCATCAGCGCTTTGTCTTTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGAGCAGGCAGATCGCCTAAATGCCCGGCACCTTCCGGCCCTTCATGTTTACCGGTa                                                                                                                             >  2:1265218/1‑150 (MQ=255)
                     gACCATCAGCGCTTTGTCTTTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCCGGCCCTTCATGTTTACCGGTAttt                                                                                                                          <  2:885725/151‑1 (MQ=255)
                         aTCAGCGCTTTGTCTTTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCCGGCCCTTCATGTTTACCGGTATTTtgtg                                                                                                                      <  2:730919/151‑1 (MQ=255)
                           cAGCGCTTTGTCTTTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCCGGCCCTTCATGTTTACCGGTATTTTGTgg                                                                                                                     >  1:580962/1‑150 (MQ=255)
                            aGCGCTTTGTCTTTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCCGGCCCTTCATGTTTACCGGTATTTTGTGGAt                                                                                                                   >  2:238978/1‑151 (MQ=255)
                              cgcTTTGTCTTTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCCGGCCCTTCATGTTTACCGGTATTTTGTGGATc                                                                                                                  >  2:168117/1‑150 (MQ=255)
                              cgcTTTGTCTTTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCCGGCCCTTCATGTTTACCGGTATTTTGTGGATCa                                                                                                                 <  2:21887/151‑1 (MQ=255)
                                 tttGTCTTTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTCTTGACGACCAGTGCAGGCAGATCGCCTACATGCCCGGCACCTTCCGGACCTTCCTGTTTTACGGTATTTTGTGGCTCaa                                                                                                                >  2:474899/1‑149 (MQ=255)
                                  ttGTCTTTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCCGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAg                                                                                                               <  1:353519/149‑1 (MQ=255)
                                  ttGTCTTTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCCGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGa                                                                                                              <  1:538906/150‑1 (MQ=255)
                                   tGTCTTTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCCGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGAt                                                                                                             <  1:776858/150‑1 (MQ=255)
                                   tGTCTTTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCCGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGAt                                                                                                             <  2:1017595/150‑1 (MQ=255)
                                     tCTTTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCCGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGcc                                                                                                          >  1:698956/1‑151 (MQ=255)
                                      cTTTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCTTTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCCGGCCATTCATGTTTACCGGTATTTTGTGGATCAAGATGccc                                                                                                         >  2:43040/1‑151 (MQ=255)
                                        ttGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCCGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGc                                                                                                       >  1:1101909/1‑151 (MQ=255)
                                          gATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCCGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGcc                                                                                                      <  2:707703/150‑1 (MQ=255)
                                          gATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCCGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCt                                                                                                     <  2:1703053/151‑1 (MQ=255)
                                           aTTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCCGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTg                                                                                                    <  2:1039557/151‑1 (MQ=255)
                                           aTTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCCGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTg                                                                                                    >  2:62695/1‑151 (MQ=255)
                                                    aGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCCGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCcg                                                                                           <  2:224907/151‑1 (MQ=255)
                                                        aTTTCAGACGAGGTGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCTCCTTCCGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCgcgg                                                                                        <  1:1134164/150‑1 (MQ=255)
                                                          ttCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCCGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCGCGgcgc                                                                                     >  1:1116962/1‑151 (MQ=255)
                                                          ttCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCCGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCGCGgcgc                                                                                     >  1:127782/1‑151 (MQ=255)
                                                            cAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGCCGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCCGGCCCTGCATGTTTCCCGGTATTTTGTCGATCAAGATGCCCGCCTGCGGATTCCGCGGCGCt                                                                                    <  1:1275090/150‑1 (MQ=255)
                                                            cAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCCGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCGCGGCGCTg                                                                                   >  1:12222/1‑151 (MQ=255)
                                                               aCGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCCGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCGCGGCGCTGGc                                                                                 >  1:715064/1‑150 (MQ=255)
                                                                cGAGGCGCGATGACGGCATCGTTAGCTTTGCCGTCATTATTGCCGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCCGGCCCTTCATGTTTACCGGTATTTTGTGGAGCAAGATGCCCGCCTGCGGATTCCGCGGCGCTGGCt                                                                                <  2:1709864/150‑1 (MQ=255)
                                                                cGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCCGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCGCGGCGCTGGCtt                                                                               >  2:1525872/1‑151 (MQ=255)
                                                                cGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCCGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCGCGGCGCTGGCtt                                                                               <  2:971979/151‑1 (MQ=255)
                                                                 gAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCCGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCGCGGCGCTGGCtt                                                                               <  1:92856/150‑1 (MQ=255)
                                                                   ggCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCCGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCGCGGCGCTGGCTTTg                                                                             <  1:1612868/150‑1 (MQ=255)
                                                                   ggCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCCGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCGCGGCGCTGGCTTTGc                                                                            <  1:846908/151‑1 (MQ=255)
                                                                    gcgcGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCCGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCGCGGCGCTGGCTTTGcc                                                                           >  2:1250697/1‑151 (MQ=255)
                                                                       cgATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGCCCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCCGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCGCGGCGCTGGCTTTGCCATc                                                                        <  2:852535/151‑1 (MQ=255)
                                                                        gATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCCGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCGCGGCGCTGGCTTTGCCATCt                                                                       <  1:1303302/151‑1 (MQ=255)
                                                                           gACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCCGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCGCGGCGCTGGCTTTGCCATCttt                                                                     <  2:28854/150‑1 (MQ=255)
                                                                           gACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCCGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCGCGGCGCTGGCTTTGCCATCttt                                                                     <  1:925524/150‑1 (MQ=255)
                                                                            aCGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCCGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTtgg                                                                   >  2:1347744/1‑151 (MQ=255)
                                                                              ggCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCCGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTtggtg                                                                 >  2:1235556/1‑151 (MQ=255)
                                                                               gCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCCGGCCCTTCATGTTTCCCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTtggtgg                                                                <  2:626482/151‑1 (MQ=255)
                                                                                cATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCCGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTTGGtggc                                                               >  2:1629780/1‑151 (MQ=255)
                                                                                      tAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCCGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTTGGtggctggct                                                          <  2:252442/150‑1 (MQ=255)
                                                                                       aGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCCGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTTGGtggctggctgg                                                        >  1:1616555/1‑151 (MQ=255)
                                                                                          tttGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCTCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTTGGtggctggctggc                                                       >  2:1500294/1‑149 (MQ=255)
                                                                                          tttGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCCGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCGCGGGGCTGGCATTGCCATCTTTGGtgcctggctggc                                                       >  1:1455607/1‑149 (MQ=255)
                                                                                          tttGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCCGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTTGGtggctggctggc                                                       >  1:287543/1‑149 (MQ=255)
                                                                                                 tCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCCGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTTGGTGGCTGGCTGGCAGCTTCCtt                                              >  2:1428679/1‑151 (MQ=255)
                                                                                                 tCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCCGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTTGGTGGCTGGCTGGCAGCTTCCt                                               <  1:1043139/150‑1 (MQ=255)
                                                                                                   attattGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCCGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTTGGTGGCTGGCTGGCAGCTTCCTTTg                                            >  1:182608/1‑151 (MQ=255)
                                                                                                         gacgacCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCCGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTTGGTGGCTGGCTGGCAGCTTCCTTTGGCATg                                       >  1:347919/1‑150 (MQ=255)
                                                                                                         gacgacCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCCGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTTGGTGGCTGGCTGGCAGCTTCCTTTGGCATGa                                      <  1:394395/151‑1 (MQ=255)
                                                                                                           cgacCAGTGCAGGCAGATCGCCTAAATGCCCGGCTCCTCCCGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTTGGTGGCTGGCTGGCAGCTTCCTTTGGCATGaa                                     <  2:1101994/150‑1 (MQ=255)
                                                                                                           cgacCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCCGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTTGGTGGCTGGCTGGCAGCTTCCTTTGGCATGAat                                    >  1:424577/1‑151 (MQ=255)
                                                                                                            gacCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCCGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTTGGTGGCTGGCTGGCAGCTTCCTTTGGCATGAat                                    >  1:1031246/1‑150 (MQ=255)
                                                                                                            gacCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCCGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTTGGTGGCTGGCTGGCAGCTTCCTTTGGCATGAat                                    >  2:165656/1‑150 (MQ=255)
                                                                                                             acCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCCGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTTGGTGGCTGGCTGGCAGCTTCCTTTGGCATGAatat                                  >  2:1646527/1‑151 (MQ=255)
                                                                                                              ccAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCCGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTTGGTGGCTGGCTGGCAGCTTCCTTTGGCATGAATATg                                 >  2:253071/1‑151 (MQ=255)
                                                                                                                 gTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCCGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTTGGTGGCTGGCTGGCAGCTTCCTTTGGCATGAATATGGaa                              <  1:208580/151‑1 (MQ=255)
                                                                                                                 gTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCCGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTTGGTGGCTGGCTGGCAGCTTCCTTTGGCATGAATATGGa                               <  1:844774/150‑1 (MQ=255)
                                                                                                                   gcaggcagATCGCCTAAATGCCCGGCACCTTCCGGCCTTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTTGGTGGCTGGCTGGCAGCTTCCTTTGGCATGAATATGGAAGc                            >  2:1596111/1‑151 (MQ=255)
                                                                                                                    caggcagATCGCCTAAATGCCCGGCACCTTCCGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTTGGTGGCTGGCTGGCAGCTTCCTTTGGCATGAATATGGAAGcc                           >  2:975181/1‑151 (MQ=255)
                                                                                                                     aggcagATCGCCTAAATGCCCGGCACCTTCCGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTTGGTGGCTGGCTGGCAGCTTCCTTTGGCATGAATATGGAAGcc                           >  1:644868/1‑150 (MQ=255)
                                                                                                                       gcagATCGCCTAAATGCCCGGCACCTTCCGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTTGGTGGCTGGCTGGCAGCTTCCTTTGGCATGAATATGGAAGCCATg                        >  1:507449/1‑151 (MQ=255)
                                                                                                                       gcagATCGCCTAAATGCCCGGCACCTTCCGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTTGGTGGCTGGCTGGCAGCTTCCTTTGGCATGAATATGGAAGCCATg                        <  1:836682/151‑1 (MQ=255)
                                                                                                                        cagATCGCCTAAATGCCCGGCACCTTCCGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTTGGTGGCTGGCTGGCAGCTTCCTTTGGCATGAATATGGAAGCCATg                        <  2:1079951/150‑1 (MQ=255)
                                                                                                                           accGCCTAAATGCCCGGCCCCTTCCGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTTGGTGGCTGGCTGGCAGCTTCCTTTGGCATGAATATGGAAGCCATGTTc                     <  2:750445/148‑1 (MQ=255)
                                                                                                                             cGCCTAAATGCCCGGCACCTTCCGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTTGGTGGCTGGCTGGCAGCTTCCTTTGGCATGAATATGGAAGCCATGTTCAcc                  >  1:40747/1‑151 (MQ=255)
                                                                                                                              gCCTAAATGCCCGGCACCTTCCGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTTGGTGGCTGGCTGGCAGCTTCCTTTGGCATGAATATGGAAGCCATGTTCACCg                 >  1:192364/1‑151 (MQ=255)
                                                                                                                               ccTAAATGCCCGGCCCCTTCCGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTTGGTGGCTGGCTGGCAGCTTCCTTTGGCATGAATATGGAAGCCATGTTCACCgg                <  1:774879/151‑1 (MQ=255)
                                                                                                                               ccTAAATGCCCGGCACCTTCCGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTTGGTGGCTGGCTGGCAGCTTCCTTTGGCATGAATATGGAAGCCATGTTCACCgg                <  2:12225/151‑1 (MQ=255)
                                                                                                                               ccTAAATGCCCGGCACCTTCCGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTTGGTGGCTGGCTGGCAGCTTCCTTTGGCATGAATATGGAAGCCATGTTCACCg                 <  2:121624/150‑1 (MQ=255)
                                                                                                                                       cccGGCCCCTTCCGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTTGGTGGCTGGCTGGCAGCTTCCTTTGGCATGAATATGGAAGCCATGTTCACCGGGGGGTa          <  2:1115525/149‑1 (MQ=255)
                                                                                                                                       cccGGCACCTTCCGGCGCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTTGGTGGCTGGCTGGCAGCTTCCTTTGGCATGAATATGGAAGCCATGTTCACCGGGGGGTAAt        >  2:287836/1‑151 (MQ=255)
                                                                                                                                       cccGGCACCTTCCGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTTGGTGGCTGGCTGGCAGCTTCCTTTGGCATGAATATGGAAGCCATGTTCACCGGGGGGTa          <  2:1023618/149‑1 (MQ=255)
                                                                                                                                        ccGGCCCCTTCCGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTTGGTGGCTGGCTGGCAGCTTCCTTTGGCATGAATATGGAAGCCATGTTCACCGGGGGGTAATg       <  2:294900/151‑1 (MQ=255)
                                                                                                                                        ccGGCACCTTCCGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTTGGTGGCTGGCTGGCAGCTTCCTTTGGCATGAATATGGAAGCCATGTTCACCGGGGGGTAATg       <  2:1074080/151‑1 (MQ=255)
                                                                                                                                         cGGCACCTTCCGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTTGGTGGCTGGCTGGCAGCTTCCTTTGGCATGAATATGGAAGCCATGTTCACCGGGGGGTAATGc      >  2:693537/1‑151 (MQ=255)
                                                                                                                                          ggCACCTTCCGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTTGGTGGCTGGCTGGCAGCTTCCTTTGGCATGAATATGGAAGCCATGTTCACCGGGGGGTAATGCt     <  1:843576/151‑1 (MQ=255)
                                                                                                                                           gCACTTTCCGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTTGGTGGCTGGCTGGCAGCTTCCTTTGGCATGAATATGGAAGCCATGTTCACCGGGGGGTAATGCt     <  2:1326485/150‑1 (MQ=255)
                                                                                                                                             aCCTTCCGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTTGGTGGCTGGCTGGCAGCTTCCTTTGGCATGAATATGGAAGCCATGTTCACCGGGGGGTAATGCTTTc  >  1:810095/1‑151 (MQ=255)
                                                                                                                                              ccTTCCGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTTGGTGGCTGGCTGGCAGCTTCCTTTGGCATGAATATGGAAGCCATGTTCACCGGGGGGTAATGCTTTc  >  2:1048769/1‑150 (MQ=255)
                                                                                                                                              ccTTCCGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTTGGTGGCTGGCTGGCAGCTTCCTTTGGCATGAATATGGAAGCCATGTTCACCGGGGGGTAATGCTTTc  >  2:165990/1‑150 (MQ=255)
                                                                                                                                               cTTCCGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTTGGTGGCTGGCTGGCAGCTTCCTTTGGCATGAATATGGAAGCCATGTTCACCGGGGGGTAATGCTTTc  >  2:1455075/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                
CATATTATCGCCGCCAACGTGGACCATCAGCGCTTTGTCTTTGATTTCATCCAGTGATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTTGGTGGCTGGCTGGCAGCTTCCTTTGGCATGAATATGGAAGCCATGTTCACCGGGGGGTAATGCTTTC  >  NC_000913/1724197‑1724488

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: