Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,931,721 T→C S281P (TCC→CCC)  purT → phosphoribosylglycinamide formyltransferase 2

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,931,7210TC97.1% 236.1 / ‑1.0 71S281P (TCC→CCC) purTphosphoribosylglycinamide formyltransferase 2
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base C (27/41);  total (28/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 = 5.17e-01

TGAATCCTGGCAACCACAGCAAATGAGCCCGCTTGCCCTTGAACGTGCGCAGGAGATTGCCCGTAAAGTGGTGCTGGCACTGGGCGGTTATGGGTTGTTTGGTGTCGAGCTATTTGTCTGTGGTGATGAGGTGATTTTCAGTGAGGTCTCCCCTCGTCCACATGATACCGGGATGGTGACGTTAATTTCTCAAGATCTCTCAGAGTTTGCCCTGCATGTACGTGCCTTCCTCGGACTTCCGGTTGGCGGGATCCGTCAGTATGGTCCTGCAGCTTCTGCCGTTATTCTGCCACA  >  NC_000913/1931573‑1931866
                                                                                                                                                    |                                                                                                                                                 
tGAATCCTGGCAACCACAGCAAATGAGCCCGCTTGCCCTTGAACGTGCGCAGGAGATTGCCCGTAAAGTGGTGCTGGCACTGGGCGGTTATGGGTTGTTTGGTGTCGAGCTATTTGTCTGTGGTGATGAGGTGATTTTCAGTGAGGTcccc                                                                                                                                                 >  2:828430/1‑151 (MQ=255)
tGAATCCTGGCAACCACAGCAAATGAGCCCGCTTGCCCTTGAACGTGCGCAGGAGATTGCCCGTAAAGTGGTGCTGGCACTGGGCGGTTATGGGTTGTTTGGTGTCGAGCTATTTGTCTGTGGTGATGAGGTGATTTTCAGTGAGGTcccc                                                                                                                                                 >  1:950064/1‑151 (MQ=255)
  aaTCCTGGCAACCACAGCAAATGAGCCCGCTTGCCCTTGAACGTGCGCAGGAGATTGCCCGTAAAGTGGTGCTGGCACTGGGCGGTTATGGGTTGTTTGGTGTCGAGCTATTTGTCTGTGGTGATGAGGTGATTTTCAGTGAGGTcccccc                                                                                                                                               <  2:922160/151‑1 (MQ=255)
  aaTCCTGGCAACCACAGCAAATGAGCCCGCTTGCCCTTGAACGTGCGCAGGAGATTGCCCGTAAAGTGGTGCTGGCACTGGGCGGTTATGGGTTGTTTGGTGTCGAGCTATTTGTCTGTGGTGATGAGGTGATTTTCAGTGAGGTcccccc                                                                                                                                               >  1:312200/1‑151 (MQ=255)
  aaTCCTGGCAACCACAGCAAATGAGCCCGCTTGCCCTTGAACGTGCGCAGGAGATTGCCCGTAAAGTGGTGCTGGCACTGGGCGGTTATGGGTTGTTTGGTGTCGAGCTATTTGTCTGTGGTGATGAGGTGATTTTCAGTGAGGTcccccc                                                                                                                                               >  1:310493/1‑151 (MQ=255)
  aaTCCTGGCAACCACAGCAAATGAGCCCGCTTGCCCTTGAACGTGCGCAGGAGATTGCCCGTAAAGTGGTGCTGGCACTGGGCGGTTATGGGTTGTTTGGTGTCGAGCTATTTGTCTGTGGTGATGAGGTGATTTTCAGTGAGGTcccccc                                                                                                                                               >  1:1286971/1‑151 (MQ=255)
   aTCCTGGCAACCACAGCAAATGAGCCCGCTTGCCCTTGAACGTGCGCAGGAGATTGCCCGTAAAGTGGTGCTGGCACTGGGCGGTTATGGGTTGTTTGGTGTCGAGCTATTTGTCTGTGGTGATGAGGTGATTTTCAGTGAGGTCCCCCCt                                                                                                                                              <  2:900978/151‑1 (MQ=255)
        ggCAACCACAGCAAATGAGCCCGCTTGCCCTTGAACGTGCGCAGGAGATTGCCCGTAAAGTGGTGCTGGCACTTGGCGGTTATGGGTTGTTTGGTGTCGAGCTATTTGTCTGTGGTGATGCGGTGATTTTCAGTGAGGTCCCCCCTCGTcc                                                                                                                                         <  1:954686/151‑1 (MQ=255)
        ggCAACCACAGCAAATGAGCCCGCTTGCCCTTGAACGTGCGCAGGAGATTGCCCGTAAAGTGGTGCTGGCACTGGGCGGTTATGGGTTGTTTGGTGTCGAGCTATTTGTCTGTGGTGATGAGGTGATTTTCAGTGAGGTCCCCCCTCGTcc                                                                                                                                         <  2:1180977/151‑1 (MQ=255)
            aCCACAGCAAATGAGCCCGCTTGCCCTTGAACGTGCGCAGGAGATTGCCCGTAAAGTGGTGCTGGCACTGGGCGGTTATGGGTTGTTTGGTGTCGAGCTATTTGTCTGTGGTGATGAGGTGATTTTCAGTGAGGTCCCCCCTCGTCCACAt                                                                                                                                     <  1:1189255/151‑1 (MQ=255)
               acaGCAAATGAGCCCGCTTGCCCTTGAACGTGCGCAGGAGATTGCCCGTAAAGTGGTGCTGGCACTGGGCGGTTATGGGTTGTTTGGTGTCGAGCTATTTGTCTGTGGTGATGAGGTGATTTTCAGTGAGGTCCCCCCTCGTCCACATGAt                                                                                                                                  <  2:588769/151‑1 (MQ=255)
                  gCAAATGAGCCCGCTTGCCCTTGAACGTGCGCAGGAGATTGCCCGTAAAGTGGTGCTGGCACTGGGCGGTTATGGGTTGTTTGGTGTCGAGCTATTTGTCTGTGGTGATGAGGTGATTTTCAGTGAGGTCCCCCCTCGTCCACATGATAcc                                                                                                                               >  1:1006899/1‑151 (MQ=255)
                  gCAAATGAGCCCGCTTGCCCTTGAACGTGCGCAGGAGATTGCCCGTAAAGTGGTGCTGGCACTGGGCGGTTATGGGTTGTTTGGTGTCGAGCTATTTGTCTGTGGTGATGAGGTGATTTTCAGTGAGGTCCCCCCTCGTCCACATGATAcc                                                                                                                               >  2:1079828/1‑151 (MQ=255)
                      aTGAGCCCGCTTGCCCTTGAACGTGCGCAGGAGATTGCCCGTAAAGTGGTGCTGGCACTGGGCGGTTATGGGTTGTTTGGTGTCGAGCTATTTGTCTGTGGTGATGAGGTGATTTTCAGTGAGGTCCCCCCTCGTCCACATGATACCGGGa                                                                                                                           <  2:265143/151‑1 (MQ=255)
                      aTGAGCCCGCTTGCCCTTGAACGTGCGCAGGAGATTGCCCGTAAAGTGGTGCTGGCACTGGGCGGTTATGGGTTGTTTGGTGTCGAGCTATTTGTCTGTGGTGATGAGGTGATTTTCAGTGAGGTCCCCCCTCGTCCACATGATACCGGGa                                                                                                                           <  1:224386/151‑1 (MQ=255)
                      aTGAGCCCGCTTGCCCTTGAACGTGCGCAGGAGATTGCCCGTAAAGTGGCGCTGGCACTGGGCGGTTATGGGTTGTTTGGTGTCGCGCTATTTGTCTGTGGTGATGAGGTGATTTTCAGTGAGGTCCCCCCTCGTCCACATGATAACggg                                                                                                                            >  1:612319/1‑150 (MQ=255)
                        gAGCCCGCTTGCCCTTGAACGTGCGCAGGAGATTGCCCGTAAAGTGGTGCTGGCACTGGGCGGTTATGGGTTGTTTGGTGTCGAGCTATTTGTCTGTGGTGATGAGGTGATTTTCAGTGAGGTCCCCCCTCGTCCACATGATACCGGGAt                                                                                                                          <  1:1198098/150‑1 (MQ=255)
                        gAGCCCGCTTGCCCTTGAACGTGCGCAGGAGATTGCCCGTAAAGTGGTGCTGGCACTGGGCGGTTATGGGTTGTTTGGTGTCGAGCTATTTGTCTGTGGTGATGAGGTGATTTTCAGTGAGGTCCCCCCTCGTCCACATGATACCGGGATg                                                                                                                         >  2:144707/1‑151 (MQ=255)
                           cccGCTTGCCCTTGAACGTGCGCAGGAGATTGCCCGTAAAGTGGTGCTGGCACTGGGCGGTTATGGGTTGTTTGGTGTCGAGCTATTTGTCTGTGGTGATGAGGTGATTTTCAGTGAGGTCCCCCCTCGTCCACATGATACCGGGATGGTg                                                                                                                      >  1:1210185/1‑151 (MQ=255)
                             cGCTTGCCCTTGAACGTGCGCAGGAGATTGCCCGTAAAGTGGTGCTGGCACTGGGCGGTTATGGGTTGTTTGGTGTCGAGCTATTTGTCTGTGGTGATGAGGTGATTTTCAGTGAGGTCCCCCCTCGTCCACATGATACCGGGATGGTGAc                                                                                                                    <  1:1194195/151‑1 (MQ=255)
                              gCTTGCCCTTGAACGTGCGCAGGAGATTGCCCGTAAAGTGGTGCTGGCACTGGGCGGTTATGGGTTGTTTGGTGTCGAGCTATTTGTCTGTGGTGATGAGGTGATTTTCAGTGAGGTCCCCCCTCGTCCACATTATACCGGGATGGTGACg                                                                                                                   >  1:681651/1‑151 (MQ=255)
                               cTTGCCCTTGAACGTGCGCAGGAGATTGCCCGTAAAGTGGTGCTGGCACTGGGCGGTTATGGGTTGTTTGGTGTCGAGCTATTTGTCTGTGGTGATGAGGTGATTTTCAGTGAGGTCCCCCCTCGTCCACATGATACCGGGATGGTGACg                                                                                                                   <  1:1295758/150‑1 (MQ=255)
                                ttGCCCTTGAACGTGCGCAGGAGATTGCCCGTAAAGTGGTGCTGGCACTGGGCGGTTATGGGTTGTTTGGTGTCGAGCTATTTGTCTGTGGTGATGAGGTGATTTTCAGTGAGGTCCCCCCTCGTCCACATGATACCGGGATGGTGACGtt                                                                                                                 >  2:42521/1‑151 (MQ=255)
                                  gCCCTTGAACGTGCGCAGGAGATTGCCCGTAAAGTGGTGCTGGCACTGGGCGGTTATGGGTTGTTTGGTGTCGAGCTATTTGTCTGTGGTGATGAGGTGATTTTCAGTGAGGTCCCCCCTCGTCCACATGATACCGGGATGGTGACGTTa                                                                                                                >  1:1106720/1‑150 (MQ=255)
                                    ccTTGAACGTGCGCAGGAGATTGCCCGTAAAGTGGTGCTGGCACTGGGCGGTTATGGGTTGTTTGGTGTCGAGCTATTTGTCTGTGGTGATGAGGTGATTTTCAGTGAGGTCCCCCCTCGTCCACATGATACCGGGATGGTGACGTTAAtt                                                                                                             >  2:354224/1‑151 (MQ=255)
                                       tGAACGTGCGCAGGAGATTGCCCGTAAAGTGGTGCTGGCACTGGGCGGTTATGGGTTGTTTGGTGTCGAGCTATTTGTCTGTGGTGATGAGGTGATTTTCAGTGAGGTCCCCCCTCGTCCACATGATACCGGGATGGTGACGTTAATTtc                                                                                                           <  1:93971/150‑1 (MQ=255)
                                        gAACGTGCGCAGGAGATTGCCCGTAAAGTGGTGCTGGCACTGGGCGGTTATGGGTTGTTTGGTGTCGAGCTATTTGTCTGTGGTGATGAGGTGATTTTCAGTGAGGTCCCCCCTCGTCCACATGATACCGGGATGGTGACGTTAATTtct                                                                                                          <  1:108784/150‑1 (MQ=255)
                                        gAACGTGCGCAGGAGATTGCCCGTAAAGTGGTGCTGGCACTGGGCGGTTATGGGTTGTTTGGTGTCGAGCTATTTGTCTGTGGTGATGAGGTGATTTTCAGTGAGGTCCCCCCTCGTCCACATGATACCGGGATGGTGACGTTAATTtct                                                                                                          <  2:95270/150‑1 (MQ=255)
                                         aaCGTGCGCAGGAGATTGCCCGTAAAGTGGTGCTGGCACTGGGCGGTTATGGGTTGTTTGGTGTCGAGCTATTTGTCTGTGGTGATGAGGTGATTTTCAGTGAGGTCCCCCCTCGTCCACATGATACCGGGATGGTGACGTTAATTTCTCa                                                                                                        <  2:220398/151‑1 (MQ=255)
                                            gttcgcAGGAGATTGCCCGTAAAGTGGTGTTGGCACTGTGCGGTTATGGGTTGTTTGGTGTCGAGCTATTTGTCTGTGGTGATGAGGTGATTTTCAGTGAGGTCCCCCCTCGTCCACATGATACCGGGATGGTGACGTTAATTTCTCAAg                                                                                                      <  1:259451/147‑1 (MQ=255)
                                            gTGCGCAGGAGATTGCCCGTAAAGTGGTGCTGGCACTGGGCGGTTATGGGTTGTTTGGTGTCGAGCTATTTGTCTGTGGTGATGAGGTGATTTTCAGTGAGGTCCCCCCTCGTCCACATGATACCGGGATGGTGACGTTAATTTCTCa                                                                                                        >  2:551457/1‑148 (MQ=255)
                                            gTGCGCAGGAGATTGCCCGTAAAGTGGTGCTGGCACTGGGCGGTTATGGGTTGTTTGGTGTCGAGCTATTTGTCTGTGGTGATGAGGTGATTTTCAGTGAGGTCCCCCCTCGTCCACATGATACCGGGATGGTGACGTTAATTTCTCa                                                                                                        >  2:686599/1‑148 (MQ=255)
                                            gTGCGCAGGAGATTGCCCGTAAAGTGGTGCTGGCACTGGGCGGTTATGGGTTGTTTGGTGTCGAGCTATTTGTCTGTGGTGATGAGGTGATTTTCAGTGAGGTCCCCCCTCGTCCACATGATACCGGGATGGTGACGTTAATTTCTCa                                                                                                        >  2:1177978/1‑148 (MQ=255)
                                                gcAGGAGATTGCCCGTAAAGTGGTGCTGGCACTGGGCGGTTATGGGTTGTTTGGTGTCGAGCTATTTGTCTGTGGTGATGAGGTGATTTTCAGTGAGGTCCCCCCTCGTCCACATGATACCGGGATGGTGACGTTAATTTCTCAAGAtctc                                                                                                 >  1:1210876/1‑151 (MQ=255)
                                                gcAGGAGATTGCCCGTAAAGTAGTGCAGGCACTGAGCGGTTATGGGTTGTTTTGTGCCGTGCTATTTGTCTGGGGTGATGAGGGGCTTATCAGTGAGGTCCCCCCTCGTCCACATGATACCGGGATGGAGACGTTAATTTCTCATGAtctc                                                                                                 >  1:1044624/1‑151 (MQ=255)
                                                  aGGAGATTGCCCGTAAAGTGGTGCTGGCACTGGGCGGTTATGGGTTGTTTGGTGTCGAGCTATTTGTCTGTGGTGATGAGGTGATTTTCAGTGAGGTCCCCCCTCGTCCACATGATACCGGGATGGTGACGTTAATTTCTCAAGAtctctc                                                                                               >  1:747469/1‑151 (MQ=255)
                                                    gagaTTGCCCGTAAAGTGGTGCTGGCACTGGGCGGTTATGGGTTGTTTGGTGTCGAGCTATTTGTCTGTGGTGATGAGGTGATTTTCAGTGAGGTCCCCCCTCGTCCACATGATACCGGGATGGTGACGTTAATTTCTCAAGATCTCTCag                                                                                             <  2:365467/151‑1 (MQ=255)
                                                                  aGTGGTGCTGGCACTGGGCGGTTATGGGTTGTTTGGTGTCGAGCTATTTGTCTGTGGTGATGAGGTGATTTTCAGTGAGGTCCCCCCTCGTCCACATGATACCGGGATGGTGACGTTAATTTCTCAAGATCTCTCAGAGTTTGCCCTGCAt                                                                               >  2:579371/1‑151 (MQ=255)
                                                                     ggtgCTGGCACTGGGCGGTTATGGGTTGTTTGGTGTCGAGCTATTTGTCTGTGGTGATGAGGTGATTTTCAGTGAGGTCCCCCCTCGTCCACATGATACCGGGATGGTGACGTTAATTTCTCAAGATCTCTCAGAGTTTGCCCTGCATGt                                                                             <  2:887645/150‑1 (MQ=255)
                                                                     ggtgCTGGCACTGGGCGGTTATGGGTTGTTTGGTGTCGAGCTATTTGTCTGTGGTGATGAGGTGATTTTCAGTGAGGTCCCCCCTCGTCCACATGATACCGGGATGGTGACGTTAATTTCTCAAGATCTCTCAGAGTTTGCCCTGCATGTa                                                                            <  2:168328/151‑1 (MQ=255)
                                                                     ggtgCTGGCACTGGGCGGTTATGGGTTGTTTGGTGTCGAGCTATTTGTCTGTGGTGATGAGGTGATTTTCAGTGAGGTCCCCCCTCGTCCACATGATACCGGGATGGTGACGTTAATTTCTCAAGATCTCTCAGAGTTTGCCCTGCATGTa                                                                            <  2:1042944/151‑1 (MQ=255)
                                                                       tgCTGGCACTGGGCGGTTATGGGTTGTTTGGTGTCGAGCTATTTGTCTGTGGTGATGAGGTGATTTTCAGTGAGGTCCCCCCTCGTCCACATGATACCGGGATGGGGAAGTTAATTTCTCAAGATCTCTCAGAGTTTGCCCTGCATGTACg                                                                          >  1:1236292/1‑151 (MQ=255)
                                                                        gCTGGCACTGGGCGGTTATGGGTTGTTTGGTGTCGAGCTATTTGTCTGTGGTGATGAGGTGATTTTCAGTGAGGTCCCCCCTCGTCCACATGATACCGGGATGGTGACGTTGATTTCTCAAGATCTCTCAGAGTTTGCCCTGCATGTACGt                                                                         <  2:6160/151‑1 (MQ=255)
                                                                        gCTGGCACTGGGCGGTTATGGGTTGTTTGGTGTCGAGCTATTTGTCTGTGGTGATGAGGTGATTTTCAGTGAGGTCCCCCCTCGTCCACATGATACCGGGATGGTGACGTTAATTTCTCAAGATCTCTCAGAGTTTGCCCTGCATGTACGt                                                                         <  2:511936/151‑1 (MQ=255)
                                                                          tGGCACTGGGCGGTTATGGGTTGTTTGGTGTCGAGCTATTTGTCTGTGGTGATGAGGTGATTTTCAGTGAGGTCCCCCCTCGTCCACATGATACCGGGATGGTGACGTTAATTTCTCAAGATCTCTCAGAGTTTGCCCTGCATGTACGTGc                                                                       >  2:692689/1‑151 (MQ=255)
                                                                           ggCACTGGGCGGTTATGGGTTGTTTGGTGTCGAGCTATTTGTCTGTGGTGATGAGGTGATTTTCAGTGAGGTCCCCCCTCGTCCACATGATACCGGGATGGTGACGTTAATTTCTCAAGATCTCTCAGAGTTTGCCCTGCATGTACGTGc                                                                       <  1:328775/150‑1 (MQ=255)
                                                                                 gggCGGTTATGGGTTGTTTGGTGTCGAGCTATTTGTCTGTGGTGATGAGGTGATTTTCAGTGAGGTCCCCCCTCGTCCACATGATACCGGGATGGTGACGTTAATTTCTCAAGATCTCTCAGAGTTTGCCCTGCATGTACGTGCCTTcc                                                                  <  2:440104/149‑1 (MQ=255)
                                                                                 gggCGGTTATGGGTTGTTTGGTGTCGAGCTATTTGTCTGTGGTGATGAGGTGATTTTCAGTGAGGTCCCCCCTCGTCCACATGATACCGGGATGGTGACGTTAATTTCTCAAGATCTCTCAGAGTTTGCCCTGCATGTACGTGCCTTCCTc                                                                >  1:7937/1‑151 (MQ=255)
                                                                                 gggCGGTTATGGGTTGTTTGGTGTCGAGCTATTTGTCTGTGGTGATGAGGTGATTTTCAGTGAGGTCCCCCCTCGTCCACATGATACCGGGATGGTGACGTTAATTTCTCAAGATCTCTCAGAGTTTGCCCTGCATGTACGTGCCTTCCTc                                                                <  2:906997/151‑1 (MQ=255)
                                                                                    cGGTTATGGGTTGTTTGGTGTCGAGCTATTTGTCTGTGGTGATGAGGTGATTTTCAGTGAGGTCCCCCTCGTCCACATGATACCGGGATGGTGACGTTAATTTCTCAAGATCTCTCAGAGTTTGCCCTGCATGTACGTGCCTTCCTCGGAc                                                            <  1:677121/151‑1 (MQ=255)
                                                                                     ggTTATGGGTTGTTTGGTGTCGAGCTATTTGTCTGTGGTGATGAGGTGATTTTCAGTGAGGTCTCCCCTCGTCCACATGATACCGGGATGGTGACGTTAATTTCTCAAGATCTCTCAGAGTTTGCCCTGCATGTACGTGCCTTCCTCGGAc                                                            >  2:578000/1‑151 (MQ=255)
                                                                                        tATGGGTTGTTTGGTGTCGAGCTATTTGTCTGTGGTGATGAGGTGATTTTCAGTGAGGTCCCCCCTCGTCCACATGATACCGGGATGGTGACGTTAATTTCTCAAGATCTCTCAGAGTTTGCCCTGCATGTACGTGCCTTCCTCGGACtt                                                          <  1:708741/150‑1 (MQ=255)
                                                                                           gggTTGTTTGGTGTCGAGCTATTTGTCTGTGGTGATGAGGTGATTTTCAGTGAGGTCCCCCCTCGTCCACATGATACCGGGATGGTGACGTTAATTTCTCAAGATCTCTCAGAGTTTGCCCTGCATGTACGTGCCTTCCTCGGACTTCCgg                                                      >  2:93636/1‑151 (MQ=255)
                                                                                            ggTTGTTTGGTGTCGAGCTATTTGTCTGTGGTGATGAGGTGATTTTCAGTGAGGTCCCCCCTCGTCCACATGATACCGGGATGGTGACGTTAATTTCTCAAGATCTCTCAGAGTTTGCCCTGCATGTACGTGCCTTCCTCGGACTTCCGGt                                                     <  2:337032/151‑1 (MQ=255)
                                                                                             gttgttTGGTGTCGAGCTATTTGTCTGTGGTGATGAGGTGATTTTCAGTGAGGTCCCCCCTCGTCCACATGATACCGGGATGGTGACGTTAATTTCTCAAGATCTCTCAGAGTTTGCCCTGCATGTACGTGCCTTCCTCGGACTTCCGGtt                                                    <  1:1147522/151‑1 (MQ=255)
                                                                                                  ttGGTGTCGAGCTATTTGTCTGTGGTGATGAGGTGATTTTCAGTGAGGTCCCCCCTCGTCCACATGATACCGGGATGGTGACGTTAATTTCTCAAGATCTCTCAGAGTTTGCCCTGCATGTACGTGCCTTCCTCGGACTTCCGGTTGGCg                                                <  1:1250534/150‑1 (MQ=255)
                                                                                                    ggTGTCGAGCTATTTGTCTGTGGTGATGAGGTGATTTTCAGTGAGGTCCCCCCTCGTCCACATGATACCGGGATGGTGACGTTAATTTCTCAAGATCTCTCAGAGTTTGCCCTGCATGTACGTGCCTTCCTCGGACTTCCGGTTGGCggg                                              <  2:401069/150‑1 (MQ=255)
                                                                                                       gtCGAGCTATTTGTCTGTGGTGATGAGGTGATTTTCAGTGAGGTCCCCCCTCGTCCACATGATACCGGGATGGTGACGTTAATTTCTCAAGATCTCTCAGAGTTTGCCCTGCATGTACGTGCCTTCCTCGGACTTCCGGTTGGCGGGATcc                                          <  1:108753/151‑1 (MQ=255)
                                                                                                       gtCGAGCTATTTGTCTGTGGTGATGAGGTGATTTTCAGTGAGGTCCCCCCTCGTCCACATGATACCGGGATGGTGACGTTAATTTCTCAAGATCTCTCAGAGTTTGCCCTGCATGTACGTGCCTTCCTCGGACTTCCGGTTGGCGGGATc                                           <  1:841554/150‑1 (MQ=255)
                                                                                                        tCGAGCTATTTGTCTGTGGTGATGAGGTGATTTTCAGTGAGGTCCCCCCTCGTCCACATGATACCGGGATGGTGACGTTAATTTCTCAAGATCTCTCAGAGTTTGCCCTGCATGTACGTGCCTTCCTCGGACTTCCGGTTGGCGGGATCCg                                         >  2:1166843/1‑151 (MQ=255)
                                                                                                           aGCTATTTGTCTGTGGTGATGAGGTGATTTTCAGTGAGGTCCCCCCTCGTCCACATGATACCGGGATGGTGACGTTAATTTCTCAAGATCTCTCAGAGTTTGCCCTGCATGTACGTGCCTTCCTCGGACTTCCGGTTGGCGGGATCCGTCa                                      <  2:231133/151‑1 (MQ=255)
                                                                                                           aGCTATTTGTCTGTGGTGATGAGGTGATTTTCAGTGAGGTCCCCCCTCGTCCACATGATACCGGGATGGTGACGTTAATTTCTCAAGATCTCTCAGAGTTTGCCCTGCATGTACGTGCCTTCCTCGGACTTCCGGTTGGCGGGATCCGTCa                                      <  1:44154/151‑1 (MQ=255)
                                                                                                           aGCTATTTGTCTGTGGTGATGAGGTGATTTTCAGTGAGGTCCCCCCTCGTCCACATGATACCGGGATGGTGACGTTAATTTCTCAAGATCTCTCAGAGTGTGCCCTGCATGTACGTGCCTTCCTCGGACTTCCGGTTGGCGGGATCCGTCa                                      <  1:82086/151‑1 (MQ=255)
                                                                                                              tATTTGTCTGTGGTGATGAGGTGATTTTCAGTGAGGTCCCCCCTCGTCCACATGATACCGGGATGGTGACGTTAATTTCTCAAGATCTCTCAGAGTTTGCCCTGCATGTACGTGCCTTCCTCGGACTTCCGGTTGGCGGGATCCGTCAGt                                    >  2:1022420/1‑150 (MQ=255)
                                                                                                                 ttGTCTGTGGTGATGAGGTGATTTTCAGTGAGGTCCCCCCTCGTCCACATGATACCGGGATGGTGACGTTAATTTCTCAAGATCTCTCAGAGTTTGCCCTGCATGTACGTGCCTTCCTCGGACTTCCGGTTGGCGGGATCCGTCAGTATg                                 <  1:835835/150‑1 (MQ=255)
                                                                                                                  tGTCTGTGGTGATGAGGTGATTTTCAGTGAGGTCCCCCCTCGTCCACATGATACCGGGATGGTGACGTTAATTTCTCAAGATCTCTCAGAGTTTGCCCTGCATGTACGTGCCTTCCTCGGACTTCCGGTTGGCGGGATCCGTCAGTATgg                                <  1:1000709/150‑1 (MQ=255)
                                                                                                                    tCTGTGGTGATGAGGTGATTTTCAGTGAGGTCCCCCCTCGTCCACATGATACCGGGATGGTGACGTTAATTTCTCAAGATCTCTCAGAGTTTGCCCTGCATGTACGTGCCTTCCTCGGACTTCCGGTTGGCGGGATCCGTCAGTATGGTcc                             >  1:1133135/1‑151 (MQ=255)
                                                                                                                         ggtgATGAGGTGATTTTCAGTGAGGTCCCCCCTCGTCCACATGATACCGGGATGGTGACGTTAATTTCTCAAGATCTCTCAGAGTTTGCCCTGCATGTACGTGCCTTCCTCGGACTTCCGGTTGGCGGGATCCGTCAGTATGGTCCTGCa                         <  2:73142/150‑1 (MQ=255)
                                                                                                                                aGGTGATTTTCAGTGAGGTCCCCCCTCGTCCACATGATACCGGGATGGTGACGTTAATTTCTCAAGATCTCTCAGAGTTTGCCCTGCATGTACGTGCCTTCCTCGGACTTCCGGTTGGCGGGATCCGTCAGTATGGTCCTGCAGCTTCTGc                 >  1:962973/1‑151 (MQ=255)
                                                                                                                                  gTGATTTTCAGTGAGGTCCCCCCTCGTCCACATGATACCGGGATGGTGACGTTAATTTCTCAAGATCTCTCAGAGTTTGCCCTGCATGTACGTGCCTTCCTCGGACTTCCGGTTGGCGGGATCCGTCAGTATGGTCCTGCAGCTTCTGCCg               <  1:1177944/151‑1 (MQ=255)
                                                                                                                                     aTTTTCAGTGAGGTCCCCCCTCGTCCACATGATACCGGGATGGTGACGTTAATTTCTCAAGATCTCTCAGAGTTTGCCCTGCATGTACGTGCCTTCCTCGGACTTCCGGTTGGCGGGATCCGTCAGTATGGTCCTGCAGCTTCTGCCGtt             >  1:137779/1‑150 (MQ=255)
                                                                                                                                         tCAGTGAGGTCCCCCCTCGTCCACATGATACCGGGATGGTGACGTTAATTTCTCAAGATCTCTCAGAGTTTGCCCTGCATGTACGTGCCTTCCTCGGACTTCCGGTTGGCGGGATCCGTCAGTATGGTCCTGCAGCTTCTGCCGTTATTCt        >  1:634396/1‑151 (MQ=255)
                                                                                                                                          cAGTGAGGTCTCCCCTCGTCCACATGATACCGGGATGGTGACGTTAATTTCTCAAGATCTCTCAGAGTTTGCCCTGCATGTACGTGCCTTCCTCGGACTTCCGGTTGGCGGGATCCGTCAGTATGGTCCTGCAGCTTCTGCCGTTATTCTg       <  1:577968/151‑1 (MQ=255)
                                                                                                                                          cAGTGAGGTCCCCCCTCGTCCACATGATACCGGGATGGTGACGTTAATTTCTCAAGATCTCTCAGAGTTTGCCCTGCATGTACGTGCCTTCCTCGGACTTCCGGTTGGCGGGATCCGTCAGTATGGTCCTGCAGCTTCTGCCGTTATTCt        <  2:477901/150‑1 (MQ=255)
                                                                                                                                          cAGTGAGGTCCCCCCTCGTCCACATGATACCGGGATGGTGACGTTAATTTCTCAAGATCTCTCAGAGTTTGCCCTGCATGTACGTGCCTTCCTCGGACTTCCGGTTGGCGGGATCCGTCAGTATGGTCCTGCAGCTTCTGCCGTTATTCTg       <  2:957177/151‑1 (MQ=255)
                                                                                                                                               aGGTCCCCCCTCGTCCACATGATACCGGGATGGTGACGTTAATTTCTCAAGATCTCTCAGAGTTTGCCCTGCATGTACGTGCCTTCCTCGGACTTCCGGTTGGCGGGATCCGTCAGTATGGTCCTGCAGCTTCTGCCGTTATTCTGCcaca  <  1:788570/151‑1 (MQ=255)
                                                                                                                                                ggTCCCCCCTCGTCCACATGATACCGGGATGGTGACGTTAATTTCTCAAGATCTCTCAGAGTTTGCCCTGCATGTACGTGCCTTCCTCGGACTTCCGGTTGGCGGGATCCGTCAGTATGGTCCTGCAGCTTCTGCCGTTATTCTGCcaca  >  2:1191581/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
TGAATCCTGGCAACCACAGCAAATGAGCCCGCTTGCCCTTGAACGTGCGCAGGAGATTGCCCGTAAAGTGGTGCTGGCACTGGGCGGTTATGGGTTGTTTGGTGTCGAGCTATTTGTCTGTGGTGATGAGGTGATTTTCAGTGAGGTCTCCCCTCGTCCACATGATACCGGGATGGTGACGTTAATTTCTCAAGATCTCTCAGAGTTTGCCCTGCATGTACGTGCCTTCCTCGGACTTCCGGTTGGCGGGATCCGTCAGTATGGTCCTGCAGCTTCTGCCGTTATTCTGCCACA  >  NC_000913/1931573‑1931866

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