Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,228,558 T→C pseudogene (201/1436 nt) mdtQ ← putative multidrug resistance outer membrane protein MdtQ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,228,5580TC98.7% 269.7 / ‑4.8 78pseudogene (201/1436 nt)mdtQputative multidrug resistance outer membrane protein MdtQ
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (44/33);  total (45/33)
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.99e-01

CTCCATATCCGCCGAAAAGTCGATCTGCGGACCATCCTGAGTGGCAACTGCCTTCGCCTGGGCTTCAGCTAGCTGAATGCGCTGTTCAGCCACCTGCATATCCGGTGCGTTCTGTAGTGCATTGTTAATTAAGGAAGTGAGTTGATTATCGTGATACTCCAGCCACCATTGGCTGTCTGGCCAACCATTTTCAGCGCCGTGGGTAATGCGGTGTCAACTTGTGCAGCGGGCGTTTGCTGGCTTAACGCCTGGCGGGTTTCATGCATAGGCGCAC  >  NC_000913/2228421‑2228694
                                                                                                                                         |                                                                                                                                        
cTCCATATCCGCCGAAAAGTCGATCTGCGGACCATCCTGAGTGGCAACTGCCTTCGCCTGGGCTTCAGCTAGCTGAATGCGCTGTTCAGCCACCTGCATATCCGGTGCGTTCTGTAGTGCATTGTTAATTAAGGAAGCGa                                                                                                                                        <  1:720620/140‑1 (MQ=255)
   cATATCCGCCGAAAAGTCGATCTGCGGACCATCCTGAGTGGCAACTGCCTTCGCCTGGGCTTCAGCTAGCTGAATGCGCTGTTCAGCCACCTGCATATCCGGTGCGTTCTGTAGTGCATTGTTAATTAAGGAAGCGAGtt                                                                                                                                     >  1:755015/1‑140 (MQ=255)
             gAAAAGTCGATCTGCGGACCATCCTGAGTGGCAACTGCCTTCGCCTGGGCTTCAGCTAGCTGAATGCGCTGTTCAGCCACCTGCATATCCGGTGCGTTCTGTAGTGCATTGTTAATTAAGGAAGCGAGTTGATTATCGTg                                                                                                                           >  2:589759/1‑140 (MQ=255)
                  gTCGATCTGCGGACCATCCTGAGTGGCAACTGCCTTCGCCTGGGCTTCAGCTAGCTGAATGCGCTGTTCAGCCACCTGCATATCCGGTGCGTTCTGTAGTGCATTGTTAATTAAGGAAGCGAGTTGATTATCGTGATACt                                                                                                                      <  2:717684/140‑1 (MQ=255)
                      aTCTGCGGACCATCCTGAGTGGCAACTGCCTTCGCCTGGGCTTCAGCTAGCTGAATGCGCTGTTCAGCCACCTGCATATCCGGTGCGTTCTGTAGTGCATTGTTAATTAAGGAAGCGAGTTGATTATCGTGATACTCCAg                                                                                                                  <  2:755015/140‑1 (MQ=255)
                             gACCATCCTGAGTGGCAACTGCCTTCGCCTGGGCTTCAGCTAGCTGAATGCGCTGTTCAGCCACCTGCATATCCGGTGCGTTCTGTAGTGCATTGTTAATTAAGGAAGCGAGTTGATTATCGTGATACTCCAGCCACCAt                                                                                                           >  2:801368/1‑140 (MQ=255)
                                cATCCTGAGTGGCAACTGCCTTCGCCTGGGCTTCAGCTAGCTGAATGCGCTGTTCAGCCACCTGCATATCCGGTGCGTTCTGTAGTGCATTGTTAATTAAGGAAGCGAGTTGATTATCGTGATACTCCAGCCACCATTgg                                                                                                        <  2:148631/140‑1 (MQ=255)
                                 aTCCTGAGTGGCAACTGCCTTCGCCTGGGCTTCAGCTAGCTGAATGCGCTGTTCAGCCACCTGCATATCCGGTGCGTTCTGTAGTGCATTGTTAATTAAGGAAGTGAGTTGATTATCGTGATACTCCAGCCACCATTGGc                                                                                                       >  2:395215/1‑140 (MQ=255)
                                  tCCTGAGTGGCAACTGCCTTCGCCTGGGCTTCAGCTAGCTGAATGCGCTGTTCAGCCACCTGCATATCCGGTGCGTTCTGTAGTGCATTGTTAATTAAGGAAGCGAGTTGATTATCGTGATACTCCa                                                                                                                   >  2:378646/1‑127 (MQ=255)
                                  tCCTGAGTGGCAACTGCCTTCGCCTGGGCTTCAGCTAGCTGAATGCGCTGTTCAGCCACCTGCATATCCGGTGCGTTCTGTAGTGCATTGTTAATTAAGGAAGCGAGTTGATTATCGTGATACTCCa                                                                                                                   <  1:378646/127‑1 (MQ=255)
                                  tCCTGAGTGGCAACTGCCTTCGCCTGGGCTTCAGCTAGCTGAATGCGCTGTTCAGCCACCTGCATATCCGGTGCGTTCTGTAGTGCATTGTTAATTAAGGAAGCGAGTTGATTATCGTGATACTCCAGCCACCATTGGCt                                                                                                      >  2:1056715/1‑140 (MQ=255)
                                   ccTGAGTGGCAACTGCCTTCGCCTGGGCTTCAGCTAGCTGAATGCGCTGTTCAGCCACCTGCATATCCGGTGCGTTCTGTAGTGCATTGTTAATTAAGGAAGCGAGTTGATTATCGTGATACTCCAGCCACCATTGGCTg                                                                                                     >  1:553890/1‑140 (MQ=255)
                                    cTGAGTGGCAACTGCCTTCGCCTGGGCTTCAGCTAGCTGAATGCGCTGTTCAGCCACCTGCATATCCGGTGCGTTCTGTAGTGCATTGTTAATTAAGGAAGCGAGTTGATTATCGTGATACTCCAGCCACCATTGGCTGt                                                                                                    <  1:801368/140‑1 (MQ=255)
                                       aGTGGCAACTGCCTTCGCCTGGGCTTCAGCTAGCTGAATGCGCTGTTCAGCCACCTGCATATCCGGTGCGTTCTGTAGTGCATTGTTAATTAAGGAAGCGAGTTGATTATCGTGATACTCCAGCCACCATTGGCTGTCTg                                                                                                 >  2:187549/1‑140 (MQ=255)
                                         tGGCAACTGCCTTCGCCTGGGCTTCAGCTAGCTGAATGCGCTGTTCAGCCACCTGCATATCCGGTGCGTTCTGTAGTGCATTGTTAATTAAGGAAGCGAGTTGATTATCGTGATACTCCAGCCACCATTGGCTGTCTGGc                                                                                               <  1:14455/140‑1 (MQ=255)
                                          ggCAACTGCCTTCGCCTGGGCTTCAGCTAGCTGAATGCGCTGTTCAGCCACCTGCATATCCGGTGCGTTCTGTAGTGCATTGTTAATTAAGGAAGCGAGTTGATTATCGTGATACTCCAGCCACCATTGGCTGTCTGGcc                                                                                              >  1:1022218/1‑140 (MQ=255)
                                            cAACTGCCTTCGCCTGGGCTTCAGCTAGCTGAATGCGCTGTTCAGCCACCTGCATATCCGGTGCGTTCTGTAGTGCATTGTTAATTAAGGAAGCGAGTTGATTATCGTGATACTCCAGCCACCATTg                                                                                                         <  2:449135/127‑1 (MQ=255)
                                            cAACTGCCTTCGCCTGGGCTTCAGCTAGCTGAATGCGCTGTTCAGCCACCTGCATATCCGGTGCGTTCTGTAGTGCATTGTTAATTAAGGAAGCGAGTTGATTATCGTGATACTCCAGCCACCATTg                                                                                                         >  1:449135/1‑127 (MQ=255)
                                                 gCCTTCGCCTGGGCTTCAGCTAGCTGAATGCGCTGTTCAGCCACCTGCATATCCGGTGCGTTCTGTAGTGCATTGTTAATTAAGGAAGCGAGTTGATTATCGTGATACTCCAGCCACCATTGGCTGTCTGGCCAACCAtt                                                                                       >  2:701681/1‑140 (MQ=255)
                                                  ccTTCGCCTGGGCTTCAGCTAGCTGAATGCGCTGTTCAGCCACCTGCATATCCGGTGCGTTCTGTAGTGCATTGTTAATTAAGGAAGCGAGTTGATTATCGTGATACTCCAGCCACCATTGGCTGTCTGGCCAACCAttt                                                                                      >  2:958243/1‑140 (MQ=255)
                                                    ttCGCCTGGGCTTCAGCTAGCTGAATGCGCTGTTCAGCCACCTGCATATCCGGTGCGTTCTGTAGTGCATTGTTAATTAAGGAAGCGAGTTGATTATCGTGATACTCCAGCCACCATTGGCTGTCTGGCCAACCATTTTc                                                                                    >  2:24063/1‑140 (MQ=255)
                                                     tCGCCTGGGCTTCAGCTAGCTGAATGCGCTGTTCAGCCACCTGCATATCCGGTGCGTTCTGTAGTGCATTGTTAATTAAGGAAGCGAGTTGATTATCGTGATACTCCAGCCACCATTGGCTGTCTGGCCAACCATTTTCa                                                                                   <  1:187549/140‑1 (MQ=255)
                                                     tCGCCTGGGCTTCAGCTAGCTGAATGCGCTGTTCAGCCACCTGCATATCCGGTGCGTTCTGTAGTGCATTGTTAATTAAGGAAGCGAGTTGATTATCGTGATACTCCAGCCACCATTGGCTGTCTGGCCAACCATTTTCa                                                                                   <  2:734423/140‑1 (MQ=255)
                                                          tGGGCTTCAGCTAGCTGAATGCGCTGTTCAGCCACCTGCATATCCGGTGCGTTCTGTAGTGCATTGTTAATTAAGGAAGCGAGTTGATTATCGTGATACTCCAGCCACCATTGGCTGTCTGGCCAACCATTTTCAGCGcc                                                                              <  1:701681/140‑1 (MQ=255)
                                                          tGGGCTTCAGCTAGCTGAATGCGCTGTTCAGCCACCTGCATATCCGGTGCGTTCTGTAGTGCATTGTTAATTAAGGAAGCGAGTTGATTATCGTGATACTCCAGCCACCATTGGCTGTCTGGCCAACCATTTTCAGCGcc                                                                              <  1:24063/140‑1 (MQ=255)
                                                          tGGGCTTCAGCTAGCTGAATGCGCTGTTCAGCCACCTGCATATCCGGTGCGTTCTGTAGTGCATTGTTAATTAAGGAAGCGAGTTGATTATCGTGATACTCCAGCCACCATTGGCTGTCTGGCCAACCATTTTCAGCGcc                                                                              <  2:957308/140‑1 (MQ=255)
                                                            ggCTTCAGCTAGCTGAATGCGCTGTTCAGCCACCTGCATATCCGGTGCGTTCTGTAGTGCATTGTTAATTAAGGAAGCGAGTTGATTATCGTGATACTCCAGCCACCATTGGCTGTCTGGCCAACCATTTTCAGCGCCGt                                                                            <  2:553890/140‑1 (MQ=255)
                                                                 cAGCTAGCTGAATGCGCTGTTCAGCCACCTGCATATCCGGTGCGTTCTGTAGTGCATTGTTAATTAAGGAAGCGAGTTGATTATCGTGATACTCCAGCCACCATTGGCTGTCTGGCCAACCATTTTCAGCGCCGTGGGTa                                                                       >  2:598991/1‑140 (MQ=255)
                                                                 cAGCTAGCTGAATGCGCTGTTCAGCCACCTGCATATCCGGTGCGTTCTGTAGTGCATTGTTAATTAAGGAAGCGAGTTGATTATCGTGATACTCCAGCCACCATTGGCTGTCTGGCCAACCATTTTCAGCGCCGTGGGTa                                                                       >  1:1016029/1‑140 (MQ=255)
                                                                 cAGCTAGCTGAATGCGCTGTTCAGCCACCTGCATATCCGGTGCGTTCTGTAGTGCATTGTTAATTAAGGAAGCGAGTTGATTATCGTGATACTCCAGCCACCATTGGCTGTCTGGCCAACCATTTTCAGCGCCGTGGGTa                                                                       >  1:509432/1‑140 (MQ=255)
                                                                    ctagctGAATGCGCTGTTCAGCCACCTGCATATCCGGTGCGTTCTGTAGTGCATTGTTAATTAAGGAAGCGAGTTGATTATCGTGATACTCCAGCCACCATTGGCTGTCTGGCCAACCATTTTCAGCGCCGTGGGTAATg                                                                    <  2:1016029/140‑1 (MQ=255)
                                                                    ctagctGAATGCGCTGTTCAGCCACCTGCATATCCGGTGCGTTCTGTAGTGCATTGTTAATTAAGGAAGCGAGTTGATTATCGTGATACTCCAGCCACCATTGGCTGTCTGGCCAACCATTTTCAGCGCCGTGGGTAATg                                                                    <  2:509432/140‑1 (MQ=255)
                                                                      agctGAATGCGCTGTTCAGCCACCTGCATATCCGGTGCGTTCTGTAGTGCATTGTTAATTAAGGAAGCGAGTTGATTATCGTGATACTCCAGCCACCATTGGCTGTCTGGCCAACCATTTTCAGCGCCGTGGGTAATGCg                                                                  >  1:788429/1‑140 (MQ=255)
                                                                       gctGAATGCGCTGTTCAGCCACCTGCATATCCGGTGCGTTCTGTAGTGCATTGTTAATTAAGGAAGCGAGTTGATTATCGTGATACTCCAGCCACCATTGGCTGTCTGGCCAACCATTTTCAGCGCCGTGGGTAATg                                                                    <  2:1253043/137‑1 (MQ=255)
                                                                       gctGAATGCGCTGTTCAGCCACCTGCATATCCGGTGCGTTCTGTAGTGCATTGTTAATTAAGGAAGCGAGTTGATTATCGTGATACTCCAGCCACCATTGGCTGTCTGGCCAACCATTTTCAGCGCCGTGGGTAATg                                                                    >  1:1253043/1‑137 (MQ=255)
                                                                         tGAATGCGCTGTTCAGCCACCTGCATATCCGGTGCGTTCTGTAGTGCATTGTTAATTAAGGAAGCGAGTTGATTATCGTGATACTCCAGCCACCATTGGCTGTCTGGCCAACCATTTTCAGCGCCGTGGGTAATGCGgtg                                                               <  1:364675/140‑1 (MQ=255)
                                                                           aaTGCGCTGTTCAGCCACCTGCATATCCGGTGCGTTCTGTAGTGCATTGTTAATTAAGGAAGCGAGTTGATTATCGTGATACTCCAGCCACCATTGGCTGTCTGGCCAACCATTTTCAGCGCCGTGGGTAATGCGGTGTc                                                             <  1:589759/140‑1 (MQ=255)
                                                                            aTGCGCTGTTCAGCCACCTGCATATCCGGTGCGTTCTGTAGTGCATTGTTAATTAAGGAAGCGAGTTGATTATCGTGATACTCCAGCCACCATTGGCTGTCTGGCCAACCATTTTCAGCGCCGTGGGTAATGCgg                                                                 >  2:6320/1‑135 (MQ=255)
                                                                            aTGCGCTGTTCAGCCACCTGCATATCCGGTGCGTTCTGTAGTGCATTGTTAATTAAGGAAGCGAGTTGATTATCGTGATACTCCAGCCACCATTGGCTGTCTGGCCAACCATTTTCAGCGCCGTGGGTAATGCgg                                                                 <  1:6320/135‑1 (MQ=255)
                                                                                  tGTTCAGCCACCTGCATATCCGGTGCGTTCTGTAGTGCATTGTTAATTAAGGAAGCGAGTTGATTATCGTGATACTCCAGCCACCATTGGCTGTCTGGCCAACCATTTTCAGCGCCGTGGGTAATGCGGTGTCAACTtgt                                                      >  1:428457/1‑140 (MQ=255)
                                                                                     tCAGCCACCTGCATATCCGGTGCGTTCTGTAGTGCATTGTTAATTAAGGAAGCGAGTTGATTATCGTGATACTCCAGCCACCATTGGCTGTCTGGCCAACCATTTTCAGCGCCGTGGGTAATGCGgtg                                                               >  1:1146057/1‑128 (MQ=255)
                                                                                     tCAGCCACCTGCATATCCGGTGCGTTCTGTAGTGCATTGTTAATTAAGGAAGCGAGTTGATTATCGTGATACTCCAGCCACCATTGGCTGTCTGGCCAACCATTTTCAGCGCCGTGGGTAATGCGgtg                                                               <  2:1146057/128‑1 (MQ=255)
                                                                                          cACCTGCATATCCGGTGCGTTCTGTAGTGCATTGTTAATTAAGGAAGCGAGTTGATTATCGTGATACTCCAGCCACCATTGGCTGTCTGGCCAACCATTTTCAGCGCCGTGGGTAATGCGGTGTCAACTTGTGCAGCggg                                              <  2:396183/140‑1 (MQ=255)
                                                                                              tGCATATCCGGTGCGTTCTGTAGTGCATTGTTAATTAAGGAAGCGAGTTGATTATCGTGATACTCCAGCCACCATTGGCTGTCTGGCCAACCATTTTCAGCGCCGTGGGTAATGCGGTGTCAACTTGTGCAGCGGGCGtt                                          >  1:65619/1‑140 (MQ=255)
                                                                                              tGCATATCCGGTGCGTTCTGTAGTGCATTGTTAATTAAGGAAGCGAGTTGATTATCGTGATACTCCAGCCACCATTGGCTGTCTGGCCAACCATTTTCAGCGCCGTGGGTAATGCGGTGTCAACTTGTGCAGCGGGCGtt                                          >  1:317237/1‑140 (MQ=255)
                                                                                              tGCATATCCGGTGCGTTCTGTAGTGCATTGTTAATTAAGGAAGCGAGTTGATTATCGTGATACTCCAGCCACCATTGGCTGTCTGGCCAACCATTTTCAGCGCCGTGGGTAATGCGGTGTCAACTTGTGCAGCGGGCGtt                                          >  2:82015/1‑140 (MQ=255)
                                                                                                cATATCCGGTGCGTTCTGTAGTGCATTGTTAATTAAGGAAGCGAGTTGATTATCGTGATACTCCAGCCACCATTGGCTGTCTGGCCAACCATTTTCAGCGCCGTGGGTAATGCGGTGTCAACTTGTGCAGCGGGCGTTTg                                        <  1:1056715/140‑1 (MQ=255)
                                                                                                  tatCCGGTGCGTTCTGTAGTGCATTGTTAATTAAGGAAGCGAGTTGATTATCGTGATACTCCAGCCACCATTGGCTGTCTGGCCAACCATTTTCAGCGCCGTGGGTAATGCGGTGTCAACTTGTGCAGCGGGCGTTTGCt                                      <  2:344049/140‑1 (MQ=255)
                                                                                                       ggTGCGTTCTGTAGTGCATTGTTAATTAAGGAAGCGAGTTGATTATCGTGATACTCCAGCCACCATTGGCTGTCTGGCCAACCATTTTCAGCGCCGTGGGTAATGCGGTGTCAACTTGTGCAGCGGGCGTTTGCt                                      <  2:709759/135‑1 (MQ=255)
                                                                                                       ggTGCGTTCTGTAGTGCATTGTTAATTAAGGAAGCGAGTTGATTATCGTGATACTCCAGCCACCATTGGCTGTCTGGCCAACCATTTTCAGCGCCGTGGGTAATGCGGTGTCAACTTGTGCAGCGGGCGTTTGCt                                      >  1:709759/1‑135 (MQ=255)
                                                                                                       ggTGCGTTCTGTAGTGCATTGTTAATTAAGGAAGCGAGTTGATTATCGTGATACTCCAGCCACCATTGGCTGTCTGGCCAACCATTTTCAGCGCCGTGGGTAATGCGGTGTCAACTTGTGCAGCGGGCGTTTGCTGGCtt                                 <  2:788429/140‑1 (MQ=255)
                                                                                                        gTGCGTTCTGTAGTGCATTGTTAATTAAGGAAGCGAGTTGATTATCGTGATACTCCAGCCACCATTGGCTGTCTGGCCAACCATTTTCAGCGCCGTGGGTAATGCGGTGTCAACTTGTGCAGCGGGCGTTTGCTGGCTTa                                >  1:625180/1‑140 (MQ=255)
                                                                                                        gTGCGTTCTGTAGTGCATTGTTAATTAAGGAAGCGAGTTGATTATCGTGATACTCCAGCCACCATTGGCTGTCTGGCCAACCATTTTCAGCGCCGTGGGTAATGCGGTGTCAACTTGTGCAGCGGGCGTTTGCTGGCTTa                                >  2:608678/1‑140 (MQ=255)
                                                                                                        gTGCGTTCTGTAGTGCATTGTTAATTAAGGAAGCGAGTTGATTATCGTGATACTCCAGCCACCATTGGCTGTCTGGCCAACCATTTTCAGCGCCGTGGGTAATGCGGTGTCAACTTGTGCAGCGGGCGTTTGCTGGCTTa                                <  2:695834/140‑1 (MQ=255)
                                                                                                        gTGCGTTCTGTAGTGCATTGTTAATTAAGGAAGCGAGTTGATTATCGTGATACTCCAGCCACCATTGGCTGTCTGGCCAACCATTTTCAGCGCCGTGGGTAATGCGGTGTCAACTTGTGCAGCGGGCGTTTGCTGGCTTa                                >  1:695834/1‑140 (MQ=255)
                                                                                                         tGCGTTCTGTAGTGCATTGTTAATTAAGGAAGCGAGTTGATTATCGTGATACTCCAGCCACCATTGGCTGTCTGGCCAACCATTTTCAGCGCCGTGGGTAATGCGGTGTCAACTTGTGCAGCGGGCGTTTGCTGGCtt                                 <  2:137079/138‑1 (MQ=255)
                                                                                                         tGCGTTCTGTAGTGCATTGTTAATTAAGGAAGCGAGTTGATTATCGTGATACTCCAGCCACCATTGGCTGTCTGGCCAACCATTTTCAGCGCCGTGGGTAATGCGGTGTCAACTTGTGCAGCGGGCGTTTGCTGGCtt                                 >  1:137079/1‑138 (MQ=255)
                                                                                                          gCGTTCTGTAGTGCATTGTTAATTAAGGAAGCGAGTTGATTATCGTGATACTCCAGCCACCATTGGCTGTCTGGCCAACCATTTTCAGCGCCGTGGGTAATGCGGTGTCAACTTGTGCAGCGGGCGTTTGCTGGCTTAAc                              >  1:1089861/1‑140 (MQ=255)
                                                                                                          gCGTTCTGTAGTGCATTGTTAATTAAGGAAGCGAGTTGATTATCGTGATACTCCAGCCACCATTGGCTGTCTGGCCAACCATTTTCAGCGCCGTGGGTAATGCGGTGTCAACTTGTGCAGCGGGCGTTTGCTGGCTTAAc                              >  1:1043915/1‑140 (MQ=255)
                                                                                                                  tAGTGCATTGTTAATTAAGGAAGCGAGTTGATTATCGTGATACTCCAGCCACCATTGGCTGTCTGGCCAACCATTTTCAGCGCCGTGGGTAATGCGGTGTCAACTTGTGCAGCGGGCGTTTGCTGGCTTAACGCCTGGCg                      >  1:127566/1‑140 (MQ=255)
                                                                                                                  tAGTGCATTGTTAATTAAGGAAGCGAGTTGATTATCGTGATACTCCAGCCACCATTGGCTGTCTGGCCAACCATTTTCAGCGCCGTGGGTAATGCGGTGTCAACTTGTGCAGCGGGCGTTTGCTGGCTTAACGCCTGGCg                      >  2:909865/1‑140 (MQ=255)
                                                                                                                     tGCATTGTTAATTAAGGAAGCGAGTTGATTATCGTGATACTCCAGCCACCATTGGCTGTCTGGCCAACCATTTTCAGCGCCGTGGGTAATGCGGTGTCAACTTGTGCAGCGGGCGTTTGCTGGCTTAACGCCTGGCGGGt                   >  1:729290/1‑140 (MQ=255)
                                                                                                                         ttGTTAATTAAGGAAGCGAGTTGATTATCGTGATACTCCAGCCACCATTGGCTGTCTGGCCAACCATTTTCAGCGCCGTGGGTAATGCGGTGTCAACTTGTGCAGCGGGCGTTTGCTGGCTTAACGCCTGGCGGGTTTCa               >  2:882725/1‑140 (MQ=255)
                                                                                                                          tGTTAATTAAGGAAGCGAGTTGATTATCGTGATACTCCAGCCACCATTGGCTGTCTGGCCAACCATTTTCAGCGCCGTGGGTAATGCGGTGTCAACTTGTGCAGCGGGCGTTTGCTGGCTTAACGCCTGGCGGGTTTCAt              >  1:570735/1‑140 (MQ=255)
                                                                                                                          tGTTAATTAAGGAAGCGAGTTGATTATCGTGATACTCCAGCCACCATTGGCTGTCTGGCCAACCATTTTCAGCGCCGTGGGTAATGCGGTGTCAACTTGTGCAGCGGGCGTTTGCTGGCTTAACGCCTGGCGGGTTTCAt              >  1:912620/1‑140 (MQ=255)
                                                                                                                            ttaattaaGGAAGCGAGTTGATTATCGTGATACTCCAGCCACCATTGGCTGTCTGGCCAACCATTTTCAGCGCCGTGGGTAATGCGGTGTCAACTTGTGCAGCGGGCGTTTGCTGGCTTAACGCCTGGCGGGTTTCATGc            >  1:128447/1‑140 (MQ=255)
                                                                                                                            ttaattaaGGAAGCGAGTTGATTATCGTGATACTCCAGCCACCATTGGCTGTCTGGCCAACCATTTTCAGCGCCGTGGGTAATGCGGTGTCAACTTGTGCAGCGGGCGTTTGCTGGCTTAACGCCTGGCGGGTTTCATGc            >  1:808564/1‑140 (MQ=255)
                                                                                                                            ttaattaaGGAAGCGAGTTGATTATCGTGATACTCCAGCCACCATTGGCTGTCTGGCCAACCATTTTCAGCGCCGTGGGTAATGCGGTGTCAACTTGTGCAGCGGGCGTTTGCTGGCTTAACGCCTGGCGGGTTTCATGc            >  1:490939/1‑140 (MQ=255)
                                                                                                                             taattaaGGAAGCGAGTTGATTATCGTGATACTCCAGCCACCATTGGCTGTCTGGCCAACCATTTTCAGCGCCGTGGGTAATGCGGTGTCAACTTGTGCAGCGGGCGTTTGCTGGCTTAACGCCTGGCGGGTTTCATGCa           <  2:808564/140‑1 (MQ=255)
                                                                                                                               attaaGGAAGCGAGTTGATTATCGTGATACTCCAGCCACCATTGGCTGTCTGGCCAACCATTTTCAGCGCCGTGGGTAATGCGGTGTCAACTTGTGCAGCGGGCGTTTGCTGGCTTAACGCCTGGCGGGTTTCATGCATa         <  1:799580/140‑1 (MQ=255)
                                                                                                                               attaaGGAAGCGAGTTGATTATCGTGATACTCCAGCCACCATTGGCTGTCTGGCCAACCATTTTCAGCGCCGTGGGTAATGCGGTGTCAACTTGTGCAGCGGGCGTTTGCTGGCTTAACGCCTGGCGGGTTTCATGCATa         >  1:383706/1‑140 (MQ=255)
                                                                                                                               attaaGGAAGCGAGTTGATTATCGTGATACTCCAGCCACCATTGGCTGTCTGGCCAACCATTTTCAGCGCCGTGGGTAATGCGGTGTCAACTTGTGCAGCGGGCGTTTGCTGGCTTAACGCCTGGCGGGTTTCATGCATa         >  2:1012987/1‑140 (MQ=255)
                                                                                                                                   aGGAAGCGAGTTGATTATCGTGATACTCCAGCCACCATTGGCTGTCTGGCCAACCATTTTCAGCGCCGTGGGTAATGCGGTGTCAACTTGTGCAGCGGGCGTTTGCTGGCTTAACGCCTGGCGGGTTTCATGCATAGgcg     <  2:383706/140‑1 (MQ=255)
                                                                                                                                   aGGAAGCGAGTTGATTATCGTGATACTCCAGCCACCATTGGCTGTCTGGCCAACCATTTTCAGCGCCGTGGGTAATGCGGTGTCAACTTGTGCAGCGGGCGTTTGCTGGCTTAACGCCTGGCGGGTTTCATGCATAGgcg     <  1:1012987/140‑1 (MQ=255)
                                                                                                                                    ggAAGCGAGTTGATTATCGTGATACTCCAGCCACCATTGGCTGTCTGGCCAACCATTTTCAGCGCCGTGGGTAATGCGGTGTCAACTTGTGCAGCGGGCGTTTGCTGGCTTAACGCCTGGCGGGTTTCATGCATAGgcgc    >  1:911023/1‑140 (MQ=255)
                                                                                                                                      aaGCGAGTTGATTATCGTGATACTCCAGCCACCATTGGCTGTCTGGCCAACCATTTTCAGCGCCGTGGGTAATGCGGTGTCAACTTGTGCAGCGGGCGTTTGCTGGCTTAACGCCTGGCGGGTTTCATGCa           <  1:154801/131‑1 (MQ=255)
                                                                                                                                      aaGCGAGTTGATTATCGTGATACTCCAGCCACCATTGGCTGTCTGGCCAACCATTTTCAGCGCCGTGGGTAATGCGGTGTCAACTTGTGCAGCGGGCGTTTGCTGGCTTAACGCCTGGCGGGTTTCATGCATAGGCGcac  >  2:585481/1‑140 (MQ=255)
                                                                                                                                      aaGCGAGTTGATTATCGTGATACTCCAGCCACCATTGGCTGTCTGGCCAACCATTTTCAGCGCCGTGGGTAATGCGGTGTCAACTTGTGCAGCGGGCGTTTGCTGGCTTAACGCCTGGCGGGTTTCATGCATAGGCGcac  >  2:442843/1‑140 (MQ=255)
                                                                                                                                         |                                                                                                                                        
CTCCATATCCGCCGAAAAGTCGATCTGCGGACCATCCTGAGTGGCAACTGCCTTCGCCTGGGCTTCAGCTAGCTGAATGCGCTGTTCAGCCACCTGCATATCCGGTGCGTTCTGTAGTGCATTGTTAATTAAGGAAGTGAGTTGATTATCGTGATACTCCAGCCACCATTGGCTGTCTGGCCAACCATTTTCAGCGCCGTGGGTAATGCGGTGTCAACTTGTGCAGCGGGCGTTTGCTGGCTTAACGCCTGGCGGGTTTCATGCATAGGCGCAC  >  NC_000913/2228421‑2228694

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