Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,596,873 Δ1 bp 11.0% intergenic (‑136/+32) ypfJ ← / ← purC uncharacterized protein YpfJ/phosphoribosylaminoimidazole‑succinocarboxamide synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,596,8730C.11.0% 251.7 / 25.4 82intergenic (‑136/+32)ypfJ/purCuncharacterized protein YpfJ/phosphoribosylaminoimidazole‑succinocarboxamide synthase
Reads supporting (aligned to +/- strand):  ref base C (32/40);  new base . (4/5);  total (37/45)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCCAACGCATAGATACCTCAACAATATATTCATTATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGGAAACGGTCTTTGTCCATTTTCTCCAGCGTTTCTTTGTCCCACAG  >  NC_000913/2596728‑2597021
                                                                                                                                                 |                                                                                                                                                     
gCCAACGCATAGATACCTCAACAATATATTCATTATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGGATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGc                                                                                                                                                   <  2:879209/150‑1 (MQ=255)
gCCAACGCATAGATACCTCAACAATATATTCATTATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGc                                                                                                                                                   <  1:318187/150‑1 (MQ=255)
 ccAACGCATAGATACCTCAACAATATATTCATTATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCt                                                                                                                                                  >  2:1469532/1‑150 (MQ=255)
    aCGCATAGATACCTCAACAATATATTCATTATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCa                                                                                                                                              <  2:8632/150‑1 (MQ=255)
     cGCATAGATACCTCAACAATATATTCATTATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCAc                                                                                                                                             >  1:1278051/1‑150 (MQ=255)
       cATAGATACCTCAACAATATATTCATTATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACgg                                                                                                                                           <  2:571738/150‑1 (MQ=255)
        aTAGATACCTCAACAATATATTCATTATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGc                                                                                                                                          <  2:1497688/150‑1 (MQ=255)
         tAGATACCTCAACAATATATTCATTATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCa                                                                                                                                         <  1:55911/150‑1 (MQ=255)
         tAGATACCTCAACAATATATTCATTATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCAGCTGGCACGGCa                                                                                                                                         <  1:617082/150‑1 (MQ=255)
          aGATACCTCAACAATATATTCATTATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCAGCTGGCACGGCaa                                                                                                                                        <  2:12374/150‑1 (MQ=255)
            aTACCTCAACAATATATTCATTATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCAGCTGGCACGGCAAga                                                                                                                                      <  2:1057035/150‑1 (MQ=255)
            aTACCTCAACAATATATTCATTATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGAGGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCAGCTGGCACGGCAAga                                                                                                                                      <  1:645523/150‑1 (MQ=255)
              aCCTCAACAATATATTCATTATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAgatg                                                                                                                                    <  2:1278051/150‑1 (MQ=255)
               ccTCAACAATATATTCATTATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAgatga                                                                                                                                   <  1:935764/150‑1 (MQ=255)
               ccTCAACAATATATTCATTATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAgatga                                                                                                                                   <  1:1040791/150‑1 (MQ=255)
                cTCAACAATATATTCATTATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAgatgat                                                                                                                                  >  2:124469/1‑150 (MQ=255)
                cTCAACAATATATTCATTATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAgatgat                                                                                                                                  >  1:29115/1‑150 (MQ=255)
                 tCAACTATATATTTATTATAGTTATGTTCGTAGGTGTTTCATAGGTTGTTTTCCACTGTAAATAGGGATGAAATGGCGGAGGAATGAGCGTTGAATTCATATATTGCGAGAGCTGTTGTCTTGCCGTGCCAGATTCACGGCAAGatggtg                                                                                                                                 <  2:783398/150‑1 (MQ=255)
                  caacaaTATATTCATTATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGAtgatga                                                                                                                                <  1:1164260/150‑1 (MQ=255)
                    acaaTATATTCATTATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAAc                                                                                                                              <  1:124469/150‑1 (MQ=255)
                          tatTCATTATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGaaaa                                                                                                                        >  1:1363607/1‑150 (MQ=255)
                           atTCATTATAGAAATGATCGTCGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTTCACGGCAAGATGATGAACAGaaaat                                                                                                                       <  2:849013/150‑2 (MQ=255)
                            ttCATTATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCATGGCAAGATGATGAACAGAAAAAt                                                                                                                      >  2:168880/1‑150 (MQ=255)
                             tCATTATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATc                                                                                                                     <  1:793315/150‑1 (MQ=255)
                              cATTATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCTCAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCGGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCa                                                                                                                    >  2:133848/1‑150 (MQ=255)
                              cATTATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAATCGGGATGAAATGTCGGAGGAATGAGCGTTCAATGAATATATTCCGAGAGAGGTTGTCTTTCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCa                                                                                                                    >  1:610901/1‑150 (MQ=255)
                              cATTATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGaa                                                                                                                               >  2:780592/1‑139 (MQ=255)
                              cATTATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCa                                                                                                                    >  2:108025/1‑150 (MQ=255)
                              cATTATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAAGGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGaa                                                                                                                               <  1:780592/139‑1 (MQ=255)
                                   taGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCa                                                                                                               >  1:482755/1‑150 (MQ=255)
                                   taGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCa                                                                                                               >  2:1082379/1‑150 (MQ=255)
                                    aGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAg                                                                                                              <  1:1469532/150‑1 (MQ=255)
                                     gAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGc                                                                                                             <  1:993088/150‑1 (MQ=255)
                                          gATCGTAGGTTGTTCAGAGGAAGATTACCACAGTAAAGCGGGATGAAATGGCGGAGTAATGAGCGTGCACTGCATATATTCCGAGAGCGGTGGTCTTGCCGGGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTac                                                                                                        <  2:520875/150‑1 (MQ=255)
                                              gTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACAccc                                                                                                    <  2:1372714/150‑1 (MQ=255)
                                                  gtggtTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAggc                                                                                                >  2:79056/1‑150 (MQ=255)
                                                  gtggtTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAggc                                                                                                <  1:1156606/150‑1 (MQ=255)
                                                     gtTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAggcggc                                                                                             >  1:1428761/1‑150 (MQ=255)
                                                     gtTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAggcggc                                                                                             <  1:108025/150‑1 (MQ=255)
                                                       tCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGcggcgg                                                                                           >  1:815315/1‑150 (MQ=255)
                                                        cAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCACCTGCCCGGCAAGATGATGAACAGAAAAATCAGTCCAGCCGTACACCCAGGCggcggg                                                                                          >  2:417515/1‑150 (MQ=255)
                                                           agGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCggcgggcgc                                                                                       >  1:612429/1‑149 (MQ=255)
                                                            ggAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTTCCGTTCCAGCTGGCCGGCAAGAATGTTAAAAGAAAAAAAAGTGCAGCTGTAAACCCAGGCcgcgggcgcc                                                                                      >  1:406306/1‑148 (MQ=255)
                                                            ggAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGAc                                                                                      >  2:127618/1‑150 (MQ=255)
                                                            ggAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGGCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGAc                                                                                      <  1:1069157/150‑1 (MQ=255)
                                                              aaGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACgg                                                                                    >  2:518277/1‑150 (MQ=255)
                                                                      cacaGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGaaa                                                                                                                         <  2:1202406/105‑1 (MQ=255)
                                                                      cacaGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGaaa                                                                                                                         >  1:1202406/1‑105 (MQ=255)
                                                                       acaGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGGCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAgg                                                                           <  1:748251/150‑1 (MQ=255)
                                                                        caGGAAAGCGGGATGAAATGGCGGAGGACTGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGGGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCTGGCGACGGCTTCATAGGc                                                                          <  1:1538961/150‑1 (MQ=255)
                                                                        caGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGc                                                                          >  1:538849/1‑150 (MQ=255)
                                                                         aGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCt                                                                         <  1:727491/150‑1 (MQ=255)
                                                                            aaaGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCg                                                                      <  2:1428761/150‑1 (MQ=255)
                                                                            aaaGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCg                                                                      <  1:1622997/150‑1 (MQ=255)
                                                                                  ggATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAgg                                                                >  2:326667/1‑150 (MQ=255)
                                                                                           ggCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAgg                                                       >  1:1468127/1‑150 (MQ=255)
                                                                                            gCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTTCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACTGCTTCATAGGCTTCGATGAGGCCAGGGAGGc                                                      <  2:482755/150‑1 (MQ=255)
                                                                                             cGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGcg                                                     >  2:298406/1‑149 (MQ=255)
                                                                                                   aaTGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCg                                               <  2:538849/150‑1 (MQ=255)
                                                                                                      gAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGGaa                                            <  2:29115/150‑1 (MQ=255)
                                                                                                       aGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGGaaa                                           >  2:721704/1‑150 (MQ=255)
                                                                                                          gTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGGAAACgg                                        <  1:728070/150‑1 (MQ=255)
                                                                                                                  cATATATTCCGAGAGCGGTTGTCTTGCCGTGACAGCTGCACGGAAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTAGATCAGGCCACCCAGGCTCTGGCGGAAACGCTCTTTGta                                >  2:546677/1‑149 (MQ=255)
                                                                                                                    tatatTCCGCGAGCGGTGGGCATGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAACGCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGGAAACGGTCTTTGTcga                              <  2:1475877/150‑3 (MQ=255)
                                                                                                                    tatatTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGGAAACGGTCTTTGTCCa                              <  2:1482151/150‑1 (MQ=255)
                                                                                                                    tatatTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGGAAACGGTCTTTGTCCa                              >  2:1260590/1‑150 (MQ=255)
                                                                                                                    tatatGCCGAGAGCGGTGGTAGTGCCGGGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTCCACCCCCGCGGCGGGAGACGTCTTCATAGGCGTAGAGCAGTCCACCGAGGCTCTGGCGGAAACGGTCGTTGTCCa                              <  2:1470055/150‑1 (MQ=255)
                                                                                                                        ttCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGGAAACGGTCTTTGTCCAtttt                          <  1:1483471/150‑1 (MQ=255)
                                                                                                                           cGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGGAAACGGTCTTTGTCCATTTtctc                       >  2:335274/1‑150 (MQ=255)
                                                                                                                            gagaGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCAACGAGGCTCTGGCGGAAACGGTCTTTGTCCATTTTCTcc                      >  2:155407/1‑150 (MQ=255)
                                                                                                                                 tcatcatCTTGCCGTGCAGCTGGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGGAAACGGTCTTTGTCCATTTTCTCCAGCGt                 >  1:205047/7‑150 (MQ=255)
                                                                                                                                    ttGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGGAAACGGTCTTTGTCCATTTTCTCCAGCGTTTc              <  2:1468127/150‑1 (MQ=255)
                                                                                                                                    tcatCTTGCCGTGCAGCTGGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGGAAACGGTCTTTGTCCATTTTCTCCAGCGTTTc              >  1:225209/4‑150 (MQ=255)
                                                                                                                                     catCTTGCCGTGCAGCTGGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGGAAACGGTCTTTGTCCATTTTCTCCAGCGTTTCt             >  2:1389271/3‑150 (MQ=255)
                                                                                                                                      gTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACCCCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGGAAACGGGCTTTGTCCATTTTCTCCAGCGTTTCtt            >  1:1598875/1‑150 (MQ=255)
                                                                                                                                        cTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGGAAACGGTCTTTGTCCATTTTCTCCAGCGTTTCTTTg          <  1:155407/150‑1 (MQ=255)
                                                                                                                                         ttGCCGTGCAGCTGGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGGAAACGGTCTTTGTCCATTTTCTCCAGCGTTTCTTTGt         >  2:386496/1‑150 (MQ=255)
                                                                                                                                          tGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGGAAACGGTCTTTGTCCATTTTCTCCAGCGTTTCTTTGTc        <  1:721704/150‑1 (MQ=255)
                                                                                                                                          tGCCGTGCAGCTGGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGGAAACGGTCTTTGTCCATTTTCTCCAGCGTTTCTTTGTc        <  1:1389271/150‑1 (MQ=255)
                                                                                                                                           gCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGGAAACGGTCTTTGTCCATTTTCTCCAGCGTTTCTTTGTcc       <  1:79056/150‑1 (MQ=255)
                                                                                                                                               tGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGGAAACGGTCTTTGTCCATTTTCTCCAGCGTTTCTTTGTCCcaca   >  1:135132/1‑150 (MQ=255)
                                                                                                                                                gCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGGAAACGGTCTTTGTCCATTTTCTCCAGCGTTTCTTTGTCCCACAg  >  1:1559170/1‑150 (MQ=255)
                                                                                                                                                 |                                                                                                                                                     
GCCAACGCATAGATACCTCAACAATATATTCATTATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGGAAACGGTCTTTGTCCATTTTCTCCAGCGTTTCTTTGTCCCACAG  >  NC_000913/2596728‑2597021

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