Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,365,487 C→T 5.1% pseudogene (303/367 nt) yhcE → pseudogene fragment; interrupted by IS5

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,365,4870CT 5.1% 283.3 / 3.8 79pseudogene (303/367 nt)yhcEpseudogene fragment; interrupted by IS5
Reads supporting (aligned to +/- strand):  ref base C (40/35);  new base T (2/2);  total (42/37)
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

CTGGTAATGCAATGCGATCAGGACGCTCATATTCGTGTTCAGTTGAATACCGCTAATATTTCCAACAACGGATTTGGATTTGGCCCGAATGGTTCATTAAATCTGATCGCCAGCGATGCATTTTCAGGTAGTAACAATCTCGACCTGGCGTTAGCCAGTGGCAAAAATGACAATCCGGGCAGCACAGGAACAGCATCAATTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAG  >  NC_000913/3365353‑3365622
                                                                                                                                      |                                                                                                                                       
cTGGTAATGCAATGCGATCAGGACGCTCATATTCGTGTTCAGTTGAATACCGCTAATATTTCCAACAACGGATTTGGATTTGGCCCGAATGGTTCATTAAATCTGATCGCCAGCGATGCATTTTCAGGTAGTAACAAtct                                                                                                                                    >  2:956160/1‑140 (MQ=255)
cTGGTAATGCAATGCGATCAGGACGCTCATATTCGTGTTCAGTTGAATACCGCTAATATTTCCAACAACGGATTTGGATTTGGCCCGAATGGTTCATTAAATCTGATCGCCAGCGATGCATTTTCAGGTAGTAACAAtct                                                                                                                                    >  2:328614/1‑140 (MQ=255)
cTGGTAATGCAATGCGATCAGGACGCTCATATTCGTGTTCAGTTGAATACCGCTAATATTTCCAACAACGGATTTGGATTTGGCCCGAATGGTTCATTAAATCTGATCGCCAGCGATGCATTTTCAGGTAGTAACAAtct                                                                                                                                    >  2:328850/1‑140 (MQ=255)
           atgcGATCAGGACGCTCATATTCGTGTTCAGTTGAATACCGCTAATATTTCCAACAACGGATTTGGATTTGGCCCGAATGGTTCATTAAATCTGATCGCCAGCGATGCATTTTCAGGTAGTAACAATCTCGACCTGGCGt                                                                                                                         >  1:1195417/1‑140 (MQ=255)
             gcGATCAGGACGCTCATATTCTAGTTCAGTTGAATACCGCAAATATTACCAACAACGGATTAGGATTTGGCCCGAATGGTACAATAAATCTGATCGCCAGCGATGCATTATCAGGTAGTAACAATCTCGACCTGGCGTTa                                                                                                                       >  2:1168835/1‑140 (MQ=255)
             gcGATCAGGACGCTCATATTCGTGTTCAGTTGAATACCGCTAATATTTCCAACAACGGATTTGGATTTGGCCCAAATGGTTCATTAAATCTGATCGCCAGCGATGCATTTTCAGGTAGTAACAATCTCGACCTGGCGTTa                                                                                                                       >  1:39594/1‑140 (MQ=255)
             gcGATCAGGACGCTCATATTCGTGTTCAGTTGAATACCGCTAATATTTCCAACAACGGATTTGGATTTGGACCGAATGGTTCATTAAATCTGATCGCCAGCGATGCAGTTTCAGGTAGTAACAATCTCGACCTGGCGTTa                                                                                                                       >  2:533251/1‑140 (MQ=255)
              cGATCAGGACGCTCATATTCGTGTTCAGTTGAATACCGCTAATATTTCCAACAACGGATTTGGATTTGGCCCGAATGGTTCATTAAATCTGATCGCCAGCGATGCATTTTCAGGTAGTAACAATCTCGACCTGGCGTTAg                                                                                                                      >  2:7103/1‑140 (MQ=255)
              cGATCAGGACGCTCATATTCGTGTTCAGTTGAATACCGCTAATATTTCCAACAACGGATTTGGATTTGGCCCGAATGGTTCATTAAATCTGATCGCCAGCGATGCATTTTCAGGTAGTAACAATCTCGACCTGGCGTTAg                                                                                                                      >  2:283967/1‑140 (MQ=255)
              cGATCAGGACGCTCATATTCGTGTTCAGTTGAATACCGCTAATATTTCCAACAACGGATTTGGATTTGGCCCGAATGGTTCATTAAATCTGATCGCCAGCGATGCATTTTCAGGTAGTAACAATCTCGACCTGGCGTTAg                                                                                                                      >  2:283852/1‑140 (MQ=255)
                 tCAGGACGCTCATATTCGTGTTCAGTTGAATACCGCTAATATTTCCAACAACGGATTTGGATTTGGCCCGAATGGTTCATTAAATCTGATCGCCAGCGATGCATTTTCAGGTAGTAACAATCTCGACCTGGCGTTAGCCa                                                                                                                   >  1:1040359/1‑140 (MQ=255)
                       cGCTCATATTCGTGTTCAGTTGAATACCGCTAATATTTCCAACAACGGATTTGGATTTGGCCCGAATGGTTCATTAAATCTGATCGCCAGCGATGCATTTTCAGGTAGTAACAATCTCGACCTGGCGTTAGCCAGTGGCa                                                                                                             >  2:794859/1‑140 (MQ=255)
                       cGCTCATATTCGTGTTCAGTTGAATACCGCTAATATTTCCAACAACGGATTTGGATTTGGCCCGAATGGTTCATTAAATCTGATCGCCAGCGATGCATTTTCAGGTAGTAACAATCTCGACCTGGCGTTAGCCAGTGGCa                                                                                                             >  2:457373/1‑140 (MQ=255)
                       cGCTCATATTCGTGTTCAGTTGAATACCGATAATATTTCCAACAACGGATTTGGATTGGCCCGAATGGTTCATTAAATCTGATCGCCAGCGATGCATTTTCAGGTAGTAACAATCTCGACCTGGCGTTAGCCAGTGGCaa                                                                                                            <  2:55913/140‑1 (MQ=255)
                        gCTCATATTCGTGTTCAGTTGAATACCGCTAATATTTCCAACAACGGATTTGGATTTGGCCCGAATGGTTCATTAAATCTGATCGCCAGCGATGCATTTTCAGGTAGTAACAATCTCGACCTGGCGTTAGCCAGTGGCaa                                                                                                            <  2:310456/140‑1 (MQ=255)
                        gCTCATATTCGTGTTCAGTTGAATACCGCTAATATTTCCAACAACGGATTTGGATTTGGCCCGAATGGTTCATTAAATCTGATCGCCAGCGATGCATTTTCAGGTAGTAACAATCTCGACCTGGCGTTAGCCAGTGGCaa                                                                                                            >  1:247831/1‑140 (MQ=255)
                        gCTCATATTCGTGTTCAGTTGAATACCGCTAATATTTCCAACAACGGATTTGGATTTGGCCCGAATGGTTCATTAAATCTGATCGCCAGCGATGCATTTTCAGGTAGTAACAATCTCGACCTGGCGTTAGCCAGTGGCaa                                                                                                            <  2:310392/140‑1 (MQ=255)
                         cTCATATTCGTGTTCAGTTGAATACCGCTAATATTTCCAACAACGGATTTGGATTTGGCCCGAATGGTTCATTAAATCTGATCGCCAGCGATGCATTTTCAGGTAGTAACAATCTCGACCTGGCGTTAGCCAGTGGCaaa                                                                                                           >  1:991208/1‑140 (MQ=255)
                             tatTCGTGTTCAGTTGAATACCGCTAATATTTCCAACAACGGATTTGGATTTGGCCCGAATGGTTCATTAAATCTGATCGCCAGCGATGCATTTTCAGGTAGTAACAATCTCGACCTGGCGTTAGCCAGTGGCAAAAATg                                                                                                       <  2:1074600/140‑1 (MQ=255)
                             tatTCGTGTTCAGTTGAATACCGCTAATATTTCCAACAACGGATTTGGATTTGGCCCGAATGGTTCATTAAATCTGATCGCCAGCGATGCATTTTCAGGTAGTAACAATCTCGACCTGGCGTTAGCCAGTGGCAAAAATg                                                                                                       <  2:1072463/140‑1 (MQ=255)
                             tatTCGTGTTCAGTTGAATACCGCTAATATTTCCAACAACGGATTTGGATTTGGCCCGAATGGTTCATTAAATCTGATCGCCAGCGATGCATTTTCAGGTAGTAACAATCTCGACCTGGCGTTAGCCAGTGGCAAAAATg                                                                                                       <  2:1015682/140‑1 (MQ=255)
                               ttCGTGTTCAGTTGAATACCGCTAATATTTCCAACAACGGATTTGGATTTGGCCCGAATGGTTCATTAAATCTGATCGCCAGCGATGCATTTTCAGGTAGTAACAATCTCGACCTGGCGTTAGCCAGTGGCAAAAATGAc                                                                                                     >  2:7999/1‑140 (MQ=255)
                               ttCGTGTTCAGTTGAATACCGCTAATATTTCCAACAACGGATTTGGATTTGGCCCGAATGGTTCATTAAATCTGATCGCCAGCGATGCATTTTCAGGTAGTAACAATCTCGACCTGGCGTTAGCCAGTGGCAAAAATGAc                                                                                                     >  2:786528/1‑140 (MQ=255)
                               ttCGTGATCAGTTGAATACCGCTAATATTTCCAACAACGGATTTGGATTTGGCCCGAATGGTTCATTAAATCTGATCGCCAGCGATGCATTTTCAGGTAGTAACAATCTCGACCTGGCGTTAGCCAGTGGCAAAAATGAc                                                                                                     >  2:745962/1‑140 (MQ=255)
                                    gtTCAGTTGAATACCGCTAATATTTCCAACAACGTATTTGGATTTGGCCCGAATGGTTCATTAAATCTGATCGCCAGCGATGCATTTTCAGGTAGTCACAATCTCGACCTGGCGTTAGCCAGTGGCAAAAATGACAATcc                                                                                                >  2:1064434/1‑140 (MQ=255)
                                    gtTCAGTTGAATACCGCTAATATTTCCAACAACGGATTTGGATTTGGCCCGAATGGTTCATTAAATCTGATCGCCAGCGATGCATTTTCAGGTAGTAACAATCTCGACCTGGCGTTAGCCAGTGGCAAAAATGACAATcc                                                                                                >  2:581167/1‑140 (MQ=255)
                                    gtTCAGTTGAATACCGCTAATATTTCCAACAACGGATTTGGATTTGGCCCGAATGGTTCATTAAATCTGATCGCCAGCGATGCATTTTCAGGTAGTAACAATCTCGACCTGGCGTTAGCCAGTGGCAAAAATGACAATcc                                                                                                >  1:249441/1‑140 (MQ=255)
                                       cAGTTGAATACCGCTAATATTTCCAACAACGGATTTGGATTTGGCCCGAATGGTTCATTAAATCTGATCGCCAGCGATGCATTTTCAGGTAGTAACAATCTCGACCTGGCGTTAGCCAGTGGCAAAAATGACAATCCggg                                                                                             >  2:776470/1‑140 (MQ=255)
                                       cAGTTGAATACCGCTAATATTTCCAACAACGGATTTGGATTTGGCCCGAATGGTTCATTAAATCTGATCGCCAGCGATGCATTTTCAGGTAGTAACAATCACGACCTGGCTTTACCCAGTGGCAAAAATGACAATCCgag                                                                                             >  2:776380/1‑138 (MQ=255)
                                        aGTTGAATACCGCTAATATTTCCAACAACGGATTTGGATTTGGCCCGAATGGTTCATTAAATCTGATCGCCAGCGATGCATTTTCAGGTAGTAACAATCTCGACCTGGCGTTAGCCAGTGGCAAAACTGACAATCCGGgc                                                                                            <  1:11823/140‑1 (MQ=255)
                                        aGTTGAATACCGCTAATATTTCCAACAACGGATTTGGATTTGGCCCGAATGGTTCATTAAATCTGATCGCCAGCGATGCATTTTCAGGTAGTAACAATCTCGACCTGGCGTTAGCCAGTGGCAAAAATGACAATCCGGgc                                                                                            <  1:61803/140‑1 (MQ=255)
                                        aGTTGAATACCGCTAATATTTCCAACAACGGATTTGGATTTGGCCCGAATGGTTCATTAAATCTGATCGCCAGCGATGCATTTTCAGGTAGTAACAATCTCGACCTGGCGTTAGCCAGTGGCAAAAATGACAATCCGGgc                                                                                            <  1:61802/140‑1 (MQ=255)
                                        aGTTGAATACCGCTAATATTTCCAACAACGGATTTGGATTTGGCCCGAATGGTTCATTAAATCTGATCGCCAGCGATGCATTTTCAGGTAGTAACAATCTCGACCTGGCGTTAGCCAGTGGCAAAAATGACAATCCGGgc                                                                                            <  1:1176385/140‑1 (MQ=255)
                                        aGTTGAATACCGCTAATATTTCCAACAACGGATTTGGATTTGGCACGAATGGTTCATTAAATCTGATCGCCAGCGATGCATTTTCAGGTAGTAACAATCTCGACCTGGCGTTAGCCAGTGGCAAAAATGACAATCCGGgc                                                                                            <  1:1176396/140‑1 (MQ=255)
                                          ttGAATACCGCTAATATTTCCAACAACGGATTTGGATTTGGCCCGAATGGTTCATTAAATCTGATCGCCAGCGATGCATTTTCAGGTAGTAACAATCTCGACCTGGCGTTAGCCAGTGGCAAAAATGACAATCCGGgcag                                                                                          <  2:294685/140‑1 (MQ=255)
                                             aaTACCGCTAATATTTCCAACAACGGATTTGGATTTGGCCCGAATGGTTCATTAAATCTGATCGCCAGCGATGCATTTTCAGGTAGTAACAATCTCGACCTGGCGTTAGCCAGTGGCAAAAATGACAATCCGGGCAGcac                                                                                       <  2:1040359/140‑1 (MQ=255)
                                                  cGCTAATATTTCCAACAACGGATTTGGATTTGGCCCGAATGGTTCATTAAATCTGATCGCCAGCGATGCATTTTCAGGTAGTAACAATCTCGACCTGGCGTTAGCCAGTGGCAAAAATGACAATCCGGGCAGCACAGGaa                                                                                  >  2:1003330/1‑140 (MQ=255)
                                                  cGCTAATATTTCCAACAACGGATTTGGATTTGGCCCGAATGGTTCATTAAATCTGATCGCCAGCGATGCATTTTCAGGTAGTAACAATCTCGACCTGGCGTTAGCCAGTGGCAAAAATGACAATCCGGGCAGCACAGGaa                                                                                  >  2:1003331/1‑140 (MQ=255)
                                                    cTAATATTTCCAACAACGGATTTGGATTTGGCCCGTATGGTTCATTAAATCTGATCGCCAGCGATGCATTTTCAGGTAGTAACAATCTCGACCTGGCGTTAGCCAGTGGCAAAAATGACAATCCGGGCAGCACAGGAACa                                                                                >  1:380057/1‑140 (MQ=255)
                                                    cTAATATTTCCAACAACGGATTTGGATTTGGCCCGAATGGTTCATTAAATCTGATCGCCAGCGATGCATTTTCAGGTAGTAACAATCTCGACCTGGCGTTAGCCAGTGGCAAAAATGACAATCCGGGCAGCACAGGAACa                                                                                >  1:676761/1‑140 (MQ=255)
                                                      aaTATTTCCAACAACGGATTTGGATTTGGCCCGAATGGTTCATTAAATCTGATCGCCAGCGATGCATTTTCAGGTAGTAACAATCTCGACCTGGCGTTAGCCAGTGGCAAAAATGACAATCCGGGCAGCACAGGAACAGc                                                                              <  1:939540/140‑1 (MQ=255)
                                                      aaTATTTCCAACAACGGATTTGGATTTGGCCCGAATGGTTCATTAAATCTGATCGCCAGCGATGCATTTTCAGGTAGTAACAATCTCGACCTGGCGTTAGCCAGTGGCAAAAATGACAATCCGGGCAGCACAGGAACAGc                                                                              <  1:939429/140‑1 (MQ=255)
                                                      aaTATTTCCAACAACGGATTTGGATTTGGCCCGAATGGTTCATTAAATCTGATCGCCAGCGATGCATTTTCAGGTAGTAACAATCTCGACCTGGCGTTAGCCAGTGGCAAAAATGACAATCCGGGCAGCACAGGAACAGc                                                                              <  1:189162/140‑1 (MQ=255)
                                                       atatTTCCAACAACGGATTTGGATTTGGCCCGAATGGTTCATTAAATCTGATCGCCAGCGATGCATTTTCGGGTAGTAACAATCTCGACCTGGCGTTAGCCAGTGGCAAAAATGACAATCCGGGCAGCACAGGAACAGCa                                                                             <  2:224486/140‑1 (MQ=255)
                                                       atatTTCCAACAACGGATTTGGATTTGGCCCGAATGGTTCATTAAATCTGATCGCCAGCGATGCATTTTCAGGTAGTAACAATCTCGACCTGGCGTTAGCCAGTGGCAAAAATGACAATCCGGGCAGCACAGGAACAGCa                                                                             <  2:990072/140‑1 (MQ=255)
                                                       atatTTCCAACAACGGATTTGGATTTGGCCCGAATGGTTCATTAAATCTGATCGCCAGCGATGCATTTTCAGGTAGTAACAATCTCGACCTGGCGTTAGCCAGTGGCAAAAATGACAATCCGGGCAGCACAGGAACAGCa                                                                             <  2:224395/140‑1 (MQ=255)
                                                        tatTTCCAACAACGGATTTGGATTTGGCCCGAATGGTTCATTAAATCTGATCGCCAGCGATGCATTTTCAGGTAGTAACAATCTCGACCTGGCGTTAGCCAGTGGCAAAAATGACAATCCGGGCAGCACAGGAACAGCAt                                                                            <  1:33164/140‑1 (MQ=255)
                                                            tCCAACAACGGATTTGGATTTGGCCCGAATGGTTCATTAAATCTGATCGCCAGCGATGCATTTTCAGGTAGTAATAATCTCGAcc                                                                                                                               <  1:1099719/85‑1 (MQ=255)
                                                            tCCAACAACGGATTTGGATTTGGCCCGAATGGTTCATTAAATCTGATCGCCAGCGATGCATTTTCAGGTAGTAATAATCTCGAcc                                                                                                                               <  1:1100013/85‑1 (MQ=255)
                                                            tCCAACAACGGATTTGGATTTGGCCCGAATGGTTCATTAAATCTGATCGCCAGCGATGCATTTTCAGGTAGTAATAATCTCGAcc                                                                                                                               >  2:1100013/1‑85 (MQ=255)
                                                            tCCAACAACGGATTTGGATTTGGCCCGAATGGTTCATTAAATCTGATCGCCAGCGATGCATTTTCAGGTAGTAATAATCTCGAcc                                                                                                                               >  2:1099719/1‑85 (MQ=255)
                                                            tCCAACAACGGATTTGGATTTGGCCCGAATGGTTCATTAAATCTGATCGCCAGCGATGCATTTTCAGGTAGTAACAATCTCGACCTGGCGTTAGCCAGTGGCAAAAATGACAATCCGGGCAGCACAGGAACAGCATCAAt                                                                        <  2:25527/140‑1 (MQ=255)
                                                            tCCAACAACGGATTTGGATTTGGCCCGAATGGTTCATTAAATCTGATCGCCAGCGATGCATTTTCAGGTAGTAACAATCTCGACCTGGCGTTAGCCAGTGGCAAAAATGACAATCCGGGCAGCACAGGAACAGCATCAAt                                                                        >  2:74695/1‑140 (MQ=255)
                                                            tCCAACAACGGATTTGGATTTGGCCCGAATGGTTCATTAAATCTGATCGCCAGCGATGCATTTTCAGGTAGTAACAATCTCGACCTGGCGTTAGCCAGTGGCAAAAATGACAATCCGGGCAGCACAGGAACAGCATCAAt                                                                        >  1:559127/1‑140 (MQ=255)
                                                                  aacGGATTTGGATTTGGCCCGAATGGTTCATTAAATCTGATCGCCAGCGATGCATTTTCAGGTAGTAACAATCTCGACCTGGCGTTAGCCAGTGGCAAAAATGACAATCCGGGCAGCACAGGAACAGCATCAATTAGGGa                                                                  <  1:590066/140‑1 (MQ=255)
                                                                   acGGATTTGGATTTGGCCCGAATGGTTCATTAAATCTGATCGCCAGCGATGCATTTTCAGGTAGTAACAATCTCGACCTGGCGTTAGCCAGTGGCAAAAATGACAATCCGGGCAGCACAGGAACAGCATCAATTAGGGat                                                                 >  1:848203/1‑139 (MQ=255)
                                                                   acGGATTTGGATTTGGCCCGAATGGTTCATTAAATCTGATCGCCAGCGATGCATTTTCAGGTAGTAACAATCTCGACCTGGCGTTAGCCAGTGGCAAAAATGACAATCCGGGCAGCACAGGAACAGCATCAATTAGGGaa                                                                 >  1:848998/1‑140 (MQ=255)
                                                                    cGGATTTGGATTTGGCCCGAATGGTTCATTAAATCTGATCGCCAGCGATGCATTTTCAGGTAGTAACAATCTCGACCTGGCGTTAGCCAGTGGCAAAAATGACAATCCGGGCAGCACAGGAACAGCATCAATTAGGGAAg                                                                <  1:1006852/140‑1 (MQ=255)
                                                                    cGGATTTGGATTTGGCCCGAATGGTTCATTAAATCTGATCGCCAGCGATGCATTTTCAGGTAGTAACAATCTCGACCTGGCGTTAGCCAGTGGCAAAAATGACAATCCGGGCAGCACAGGAACAGCATCAATTAGGGAAg                                                                <  1:1006885/140‑1 (MQ=255)
                                                                    cGGATTTGGATTTGGCCCGAATGGTTCATTAAATCTGATCGCCAGCGATGCATTTTCAGGTAGTAACAATCTCGACCTGGCGTTAGCCAGTGGCAAAAATGACAATCCGGGCAGCACAGGAACAGCATCAATTAGGGAAg                                                                <  1:1182786/140‑1 (MQ=255)
                                                                               ttggCCCGACTGGTTCATTATATCTGATCTCCAGCGATGCTTTTTCAGGTAGTAACAATCTCGACCTGGCGTTAGCCAGTGGCAAAAATGACAATCCGGGCAGCACAGGAACAGCATCAATTAGGGAAGGTGCGAATAAg                                                     <  2:320214/140‑1 (MQ=255)
                                                                               ttggCCCGAATGGTTCATTAAATCTGATCGCCAGCGTTGCATTTTCAGGTAGTAACAATCTCGACCTGGCGTTAGCCAGTGGCAAAAATGACATTCCGGGCAGCACAGGAACAGCATCAATTAGGGAAGGTGCGAATAAg                                                     <  2:126588/140‑1 (MQ=255)
                                                                               ttggCCCGAATGGTTCATTAAATCTGATCGCCAGCGATGCATTTTCAGGTAGTAACAATCTCGACCTGGCGTTAGCCAGTGGCAAAAATGACAATCCGGGCAGCACAGGAACAGCATCAATTAGGGAAGGTGCGAATAAg                                                     <  2:126612/140‑1 (MQ=255)
                                                                               ttggCCCGAATGGTTCATTAAATCTGATCGCCAGCGATGCATTTTCAGGTAGTAACAATCTCGACCTGGCGTTAGCCAGTGGCAAAAATGACAATCCGGGCAGCACAGGAACAGCATCAATTAGGGAAGGTGCGAATAAg                                                     <  2:126583/140‑1 (MQ=255)
                                                                               ttggCCCGAATGGTTCATTAAATCTGATCGCCAGCGATGCATTTTCAGGTAGTAACAATCTCGACCTGGCGTTAGCCAGTGGCAAAAATGACAATCCGGGCAGCACAGGAACAGCATCAATTAGGGAAGGTGCGAATAAg                                                     <  2:122997/140‑1 (MQ=255)
                                                                                  gCCCGAATGGTTCATTAAATCTGATCGCCAGCGATGCATTTTCAGGTAGTAACAATCTCGACCTGGCGTTAGCCAGTGGCAAAAATGACAATCCGGGCAGCACAGGAACAGCATCAATTAGGGAAGGTGCGAATAAGCgg                                                  <  1:922531/140‑1 (MQ=255)
                                                                                            ttCATTAAATCTGATCGCCAGCGATGCATTTTCAGGTAGTAACAATCTCGACCTGGCGTTAGCCAGTGGCAAAAATGACAATCCGGGCAGCACAGGAACAGCATCAATTAGGGAAGGTGCGAATAAGCGGGGAAAttctt                                        >  2:579231/1‑140 (MQ=255)
                                                                                                 ttaaTCTGATCGCCAGCGATGCCTTTTCAGGTAGTAACTATCTCGTCCTGGCGTTAGCCAGTGGCCATAATGACAATCCGGGCAGCACAGGAACAGCATCAATTAGGTAAGGTGCGATAATGAGGCGAAATAGTTCTCgg                                   >  2:1042159/3‑140 (MQ=255)
                                                                                                 tAAATCTGATCGCCAGCGATGCATTTTCAGGTAGTAACTATCTCGACCTGGCGTTAGCCAGTGTCAATAATGTCAATCCGGGCAGCACAGGAACAGCCTCTATTAGGGAAGGTGCGTATTAGCGGGGAGATTATTCTCgg                                   >  2:3245/1‑140 (MQ=255)
                                                                                                 tAAATCTGATCGCCAGCGATGCATTTTCAGGTAGTAACACTCTCGACCTGGCGTTAGCCAGTGGCAAAAATGACAATCCGGGCAGCACAGGAACAGCATCAATTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCgg                                   <  1:770735/140‑1 (MQ=255)
                                                                                                 tAAATCTGATCGCCAGCGATGCATTTTCAGGTAGTAACAATCTCGACCTGGCGTTAGCCAGTGGCAAAAATGACAATCCGGGCAGCACAGGAACAGCATCAATTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCgg                                   >  2:586009/1‑140 (MQ=255)
                                                                                                     tCTGATCGCCAGCGATGCATTTTCAGGTAGTAACAATCTCGACCTGGCGTTAGCCAGTGGCAAAAATGACAATCCGGGCAGCACAGGAACAGCATCAATTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGa                               >  2:765988/1‑140 (MQ=255)
                                                                                                     tCTGATCGCCAGCGATGCATTTTCAGGTAGTAACAATCTCGACCTGGCGTTAGCCAGTGGCAAAAATGACAATCCGGGCAGCACAGGAACAGCATCAATTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGa                               >  1:1176391/1‑140 (MQ=255)
                                                                                                      cTGATCGCCAGCGATGCATTTTCAGGTAGTAACAATCTCGACCTGGCGTTAGCCAGTGGCAAAAATGACAATCCGGGCAGCACAGGAACAGCATCAATTAGGGAAGGTGCGACTAAGCGGGGAAATTCTTCTCGGCTGAc                              >  2:753208/1‑140 (MQ=255)
                                                                                                      cTGATCGCCAGCGATGCATTTTCAGGTAGTAACAATCTCGACCTGGCGTTAGCCAGTGGCAAAAATGACAATCCGGGCAGCACAGGAACAGCATCAATTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGAc                              >  2:753163/1‑140 (MQ=255)
                                                                                                           cGCCAGCGATGCATTTTCAGGTAGTAACAATCTCGACCTGGCGTTAGCCAGTGGCAAAAATGACAATCCGGGCAGCACAGGAACAGCATCAATTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGt                         <  2:1064277/140‑1 (MQ=255)
                                                                                                                        ttttCAGGTAGTAACAATCTCGACCTGGCGTTAGCCAGTGGCAAAAATGACAATCCGGGCAGCACAGGAACAGCATCAATTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATttctt            <  1:1218866/140‑1 (MQ=255)
                                                                                                                        ttttCAGGTAGTAACAATCTCGACCTGGCGTTAGCCAGTGGCAAAAATGACAATCCGGGCAGCACAGGAACAGCATCAATTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATttctt            >  1:667244/1‑140 (MQ=255)
                                                                                                                                  gtaACAATCTCGACCTGGCGTTAGCCAGTGGCAAAAATGACAATCCGGGCAGCACAGGAACAGCATCAATTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAg  <  1:509375/140‑1 (MQ=255)
                                                                                                                                      |                                                                                                                                       
CTGGTAATGCAATGCGATCAGGACGCTCATATTCGTGTTCAGTTGAATACCGCTAATATTTCCAACAACGGATTTGGATTTGGCCCGAATGGTTCATTAAATCTGATCGCCAGCGATGCATTTTCAGGTAGTAACAATCTCGACCTGGCGTTAGCCAGTGGCAAAAATGACAATCCGGGCAGCACAGGAACAGCATCAATTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAG  >  NC_000913/3365353‑3365622

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