Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,329,296 A→T 7.1% K112* (AAA→TAA)  dacB → D‑alanyl‑D‑alanine carboxypeptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,329,2960AT 7.1% 296.2 / 10.9 84K112* (AAA→TAA) dacBD‑alanyl‑D‑alanine carboxypeptidase
Reads supporting (aligned to +/- strand):  ref base A (39/39);  new base T (3/3);  total (42/42)
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.81e-01

CAGAAAGTGATTACTGCGCTGGCGGCGTTGATTCAACTCGGCCCCGATTTTCGTTTTACCACGACGCTTGAAACCAAAGGCAATGTGGAAAACGGCGTACTTAAGGGTGACTTAGTGGCGCGATTTGGTGCCGATCCGACGTTAAAACGTCAGGATATTCGCAATATGGTCGCGACTTTGAAAAAATCTGGCGTCAACCAAATCGATGGCAATGTGTTGATAGATACCTCCATTTTCGCCAGCCACGATAAAGCCCCCGGCTGGCCATGGAATGACATGACACAA  >  NC_000913/3329152‑3329436
                                                                                                                                                |                                                                                                                                            
cAGAAAGTGATTACTGCGCTGGCGGCGTTGATTCAACTCGGCCCCGATTTTCGTTTTACCACGACGCTTGAAACCAAAGGCAATGTGGAAAACGGCGTACTTAAGGGTGACTTAGTGGCGCGATTTGGTGCCGATCCGACGTTaaa                                                                                                                                             <  3:37026/146‑1 (MQ=255)
cAGAAAGTGATTACTGCGCTGGCGGCGTTGATTCAACTCGGCCCCGATTTTCGTTTTACCACGACGCTTGAAACCAAAGGCAATGTGGAAAACGGCGTACTTAAGGGTGACTTAGTGGCGCGATTTGGTGCCGATCCGACGTTaaa                                                                                                                                             >  4:37026/1‑146 (MQ=255)
  gAAAGTGATTACTGCGCTGGCGGCGTTGATTCAACTCGGCCCCGATTTTCGTTTTACCACGACGCTTGAAACCAAAGGCAATGTGGAAAACGGCGTACTTAAGGGTGACTTAGTGGCGCGATTTGGTGCCGATCCGACGTTAAAACGt                                                                                                                                         >  4:671334/1‑148 (MQ=255)
   aaaGTGATTACTGCGCTGGCGGCGTTGATTCAACTCGGCCCCGATTTTCGTTTTACCACGACGCTTGAAACCAAAGGCAATGTGGAAAACGGCGTACTTAAGGGTGACTTAGTGGCGCGATTTGGTGCCGATCCGACGTTAAAACGTc                                                                                                                                        <  3:671334/148‑1 (MQ=255)
                        gcgTTGATTCAACTCGGCCCCGATTTTCGTTTTACCACGACGCTTGAAACCAAAGGCAATGTGGAAAACGGCGTACTTAAGGGTGACTTAGTGGCGCGATTTGGTGCCGATCCGACGTTAAAACGTCAGGATATTCGCAATATGGTcg                                                                                                                   <  1:49176/148‑1 (MQ=255)
                         cgTTGATTCAACTCGGCCCCGATTTTCGTTTTACCACGACGCTTGAAACCAAAGGCAATGTGGAAAACGGCGTACTTAAGGGTGACTTAGTGGCGCGATTTGGTGCCGATCCGACGTTAAAACGTCa                                                                                                                                       <  2:1123098/127‑1 (MQ=255)
                         cgTTGATTCAACTCGGCCCCGATTTTCGTTTTACCACGACGCTTGAAACCAAAGGCAATGTGGAAAACGGCGTACTTAAGGGTGACTTAGTGGCGCGATTTGGTGCCGATCCGACGTTAAAACGTCa                                                                                                                                       >  1:1123098/1‑127 (MQ=255)
                                                 ttCGTTTTACCACGACGCTTGAAACCAAAGGCAATGTGGAAAACGGCGTACTTAAGGGTGACTTAGTGGCGCGATTTGGTGCCGATCCGACGTTAAAACGTCAGGATATTCGCAATATGGTCGCGACTTTGAAAAAATCTGGCGTCaa                                                                                          >  1:766832/1‑148 (MQ=255)
                                                   cGTTTTACCACGACGCTTGAAACCAAAGGCAATGTGGAAAACGGCGTACTTAAGGGTGACTTAGTGGCGCGATTTGGTGCCGATCCGACGTTAAAACGt                                                                                                                                         >  1:750094/1‑99 (MQ=255)
                                                   cGTTTTACCACGACGCTTGAAACCAAAGGCAATGTGGAAAACGGCGTACTTAAGGGTGACTTAGTGGCGCGATTTGGTGCCGATCCGACGTTAAAACGt                                                                                                                                         <  2:750094/99‑1 (MQ=255)
                                                     ttttACCACGACGCTTGAAACCAAAGGCAATGTGGAAAACGGCGTACTTAAGGGTGACTTAGTGGCGCGATTTGGTGCCGATCCGACGTTAAAACGTc                                                                                                                                        <  3:98142/98‑1 (MQ=255)
                                                     ttttACCACGACGCTTGAAACCAAAGGCAATGTGGAAAACGGCGTACTTAAGGGTGACTTAGTGGCGCGATTTGGTGCCGATCCGACGTTAAAACGTc                                                                                                                                        >  4:98142/1‑98 (MQ=255)
                                                        tACCACGACGCTTGAAACCAAAGGCAATGTGGAAAACGGCGTACTTAAGGGTGACTTAGTGGCGCGATTTGGTGCCGATCCGACGTTTTAACGTcggat                                                                                                                                    >  2:715302/1‑95 (MQ=255)
                                                        tACCACGACGCTTGAAACCAAAGGCAATGTGGAAAACGGCGTACTTAAGGGTGACTTAGTGGCGCGATTTGGTGCCGATCCGACGTTTTAACGTcggat                                                                                                                                    <  1:715302/99‑5 (MQ=255)
                                                        tACCACGACGCTTGAAACCAAAGGCAATGTGGAAAACGGCGTACTTAAGGGTGACTTAGTGGCGCGATTTGGTGCCGATCCGACGTTAAAACGTCAGGATATTCGCAATATGGTCGCGACttt                                                                                                            <  3:201918/123‑1 (MQ=255)
                                                        tACCACGACGCTTGAAACCAAAGGCAATGTGGAAAACGGCGTACTTAAGGGTGACTTAGTGGCGCGATTTGGTGCCGATCCGACGTTAAAACGTCAGGATATTCGCAATATGGTCGCGACttt                                                                                                            >  4:201918/1‑123 (MQ=255)
                                                            acgacgCTTGAAACCAAAGGCAATGTGGAAAACGGCGTACTTAAGGGTGACTTAGTGGCGCGATTTGGTGCCGATCCGACGTTAAAACGTCAGGATATTCGCAATATGGTCGCGACTTTGAAAAAATCTGGCGTCAACCAAATCGATg                                                                               <  2:766832/148‑1 (MQ=255)
                                                            acgacgCCTGAAACCAAAGGCAATGTGGAAAACGGCGTACTTAAGGGTGACTTAGTGGCGCGATTTGGTGCCGATCCGACGTTAAAACGt                                                                                                                                         <  4:35538/90‑1 (MQ=255)
                                                            acgacgCCTGAAACCAAAGGCAATGTGGAAAACGGCGTACTTAAGGGTGACTTAGTGGCGCGATTTGGTGCCGATCCGACGTTAAAACGt                                                                                                                                         >  3:35538/1‑90 (MQ=255)
                                                               acgCTTGAAACCAAAGGCAATGTGGAAAACGGCGTACTTAAGGGTGACTTAGTGGCGCGATTTGGTGCCGATCCGACGTTAAAACGTCAGGATATTCGCAATATGGTCGCGACTTTGAAAAAATCTGGCGTCAACCAAATCGATGGCa                                                                            >  2:423553/1‑148 (MQ=255)
                                                                 gCTTGAAACCAAAGGCAATGTGGAAAACGGCGTACTTAAGGGTGACTTAGTGGCGCGATTTGGTGCCGATCCGACGTTaaa                                                                                                                                             >  1:1105138/1‑81 (MQ=255)
                                                                 gCTTGAAACCAAAGGCAATGTGGAAAACGGCGTACTTAAGGGTGACTTAGTGGCGCGATTTGGTGCCGATCCGACGTTaaa                                                                                                                                             <  2:1105138/81‑1 (MQ=255)
                                                                    tGAAACCAAAGGCAATGTGGAAAACGGCGTACTTAAGGGTGACTTAGTGGCGCGATTTGGTGCCGATCCGACGTTAAAACg                                                                                                                                          >  1:687988/1‑81 (MQ=255)
                                                                    tGAAACCAAAGGCAATGTGGAAAACGGCGTACTTAAGGGTGACTTAGTGGCGCGATTTGGTGCCGATCCGACGTTAAAACg                                                                                                                                          <  2:687988/81‑1 (MQ=255)
                                                                    tGAAACCAAAGGCAATGTGGAAAACGGCGTACTTAAGGGTGACTTAGTGGCGCGATTTGGTGCCGATCCGACGTTAAAACGTc                                                                                                                                        <  4:150299/83‑1 (MQ=255)
                                                                    tGAAACCAAAGGCAATGTGGAAAACGGCGTACTTAAGGGTGACTTAGTGGCGCGATTTGGTGCCGATCCGACGTTAAAACGTc                                                                                                                                        >  3:150299/1‑83 (MQ=255)
                                                                    tGAAACCAAAGGCAATGTGGAAAACGGCGTACTTAAGGGTGACTTAGTGGCGCGATTTGGTGCCGATCCGACGTTAAAACGTCAGGATAtt                                                                                                                                <  3:805152/91‑1 (MQ=255)
                                                                    tGAAACCAAAGGCAATGTGGAAAACGGCGTACTTAAGGGTGACTTAGTGGCGCGATTTGGTGCCGATCCGACGTTAAAACGTCAGGATAtt                                                                                                                                >  4:805152/1‑91 (MQ=255)
                                                                    tGAAACCAAAGGCAATGTGGAAAACGGCGTACTTAAGGGTGACTTAGTGGCGCGATTTGGTGCCGATCCGACGTTAAAACGTCAGGATATTCGCAATATGGTCGCGACTTTGAAAAAATCTGGCGTCAACCAAATCGATGGCAAtgtg                                                                       <  1:423553/148‑1 (MQ=255)
                                                                      aaaCCAAAGGCAATGTGGAAAACGGCGTACTTAAGGGTGACTTAGTGGCGCGATTTGGTGCCGATCCGACGTTAAAACGTCAGGATATTCGCAATATGGTCGCGACTTTGAAAAAATCTGGCGTCAACCAAATCGa                                                                                 >  3:686782/1‑136 (MQ=255)
                                                                      aaaCCAAAGGCAATGTGGAAAACGGCGTACTTAAGGGTGACTTAGTGGCGCGATTTGGTGCCGATCCGACGTTAAAACGTCAGGATATTCGCAATATGGTCGCGACTTTGAAAAAATCTGGCGTCAACCAAATCGa                                                                                 <  4:686782/136‑1 (MQ=255)
                                                                      aaaCCAAAGGCAATGTGGAAAACGGCGTACTTAAGGGTGACTTAGTGGCGCGATTTGGTGCCGATCCGACGTTAAAACGTCAGGATATTCGCAATATGGTCGCGACTTTGAAAAAATCTGGCGTCAACCAAATCGATGGCAATGTGtt                                                                     >  2:1003382/1‑148 (MQ=255)
                                                                       aaCCAAAGGCAATGTGGAAAACGGCGTACTTAAGGGTGACTTAGTGGCGCGATTTGGTGCCGATCCGACGTTAAAACGTCAGGATATTCGCAATATGGTCGCGACt                                                                                                              >  2:1032793/1‑106 (MQ=255)
                                                                       aaCCAAAGGCAATGTGGAAAACGGCGTACTTAAGGGTGACTTAGTGGCGCGATTTGGTGCCGATCCGACGTTAAAACGTCAGGATATTCGCAATATGGTCGCGACt                                                                                                              <  1:1032793/106‑1 (MQ=255)
                                                                         ccAAAGGCAATGTGGAAAACGGCGTACTTAAGGGTGACTTAGTGGCGCGATTTGGTGCCGATCCGACGTTAAAACGt                                                                                                                                         >  4:329036/1‑77 (MQ=255)
                                                                         ccAAAGGCAATGTGGAAAACGGCGTACTTAAGGGTGACTTAGTGGCGCGATTTGGTGCCGATCCGACGTTAAAACGt                                                                                                                                         <  3:329036/77‑1 (MQ=255)
                                                                          cAAAGGCAATGTGGAAAACGGCGTACTTAAGGGTGACTTAGTGGCGCGATTTGGTGCCGATCCGACGTTAAAACGTCAGGATATTCGCAATATGGTCGCGACTTTGAAAAAATCTGGCGTCa                                                                                           >  3:460188/1‑122 (MQ=255)
                                                                          cAAAGGCAATGTGGAAAACGGCGTACTTAAGGGTGACTTAGTGGCGCGATTTGGTGCCGATCCGACGTTAAAACGTCAGGATATTCGCAATATGGTCGCGACTTTGAAAAAATCTGGCGTCa                                                                                           <  4:460188/122‑1 (MQ=255)
                                                                           aaaGGCAATGTGGAAAACGGCGTACTTAAGGGTGACTTAGTGGCGCGATTTGGTGCCGATCCGACGTTAAAACGTCAGGATATTCGCAATATGGTCGCGACTTTGAAAAAATCTGGCGTCAACCAAATCGATGGCAAt                                                                          <  3:368700/138‑1 (MQ=255)
                                                                           aaaGGCAATGTGGAAAACGGCGTACTTAAGGGTGACTTAGTGGCGCGATTTGGTGCCGATCCGACGTTAAAACGTCAGGATATTCGCAATATGGTCGCGACTTTGAAAAAATCTGGCGTCAACCAAATCGATGGCAAt                                                                          >  4:368700/1‑138 (MQ=255)
                                                                           aaaGGCAATGTGGAAAACGGCGTACTTAAGGGTGACTTAGTGGCGCGATTTGGTGCCGATCCGACGTTAAAACGTCAGGATATTCGCAATATGGTCGCGACTTTGAAAAAATCTGGCGTCAACCAAATCGATGGCAATGTGTTgatag                                                                <  1:1003382/148‑1 (MQ=255)
                                                                            aaGGCAATGTGGAAAACGGCGTACTTAAGGGTGACTTAGTGGCGCGATTTGGTGCCGATCCGACGTTaaa                                                                                                                                             >  1:1005804/1‑70 (MQ=255)
                                                                            aaGGCAATGTGGAAAACGGCGTACTTAAGGGTGACTTAGTGGCGCGATTTGGTGCCGATCCGACGTTaaa                                                                                                                                             <  2:1005804/70‑1 (MQ=255)
                                                                              ggCAATGTGGAAAACGGCGTACTTAAGGGTGACTTAGTGGCGCGATTTGGTGCCGATCCGACGTTAAAACGTCAGGATATTCGCAatat                                                                                                                        <  1:412068/89‑1 (MQ=255)
                                                                              ggCAATGTGGAAAACGGCGTACTTAAGGGTGACTTAGTGGCGCGATTTGGTGCCGATCCGACGTTAAAACGTCAGGATATTCGCAatat                                                                                                                        >  2:412068/1‑89 (MQ=255)
                                                                               gCAATGTGGAAAACGGCGTACTTAAGGGTGACTTAGTGGCGCGATTTGGTGCCGATCCGACGTTAAAACGTCAGGATATTCGCAATATGGTCGCGACTTTGAAAAAATCTGGCGt                                                                                             <  1:377115/115‑1 (MQ=255)
                                                                               gCAATGTGGAAAACGGCGTACTTAAGGGTGACTTAGTGGCGCGATTTGGTGCCGATCCGACGTTAAAACGTCAGGATATTCGCAATATGGTCGCGACTTTGAAAAAATCTGGCGt                                                                                             >  2:377115/1‑115 (MQ=255)
                                                                                  aTGTGGAAAACGGCGTACTTAAGGGTGACTTAGTGGCGCGATTTGGTGCCGATCCGACGTTAAAACGTCAGGATATTCGCAATATGGTcgcg                                                                                                                 >  4:665954/1‑92 (MQ=255)
                                                                                  aTGTGGAAAACGGCGTACTTAAGGGTGACTTAGTGGCGCGATTTGGTGCCGATCCGACGTTAAAACGTCAGGATATTCGCAATATGGTcgcg                                                                                                                 <  3:665954/92‑1 (MQ=255)
                                                                                        aaaaCGGCGTACTTAAGGGTGACTTAGTGGCGCGATTTGGTGCCGATCCGACGTTAAAACGTCAGGATATTCGCAATATGGTCGCGACTTTGAAAAAATCTGGCGt                                                                                             <  2:82998/106‑1 (MQ=255)
                                                                                        aaaaCGGCGTACTTAAGGGTGACTTAGTGGCGCGATTTGGTGCCGATCCGACGTTAAAACGTCAGGATATTCGCAATATGGTCGCGACTTTGAAAAAATCTGGCGt                                                                                             >  1:82998/1‑106 (MQ=255)
                                                                                         aaaCGGCGTACTTAAGGGTGACTTAGTGGCGCGATTTGGTGCCGATCCGACGTTAAAACGTCAGGATATTCGCAATATGGTCGCGACTTTGAAAAAATCTGGCGTCAACCAAATCGATGGCAATGTGTTgat                                                                  >  2:661037/1‑132 (MQ=255)
                                                                                         aaaCGGCGTACTTAAGGGTGACTTAGTGGCGCGATTTGGTGCCGATCCGACGTTAAAACGTCAGGATATTCGCAATATGGTCGCGACTTTGAAAAAATCTGGCGTCAACCAAATCGATGGCAATGTGTTgat                                                                  <  1:661037/132‑1 (MQ=255)
                                                                                                    ttAAGGGTGACTTAGTGGCGCGATTTGGTGCCGATCCGACGTTAAAACGTCAGGATATTCGCAATATGGTCGCGACTTTGAAAAAATCTGGCGTCAACCAAATCGATGGCAATGTGTTgat                                                                  >  4:58549/1‑121 (MQ=255)
                                                                                                    ttAAGGGTGACTTAGTGGCGCGATTTGGTGCCGATCCGACGTTAAAACGTCAGGATATTCGCAATATGGTCGCGACTTTGAAAAAATCTGGCGTCAACCAAATCGATGGCAATGTGTTgat                                                                  <  3:58549/121‑1 (MQ=255)
                                                                                                     tAAGGGTGACTTAGTGGCGCGATTTGGTGCCGATCCGACGTTAAAACGTCAGGATATTCGCAATATGGTCGCGACTTTGAAAAAATCTGGCGTCAACCAAATCGATGGCAATGTGTTgat                                                                  <  3:475178/120‑1 (MQ=255)
                                                                                                     tAAGGGTGACTTAGTGGCGCGATTTGGTGCCGATCCGACGTTAAAACGTCAGGATATTCGCAATATGGTCGCGACTTTGAAAAAATCTGGCGTCAACCAAATCGATGGCAATGTGTTgat                                                                  >  4:475178/1‑120 (MQ=255)
                                                                                                         ggTGACTTAGTGGCGCGATTTGGTGCCGATCCGACGTTAAAACGTCAGGATATTCGCAATATGGTCGCGACTTTGAAAAAATCTGGCGTCAACCAAATCGATGGCAATGTGTTGATAGATACCTCCAtttt                                                   >  2:911494/1‑131 (MQ=255)
                                                                                                         ggTGACTTAGTGGCGCGATTTGGTGCCGATCCGACGTTAAAACGTCAGGATATTCGCAATATGGTCGCGACTTTGAAAAAATCTGGCGTCAACCAAATCGATGGCAATGTGTTGATAGATACCTCCAtttt                                                   <  1:911494/131‑1 (MQ=255)
                                                                                                          gTGACTTAGTGGCGCGATTTGGTGCCGATCCGACGTTAAAACGTCAGGATATTCGCAATATGGt                                                                                                                     <  2:1106311/64‑1 (MQ=255)
                                                                                                          gTGACTTAGTGGCGCGATTTGGTGCCGATCCGACGTTAAAACGTCAGGATATTCGCAATATGGt                                                                                                                     >  1:1106311/1‑64 (MQ=255)
                                                                                                          gTGACTTAGTGGCGCGATTTGGTGCCGATCCGACGTTAAAACGTCAGGATATTCGCAATATGGTCGCGACTTTGAAAAAATCTGGCGTCAACCAAATCGATGGCAATGTGTTGATAGATACCTCCATTTTCGCCAGCCACGATaa                                    <  1:738388/145‑1 (MQ=255)
                                                                                                          gTGACTTAGTGGCGCGATTTGGTGCCGATCCGACGTTAAAACGTCAGGATATTCGCAATATGGTCGCGACTTTGAAAAAATCTGGCGTCAACCAAATCGATGGCAATGTGTTGATAGATACCTCCATTTTCGCCAGCCACGATaa                                    >  2:738388/1‑145 (MQ=255)
                                                                                                                tAGTGGCGCGATTTGGTGCCGATCCGACGTTAAAACGTCAGGATATTCGCAATATGGTCGCGACTTTGAAAAAATCTGGCGTCAACCAAATCGATGGCAATGTGTTga                                                                   <  3:440213/108‑1 (MQ=255)
                                                                                                                tAGTGGCGCGATTTGGTGCCGATCCGACGTTAAAACGTCAGGATATTCGCAATATGGTCGCGACTTTGAAAAAATCTGGCGTCAACCAAATCGATGGCAATGTGTTga                                                                   >  4:440213/1‑108 (MQ=255)
                                                                                                                   tGGCGCGATTTGGTGCCGATCCGACGTTAAAACGTCAGGATATTCGCAATATGGTCGCGACTTTGAAAAAATCTgg                                                                                                >  1:504403/1‑76 (MQ=255)
                                                                                                                   tGGCGCGATTTGGTGCCGATCCGACGTTAAAACGTCAGGATATTCGCAATATGGTCGCGACTTTGAAAAAATCTgg                                                                                                <  2:504403/76‑1 (MQ=255)
                                                                                                                      cgcgATTTGGTGCCGATCCGACGTTAAAACGTCAGGATATTCGCAATATGGTCGCGACTTTGAAAAAATCTGGCGTCAACCAAATCGATGGCAATGTGTTga                                                                   >  1:214993/1‑102 (MQ=255)
                                                                                                                      cgcgATTTGGTGCCGATCCGACGTTAAAACGTCAGGATATTCGCAATATGGTCGCGACTTTGAAAAAATCTGGCGTCAACCAAATCGATGGCAATGTGTTga                                                                   <  2:214993/102‑1 (MQ=255)
                                                                                                                      cgcgATTTGGTGCCGATCCGACGTTAAAACGTCAGGATATTCGCAATATGGTCGCGACTTTGAAAAAATCTGGCGTCAACCAAATCGATGGCAATGTGTTGATAGATACCTCCATTTTCGCCAGCCACGATAAAGCCCCCGGCTGGcc                     >  2:482108/1‑148 (MQ=255)
                                                                                                                        cgATTTGGTGCCGATCCGACGTTAAAACGTCAGGATATTCGCAATATGGTCGCGACttt                                                                                                            <  3:877573/59‑1 (MQ=255)
                                                                                                                        cgATTTGGTGCCGATCCGACGTTAAAACGTCAGGATATTCGCAATATGGTCGCGACttt                                                                                                            >  4:877573/1‑59 (MQ=255)
                                                                                                                        cgATTTGGTGCCGATCCGACGTTAAAACGTCAGGATATTCGCAATATGGTCGCGACTTTGAAAAAATCTGGCGTCAACCaa                                                                                      >  4:387846/1‑81 (MQ=255)
                                                                                                                        cgATTTGGTGCCGATCCGACGTTAAAACGTCAGGATATTCGCAATATGGTCGCGACTTTGAAAAAATCTGGCGTCAACCaa                                                                                      <  3:387846/81‑1 (MQ=255)
                                                                                                                          aTTTGGTGCCGATCCGACGTTAAAACGTGAGGATATTCGCAATATGGTCGCGACTTTGAAAAAATCTGGCGTCAACCAAATCGATGGCAATGTGTTGATAGATACCTCCATTTTCGCGAGCCACGATAAAGCCCCCGGCTGGCCATgg                 >  3:931579/1‑148 (MQ=255)
                                                                                                                                   cGATCCGACGTTAAAACGTCAGGATATTCGCAATATGGTCGCGACTTTGAAAAAATCTGGCGTCAACCAAATCGATGGCAATGTGTTGATAGATACCTCCATTTTCGCCAGCCACGa                                       >  1:614584/1‑117 (MQ=255)
                                                                                                                                   cGATCCGACGTTAAAACGTCAGGATATTCGCAATATGGTCGCGACTTTGAAAAAATCTGGCGTCAACCAAATCGATGGCAATGTGTTGATAGATACCTCCATTTTCGCCAGCCACGa                                       <  2:614584/117‑1 (MQ=255)
                                                                                                                                        tgACGTTTTAACGTCAGGATATTCGCAATATGGTCGCGACTTTGAAAAAATCTGGCGTCAACCAAATCGATGGCAATGTGTTGATAGATACCTCCATTTTCGCCAGCCACGATAAAGCCCCCGGCTGGTCATGGAatga            >  4:775619/2‑139 (MQ=255)
                                                                                                                                        tgACGTTTTAACGTCAGGATATTCGCAATATGGTCGCGACTTTGAAAAAATCTGGCGTCAACCAAATCGATGGCAATGTGTTGATAGATACCTCCATTTTCGCCAGCCACGATAAAGCCCCCGGCTGGCCATGGAatga            >  4:138622/2‑139 (MQ=255)
                                                                                                                                        tgACGTTTTAACGTCAGGATATTCGCAATATGGTCGCGACTTTGAAAAAATCTGGCGTCAACCAAATCGATGGCAATGTGTTGATAGATACCTCCATTTTCGCCAGCCACGATAAAGCCCCCGGCTGGCCATGGAatga            <  3:775619/138‑1 (MQ=255)
                                                                                                                                        tgACGTTTTAACGTCAGGATATTCGCAATATGGTCGCGACTTTGAAAAAATCTGGCGTCAACCAAATCGATGGCAATGTGTTGATAGATACCTCCATTTTCGCCAGCCACGATAAAGCCCCCGGCTGGCCATGGAatga            <  3:138622/138‑1 (MQ=255)
                                                                                                                                        cGACGTTAAAACGTCAGGATATTCGCAATATGGTCGCGACTTTGAAAAAATCTGGCGTCAACCAAATCGATGGCAATGTGTTGATAGATACCTCCAtttt                                                   >  2:1016938/1‑100 (MQ=255)
                                                                                                                                        cGACGTTAAAACGTCAGGATATTCGCAATATGGTCGCGACTTTGAAAAAATCTGGCGTCAACCAAATCGATGGCAATGTGTTGATAGATACCTCCAtttt                                                   <  1:1016938/100‑1 (MQ=255)
                                                                                                                                         gACGTTAACACGTCAGCATATTCGCAATATGGTCGCGACTTTGAAACAATCTGGCTTCAACCAATTCGATGGCAATGTGTTGATAGATTCCTCCATTTTCGTCAGCCACGATACAGCCCGCGGCTGGCCATGGAATGTCATGACACaa  <  4:931579/148‑1 (MQ=255)
                                                                                                                                          aCGTTAAAACGTCAGGATATTCGCAATATGGTCGCGACTTTGAAAACATCTGGCGTCAACCAAATCGATGGCAATGTGTTGATAGATACCt                                                          <  4:178528/91‑1 (MQ=255)
                                                                                                                                          aCGTTAAAACGTCAGGATATTCGCAATATGGTCGCGACTTTGAAAAAATCTGGCGTCAACCAAATCGATGGCAATGTGTTga                                                                   >  4:1074648/1‑82 (MQ=255)
                                                                                                                                          aCGTTAAAACGTCAGGATATTCGCAATATGGTCGCGACTTTGAAAAAATCTGGCGTCAACCAAATCGATGGCAATGTGTTga                                                                   <  3:1074648/82‑1 (MQ=255)
                                                                                                                                          aCGTTAAAACGTCAGGATATTCGCAATATGGTCGCGACTTTGAAAAAATCTGGCGTCAACCAAATCGATGGCAATGTGTTGATAGATACCt                                                          >  3:178528/1‑91 (MQ=255)
                                                                                                                                          aCGTTAAAACGTCAGGATATTCGCAATATGGTCGCGACTTTGAAAAAATCTGGCGTCAACCAAATCGATGGCAATGTGTTGATAGATACCTCCATTTTCGCCAGCCACGATAAAGCCCCCGGCTGGCCATGGAatgac           >  1:403929/1‑138 (MQ=255)
                                                                                                                                          aCGTTAAAACGTCAGGATATTCGCAATATGGTCGCGACTTTGAAAAAATCTGGCGTCAACCAAATCGATGGCAATGTGTTGATAGATACCTCCATTTTCGCCAGCCACGATAAAGCCCCCGGCTGGCCATGGAatgac           <  2:403929/138‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                            
CAGAAAGTGATTACTGCGCTGGCGGCGTTGATTCAACTCGGCCCCGATTTTCGTTTTACCACGACGCTTGAAACCAAAGGCAATGTGGAAAACGGCGTACTTAAGGGTGACTTAGTGGCGCGATTTGGTGCCGATCCGACGTTAAAACGTCAGGATATTCGCAATATGGTCGCGACTTTGAAAAAATCTGGCGTCAACCAAATCGATGGCAATGTGTTGATAGATACCTCCATTTTCGCCAGCCACGATAAAGCCCCCGGCTGGCCATGGAATGACATGACACAA  >  NC_000913/3329152‑3329436

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