Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,147,692 A→G 5.5% A24A (GCA→GCG) ‡ plsX → putative phosphate acyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,147,6920AG 5.5% 277.0 / 4.6 90A24A (GCA→GCG) ‡plsXputative phosphate acyltransferase
Reads supporting (aligned to +/- strand):  ref base A (44/41);  new base G (3/2);  total (47/43)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.88e-01

TCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGTTAGATGTCATGGGAGGGGATTTTGGCCCTTCCGTGACAGTGCCTGCAGCATTGCAGGCACTGAATTCTAATTCGCAACTCACTCTTCTTTTAGTCGGCAATTCCGACGCCATCACGCCATTACTTGCTAAAGCTGACTTTGAACAACGTTCGCGTCTGCAGATTATTCCTGCGCAGTCAGTTATCGCCAGTGATGC  >  NC_000913/1147548‑1147838
                                                                                                                                                |                                                                                                                                                  
tCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGTTAGATGTCATGGGAGGGGATTTTGGCCCTTCCGTGACAGTGCCTGCAGCATTGCa                                                                                                                                               >  2:1213787/1‑150 (MQ=255)
 cTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGTTAGATGTCATGGGAGGGGATTTTGGCCCTTCCGTGACAGTGCCTGCAGCATTGCAg                                                                                                                                              <  1:767970/150‑1 (MQ=255)
 cTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGTTAGATGTCATGGGAGGGGATTTTGGCCCTTCCGTGACAGTGCCTGCAGCATTGCAg                                                                                                                                              <  1:279403/150‑1 (MQ=255)
   gCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGTTAGATGTCATGGGAGGGGATTTTGGCCCTTCCGTGACAGTGCCTGCAGCATTGCAGGc                                                                                                                                            >  1:1324304/1‑150 (MQ=255)
   gCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGTTAGATGTCATGGGAGGGGATTTTGGCCCTTCCGTGACAGTGCCTGCAGCATTGCAGGc                                                                                                                                            <  2:155323/150‑1 (MQ=255)
    cGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGTTAGATGTCATGGGAGGGGATTTTGGCCCTTCCGTGACAGTGCCTGCAGCATTGCAGGCa                                                                                                                                           <  1:550267/150‑1 (MQ=255)
    cGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGTTAGATGTCATGGGAGGGGATTTTGGCCCTTCCGTGACAGTGCCTGCAGCATTGCAGGCa                                                                                                                                           >  1:1268710/1‑150 (MQ=255)
      tgatgaAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCTGGCGGCGACGATACCTTGACACGTCTAACCCTGTCGTTAGATGTCATGGGAGGGGATTTTGGCCCTTCCTTGACAGTGCCTGCAATGCTGCAGGCACt                                                                                                                                         <  1:596096/150‑1 (MQ=255)
       gatgaAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGTTAGATGTCATGGGAGGGGATTTTGGCCCTTCCGTGACAGTGCCTGCAATGCTGCAGGCACTg                                                                                                                                        <  1:364412/150‑1 (MQ=255)
            aggttAGTGATGATTTTCGCGAGGCTACTGGGTTAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGTTAGATGTCATGGGAGGGGATTTTGGCCCTTCCGTGACAGTGCCTGCAGCATTGCAGGCACTGAATTc                                                                                                                                   <  1:1059992/147‑1 (MQ=255)
            aGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGTTAGATGTCATGGGAGGGGATTTTGGCCCTTCCGTGACAGTGCCTGCAGCATTGCAGGc                                                                                                                                            >  2:199334/1‑141 (MQ=255)
            aGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGTTAGATGTCATGGGAGGGGATTTTGGCCCTTCCGTGACAGTGCCTGCAGCATTGCAGGc                                                                                                                                            <  1:199334/141‑1 (MQ=255)
             gCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGTTAGATGTCATGGGAGGGGATTTTGGCCCTTCCGTGACAGTGCCTGCAGCATTGCAGGCACTGAATTCt                                                                                                                                  <  1:865111/150‑1 (MQ=255)
                tAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGTTAGATGTCATGGGAGGGGATTTTGGCCCTTCCGTGACAGTGCCTGCAGCATTGCAGGCACTGAATTCTAAt                                                                                                                               >  1:1229337/1‑150 (MQ=255)
                 aGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGTTAGATGTCATGGGAGGGGATTTTGGCCCTTCCGTGACAGTGCCTGCAGCATTGCAGGCACTGAATTCTAAtt                                                                                                                              >  1:673069/1‑150 (MQ=255)
                  gTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGTTAGAAGTCATGGGAGGGGATTTTGGCCCTTCCGTGACAGTGCCTGCAGCATTGCAGGCACTGAATTCTAATTc                                                                                                                             >  2:1189716/1‑150 (MQ=255)
                   tGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGTTAGATGTCATGGGAGGGGATTTTGGCCCTTCCGTGACAGTGCCTGCAGCATTGCAGGCACTGAATTCTAATTCg                                                                                                                            <  2:598719/150‑1 (MQ=255)
                    gAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGTTAGATGTCATGGGAGGGGATTTTGGCCCTTCCGTGACAGTGCCTGCAGCATTGCAGGCACTGAATTCTAATTCGc                                                                                                                           <  1:183555/150‑1 (MQ=255)
                    gAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGTTAGATGTCATGGGAGGGGATTTTGGCCCTTCCGTGACAGTGCCTGCAGCATTGCAGGCACTGAATTCTAATTCGc                                                                                                                           >  2:1249248/1‑150 (MQ=255)
                    gAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGTTAGATGTCATGGGAGGGGATTTTGGCCCTTCCGTGACAGTGCCTGCAGCATTGCAGGCACTGAATTCTAATTCGc                                                                                                                           >  2:1060760/1‑150 (MQ=255)
                      ggATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGTTAGATGTCATGGGAGGGGATTTTGGCCCTTCCGTGACAGTGCCTGCAGCATTGCAGGCACTGAATTCTAATTCGCaa                                                                                                                         >  1:521076/1‑150 (MQ=255)
                       gATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGTTAGATGTCATGGGAGGGGATTTTGGCCCTTCCGTGACAGTGCCTGCAGCATTGCAGGCACTGAATTCTAATTCGCAac                                                                                                                        <  2:1063476/150‑1 (MQ=255)
                       gATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGTTAGATGTCATGGGAGGGGATTTTGGCCCTTCCGTGACAGTGCCTGCAGCATTGCAGGCACTGAATTCTAATTCGCAac                                                                                                                        <  1:1213787/150‑1 (MQ=255)
                           ttCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGTTAGATGTCATGGGAGGGGATTTTGGCCCTTCCGTGACAGTGCCTGCAGCATTGCAGGCACTGAATTCTAATTCGCAactcac                                                                                                                    <  1:1237096/150‑1 (MQ=255)
                            tCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGTTAGATGTCATGGGAGGGGATTTTGGCCCTTCCGTGACAGTGCCTGCAGCATTGCAGGCACTGAATTCTa                                                                                                                                 >  2:1192210/1‑136 (MQ=255)
                            tCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGTTAGATGTCATGGGAGGGGATTTTGGCCCTTCCGTGACAGTGCCTGCAGCATTGCAGGCACTGAATTCTa                                                                                                                                 <  1:1192210/136‑1 (MQ=255)
                             ccccAGGCAGCTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGTTAGATGTCATGGGAGGGGATTTTGGCCCTTCCGTGTCAGTGCCTGCAGCATTGCAGGCACTGAATTCTAATTCGCAactcactc                                                                                                                  >  2:1034542/1‑150 (MQ=255)
                              cccAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGTTAGATGTCATGGGAGGGGATTTTGGCCCTTCCGTGACAGTGCCTGCAGCATTGCAGGCACTGAATTCTAATTCGCAACTCActct                                                                                                                 >  2:383365/1‑150 (MQ=255)
                                cAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGTTAGATGTCATGGGAGGGGATTTTGGCCCTTCCGTGACAGTGCCTGCAGCATTGCAGGCACTGAATTCTAATTCGCAACTCACtcttc                                                                                                               >  2:349454/1‑150 (MQ=255)
                                 aGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGTTAGATGTCATGGGAGGGGATTTTGGCCCTTCCGTGACAGTGCCTGCAGCATTGCAGGCACTGAATTCTAATTCGCAACTCACtcttct                                                                                                              <  2:778524/150‑1 (MQ=255)
                                   gCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGTTAGATGTCATGGGAGGGGATTTTGGCCCTTCCGTGACAGTGCCTGCAGCATTGCAGGCACTGAATTCTAATTCGCAACTCACTCTTCttt                                                                                                            <  1:922501/150‑1 (MQ=255)
                                    cAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGTTAGATGTCATGGGAGGGGATTTTGGCCCTTCCGTGACAGTGCCTGCAGCATTGCAGGCACTGAATTCTAATTCGCAACTCACTCTTCtttt                                                                                                           <  1:1065443/150‑1 (MQ=255)
                                    cAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGTTAGATGTCATGGGAGGGGATTTTGGCCCTTCCGTGACAGTGCCTGCAGCATTGCAGGCACTGAATTCTAATTCGCAACTCACTCTTCtttt                                                                                                           <  1:579330/150‑1 (MQ=255)
                                    cAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGTTAGATGTCATGGGAGGGGATTTTGGCCCTTCCGTGACAGTGCCTGCAGCATTGCAGGCACTGAATTCTAATTCGCAACTCACTCTTCtttt                                                                                                           <  1:1013083/150‑1 (MQ=255)
                                     aaCTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGTTAGATGTCATGGGAGGGGATTTTGGCCCTTCCGTGACAGTGCCTGCAGCATTGCAGGCACTGAATTCTAATTCGCAACTCACTCTTCTTTTa                                                                                                          <  2:673069/150‑1 (MQ=255)
                                        tGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGTTAGATGTCATGGGAGGGGATTTTGGCCCTTCCGTGACAGTGCCTGCAGCATTGCAGGCACTGAATTCTAATTCGCAACTCACTCTTCTTTTAGTc                                                                                                       >  1:1259733/1‑150 (MQ=255)
                                        tGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGTTAGATGTCATGGGAGGGGATTTTGGCCCTTCCGTGACAGTGCCTGCAGCATTGCAGGCACTGAATTCTAATTCGCAACTCACTCTTCTTTTAGTc                                                                                                       <  1:280552/150‑1 (MQ=255)
                                         ggggAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGTTAGATGTCATGGGAGGGGATTTTGGCCCTTCCGTGACAGTGCCTGCAGCATTGCAGGCACTGAATTCTAATTCGCAACTCACTCTTCTTTTAGTCg                                                                                                      >  1:485761/1‑150 (MQ=255)
                                           ggAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGTTAGATGTCATGGGAGGGGATTTTGGCCCTTCCGTGTCAGTGCCTGCAGCATTGCAGGCACTGAATTCTAATTCGCAACTCACTCTTCTTTTAGTCGGc                                                                                                    <  1:349454/150‑1 (MQ=255)
                                                           gcggcgACGATACCTTGACACGTCTAACCCTGGCGTTAGATGTCATGGGAGGGGATTTTGGCCCTTCCGTGACAGTGCCTGCAGCATTGCAGGCACTGAATTCTAATTCGCAACTCACTCTTCTTTTAGTCGGCAATTCCGACGCCATCa                                                                                    <  1:528621/150‑1 (MQ=255)
                                                           gcggcgACGATACCTTGACACGTCTAACCCTGGCGTTAGATGTCATGGGAGGGGATTTTGGCCCTTCCGTGACAGTGCCTGCAGCATTGCAGGCACTGAATTCTAATTCGCAACTCACTCTTCTTTTAGTCGGCAATTCCGACGCCATCa                                                                                    >  2:1412801/1‑150 (MQ=255)
                                                             ggcgACGATACCTTGACACGTCTAACCCTGGCGTTAGATGTCATGGGAGGGGATTTTGGCCCTTCCGTGTCAGTGCCTGCAGCATTGCAGGCACTGAATTCTAATTCGCAACTCACTCTTCTTTTAGTCGGCAATTCCGACGCCATCACg                                                                                  <  1:1034542/150‑1 (MQ=255)
                                                                 acgaTACCTTGACACGTCTAACCCTGGCGTTAGATGTCATGGGAGGGGATTTTGGCCCTTCCGTGACAGTGCCTGCAGCATTGCAGGCACTGAATTCTAATTCGCAACTCACTCTTCTTTTAGTCGGCAATTCCGACGCCATCACGCCAt                                                                              <  2:1151915/150‑1 (MQ=255)
                                                                  cgaTACCTTGACACGTCTAACCCTGGCGTTAGATGTCATGGGAGGGGATTTTGGCCCTTCCGTGACGGTGCCTGCAGCATTGCAGGCACTGAATTCTAATTCGCAACACACTCTTCTTTTAGTCGGCAATTCCGACGCCATCACGCCAtt                                                                             <  2:1030886/150‑1 (MQ=255)
                                                                       ccTTGACACGTCTAACCCTGGCGTTAGATGTCATGGGAGGGGATTTTGGCCCTTCCGTGACAGTGCCTGCAGCATTGCAGGCACTGAATTCTAATTCGCAACTCACTCTTCTTTTAGTCGGCAATTCCGACGCCATCACGCCATTACTTg                                                                        >  1:1327678/1‑150 (MQ=255)
                                                                       ccTTGACACGTCTAACCCTGGCGTTAGATGTCATGGGAGGGGATTTTGGCCCTTCCGTGACAGTGCCTGCAGCATTGCAGGCACTGAATTCTAATTCGCAACTCACTCTTCTTTTAGTCGGCAATTCCGACGCCATCACGCCATTACTTg                                                                        <  1:1060760/150‑1 (MQ=255)
                                                                       ccTTGACACGTCTAACCCTGGCGTTAGATGTCATGGGAGGGGATTTTGGCCCGTCCGTGACAGTGCCTGCAGCATTGCAGGCACTGAAGTCTAATTCGCAACTCACTCGTCTTTTAGTCGGCAATTCCGCCGCCAGCACGCCCTTACtgg                                                                        >  2:1278929/1‑148 (MQ=255)
                                                                          tGACACGTCTAACCCTGGCGTTAGATGTCATGGGAGGGGATTTTGGCCCTTCCGTGACAGTGCCTGCAGCATTGCAGGCACTGAATTCTAATTCGCAACTCACTCTTCTTTTAGTCGGCAATTCCGACGCCATCACGCCATTACTTGCTa                                                                     <  2:485761/150‑1 (MQ=255)
                                                                             cacGTCTAACCCTGGCGTTAGATGTCATGGGAGGGGATTTTGGCCCTTCCGTGACAGTGCCTGCAGCATTGCAGGCACTGAATTCTAATTCGCAACTCACTCTTCTTTTAGTCGGCAATTCCGACGCCATCACGCCATTACTTGCTAAAg                                                                  <  1:1412801/150‑1 (MQ=255)
                                                                              acGTCTAACCCTGGCGTTAGATGTCATGGGAGGGGATTTTGGCCCTTCCGTGACAGTGCCTGCAGCATTGCAGGCACTGAATTCTAATTCGCAACTCACTCTTCTTTTAGTCGGCAATTCCGACGCCATCACGCCATTACTTGCTAAAGc                                                                 <  1:383365/150‑1 (MQ=255)
                                                                                  cTAACCCTGGCGTTAGATGTCATGGGAGGGGATTTTGGCCCTTCCGTGACAGTGCCTGCAGCATTGCAGGCACTGAATTCTAATTCGCAACTCACTCTTCTTTTAGTCGGCAATTCCGACGCCATCACGCCATTACTTGCTAAAGCTGAc                                                             >  1:1086562/1‑150 (MQ=255)
                                                                                   tAACCCTGGCGTTAGATGTCATGGGAGGGGATTTTGGCCCTTCCGTGACAGTGCCTGCAGCATTGCAGGCACTGAATTCTAATTCGCAACTCACTCTTCTTTTAGTCGGCAATTCCGACGCCATCACGCCATTACTTGCTAAAGCTGACt                                                            >  1:1109187/1‑150 (MQ=255)
                                                                                     aCCCTGGCGTTAGATGTCATGGGAGGGGATTTTGGCCCTTCCGTGACAGTGCCTGCAGCATTGCAGGCACTGAATTCTAATTCGCAACTCACTCTTCTTTTAGTCGGCAATTCCGACGCCATCACGCCATTACTTGCTAAAGCTGACttt                                                          >  2:799575/1‑150 (MQ=255)
                                                                                      cccTGGCGTTAGATGTCATGGGAGGGGATTTTGGCCCTTCCGTGACAGTGCCTGCAGCATTGCAGGCACTGAATTCTAATTCGCAACTCACTCTTCTTTTAGTCGGCAATTCCGACGCCATCACGCCATTACTTGCTAAAGCTGACTTTg                                                         >  1:1265645/1‑150 (MQ=255)
                                                                                       ccTGGCGTTAGATGTCATGGGAGGGGATTTTGGCCCTTCCGTGACAGTGCCTGCAGCATTGCAGGCACTGAATTCTAATTCGCAACTCACTCTTCTTTTAGTCGGCAATTCCGACGCCATCACGCCATTACTTGCTAAAGCTGACTTTGa                                                        <  2:875635/150‑1 (MQ=255)
                                                                                        cTGGCGTTAGATGTCATGGGAGGGGATTTTGGCCCTTCCGTGACAGTGCCTGCAGCATTGCAGGCACTGAATTCTAATTCGCAACTCACTCTTCTTTTAGTCGGCAATTCCGACGCCATCACGCCATTACTTGCTAAAGCTGACTTTGaa                                                       <  2:1327678/150‑1 (MQ=255)
                                                                                         tGGCGTTAGATGTCATGGGAGGGGATTTTGGCCCTTCCGTGACAGTGCCTGCAGCATTGCAGGCACTGAATTCTAATTCGCAACTCACTCTTCTTTTAGTCGGCAATTCCGACGCCATCACGCCATTACTTGCTAAAGCTGACTTTGaac                                                      <  1:851413/150‑1 (MQ=255)
                                                                                           gCGTTAGATGTCATGGGAGGGGATTTTGGCCCTTCCGTGACAGTGCCTGCAGCATTGCAGGCACTGAATTCTAATTCGCAACTCACTCTTCTTTTAGTCGGCAATTCCGACGCCATCACGCCATTACTTGCTAAAGCTGACTTTGa                                                        >  1:1238173/1‑146 (MQ=255)
                                                                                           gCGTTAGATGTCATGGGAGGGGATTTTGGCCCTTCCGTGACAGTGCCTGCAGCATTGCAGGCACTGAATTCTAATTCGCAACTCACTCTTCTTTTAGTCGGCAATTCCGACGCCATCACGCCATTACTTGCTAAAGCTGACTTTGa                                                        <  2:1238173/146‑1 (MQ=255)
                                                                                             gTTAGATGTCATGGGAGGGGATTTTGGCCCTTCCGTGACAGTGCCTGCAGCATTGCAGGCACTGAATTCTAATTCGCAACTCACTCTTCTTTTAGTCGGCAATTCCGACGCCATCACGCCATTACTTGCTAAAGCTGACTTTGAACAACg                                                  >  1:1408020/1‑150 (MQ=255)
                                                                                               tAGATGTCATGGGAGGGGATTTTGGCCCTTCCGTGACAGTGCCTGCAGCATTGCAGGCACTGAATTCTAATTCGCAACTCACTCTTCTTTTAGTCGGc                                                                                                    <  1:638480/98‑1 (MQ=255)
                                                                                               tAGATGTCATGGGAGGGGATTTTGGCCCTTCCGTGACAGTGCCTGCAGCATTGCAGGCACTGAATTCTAATTCGCAACTCACTCTTCTTTTAGTCGGc                                                                                                    >  2:638480/1‑98 (MQ=255)
                                                                                                 gATGTCATGGGAGGGGATTTTGGCCCTTCCGTGACAGTGCCTGCAGCATTGCAGGCACTGAATTCTAATTCGCAACTCACTCTTCTTTTAGTCGGCAATTCCGACGCCATCACGCCATTACTTGCTAAAGCTGACTTTGAACAACGTTcg                                              >  1:979774/1‑150 (MQ=255)
                                                                                                   tGTCATGGGAGGGGATTTTGGCCCTTCCGTGACAGTGCCTGCAGCATTGCAGGCACTGAATTCTAATTCGCAACTCACTCTTCTTTTAGTCGGCAATTCCGACGCCATCACGCCATTACTTGCTAAAGCTGACTTTGAACAACGTTcgcg                                            >  1:414778/1‑150 (MQ=255)
                                                                                                      cATGGGAGGGGATTTTGGCCCTTCCGTGACAGTGCCTGCAGCATTGCAGGCACTGAATTCTAATTCGCAACTCACTCTTCTTTTAGTCGGCAATTCCGACGCCATCACGCCATTACTTGCTAAAGCTGACTTTGAACAACGTTCGCGTCt                                         >  1:672799/1‑150 (MQ=255)
                                                                                                      cATGGGAGGGGATTTTGGCCCTTCCGTGACAGTGCCTGCAGCATTGCAGGCACTGAATTCTAATTCGCAACTCACTCTTCTTTTAGTCGGCAATTCCGACGCCATCACGCCATTACTTGCTAAAGCTGACTTTGAACAACGTTCGCGTCt                                         <  1:628262/150‑1 (MQ=255)
                                                                                                       aTGGGAGGGGATTTTGGCCCTTCCGTGACAGTGCCTGCAGCATTGCAGGCACTGAATTCTAATTCGCAACTCACTCTTCTTTTAGTCGGCAATTCCGACGCCATCACGCCATTACTTGCTAAAGCTGACTTTGAACAACGTTCACGTCTg                                        >  2:193637/1‑150 (MQ=255)
                                                                                                         gggAGGGGATTTTGGCCCTTCCGTGACAGTGCCTGCAGCATTGCAGGCACTGAATTCTAATTCGCAACTCACTCTTCTTTTAGTCGGCAATTCCGACGCCATCACGCCATTACTTGCTAAAGCTGACTTTGAACAACGTTCGCGTCTGCa                                      <  1:1314686/150‑1 (MQ=255)
                                                                                                          ggAGGGGATTTTGGCCCTTCCGTGACAGTGCCTGCAGCATTGCAGGCACTGAATTCTAATTCGCAACTCACTCTTCTTTTAGTCGGCAATTCCGACGCCATCACGCCATTACTTGCTAAAGc                                                                 <  1:344683/122‑1 (MQ=255)
                                                                                                          ggAGGGGATTTTGGCCCTTCCGTGACAGTGCCTGCAGCATTGCAGGCACTGAATTCTAATTCGCAACTCACTCTTCTTTTAGTCGGCAATTCCGACGCCATCACGCCATTACTTGCTAAAGc                                                                 >  2:344683/1‑122 (MQ=255)
                                                                                                          ggAGGGGATTTTGGCCCTTCCGTGACAGTGCCTGCAGCATTGCAGGCACTGAATTCTAATTCGCAACTCACTCTTCTTTTAGTCGGCAAATCCGACGCCATCACGCCATTACTTGCTAAAGCTGACTTTGAACAACGTTCGCGTCTGCAg                                     >  1:1243314/1‑150 (MQ=255)
                                                                                                             ggggATTTTGGCCCTTCCGTGACAGTGCCTGCAGCATTGCAGGCACTGAATTCTAATTCGCAACTCACTCTTCTTTTAGTCGGCAATTCCGACGCCATCACGCCATTACTTGCTAAAGCTGACTTTGAACAACGTTCGCGTCTGCAGatt                                  >  2:534182/1‑150 (MQ=255)
                                                                                                              gggATTTTGGCCCTTCCGTGACAGTGCCTGCAGCATTGCAGGCACTGAATTCTAATTCGCAACTCACTCTTCTTTTAGTCGGCAATTCCGACGCCATCACGCCATTACTTGCTAAAGCTGACTTTGAACAACGTTCGCGTCTGCAGatta                                 >  2:405120/1‑150 (MQ=255)
                                                                                                              gggATTTTGGCCCTTCCGTGACAGTGCCTGCAGCATTGCAGGCACTGAATTCTAATTCGCAACTCACTCTTCTTTTAGTCGGCAATTCCGACGCCATCACGCCATTACTTGCTAAAGCTGACTTTGAACAACGTTCGCGTCTGCAGatta                                 >  1:1128853/1‑150 (MQ=255)
                                                                                                                    ttGGCCCTTCCGTGACAGTGCCTGCAGCATTGCAGGCACTGAATTCTAATTCGCAACTCACTCTTCTTTTAGTCGGCAATTCCGACGCCATCACGCCATTACTTGCTAAAGCTGACTTTGAACAACGTTCGCGTCTGCAGATTATTCCTg                           >  2:449637/1‑150 (MQ=255)
                                                                                                                      ggCCCTTCCGTGACAGTGCCTGCAGCATTGCAGGCACTGAATTCTAATTCGCAACTCACTCTTCTTTTAGTCGGCAATTCCGACGCCATCACGCCATTACTTGCTAAAGCTGACTTTGAACAACGTTCGCGTCTGCAGATTATTCCTgcg                         >  2:1073932/1‑150 (MQ=255)
                                                                                                                      ggCCCTTCCGTGACAGTGCCTGCAGCATTGCAGGCACTGAATTCTAATTCGCAACTCACTCTTCTTTTAGTCGGCAATTCCGACGCCATCACGCCATTACTTGCTAAAGCTGACTTTGAACAACGTTCGCGTCTGCAGATTATTCCTgcg                         <  2:218450/150‑1 (MQ=255)
                                                                                                                        cccTTCCGTGACAGTGCCTGCAGCATTGCAGGCACTGAATTCTAATTCGCAACTCACTCTTCTTTTAGTCGGCAATTCCGACGCCATCACGCCATTACTTGCTAAAGCTGACTTTGAACAACGTTCGCGTCTGCAGATTATTCCTGCGca                       >  2:183278/1‑150 (MQ=255)
                                                                                                                            tCCGTGACAGTGCCTGCAGCATTGCAGGCACTGAATTCTAATTCGCAACTCACTCTTCTTTTAGTCGGCAATTCCGACGCCATCACGCCATTACTTGCTAAAGCTGACTTTGAACAACGTTCGCGTCTGCAGATTATTCCTGCGcagtca                   <  1:351547/150‑1 (MQ=255)
                                                                                                                              gaattcAGTGCCTGCAATGCTGCAGGCACTGAATTCTAATTCGCAACTCACTCTTCTTTTAGTCGGCAATTCCGACGCCATCACGCCATTACTTGCTAAAGCTGACTTTGAACAACGTTCGCGTCTGCAGATTATTCCTGCGcagtcagt                 >  1:600402/6‑150 (MQ=255)
                                                                                                                                tGACAGTGCCTGCAGCATTGCAGGCACTGAATTCTAATTCGCAACTCACTCTTCTTTTAGTCGGCAATTCCGACGCCATCACGCCATTACTTGCTAAAGCTGACTTTGAACAACGTTCGCGTCTGCAGATTATTCCTGCGCAGTCAGTTa               >  1:993941/1‑150 (MQ=255)
                                                                                                                                attcAGTGCCTGCAATGCTGCAGGCACTGAATTCTAATTCGCAACTCACTCTTCTTTTAGTCGGCAATTCCGACGCCATCACGCCATTACTTGCTAAAGCTGACTTTGAACAACGTTCGCGTCTGCAGATTATTCCTGCGCAGTCAGTTa               >  1:712419/4‑150 (MQ=255)
                                                                                                                                 gACAGTGCCTGCAGCATTGCAGGCACTGAATTCTAATTCGCAACTCACTCTTCTTTTAGTCGGCAATTCCGACGCCATCACGCCATTACTTGCTAAAGCTGACTTTGAACAACGTTCGTGTCTGCAGATTATTCCTGCGCAGTCAGTTAt              >  1:614845/1‑150 (MQ=255)
                                                                                                                                    aGTGCCTGCAATGCTGCAGGCACTGAATTCTAATTCGCAACTCACTCTTCTTTTAGTCGGCAATTCCGACGCCATCACGCCATTACTTGCTAAAGCTGACTTTGAACAACGTTCGCGTCTGCAGATTATTCCTGCGCAGTCAGTTATCGc           >  1:594850/1‑150 (MQ=255)
                                                                                                                                     gTGCCTGCAGCATTGCAGGCACTGAATTCTAATTCGCAACTCACTCTTCTTTTAGTCGGCAATTCCGACGCCATCACGCCATTACTTGCTAAAGCTGACTTTGAACAACGTTCGCGTCTGCAGATTATTCCTGCGCAGTCAGTTATCGcc          <  2:1408020/150‑1 (MQ=255)
                                                                                                                                     gTGCCTGCAGCATTGCAGGCACTGAATTCTAATTCGCAACTCACTCTTCTTTTAGTCGGCAATTCCGACGCCATCACGCCATTACTTGCTAAAGCTGACTTTGAACAACGTTCGCGTCTGCAGATTATTCCTGCGCAGTCAGTTATCGcc          >  1:616410/1‑150 (MQ=255)
                                                                                                                                     gTGCCTGCAGCATTGCAGGCACTGAATTCTAATTCGCAACTCACTCTTCTTTTAGTCGGCAATTCCGACGCCATCACGCCATTACTTGCTAAAGCTGACTTTGAACAACGTTCGCGTCTGCAGATTATTCCTGCGCAGTCAGTTATCGcc          <  2:1401713/150‑1 (MQ=255)
                                                                                                                                       gCCTGCAGCATTGCAGGCACTGAATTCTAATTCGCAACTCACTCTTCTTTTAGTCGGCAATTCCGACGCCATCACGCCATTACTTGCTAAAGCTGACTTTGAACAACGTTCGCGTCTGCAGATTATTCCTGCGCAGTCAGTTATCGCCAg        >  2:56165/1‑150 (MQ=255)
                                                                                                                                         cTGCAGCATTGCAGGCACTGAATTCTAATTCGCAACTCACTCTTCTTTTAGTCGGCAATTCCGACGCCATCACGCCATTACTTGCTAAAGCTGACTTTGAACAACGTTCGCGTCTGCAGATTATTCCTGCGCAGTCAGTTATCGCCAGTg      >  2:669576/1‑150 (MQ=255)
                                                                                                                                         cTGCAGCATTGCAGGCACTGAATTCTAATTCGCAACTCACTCTTCTTTTAGTCGGCAATTCCGACGCCATCACGCCATTACTTGCTAAAGCTGACTTTGAACAACGTTCGCGTCTGCAGATTATTCCTGCGCAGTCAGTTATCGCCAGTg      >  2:329517/1‑150 (MQ=255)
                                                                                                                                             agcaTTGCAGGCACTGAATTCTAATTCGCAACTCACTCTTCTTTTAGTCGGCAATTCCGACGCCATCACGCCATTACTTGCTAAAGCTGACTTTGAACAACGTTCGCGTCTGCAGATTATTCCTGCGCAGTCAGTTATCGCCAGTGATGc  >  1:501713/1‑150 (MQ=255)
                                                                                                                                                |                                                                                                                                                  
TCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGTTAGATGTCATGGGAGGGGATTTTGGCCCTTCCGTGACAGTGCCTGCAGCATTGCAGGCACTGAATTCTAATTCGCAACTCACTCTTCTTTTAGTCGGCAATTCCGACGCCATCACGCCATTACTTGCTAAAGCTGACTTTGAACAACGTTCGCGTCTGCAGATTATTCCTGCGCAGTCAGTTATCGCCAGTGATGC  >  NC_000913/1147548‑1147838

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