Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,287,682 C→A 5.3% pseudogene (1233/2525 nt) yejO ← pseudogene, autotransporter outer membrane homology;putative transport; Not classified; putative ATP‑binding component of a transport system

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,287,6820CA 5.3% 267.2 / 3.2 75pseudogene (1233/2525 nt)yejOpseudogene, autotransporter outer membrane homology;putative transport; Not classified; putative ATP‑binding component of a transport system
Reads supporting (aligned to +/- strand):  ref base C (32/39);  new base A (2/2);  total (34/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.73e-01

GATGATTAGCCATATCGGTACGCATATAAAAATTGCCGTTTCCCGACAGGGTATTGGTGGTCAACGTAATGAAATTTTCCGCACTTGCCGTCACTGATGAACGAGTCACTGGCTCATCCATCATAATCACATTTGCGTCAGACATCGCAAGGTTGGCTAAGTTGAGGAGGGTGTTTTGGTTCGTCGTAGTGATATCGGCAATGATGGTCCATACGCTATTTTCGAGCGAAACGTCCGCTTTGCTGGTATCCACGGCACCATGCGTTCTAAAAGAAGCGCCCTCAGAAACCAC  >  NC_000913/2287535‑2287826
                                                                                                                                                   |                                                                                                                                                
gatgatTAGCCATATCGGTACGCATATAAAAATTGCCGTTTCCCGACAGGGTATTGGTGGTCAACGTAATGAAATTTTCCGCACTTGCCGTCACTGATGAACGAGTCACTGGCTCATCCATCATAATCACATTTGCGTCAGACATCGCa                                                                                                                                                 >  3:25321/1‑149 (MQ=255)
   gatTAGCCATATCGGTACGCATATAAAAATTGCCGTTTCCCGACAGGGTATTGGTGGTCAACGTAATGAAATTTTCCGCACTTGCCGTCACTGATGAACGAGTCACTGGCTCATCCATCATAATCACATTTGCGTCAGACATCGCAAgg                                                                                                                                              >  4:523763/1‑149 (MQ=255)
    atTAGCCATATCGGTACGCATATAAAAATTGCCGTTTCCCGACAGGGTATTGGTGGTCAACGTAATGAAATTTTCCGCACTTGCCGTCACTGATGAACGAGTCACTGGCTCATCCATCATAATCACATTTGCGTCAGACATCGCAAGGt                                                                                                                                             <  2:320409/149‑1 (MQ=255)
     ttAGCCATATCGGTACGCATATAAAAATTGCCGTTTCCCGACAGGGTATTGGTGGTCAACGTAATGAAATTTTCCGCACTTGCCGTCACTGATGAACGAGTCACTGGCTCATCCATCATAATCACATTTGCGTCAGACATCGCAAGGtt                                                                                                                                            <  1:226947/149‑1 (MQ=255)
           atatCGGTACGCATATAAAAATTGCCGTTTCCCGACAGGGTATTGGTGGTCAACGTAATGAAATTTTCCGCACTTGCCGTCACTGATGAACGAGTCACTGGCTCATCCATCATAATCACATTTGCGTCAGACATCGCAAGGTTGGCTaa                                                                                                                                      >  1:247572/1‑149 (MQ=255)
            tatCGGTACGCATATAAAAATTGCCGTTTCCCGACAGGGTATTGGTGGTCAACGTAATGAAATTTTCCGCACTTGCCGTCACTGATGAACGAGTCACTGGCTCATCCATCATAATCACATTTGCGTCAGACATCGCAAGGTTGGCTAAg                                                                                                                                     <  3:429393/149‑1 (MQ=255)
            tatCGGTACGCATATAAAAATTGCCGTTTCCCGACAGGGTATTGGTGGTCAACGTAATGAAATTTTCCGCACTTGCCGTCACTGATGAACGAGTCACTGGCTCATCCATCATAATCACATTTGCGTCAGACATCGCAAGGTTGGCTAAg                                                                                                                                     <  1:456311/149‑1 (MQ=255)
            tatCGGTACGCATATAAAAATTGCCGTTTCCCGACAGGGTATTGGTGGTCAACGTAATGAAATTTTCCGCACTTGCCGTCACTGATGAACGAGTCACTGGCTCATCCATCATAATCACATTTGCGTCAGACATCGCAAGGTTGGCTAAg                                                                                                                                     <  3:337114/149‑1 (MQ=255)
                  tACGCATATAAAAATTGCCGTTTCCCGACAGGGTATTGGTGGTCAACGTAATGAAATTTTCCGCACTTGCCGTCACTGATGAACGAGTCACTGGCTCATCCATCATAATCACATTTGCGTCAGACATCGCAAGGTTGGCTAAGTTgagg                                                                                                                               >  1:219574/1‑149 (MQ=255)
                   aCGCATATAAAAATTGCCGTTTCCCGACAGGGTATTGGTGGTCAACGTAATGAAATTTTCCGCACTTGCCGTCACTGATGAACGAGTCACTGGCTCATCCATCATAATCACATTTGCGTCAGACATCGCAAGGTTGGCTAAGTTgagga                                                                                                                              <  1:60793/149‑1 (MQ=255)
                     gCATATAAAAATTGCCGTTTCCCGACAGGGTATTGGTGGTCAACGTAATGAAATTTTCCGCACTTGCCGTCACTGATGAACGAGTCACTGGCTCATCCATCATAATCACATTTGCGTCAGACATCGCAAGGTTGGCTAAGTTGaggagg                                                                                                                            <  3:229819/149‑1 (MQ=255)
                     gCATATAAAAATTGCCGTTTCCCGACAGGGTATTGGTGGTCAACGTAATGAAATTTTCCGCACTTGCCGTCACTGATGAACGAGTCACTGGCTCATCCATCATAATCACATTTGCGTCAGACATCGCAAGGTTGGCTAAGTTGaggagg                                                                                                                            <  3:528394/149‑1 (MQ=255)
                           aaaaaTTGCCGTTTCCCGACAGGGTATTGGTGGTCAACGTAATGAAATTTTCCGCACTTGCCGTCACTGATGAACGAGTCACTGGCTCATCCATCATAATCACATTTGCGTCAGACATCGCAAGGTTGGCTAAGTTGAGGAGGGTGttt                                                                                                                      >  2:397024/1‑149 (MQ=255)
                            aaaaTTGCCGTTTCCCGACAGGGTATTGGTGGTCAACGTAATGAAATTTTCCGCACTTGCCGTCACTGATGAACGAGTCACTGGCTCATCCATCATAATCACATTTGCGTCAGACATCGCAAGGTTGGCTAAGTTGAGGAGGGTGtttt                                                                                                                     >  2:88810/1‑149 (MQ=255)
                             aaaTTCCCGTTTCCCGACAGGGTATTGGTGGTCAACGTAATGAAATTTTCCGCACTTGCCGTCACTGATGAACGAGTCACTGGCTCATCCATCATAATCACATTTGCGTCAGACATCGCAAGGTTGGCTAAGTTGAGGAGGGTGTTTTg                                                                                                                    >  4:190212/1‑149 (MQ=255)
                                       ttCCCGACAGGGTATTGGTGGTCAACGTAATGAAATTTTCCGCACTTGCCGTCACTGATGAACGAGTCACTGGCTCATCCATCATAATCACATTTGCGTCAGACATCGCAAGGTTGGCTAAGTTGAGGAGGGTGTTTTGGTTCGTCGTa                                                                                                          <  1:88810/149‑1 (MQ=255)
                                            gACAGGGTATTGGTGGTCAACGTAATGAAATTTTCCGCACTTGCCGTCACTGATGAACGAGTCACTGGCTCATCCATCATAATCACATTTGCGTCAGACATCGCAAGGTTGGCTAAGTTGAGGAGGGTGTTTTGGTTCGTCGTAGTGat                                                                                                     >  3:92094/1‑149 (MQ=255)
                                                gggTATTGGTGGTCAACGTAATGAAATTTTCCGCACTTGCCGTCACTGATGAACGAGTCACTGGCTCATCCATCATAATCACATTTGCGTCAGACATCGCAAGGTTGGCTAAGTTGAGGAGGGTGTTTTGGTTCGTCGTAGTGATATCg                                                                                                 <  1:241271/149‑1 (MQ=255)
                                                 ggTATTGGTGGTCAACGTAATGAAATTTTCCGCACTTGCCGTCACTGATGAACGAGTCACTGGCTCATCCATCATAATCACATTTGCGTCAGACATCGCAAGGTTGGCTAAGTTGAGGAGGGTGTTTTGGTTCGTCGTAGTGATATCgg                                                                                                >  1:45313/1‑149 (MQ=255)
                                                   tATTGGTGGTCAACGTAATGAAATTTTCCGCACTTGCCGTCACTGATGAACGAGTCACTGGCTCATCCATCATAATCACATTTGCGTCAGACATCGCAAGGTTGGCTAAGTTgagga                                                                                                                              >  2:267097/1‑117 (MQ=255)
                                                   tATTGGTGGTCAACGTAATGAAATTTTCCGCACTTGCCGTCACTGATGAACGAGTCACTGGCTCATCCATCATAATCACATTTGCGTCAGACATCGCAAGGTTGGCTAAGTTgagga                                                                                                                              <  1:267097/117‑1 (MQ=255)
                                                       ggtggtCAACGTAATGAAATTTTCCGCACTTGCCGTCACTGATGAACGAGTCACTGGCTCATCCATCATAATCACATTTGCGTCAGACATCGCAAGGTTGGCTAAGTTGAGGAGGGTGTTTTGGTTCGTCGTAGTGATATCGGCAatga                                                                                          <  1:222741/149‑1 (MQ=255)
                                                        gtggtCAACGTAATGAAATTTTCCGCACTTGCCGTCACTGATGAACGAGTCACTGGCTCATCCATCATAATCACATTTGCGTCAGACATCGCAAGGTTGGCTAAGTTGAGGAGGgtg                                                                                                                         >  3:456600/1‑117 (MQ=255)
                                                        gtggtCAACGTAATGAAATTTTCCGCACTTGCCGTCACTGATGAACGAGTCACTGGCTCATCCATCATAATCACATTTGCGTCAGACATCGCAAGGTTGGCTAAGTTGAGGAGGgtg                                                                                                                         <  4:151587/117‑1 (MQ=255)
                                                        gtggtCAACGTAATGAAATTTTCCGCACTTGCCGTCACTGATGAACGAGTCACTGGCTCATCCATCATAATCACATTTGCGTCAGACATCGCAAGGTTGGCTAAGTTGAGGAGGgtg                                                                                                                         <  4:456600/117‑1 (MQ=255)
                                                        gtggtCAACGTAATGAAATTTTCCGCACTTGCCGTCACTGATGAACGAGTCACTGGCTCATCCATCATAATCACATTTGCGTCAGACATCGCAAGGTTGGCTAAGTTGAGGAGGgtg                                                                                                                         >  3:151587/1‑117 (MQ=255)
                                                           gtCAACGTAATGAAATTTTCCGCACTTGCCGTCACTGATGAACGAGTCACTGGCTCATCCATCATAATCACATTTGCGTCAGACATCGCAAGGTTGGCTAAGTTGAGGAGGGTGTTTTGGTTCGTCGTAGTGATATCGGCAATGATGGt                                                                                      <  4:25321/149‑1 (MQ=255)
                                                               aCGTAATGAAATTTTCCGCACTTGCCGTCACTGATGAACGAGTCACTGGCTCATCCATCATAATCACATTTGCGTCAGACATCGCAAGGTTGGCTAAGTTGAGGAGGGTGTTTTGGTTCGTCGTAGTGATATCGGCAATGATGGTCCAt                                                                                  <  2:45308/149‑1 (MQ=255)
                                                                   aaTGAAATTTTCCGCACTTGCCGTCACTGATGAACGAGTCACTGGCTCATCCATCATAATCACATTTGCGTCAGACATCGCAAGGTTGGCTAAGTTGAGGAGGGTGTTTTGGTTCGTCGTAGTGATATCGGCAATGATGGTCCATACGc                                                                              <  2:445524/149‑1 (MQ=255)
                                                                   aaTGAAATTTTCCGCACTTGCCGTCACTGATGAACGAGTCACTGGCTCATCCATCATAATCACATTTGCGTCAGACATCGCAAGGTTGGCTAAGTTGAGGAGGGTGTTTTGGTTCGTCGTAGTGATATCGGCAATGATGGTCCATACGc                                                                              <  4:139497/149‑1 (MQ=255)
                                                                   aaTGAAATTTTCCGCACTTGCCGTCACTGATGAACGAGTCACTGGCTCATCCATCATAATCACATTTGCGTCAGACATCGCAAGGTTGGCTAAGTTGAGGAGGGTGTTTTGGTTCGTCGTAGTGATATCGGCAATGATGGTCCATACGc                                                                              <  4:223191/149‑1 (MQ=255)
                                                                            ttCCGCACTTGCCGTCACTGATGAACGAGTCACTGGCTCATCCATCATAATCACATTTGCGTCAGACATCGCAAGGTTGGCTAAGTTGAGGAGGGTGTTTTGGTTCGTCGTAGTGATATCGGCAATGATGGTCCATACGCTATTTTCGa                                                                     <  4:301440/149‑1 (MQ=255)
                                                                                gCACTTGCCGTCACTGATGAACGAGTCACTGGCTCATCCATCATAATCACATTTGCGTCATACATCGCAAGGTTGGCTAAGTTGAGGAGGGTGTTTTGGTTCGTCTTAGTGATATCGGCAATGATGGTCCATACGCTATTTTCGAGCg                                                                  <  2:118716/148‑1 (MQ=255)
                                                                                gCACTTGCCGTCACTGATGAACGAGTCACTGGCTCATCCATCATAATCACATTTGCGTCAGACATCGCAAGGTTGGCTAAGTTGAGGAGGGTGTTTTGGTTCGTCGTAGTGATATCGGCAATGATGGTCCATACGCTATTTTCGAGCg                                                                  >  1:118716/1‑148 (MQ=255)
                                                                                gCACTTGCCGTCACTGATGAACGAGTCACTGGCTCATCCATCATAATCACATTTGCGTCAGACATCGCAAGGTTGGCTAAGTTGAGGAGGGTGTTTTGGTTCGTCGTAGTGATATCGGCAATGATGGTCCATACGCTATTTTCGAGCGa                                                                 >  3:285024/1‑149 (MQ=255)
                                                                                   cTTGCCGTCACTGATGAACGAGTCACTGGCTCATCCATCATAATCACATTTGCGTCAGACATCGCAAGGTTGGCTAAGTTGAGGAGGGTGTTTTGGTTCGTCGTAGTGATATCGGCAATGATGGTCCATACGCTATTTTCGAGCGAAAc                                                              <  3:457609/149‑1 (MQ=255)
                                                                                    ttGCCGTCACTGATGAACGAGTCACTGGCTCATCCATCATAATCACATTTGCGTCAGACATCGCAAGGTTGGCTAAGTTGAGGAGGGTGTTTTGGTTCGTCGTAGTGATATCGGCAATGATGGTCCATACGCTATTTTCGAGCGAAACg                                                             >  1:344105/1‑149 (MQ=255)
                                                                                          tCACTGATGAACGAGTCACTGGCTCATCCATCATAATCACATTTGCGTCAGACATCGCAAGGTTGGCTAAGTTGAGGAGGGTGTTTTGGTTCGTCGTAGTGATATCGGCAATGATGGTCCATACGCTATTTTCGAGCGAAACGTCCGCt                                                       >  1:522133/1‑149 (MQ=255)
                                                                                             cTGATGAACGAGTCACTGGCTCATCCATCATAATCACATTTGCGTCAGACATCGCAAGGTTGGCTAAGTTGAGGAGGGTGTTTTGGTTCGTCGTAGTGATATCGGCAATGATGGTCCATACGCTATTTTCGAGCGAAACGTCCGCTTTg                                                    <  2:62578/149‑1 (MQ=255)
                                                                                             cTGATGAACGAGTCACTGGCTCATCCATCATAATCACATTTGCGTCAGACATCGCAAGGTTGGCTAAGTTGAGGAGGGTGTTTTGGTTCGTCGTAGTGATATCGGCAATGATGGTCCATACGCTATTTTCGAGCGAAACGTCCGCTTTg                                                    <  4:59685/149‑1 (MQ=255)
                                                                                              tgatgaACGAGTCACTGGCTCATCCATCATAATCACATTTGCGTCAGACATCGCAAGGTTGGCTAAGTTGAGGAGGGTGTTTTGGTTCGTCGTAGTGATATCGGCAATGATGGTCCATACGCTATTTTCGAGCGAAACGTCCGCTTTGc                                                   >  1:448381/1‑149 (MQ=255)
                                                                                              tgatgaACGAGTCACTGGCTCATCCATCATAATCACATTTGCGTCAGACATCGCAAGGTTGGCTAAGTTGAGGAGGGTGTTTTGGTTCGTCGTAGTGATATCGGCAATGATGGTCCATACGCTATTTTCGAGCGAAACGTCCGCTTTGc                                                   >  4:230342/1‑149 (MQ=255)
                                                                                               gatgaACGAGTCACTGGCTCATCCATCATAATCACATTTGCGTCAGACATCGCAAGGTTGGCTAAGTTGAGGAGGGTGTTTTGGTTCGTCGTAGTGATATCGGCAATGATGGTCCATACGCTATTTTCGAGCGAAACGTCCGCTTTGCt                                                  <  4:550403/149‑1 (MQ=255)
                                                                                                atgaACGAGTCACTGGCTCATCCATCATAATCACATTTGCGTCAGACATCGCAAGGTTGGCTAAGTTGAGGAGGGTGTTTTGGTTCGTCGTAGTGATATCGGCAATGATGGTCCATACGCTATTTTCGAGCGAAACGTCCGCTTTGCTg                                                 <  1:188692/149‑1 (MQ=255)
                                                                                                       aGTCACTGGCTCATCCATCATAATCACATTTGCGTCAGACATCGCAAGGTTGGCTAAGTTGAGGAGGGTGTTTTGGTTCGTCGTAGTGATATCGGCAATGATGGTCCATACGCTATTTTCGAGCGAAACGTCCGCTTTGCTGGTATCCa                                          <  4:529198/149‑1 (MQ=255)
                                                                                                         gcACTGGCTCATCCATCATAATCACATTTGCGTCAGACATCGCAAGGTTGGCTAAGTTGAGGAGGGTGTTTTGGTTCGTCGTAGTGATATCGGCAATGATGGTCCATACGCTATTTTCGAGCGAAACGTCCGCTTTGCTGGTATCCACg                                        >  2:502286/2‑149 (MQ=255)
                                                                                                              ggCTCATCCATCATAATCACATTTGCGTCAGACATCGCAAGGTTGGCTAAGTTGAGGAGGGTGTTTTGGTTCGTCGTAGTGATATCGGCAATGATGGTCCATACGCTATTTTCGAGCGAAACGTCCGCTTTGCTGGTATCCACGGCAcc                                   >  1:75850/1‑149 (MQ=255)
                                                                                                                 tCATCCATCATAATCACATTTGCGTCAGACATCGCAAGGTTGGCTAAGTTGAGGAGGGTGTTTTGGTTCGTCGTAGTGATATCGGCAATGATGGTCCATACGCTATTTTCGAGCGAAACGTCCGCTTTGCTGGTATCCACGGCACCATg                                >  1:111153/1‑149 (MQ=255)
                                                                                                                 tCATCCATCATAATCACATTTGCGTCAGACATCGCAAGGTTGGCTAAGTTGAGGAGGGTGTTTTGGTTCGTCGTAGTGATATCGGCAATGATGGTCCATACGCTATTTTCGAGCGAAACGTCCGCTTTGCTGGTATCCACGGCACCATg                                >  1:98183/1‑149 (MQ=255)
                                                                                                                 tCATCCATCATAATCACATTTGCGTCAGACATCGCAAGGTTGGCTAAGTTGAGGAGGGTGTTTTGGTTCGTCGTAGTGATATCGGCAATGATGGTCCATACGCTATTTTCGAGCGAAACGTCCGCTTTGCTGGTATCCACGGCACCATg                                >  1:55416/1‑149 (MQ=255)
                                                                                                                   atccatcATAATCACATTTGCGTCAGACATCGCAAGGTTGGCTAAGTTGAGGAGGGTGTTTTGGTTCGTCGTAGTGATATCGGCAATGATGGTCCATACGCTATTTTCGAGCGAAACGTCCGCTTTGCTGGTATCCACGGCACCATGCg                              >  3:285442/1‑149 (MQ=255)
                                                                                                                   atccatcATAATCACATTTGCGTCAGACATCGCAAGGTTGGCTAAGTTGAGGAGGGTGTTTTGGTTCGTCGTAGTGATATCGGCAATGATGGTCCATACGCTATTTTCGAGCGAAACGTCCGCTTTGCTGGTATCCACGGCACCATGCg                              >  4:212857/1‑149 (MQ=255)
                                                                                                                    tccatcATAATCACATTTGCGTCAGACATCGCAAGGTTGGCTAAGTTGAGGAGGGTGTTTTGGTTCGTCGTAGTGATATCGGCAATGATGGTCCATACGCTATTTTCGAGCGAAACGTCCGCTTTGCTGGTATCCACGGCACCATGCGt                             >  4:533599/1‑149 (MQ=255)
                                                                                                                     ccatcATAATCACATTTGCGTCAGACATCGCAAGGTTGGCTAAGTTGAGGAGGGTGTTTTGGTTCGTCGTAGTGATATCGGCACTGATGGTCCATACGCTATTTTCGAGCTAAACGTCCGCTTTGCTGGTATCCACGGCACCATGCGtt                            <  1:448820/149‑1 (MQ=255)
                                                                                                                     ccatcATAATCACATTTGCGTCAGACATCGCAAGGTTGGCTAAGTTGAGGAGGGTGCTTTGGTTCGTCGTAGTGATATCGGCAATGATGGTCCATACGCTATTTTCGAGCGAAACGTCCGCTTTGCTGGTATCCACGGCACCATGCGtt                            <  1:276472/149‑1 (MQ=255)
                                                                                                                       atcatAATCACATTTGCGTCAGACATCGAAAGGTTGGCTAAGTTGAGGAGGGTGTTTTGGTTCGTCGTAGTGATATCGGCAATGATGGTCCATACGCTATTTTCGAGCGAAACGTCCGCTTTGCTGGTATCCAc                                         >  4:218685/1‑134 (MQ=255)
                                                                                                                       atcatAATCACATTTGCGTCAGACATCGAAAGGTTGGCTAAGTTGAGGAGGGTGTTTTGGTTCGTCGTAGTGATATCGGCAATGATGGTCCATACGCTATTTTCGAGCGAAACGTCCGCTTTGCTGGTATCCAc                                         <  3:218685/134‑1 (MQ=255)
                                                                                                                       atcatAATCACATTTGCGTCAGACATCGAAAGGTTGGCTAAGTTGAGGAGGGTGTTTTGGTTCGTCGTAGTGATATCGGCAATGATGGTCCATACGCTATTTTCGAGCGAAACGTCCGCTTTGCTGGTATCCAc                                         >  4:150511/1‑134 (MQ=255)
                                                                                                                       atcatAATCACATTTGCGTCAGACATCGAAAGGTTGGCTAAGTTGAGGAGGGTGTTTTGGTTCGTCGTAGTGATATCGGCAATGATGGTCCATACGCTATTTTCGAGCGAAACGTCCGCTTTGCTGGTATCCAc                                         <  3:150511/134‑1 (MQ=255)
                                                                                                                          atAATCACATTTGCGTCAGACATCGCAAGGTTGGCTAAGTTGAGGAGGGTGTTTTGGTTCGTCGTAGTGATATCGGCAATGATGGTCCATACGCTATTTTCGAGCGAAACGTCCGCTTTGCTGGTATCCACGGCACCATGCGTTCTaaa                       <  4:92094/149‑1 (MQ=255)
                                                                                                                          atAATCACATTTGCGTCAGACATCGCAAGGTTGGCTAAGTTGAGGAGGGTGTTTTGGTTCGTCGTAGTGATATCGGCAATGATGGTCCATACGCTATTTTCGAGCGAAACGTCCGCTTTGCTGGTATCCACGGCACCATGCGTTCTaaa                       >  3:410488/1‑149 (MQ=255)
                                                                                                                              tCACATTTGCGTCAGACATCGCAAGGTTGGCTAAGTTGAGGAGGGTGTTTTGGTTCGTCGTAGTGATATCGGCAATGATGGTCCATACGCTATTTTCGAGCGAAACGTCCGCTTTGCTGGTATCCACGGCACCATGCGTTCTAAaagaa                   <  4:215012/149‑1 (MQ=255)
                                                                                                                              tCACATTTGCGTCAGACATCGCAAGGTTGGCTAAGTTGAGGAGGGTGTTTTGGTTCGTCGTAGTGATATCGGCAATGATGGTCCATACGCTATTTTCGAGCGAAACGTCCGCTTTGCTGGTATCCACGGCACCATGCGTTCTAAaagaa                   <  2:111153/149‑1 (MQ=255)
                                                                                                                              tCACATTTGCGTCAGACATCGCAAGGTTGGCTAAGTTGAGGAGGGTGTTTTGGTTCGTCGTAGTGATATCGGCAATGATGGTCCATACGCTATTTTCGAGCGAAACGTCCGCTTTGCTGGTATCCACGGCACCATGCGTTCTAAaagaa                   <  2:55416/149‑1 (MQ=255)
                                                                                                                              tCACATTTGCGTCAGACATCGCAAGGTTGGCTAAGTTGAGGAGGGTGTTTTGGTTCGTCGTAGTGATATCGGCAATGATGGTCCATACGCTATTTTCGAGCGAAACGTCCGCTTTGCTGGTATCCACGGCACCATGCGTTCTAAaagaa                   <  1:397024/149‑1 (MQ=255)
                                                                                                                                  aTTTGCGTCAGACATCGCAAGGTTGGCTAAGTTGAGGAGGGTGTTTTGGTTCGTCGTAGTGATATCGGCAATGATGGTCCATACGCTATTTTCGAGCGAAACGTCCGCTTTGCTGGTATCCACGGCACCATGCGTTCTAAAAGAAgcgc               >  1:533939/1‑149 (MQ=255)
                                                                                                                                  aTTTGCGTCAGACATCGCAAGGTTGGCTAAGTTGAGGAGGGTGTTTTGGTTCGTCGTAGTGATATCGGCAATGATGGTCCATACGCTATTTTCGAGCGAAACGTCCGCTTTGCTGGTATCCACGGCACCATGCGTTCTAAAAGAAgcgc               <  4:382198/149‑1 (MQ=255)
                                                                                                                                        gTCAGACATCGCAAGGTTGGCTAAGTTGAGGAGGGTGTTTTGGTTCGTCGTAGTGATATCGGCAATGATGGTCCATACGCTATTTTCGAGCGAAACGTCCGCTTTGCTGGTATCCACGGCACCATGCGTTCTAAAAGAAGCGCCCTCAg         >  1:490660/1‑149 (MQ=255)
                                                                                                                                            gACATCGCAAGGTTGGCTAAGTTGAGGAGAGTGTTTTGGTTCGTCGTAGTGATATCGGCAATGATGGTCCATACGCTATTTTCGAGCGAAACGTCCGCTTTGCTGGTATCCACGGCACCATGCGTTCTAAAAGAAGCGCCCTCAGAAAc     <  2:533939/149‑1 (MQ=255)
                                                                                                                                             aCATCGCAAGGTTGGCTAAGTTGAGGAGGGTGTTTTGGTTCGTCGTAGTGATATCGGCAATGATGGTCCATACGCTATTTTCGAGCGAAACGTCCGCTTTGCTGGTATCCACGGCACCATGCGTTCTAAAAGAAGCGCCCTCAGAAAcc    >  1:127908/1‑149 (MQ=255)
                                                                                                                                               aTCGCAAGGTTGGCTAGGTAGAGGAGGGTGTTTTGGTTCGTCGTAGTGATATCGGCAATGATGGTCCGTACGCTATTTTCGAGCGAAACGTCCGCTTTGCTGGTATCCACGGCACCATGCGCTCTAAAAGAAGCTCCCTCAGAAACCAc  >  2:45848/1‑149 (MQ=255)
                                                                                                                                               aTCGCAAGGTTGGCTAAGTTGAGGAGGGTGTTTTGGTTCGTCGTAGTGATATCGGCAATGATGGTCCATACGCTATTTTCGAGCGAAACGTCCGCTTTGCTGGTATCCACGGCACCATGCGTTCTAAAAGAAGCGCCCTCAGAAACCAc  <  4:285442/149‑1 (MQ=255)
                                                                                                                                               aTCGCAAGGTTGGCTAAGTTGAGGAGGGTGTTTTGGTTCGTCGTAGTGATATCGGCAATGATGGTCCATACGCTATTTTCGAGCGAAACGTCCGCTTTGCTGGTATCCACGGCACCATGCGTTCTAAAAGAAGCGCCCTCAGAAACCAc  <  2:344105/149‑1 (MQ=255)
                                                                                                                                                  gCAAGGTTGGCTAAGTTGAGGAGGGTGTTTTGGTTCGTCGTAGTGATATCGGCAATGATGGTCCATACGCTATTTTCGAGCGAAACGTCCGCttt                                                     >  4:108549/1‑95 (MQ=255)
                                                                                                                                                  gCAAGGTTGGCTAAGTTGAGGAGGGTGTTTTGGTTCGTCGTAGTGATATCGGCAATGATGGTCCATACGCTATTTTCGAGCGAAACGTCCGCttt                                                     <  3:108549/95‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                
GATGATTAGCCATATCGGTACGCATATAAAAATTGCCGTTTCCCGACAGGGTATTGGTGGTCAACGTAATGAAATTTTCCGCACTTGCCGTCACTGATGAACGAGTCACTGGCTCATCCATCATAATCACATTTGCGTCAGACATCGCAAGGTTGGCTAAGTTGAGGAGGGTGTTTTGGTTCGTCGTAGTGATATCGGCAATGATGGTCCATACGCTATTTTCGAGCGAAACGTCCGCTTTGCTGGTATCCACGGCACCATGCGTTCTAAAAGAAGCGCCCTCAGAAACCAC  >  NC_000913/2287535‑2287826

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