Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,111,962 G→A 5.1% pseudogene (129/963 nt) yghE ← pseudogene, secretion pathway protein, L‑type protein homology; putative general secretion pathway for protein export (GSP)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,111,9620GA 5.1% 256.3 / 3.5 79pseudogene (129/963 nt)yghEpseudogene, secretion pathway protein, L‑type protein homology; putative general secretion pathway for protein export (GSP)
Reads supporting (aligned to +/- strand):  ref base G (36/38);  new base A (2/2);  total (39/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.80e-01

GCCAGCCCGGCGCTCCATCCCCGGCATTCATCGATGCGGCAGATAACCCGTTGCGCATCGCAAGCCATAAACCCACAAGGAATGGACATCCAGTCCGGCGCGACGGTAGCGCGGGTGATGCCGTTTTCCTGCAACCACTGCGCAATGTTGCGCATATGCTGTTGGTGAATCACTGCTACGGTTGCCAGTTGCTGGTCGATATTCAACGGGGCGAAATGCAGTTCATCGATATCCTGGTTCAGCTCTTCTTCCAGTAAGGCGGGCAGAATAGTCGGTAGCGCAGCGCGATGGAAAT  >  NC_000913/3111816‑3112110
                                                                                                                                                  |                                                                                                                                                    
gccagcCCGGCGCTCCATCCCCGGCATTCATCGATGCGGCAGATAACCCGTTGCGCATCGCAAGCCATAAACCCACAAGGAATGGACATCCAGTCCGGCGCGACGGTAGCGCGGGTGATGCCGTTTTCCTGCAACCACTGCGCAATGtt                                                                                                                                                    <  1:484863/149‑1 (MQ=255)
       cGGCGCTCCATCCCCGGCATTCATCGATGCGGCAGATAACCCGTTGCGCATCGCAAGCCATAAACCCACAAGGAATGGACATCCAGTCCGGCGCGACGGTAGCGCGGGTGATGCCGTTTTCCTGCAACCACTGCGCAACATTGCGCagt                                                                                                                                             <  2:189904/149‑3 (MQ=255)
             tCCATCCCCGGCATTCATCGATGCGGCAGATAACCCGTTGCGCATCGCAAGCCATAAACCCACAAGGAATGGACATCCAGTCCGGCGCGACGGTAGCGCGGGTGATGCCGTTTTCCTGCAACCACTGCGCAATGTTGCGCATATGCTGt                                                                                                                                       >  2:532181/1‑149 (MQ=255)
             tCCATCCCCGGCATTCATCGATGCGGCAGATAACCCGTTGCGCATCGCAAGCCATAAACCCACAAGGAATGGACATCCAGTCCGGCGCGACGGTAGCGCGGGTGATGCCGTTTTCCTGCAACCACTGCGCAACATTGCGCagtggt                                                                                                                                          <  3:326346/146‑6 (MQ=255)
             tCCATCCCCGGCATTCATCGATGCGGCAGATAACCCGTTGCGCATCGCAAGCCATAAACCCACAAGGAATGGACATCCAGTCCGGCGCGACGGTAGCGCGGGTGATGCCGTTTTCCTGCAACCACTGCGCAACATTGCGCagtggt                                                                                                                                          >  4:326346/1‑141 (MQ=255)
              ccATCCCCGGCATTCATCGATGCGGCAGATAACCCGTTGCGCATCGCAAGCCATAAACCCACAAGGAATGGACATCCAGTCCGGCGCGACGGTAGCGCGGGTGATGCCGTTTTCCTGCAACCACTGCGCAATGTTGCGCATATGCTGtt                                                                                                                                      >  2:255788/1‑149 (MQ=255)
               cATCCCCGGCATTCATCGATGCGGCAGATAACCCGTTGCGCATCGCAAGCCATAAACCCACAAGGAATGGACATCCAGTCCGGCGCGACGGTAGCGCGGGTGATGCCGTTTTCCTGCAACCACTGCGCAATGTTGCGCATATGCTGTTg                                                                                                                                     <  2:470095/149‑1 (MQ=255)
                aTCCCCGGCATTCATCGATGCGGCAGATAACCCGTTGCGCATCGCAAGCCATAAACCCACAAGGAATGGACATCCAGTCCGGCGCGACGGTAGCGCGGGTGATGCCGTTTTCCTGCAACCACTGCGCAATGTTGCGCATATGCTGTTgg                                                                                                                                    <  4:51285/149‑1 (MQ=255)
                 tCCCCGGCATTCATCGATGCGGCAGATAACCCGTTGCGCATCGCAAGCCATAAACCCACAAGGAATGGACATCCAGTCCGGCGCGACGGTAGCGCGGGTGATGCCGTTTTCCTGCAACCACTGCGCAATGTTGCGCATATGCTGTTGGt                                                                                                                                   <  3:516590/149‑1 (MQ=255)
                    ccGGCATTCATCGATGCGGCAGATAACCCGTTGCGCATCGCAAGCCATAAACCCACAAGGAATGGACATCCAGTCCGGCGCGACGGTAGCGCGGGTGATGCCGTTTTCCTGCAACCACTGCGCAATGTTGCGCATATGCTGTTGGTGaa                                                                                                                                >  1:108660/1‑149 (MQ=255)
                    ccGGCATTCATCGATGCGGCAGATAACCCGTTGCGCATCGCAAGCCATAAACCCACAAGGAATGGACATCCAGTCCGGCGCGACGGTAGCGCGGGTGATGCCGTTTTCCTGCAACCACTGCGCAATGTTGCGCATATGCTGTTGGTGaa                                                                                                                                >  1:216881/1‑149 (MQ=255)
                     cGGCATTCATCGATGCGGCAGATAACCCGTTGCGCATCGCAAGCCATAAACCCACAAGGAATGGACATCCAGTCCGGCGCGACGGTAGCGCGGGTGATGCCGTTTTCCTGCAACCACTGCGCAATGTTGCGCATATGCTGTTGGTGAAt                                                                                                                               <  3:624354/149‑1 (MQ=255)
                     cGGCATTCATCGATGCGGCAGATAACCCGTTGCGCATCGCAAGCCATAAACCCACAAGGAATGGACATCCAGTCCGGCGCGACGGTAGCGCGGGTGATGCCGTTTTCCTGCAACCACTGCGCAATGTTGCGCATATGCTATTGGTGAAt                                                                                                                               >  3:583200/1‑149 (MQ=255)
                      ggCATTCATCGATGCGGCAGATAACCCGTTGCGCATCGCAAGCCATAAACCCACAAGGAATGGACATCCAGTCCGGCGCGACGGTAGCGCGGGTGATGCCGTTTTCCTGCAACCACTGCGCAATGTTGCGCATATGCTGTTGGTGAATc                                                                                                                              <  4:385195/149‑1 (MQ=255)
                       gCATTCATCGATGCGGCAGATAACCCGTTGCGCATCGCAAGCCATAAACCCACAAGGAATGGACATCCAGTCCGGCGCGACGGTAGCGCGGGTGATGCCGTTTTCCTGCAACCACTGCGCAATGTTGCGCATATGCTGTTGGTGAATCa                                                                                                                             <  4:507643/149‑1 (MQ=255)
                         aTTCATCGATGCGGCAGATAACCCGTTGCGCATCGCAAGCCATAAACCCACAAGGAATGGACATCCAGTCCGGCGCGACGGTAGCGCGGGTGATGCCGTTTTCCTGCAACCACTGCGCAATGTTGCGCATATGCTGTTGGTGAATCACt                                                                                                                           <  1:50013/149‑1 (MQ=255)
                         aTTCATCGATGCGGCAGATAACCCGTTGCGCATCGCAAGCCATAAACCCACAAGGAATGGACATCCAGTCCGGCGCGACGGTAGCGCGGGTGATGCCGTTTTCCTGCAACCACTGCGCAATGTTGCGCATATGCTGCTGGTGAATCACt                                                                                                                           >  2:393965/1‑149 (MQ=255)
                          ttCATCGATGCGGCAGATAACCCGTTGCGCATCGCAAGCCATAAACCCACAAGGAATGGACATCCAGTCCGGCGCGACGGTAGCGCGGGTGATGCCGTTTTCCTGCAACCACTGCGCAATGTTGCGCATATGCTGTTGGTGAATCACTg                                                                                                                          <  1:289137/149‑1 (MQ=255)
                          ttCATCGATGCGGCAGATAACCCGTTGCGCATCGCAAGCCATAAACCCACAAGGAATGGACATCCAGTCCGGCGCGACGGTAGCGCGGGTGATGCCGTTTTCCTGCAACCACTGCGCAATGTTGCGCATATGCTGTTGGTGAATCACTg                                                                                                                          >  2:601104/1‑149 (MQ=255)
                             aTCGATGCGGCAGATAACCCGTTGCGCATCGCAAGCCATAAACCCACAAGGAATGGACATCCAGTCCGGCGCGACGGTAGCGCGGGTGATGCCGTTTTCCTGCAACCACTGCGCAATGTTGCGCATATGCTGTTGGTGAATCACTGCTa                                                                                                                       >  2:378162/1‑149 (MQ=255)
                             aTCGATGCGGCAGATAACCCGTTGCGCATCGCAAGCCATAAACCCACAAGGAATGGACATCCAGTCCGGCGCGACGGTAGCGCGGGTGATGCCGTTTTCCTGCAACCACTGCGCAATGTTGCGCATATGCTGTTGGTGAATCACTGCTa                                                                                                                       <  3:230364/149‑1 (MQ=255)
                                             aCCCGTTGCGCATCGCAAGCCATAAACCCACAAGGAATGGACATCCAGTCCGGCGCGACGGTAGCGCGGGTGATGCCGTTTTCCTGCAACCACTGCGCAATGTTGCGCATATGCTGTTGGTGAATCACTGCTACGGTTGCCAGTTGCTg                                                                                                       <  2:145782/149‑1 (MQ=255)
                                              cccGTTGCGCATCGCAAGCCATAAACCCACAAGGAATGGACATCCAGTCCGGCGCGACGGTAGCGCGGGTGATGCCGTTTTCCTGCAACCACTGCGCAATGTTGCGCATATGCTGTTGGTGAATCACTGCTACGGTTGCCAGTTGCTgg                                                                                                      <  1:463183/149‑1 (MQ=255)
                                               ccGTTGCGCATCGCAAGCCATAAACCCACAAGGAATGGACATCCAGTCCGGCGCGACGGTAGCGCGGGTGATGCCGTTTTCCTGCAACCACTGCGCAATGTTGCGCATATGCTGTTGGTGAATCACTGCTACGGTTGCCATTTGCTGGt                                                                                                     >  4:275734/1‑149 (MQ=255)
                                               ccGTTGCGCATCGCAAGCCATAAACCCACAAGGAATGGACATCCAGTCCGGCGCGACGGTAGCGCGGGTGATGCCGTTTTCCTGCAACCACTGCGCAATGTTGCGCATATGCTGTTGGTGAATCACTGCTACGGTTGCCAGTTGCTGGt                                                                                                     >  4:664526/1‑149 (MQ=255)
                                                cGTTGCGCATCGCAAGCCATAAACCCACAAGGAATGGACATCCAGTCCGGCGCGACGGTAGCGCGGGTGATGCCGTTTTCCTGCAACCACTGCGCAATGTTGCGCATATGCTGTTGGTGAATCACTGCTACGGTTGCCAGTTGCTGGTc                                                                                                    <  2:43790/149‑1 (MQ=255)
                                                    gcgcATCGCAAGCCATAAACCCACAAGGAATGGACATCCAGTCCGGCGCGACGGTAGCGCGGGTGATGCCGTTTTCCTGCAACCACTGCGCAATGTTGCGCATATGCTGTTGGTGAATCACTGCTACGGTTGCCAGTTGCTGGTCGata                                                                                                <  2:108660/149‑1 (MQ=255)
                                                      gcATCGCAAGCCATAAACCCACAAGGAATGGACATCCAGTCCGGCGCGACGGTAGCGCGGGTGATGCCGTTTTCCTGCAACCACTGCGCAATGTTGCGCATATGCTGTTGGTGAATCACTGCTACGGTTGCCAGTTGCTGGTCGATAtt                                                                                              >  2:79164/1‑149 (MQ=255)
                                                         tCGCAAGCCATAAACCCACAAGGAATGGACATCCAGTCCGGCGCGACGGTAGCGCGGGTGATGCCGTTTTCCTGCAACCACTGCGCAATGTTGCGCATATGCTGTTGGTGAATCACTGCTACGGTTGCCAGTTGCTGGTCGATATTCaa                                                                                           >  1:42065/1‑149 (MQ=255)
                                                            cAAGCCATAAACCCACAAGGAATGGACATCCAGTCCGGCGCGACGGTAGCGCGGGTGATGCCGTTTTCCTGCAACCACTGCGCAATGTTGCGCATATGCTGTTGGTGAATCACTGCTACGGTTGCCAGTTGc                                                                                                         <  3:396203/132‑1 (MQ=255)
                                                            cAAGCCATAAACCCACAAGGAATGGACATCCAGTCCGGCGCGACGGTAGCGCGGGTGATGCCGTTTTCCTGCAACCACTGCGCAATGTTGCGCATATGCTGTTGGTGAATCACTGCTACGGTTGCCAGTTGc                                                                                                         >  4:396203/1‑132 (MQ=255)
                                                              aGCCATAAACCCACAAGGAATGGACATCCAGTCCGGCGCGACGGTAGCGCGGGTGATGCCGTTTTCCTGCAACCACTGCGCAATGTTGCGCATATGCTGTTGGTGAATCACTGCTACGGTTGCCAGTTGCTGGTCGATATTCAACgggg                                                                                      <  1:242204/149‑1 (MQ=255)
                                                                 cATAAACCCACAAGGAATGGACATCCAGTCCGGCGCGACGGTAGCGCGGGTGATGCCGTTTTCCTGCAACCACTGCGCAATGTTGCGCATATGCTGTTGGTGAATCACTGCTACGGTTGCCAGTTGCTGGTCGATATTCAACGGGGCGa                                                                                   >  3:657582/1‑149 (MQ=255)
                                                                        ccACAAGGAATGGACATCCAGTCCGGCGCGACGGTAGCGCGGGTGATGCCGTTTTCCTGCAACCACTGCGCAATGTTGCGCATATGCTGTTGGTGAATCACTGCTACGGTTGCCAGTTGCTGGTCGATATTCAACGGGGCGAAATGCAg                                                                            >  2:313318/1‑149 (MQ=255)
                                                                        ccACAAGGAATGGACATCCAGTCCGGCGCGACGGTAGCGCGGGTGATGCCGTTTTCCTGCAACCACTGCGCAATGTTGCGCATATGCTGTTGGTGAATCACTGCTACGGTTGCCAGTTGCTGGTCGATATTCAACGGGGCGAAATGCAg                                                                            >  4:503932/1‑149 (MQ=255)
                                                                         cacaAGGAATGGACATCCAGTCCGGCGCGACGGTAGCGCGGGTGATGCCGTTTTCCTGCAACCACTGCGCAATGTTGCGCATATGCTGTTGGTGAATCACTGCTACGGTTGCCAGTTGCTGGTCGATATTCAACGGGGCGAAATGCAGt                                                                           >  3:11448/1‑149 (MQ=255)
                                                                          acaAGGAATGGACATCCAGTCCGGCGCGACGGTAGCGCGGGTGATGCCGTTTTCCTGCAACCACTGCGCAATGTTGCGCATATTCTGTTGGTGAATCACTGCTACTGTTGCCAGTTGCTGGTCGATATTCAGCGGGGCGAAATGCAGtt                                                                          <  2:385713/149‑1 (MQ=255)
                                                                          acaAGGAATGGACATCCAGTCCGGCGCGACGGTAGCGCGGGTGATGCCGTTTTCCTGCAACCACTGCGCAATGTTGCGCATATGCTGTTGGTGAATCACTGCTACGGTTGCCAGTTGCTGGTCGATATTCAACGGGGCGAAATGCAGtt                                                                          >  3:227864/1‑149 (MQ=255)
                                                                           caAGGAATGGACATCCAGTCCGGCGCGACGGTAGCGCGGGTGATGCCGTTTTCCTGCAACCACTGCGCAATGTTGCGCATATGCTGTTGGTGAATCACTGCTACGGTTGCCAGTTGCTGGTCGATATTCAACGGGGCGAAATGCAGTTc                                                                         >  4:214039/1‑149 (MQ=255)
                                                                            aaGGAATGGACATCCAGTCCGGCGCGACGGTAGCGCGGGTGATGCCGTTTTCCTGCAACCACTGCGCAATGTTGCGCATATGCTGTTGGTGAATCACTGCTACGGTTGCCAGTTGCTGGTCGATATTCAACGGGGCGAAATGCAGTTCa                                                                        >  2:466503/1‑149 (MQ=255)
                                                                             aGGAATGGACATCCAGTCCGGCGCGACGGTAGCGCGGGTGATGCCGTTTTCCTGCAACCACTGCGCAATGTTGCGCATATGCTGTTGGTGAATCACTGCTACGGTTGCCAGTTGCTGGTCGATATTCAACGGGGCGAAATGCAGTTCAt                                                                       <  4:11448/149‑1 (MQ=255)
                                                                              ggAATGGACATCCAGTCCGGCGCGACGGTAGCGCGGGTGATGCCGTTTTCCTGCAACCACTGCGCAATGTTGCGCATATGCTGTTGGTGAATCACTGCTACGGTTGCCAGTTGCTGGTCGATATTCAACGGGGCGAAATGCAGTTCATc                                                                      <  4:67300/149‑1 (MQ=255)
                                                                              ggAATGGACATCCAGTCCGGCGCGACGGTAGCGCGGGTGATGCCGTTTTCCTGCAACCACTGCGCAATGTTGCGCATATGCTGTTGGTGAATCACTGCTACGGTTGCCAGTTGCTGGTCGATATTCAACGGGGCGAAATGCAGTTCATc                                                                      <  2:322648/149‑1 (MQ=255)
                                                                                    gACATCCAGTCCGGCGCGACGGTAGCGCGGGTGATGCCGTTTTCCTGCAACCACTGCGCAATGTTGCGCATATGCTGTTGGTGAATCACTGCTACGGTTGCCAGTTGCTGGTCGATATTCAACGGGGCGAAATGCAGTTCATCGATATc                                                                >  4:446899/1‑149 (MQ=255)
                                                                                        tCCAGTCCGGCGCGACGGTAGCGCGGGTGATGCCGTTTTCCTGCAACCACTGCGCAATGTTGCGCATATGCTGTTGGTGAATCACTGCTACGGTTGCCAGTTGCTGGTCGATATTCAACGGGGCGAAATGCAGTTCATCGATATCCTgg                                                            >  3:481080/1‑149 (MQ=255)
                                                                                        tCCAGTCCGGCGCGACGGTAGCGCGGGTGATGCCGTTTTCCTGCAACCACTGCGCAATGTTGCGCATATGCTGTTGGTGAATCACTGCTACGGTTGCCAGTTGCTGGTCGATATTCAACGGGGCGAAATGCAGTTCATCGATATCCTgg                                                            >  1:120884/1‑149 (MQ=255)
                                                                                        tCCAGTCCGGCGCGACGGTAGCGCGGGTGATGCCGTTTTCCTGCAA‑CACTGCGGCATTTTGCGCATATGGTTTTGGTGGATCAATGCTTAGGTTGCCCGTTGCTGGTCGATATTCAACCGGGGGGAATGCAGTTTCTTGGTATCCTggg                                                           >  3:257205/1‑148 (MQ=255)
                                                                                         ccAGTCCGGCGCGACGGTAGCGCGGGTGATGCCGTTTTCCTGCAACCACTGCGCAATGTTGCGCATATGCTGTTGGTGAATCACTGCTACGGTTGCCAGTTGCTGGTCGATATTCAACGGGGCGAAATGCAGTTCATCGATATCCTGGt                                                           <  2:216881/149‑1 (MQ=255)
                                                                                         ccAGTCCGGCGCGACGGTAGCGCGGGTGATGCCGTTTTCCTGCAACCACTGCGCAATGTTGCGCATATGCTGTTGGTGAATCACTGCTACGGTTGCCAGTTGCTGGTCGATATTCAACGGGGCGAAATGCAGTTCATCGATATCCTGGt                                                           >  1:445755/1‑149 (MQ=255)
                                                                                         acAGTCCGGCGCGACGGTAGCGCGGGTGATGCCGTTTTCCTGCAACCACTGCGCAATGTTGCGCATATGCTGTTGGTGAATCACTGCTACGGTTGCCAGTTGCTGGTCGATATTCAACGGGGCGAAATGCAGTTCATCGATATCCTGGt                                                           >  2:586779/2‑149 (MQ=255)
                                                                                           aGTCCGGCGCGACGGTAGCGCGGGTGATGCCGTTTTCCTGCAACCACTGCGCAATGTTGCGCATATGCTGTTGGTGAATCACTGCTACGGTTGCCAGTTGCTGGTCGATATTCAACGGGGCGAAATGCAGTTCATCGATATCCTGGTTc                                                         >  4:183381/1‑149 (MQ=255)
                                                                                            gTCCGGCGCGACGGTAGCGCGGGTGATGCCGTTTTCCTGCAACCACTGCGCAATGTTGCGCATATGCTGTTGGTGAATCACTGCTACGGTTGCCAGTTGCTGGTCGATATTCAACGGGGCGAAATGCAGTTCATCGATATCCTGGTTCa                                                        <  3:183381/149‑1 (MQ=255)
                                                                                               cGGCGCGACGGTAGCGCGGGTGATGCCGTTTTCCTGCAACCACTGCGCAATGTTGCGCATATGCTGTTGGTGAATCACTGCTACGGTTGCCAGTTGCTGGTCGATATTCAACGGGGCGAAATGCAGTTCATCGATATCCTGGTTCAGct                                                     <  3:275734/149‑1 (MQ=255)
                                                                                                ggCGCGACGGTAGCGCGGGTGATGCCGTTTTCCTGCAACCACTGCGCAATGTTGCGCATATGCTGTTGGTGAATCACTGCTACGGTTGCCAGTTGCTGGTCGATATTCAACGGGGCGAAATGCAGTTCATCGATATCCTGGTTCAGctc                                                    <  3:446899/149‑1 (MQ=255)
                                                                                                   gcgACGGTAGCGCGGGTGATGCCGTTTTCCTGCAACCACTGCGCAATGTTGCGCATATGCTGTTGGTGAATCACTGCTACGGTTGCCAGTTGCTGGTCGATATTCAACGGGGCGAAATGCAGTTCATCGATATCCTGGTTCAGCtcttc                                                 <  1:255788/149‑1 (MQ=255)
                                                                                                   gcgACGGTAGCGCGGGTGATGCCGTTTTCCTGCAACCACTGCGCAATGTTGCGCATATGCTGTTGGTGAATCACTGCTACGGTTGCCAGTTGCTGGTCGATATTCAACGGGGCGAAATGCAGTTCATCGATATCCTGGTTCAGCtcttc                                                 >  1:68244/1‑149 (MQ=255)
                                                                                                      aCGGTAGCGCGGGTGATGCCGTTTTCCTGCAACCACTGCGCAATGTTGCGCATATGCTGTTGGTGAATCACTGCTACGGTTGCCAGTTGCTGGTCGATATTCAACGGGGCGAAATGCAGTTCATCGATATCCTGGTTCAGCTCttcttc                                              <  2:271145/149‑1 (MQ=255)
                                                                                                       cGGTAGCGCGGGTGATGCCGTTTTCCTGCAACCACTGCGCAATGTTGCGCATATGCTGTTGGTGAATCACTGCTACGGTTGCCAGTTGCTGGTCGATATTCAACGGGGCGAAATGCAGTTCATCGATATCCTGGTTCAGCTCTTCTTcc                                             <  1:466503/149‑1 (MQ=255)
                                                                                                        ggTAGCGCGGGTGATGCGGTTTTCCTGCAACCACTGCGCAATGTTGCGCTTATGCTGTTGGTGAATCACTGCTACGGTTGCCAGTTGCTGGTCGATATTCAACGGGGCGAAATGCAGTTCATCGATATCCTGGTTCAGCTCTTCTTCCa                                            <  4:583200/149‑1 (MQ=255)
                                                                                                        ggTAGCGCGGGTGATGCCGTTTTCCTGCAACCACTGCGCAATGTTGCGCATATGCTGTTGGTGAATCACTGCTACGGTTGCCAGTTGCTGGTCGATATTCAACGGGGCGAAATGCAGTTCATCGATATCCTGGTTCAGCTCTTCTTCCa                                            <  3:503932/149‑1 (MQ=255)
                                                                                                        ggTAGCGCGGGTGATGCCGTTTTCCTGCAACCACTGCGCAATGTTGCGCATATGCTGTTGGTGAATCACTGCTACGGTTGCCAGTTGCTGGTCGATATTCAACGGGGCGAAATGCAGTTCATCGATATCCTGGTTCAGCTCTTCTTCCa                                            <  1:601104/149‑1 (MQ=255)
                                                                                                        ggTAGCGCGGGTGATGCCGTTTTCCTGCAACCACTGCGCAATGTTGCGCATATGCTGTTGGTGAATCACTGCTACGGTTGCCAGTTGCTGGTCGATATTCAACGGGGCGAAATGCAGTTCATCGATATCCTGGTTCAGCTCTTCTTCCa                                            <  3:214039/149‑1 (MQ=255)
                                                                                                             cgcgGGTGATGCCGTTTTCCTGCAACCACTGCGCAATGTTGCGCATATGCTGTTGGTGAATCACTGCTACGGTTGCCAGTTGCTGGTCGATATTCAACGGGGCGaa                                                                                  <  2:617130/106‑1 (MQ=255)
                                                                                                             cgcgGGTGATGCCGTTTTCCTGCAACCACTGCGCAATGTTGCGCATATGCTGTTGGTGAATCACTGCTACGGTTGCCAGTTGCTGGTCGATATTCAACGGGGCGaa                                                                                  >  1:617130/1‑106 (MQ=255)
                                                                                                                gggTGATGCCGTTTTCCTGCAACCACTGCGCAATGTTGCGCATATGCTGTTGGTGAATCACTGCTACGGTTGCCAGTTGCTGGTCGATATTCAACGGGGCGAAATGCAGTTCATCGATATCCTGGTTCAGCTCTTCTTCCAGTAAGGCg                                    <  3:358096/149‑1 (MQ=255)
                                                                                                                      tGCCGTTTTCCTGCAACCACTGCGCAATGTTGCGCATATGCTGTTGGTGAATCACTGCTACGGTTGCCAGTTGCTGGTCGATATTCAACggg                                                                                       >  2:194679/1‑92 (MQ=255)
                                                                                                                      tGCCGTTTTCCTGCAACCACTGCGCAATGTTGCGCATATGCTGTTGGTGAATCACTGCTACGGTTGCCAGTTGCTGGTCGATATTCAACggg                                                                                       <  1:194679/92‑1 (MQ=255)
                                                                                                                       gCCGTTTTCCTGCAACCACTGCGCAATGTTGCGCATATGCTGTTGGTGAATCACTGCTACGGTTGCCAGTTGCTGGTc                                                                                                    >  3:645494/1‑78 (MQ=255)
                                                                                                                       gCCGTTTTCCTGCAACCACTGCGCAATGTTGCGCATATGCTGTTGGTGAATCACTGCTACGGTTGCCAGTTGCTGGTc                                                                                                    <  4:645494/78‑1 (MQ=255)
                                                                                                                       gCCGTTTTCCTGCAACCACTGCGCAATGTTGCGCATATGCTGTTGGTGAATCACTGCTACGGTTGCCAGTTGCTGGTCGATATTCAACGGGGCGAAATGCAGTTCATCGATATCCTGGTTCAGCTCTTCTTCCAGTAAGGCGGGCAGaa                             >  1:14387/1‑149 (MQ=255)
                                                                                                                          gTTTTCCTGCAACCACTGCGCAATGTTGCGCATATGCTGTTGGTGAATCACTGCTACGGTTGCCAGTTGCTGGTCGATATTCAACGGGGCGAAATGCAGTTCATCGATATCCTGGTTCAGCTCTTCTTCCAGTAAGGCGGGCAGAATAg                          >  2:460057/1‑149 (MQ=255)
                                                                                                                            tttCCTTCAACCACTGCGCAATGTTGCGCATATGCTGTTGGTGAATCACTGATACGGTTGCCAGTTGCTGGTCGATATTCAACGGGGCGAAATGCAGTTCATCGATATCCTGGTTCAGCTCTTCTTCCAGTAAGGCGGGCAGAATAGTc                        <  1:404129/149‑1 (MQ=255)
                                                                                                                                cTGCAACCACTGCGCAATGTTGCGCATATGCTGTTGGTGAATCACTGCTACGGTTGCCAGTTGCTGGTCGATATTCAACGGGGCGAAATGCAGTTCATCGATATCCTGGTTCAGCTCTTCTTCCAGTAAGGCGGGCAGAATAGTCGGTa                    >  2:483062/1‑149 (MQ=255)
                                                                                                                                cTGCAACCACTGCGCAATGTTGCGCATATGCTGTTGGTGAATCACTGCTACGGTTGCCAGTTGCTGGTCGATATTCAACGGGGCGAAATGCAGTTCATCGATATCCTGGTTCAGCTCTTCTTCCAGTAAGGCGGGCAGAATAGTCGGTa                    <  2:68244/149‑1 (MQ=255)
                                                                                                                                    aaCCACTGCGCAATGTTGCGCATATGCTGTTGGTGAATCACTGCTACGGTTGCCAGTTGCTGGTCGATATTCAACGGGGCGAAATGCAGTTCATCGATATCCTGGTTCAGCTCTTCTTCCAGTAAGGCGGGCAGAATAGTCGGTagcgc                >  2:156992/1‑149 (MQ=255)
                                                                                                                                     aCCACTGCGCAATGTTGCGCATATGCTGTTGGTGAATCACTGCTACGGTTGCCAGTTGCTGGTCGATATTCAACGGGGCGAAATGCAGTTCATCGATATCCTGGTTCAGCTCTTCTTCCAGTAAGGCGGGCAGAATAGTCGGTagcgca               >  4:66319/1‑149 (MQ=255)
                                                                                                                                       cACTGCGCAATGTTGCGCATATGCTGTTGGTGAATCACTGCTACGGTTGCCAGTTGCTGGTCGATATTCAACGGGGCGAAATGCAGTTCATCGATATCCTGGTTCAGCTCTTCTTCCAGTAAGGCGGGCAGAATAGTCGGTagcgcagc             <  1:460057/149‑1 (MQ=255)
                                                                                                                                       cACTGCGCAATGTTGCGCATATGCTGTTGGTGAATCACTGCTACGGTTGCCAGTTGCTGGTCGATATTCAACGGGGCGAAATGCAGTTCATCGATATCCTGGTTCAGCTCTTCTTCCAGTAAGGCGGGCAGAATAGTCGGTagcgcagc             >  1:174319/1‑149 (MQ=255)
                                                                                                                                        tatGCGCAACATTGCGCATATGCTGTTGGTGAATCACTGCTACGGTTGCCAGTTGCTGGTCGATATTCAACGGGGCGAAATGCAGTTCATCGATATCCTGGTTCAGCTCTTCTTCCAGTAAGGCGGGCAGAATAGTCGGTagcgcagcg            >  3:65169/3‑149 (MQ=255)
                                                                                                                                                  gTTGCGCATATGCTGTTGGTGAATCACTGCTACGGTTGCCAGTTGCTGGTCGATATTCAACGGGGCGAAATGCAGTTCATCGATATCCTGGTTCAGCTCTTCTTCCAGTAAGGCGGGCAGAATAGTCGGTAGCGCAGCGCGATGGAAAt  <  2:445755/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
GCCAGCCCGGCGCTCCATCCCCGGCATTCATCGATGCGGCAGATAACCCGTTGCGCATCGCAAGCCATAAACCCACAAGGAATGGACATCCAGTCCGGCGCGACGGTAGCGCGGGTGATGCCGTTTTCCTGCAACCACTGCGCAATGTTGCGCATATGCTGTTGGTGAATCACTGCTACGGTTGCCAGTTGCTGGTCGATATTCAACGGGGCGAAATGCAGTTCATCGATATCCTGGTTCAGCTCTTCTTCCAGTAAGGCGGGCAGAATAGTCGGTAGCGCAGCGCGATGGAAAT  >  NC_000913/3111816‑3112110

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