Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,610,720 T→A 5.1% L436Q (CTG→CAG)  prfC → peptide chain release factor RF‑3

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,610,7200TA 5.1% 728.7 / 20.9 195L436Q (CTG→CAG) prfCpeptide chain release factor RF‑3
Reads supporting (aligned to +/- strand):  ref base T (93/92);  new base A (5/5);  total (98/97)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.36e-01

Only 100 of 196 total aligned reads displayed.

CACCGGTATTCCGAACTTCGCACCAGAACTGTTCCGTCGTATCCGCCTGAAAGATCCGCTGAAGCAAAAACAGCTGCTCAAAGGGCTGGTACAGCTTTCCGAAGAGGGCGCGGTGCAGGTGTTCCGTCCAATCTCCAACAACGATCTGATCGTTGGTGCAGTTGGTGTGCTGCAGTTTGATGTGGTGGTAGCGCGCCTGAAGAGCGAATACAACGTTGAAGCAGTGTATGAGTCAGTCAACGTTGCCACTGCCCGCTGGGTAGAATGTGCAGACGCGAAGAAATTCGAAGAGTT  >  NC_000913/4610574‑4610867
                                                                                                                                                  |                                                                                                                                                   
     gTATTCCGAACTTCGCACCAGAACTGTTCCGTCGTATCCGCCTGAAAGATCCGCTGAAGCAAAAACAGCTGCTCAAAGGGCTGGTACAGCTTTCCGAAGAGGGCGCGGTGCAGGTGTTCCGTCCAATCTCCAACAACGATCTGATCGt                                                                                                                                               <  1:600157/148‑1 (MQ=255)
      tATTCCGAACTTCGCACCAGAACTGTTCCGTCGTATCCGCCTGAAAGATCCGCTGAAGCAAAAACAGCTGCTCAAAGGGCTGGTACAGCTTTCCGAAGAGGGCGCGGTGCAGGTGTTCCGTCCAATCTCCAACAACGATCTGATCGtt                                                                                                                                              <  2:996937/148‑1 (MQ=255)
      tATTCCGAACTTCGCACCAGAACTGTTCCGTCGTATCCGCCTGAAAGATCCGCTGAAGCAAAAACAGCTGCTCAAAGGGCTGGTACAGCTTTCCGAAGAGGGCGCGGTGCAGGTGTTCCGTCCAATCTCCAACAACGATCTGATCGtt                                                                                                                                              <  2:230474/148‑1 (MQ=255)
           cGAACTTCGCACCAGAACTGTTCCGTCGTATCCGCCTGAAAGATCCGCTGAAGCAAAAACAGCTGCTCAAAGGGCTGGTACAGCTCTCCGAAGAGGGCGCGGTGCAGGTGTTCCGTCCAATCTCCAACAACGATCTGATCggtggttc                                                                                                                                         >  2:60868/1‑146 (MQ=255)
                      ccAGAACTGTTCCGTCGTATCCGCCTGAAAGATCCGCTGAAGCAAAAACAGCTGCTCAAAGGGCTGGTACAGCTTTCCGAAGAGGGCGCGGTGCAGGTGTTCCGTCCAATCTCCAACAACGATCTGATCGtt                                                                                                                                              >  3:622112/1‑132 (MQ=255)
                      ccAGAACTGTTCCGTCGTATCCGCCTGAAAGATCCGCTGAAGCAAAAACAGCTGCTCAAAGGGCTGGTACAGCTTTCCGAAGAGGGCGCGGTGCAGGTGTTCCGTCCAATCTCCAACAACGATCTGATCGTTGGTGCAGTTGGTGtgc                                                                                                                              >  2:47136/1‑148 (MQ=255)
                        aGAACTGTTCCGTCGTATCCGCCTGAAAGATCCGCTGAAGCAAAAACAGCTGCTCAAAGGGCTGGTACAGCTTTCCGAAGAGGGCGCGGTGCAGGTGTTCCGTCCAATCTCCAACAACGATCTGATCGt                                                                                                                                               >  3:1012582/1‑129 (MQ=255)
                        aGAACTGTTCCGTCGTATCCGCCTGAAAGATCCGCTGAAGCAAAAACAGCTGCTCAAAGGGCTGGTACAGCTTTCCGAAGAGGGCGCGGTGCAGGTGTTCCGTCCAATCTCCAACAACGATCTGATCGt                                                                                                                                               <  4:1012582/129‑1 (MQ=255)
                                tCCGTCGTATCCGCCTGAAAGATCCGCTGAAGCAAAAACAGCTGCTCAAAGGGCTGGTACAGCTTTCCGAAGAGGGCGCGGTGCAGGTGTTCCGTCCAATCTCCAACAACGATCTGATCGTTGGTGCAGTTGGTGTGCTGCAGTTTGa                                                                                                                    >  3:661356/1‑148 (MQ=255)
                                tCCGTCGTATCCGCCTGAAAGATCCGCTGAAGCAAAAACAGCTGCTCAAAGGGCTGGTACAGCTTTCCGAAGAGGGCGCGGTGCAGGTGTTCCGTCCAATCTCCAACAACGATCTGATCGTTGGTGCAGTTGGTGTGCTGCAGTTTGa                                                                                                                    >  1:683753/1‑148 (MQ=255)
                                    tcgtATCCGCCTGAAAGATCCGCTGAAGCAAAAACAGCTGCTCAAAGGGCTGGTACAGCTTTCCGAAGAGGGCGCGGTGCAGGTGTTCCGTCCAATCTCCAACAACGATCTGATCGTTGGTGCAGTTGGTGTGCTGCAGTTTGAtgtg                                                                                                                <  2:683753/148‑1 (MQ=255)
                                    tcgtATCCGCCTGAAAGATCCGCTGAAGCAAAAACAGCTGCTCAAAGGGCTGGTACAGCTTTCCGAAGAGGGCGCGGTGCAGGTGTTCCGTCCAATCTCCAACAACGATCTGATCGTTGGTGCAGTTGGTGTGCTGCAGTTTGAtgtg                                                                                                                <  4:661356/148‑1 (MQ=255)
                                                gAAAGATCCGCTGAAGCAAAAACAGCTGCTCAAAGGGCTGGTACAGCTTTCCGAAGAGGGCGCGGTGCAGGTGTTCCGTCCAATCTCCAACAACGATCTGATCGTTGGTGCAGTTGGTGTGCTGCAGTTTGATGTGGTGGTAgcgcgc                                                                                                    >  3:406404/1‑148 (MQ=255)
                                                  aaGATCCGCTGAAGCAAAAACAGCTGCTCAAAGGGCTGGTACAGCTTTCCGAAGAGGGCGCGGTGCAGGTGTTCCGTCCAATCTCCAACAACGATCTGATCGTTGGTGCAGTTGGTGTGCTGCAGTTTGATGTGGTGGTAgcg                                                                                                       >  4:561103/1‑143 (MQ=255)
                                                    gATCCGCTGAAGCAAAAACAGCTGCTCAAAGGGCTGGTACAGCTTTCCGAAGAGGGCGCGGTGCAGGTGTTCCGTCCAATCTCCAACAACGATCTGATCGTTGGTGCAGTTgg                                                                                                                                   >  2:62247/1‑113 (MQ=255)
                                                    gATCCGCTGAAGCAAAAACAGCTGCTCAAAGGGCTGGTACAGCTTTCCGAAGAGGGCGCGGTGCAGGTGTTCCGTCCAATCTCCAACAACGATCTGATCGTTGGTGCAGTTgg                                                                                                                                   <  1:62247/113‑1 (MQ=255)
                                                      tCCGCTGAAGCAAAAACAGCTGCTCAAAGGGCTGGTACAGCTTTCCGAAGAGGGCGCGGTGCAGGTGTTCCGTCCAATCTCCAACAACGATCTGATCGTTGGTGCAGt                                                                                                                                      <  1:531221/108‑1 (MQ=255)
                                                      tCCGCTGAAGCAAAAACAGCTGCTCAAAGGGCTGGTACAGCTTTCCGAAGAGGGCGCGGTGCAGGTGTTCCGTCCAATCTCCAACAACGATCTGATCGTTGGTGCAGt                                                                                                                                      >  2:531221/1‑108 (MQ=255)
                                                      tCCGCTGAAGCAAAAACAGCTGCTCAAAGGGCTGGTACAGCTTTCCGAAGAGGGCGCGGTGCAGGTGTTCCGTCCAATCTCCAACAACGATCTGATCGTTGGTGCAGTTGGTGTGCTGCAGTTTGATGTGGTGGTAGCGCGCCTGAag                                                                                              <  2:76454/148‑1 (MQ=255)
                                                        cGCTGAAGCAAAAACAGCTGCTCAAAGGGCTGGTACAGCTTTCCGAAGAGGGCGCGGTGCAGGTGTTCCGTCCAATCTCCAACAACGATCTGATCGTTGGTGCAGTTGGTGTGCTGCAGTTTGATGTGGTGGTAGCGCGCCTGAagag                                                                                            <  4:406404/148‑1 (MQ=255)
                                                          cTGAAGCAAAAACAGCTGCTCAAAGGGCTGGTACAGCTTTCCGAAGAGGGCGCGGTGCAGGTGTTCCGTCCAATCTCCAACAACGATCTGATCGTTGGTGCAGTTGGTGTGCTGCAGTTTGATGTGGTGGTAGCGCGCCTGAAGAGCg                                                                                          >  3:774880/1‑148 (MQ=255)
                                                                cAAAAACAGCTGCTCAAAGGGCTGGTACAGCTTTCCGAAGAGGGCGCGGTGCAGGTGTTCCGTCCAATCTCCAACAACGATCTGATCGTTGGTGCAGTTGGTGTGCTGCAGTTTGATGTGGTGGTAGCGCGCCTGAAGAGCg                                                                                          <  2:1082251/142‑1 (MQ=255)
                                                                cAAAAACAGCTGCTCAAAGGGCTGGTACAGCTTTCCGAAGAGGGCGCGGTGCAGGTGTTCCGTCCAATCTCCAACAACGATCAGATCGTTg                                                                                                                                             <  1:303804/91‑1 (MQ=255)
                                                                cAAAAACAGCTGCTCAAAGGGCTGGTACAGCTTTCCGAAGAGGGCGCGGTGCAGGTGTTCCGTCCAATCTCCAACAACGATCAGATCGTTg                                                                                                                                             >  2:303804/1‑91 (MQ=255)
                                                                cAAAAACAGCTGCTCAAAGGGCTGGTACAGCTTTCCGAAGAGGGCGCGGTGCAGGTGTTCCGTCCAATCTCCAACAACGATCAGATCGTTg                                                                                                                                             <  3:179879/91‑1 (MQ=255)
                                                                  aaaaCAGCTGCTCAAAGGGCTGGTACAGCTTTCCGAAGAGGGCGCGGTGCAGGTGTTCCGTCCAATCTCCAACAACGATCTGATCGTTGGTGCAGTTGGTGTGCTGCAGTTTGATGTGGTGGTAGCGCGCCTGAAGAGCGAATACa                                                                                    <  2:452609/146‑1 (MQ=255)
                                                                  aaaaCAGCTGCTCAAAGGGCTGGTACAGCTTTCCGAAGAGGGCGCGGTGCAGGTGTTCCGTCCAATCTCCAACAACGATCTGATCGTTGGTGCAGTTGGTGTGCTGCAGTTTGATGTGGTGGTAGCGCGCCTGAAGAGCGAATACa                                                                                    >  1:452609/1‑146 (MQ=255)
                                                                   aaaCAGCTGCTCAAAGGGCTGGTACAGCTTTCCGAAGAGGGCGCGGTGCAGGTGTTCCGTCCAATCTCCAACAACGATCTGATCGTTGGTGCAGTTGGTGTGCTGCAGTTTGATGTGGTGGTAGCGCGCCTGAAGAGCGAATACAACg                                                                                 >  1:1057876/1‑148 (MQ=255)
                                                                    aaCAGCTGCTCAACGGGCTGGTACAGCTTTCCGAAGAGGGCGCGGTGCAGGTGTTCCGTCCAATCTCCAACAACGATCTGATCGTTGGTGCAGTTGGTGTGCTGCAGTTTGATGTggtggt                                                                                                           <  4:816750/121‑1 (MQ=255)
                                                                    aaCAGCTGCTCAAAGGGCTGGTACAGCTTTCCGAAGAGGGCGCGGTGCAGGTGTTCCGTCCAATCTCCAACAACGATCTGATCGTTgg                                                                                                                                            <  3:275491/88‑1 (MQ=255)
                                                                    aaCAGCTGCTCAAAGGGCTGGTACAGCTTTCCGAAGAGGGCGCGGTGCAGGTGTTCCGTCCAATCTCCAACAACGATCTGATCGTTgg                                                                                                                                            >  4:275491/1‑88 (MQ=255)
                                                                    aaCAGCTGCTCAAAGGGCTGGTACAGCTTTCCGAAGAGGGCGCGGTGCAGGTGTTCCGTCCAATCTCCAACAACGATCTGATCGTTGGTGCAGTTGGTGTGCTGCAGTTTGATGTGGTGGTAGCGCGCCTGAAGAGCGAATACAACGt                                                                                >  1:1055483/1‑148 (MQ=255)
                                                                     aCAGCTGCTCAAAGGGCTGGTACAGCTTTCCGAAGAGGGCGCGGTGCAGGTGTTCCGTCCAATCTCCAACAACGATCTGATCGTTGGTGCAGTTGGTGTGCTGCAGTTTGATGTGGTGGTAGCGCGCCTGAAGAGCGAATACa                                                                                    >  1:1087744/1‑143 (MQ=255)
                                                                       aGCTGCTCAAAGGGCTGGTACAGCTTTCCGAAGAGGGCGCGGTGCAGGTGTTCCGTCCAATCTCCAACAACGATCTGATCGTTGGTGCAGTTGGTGTGCTGCAGTTTGATGTGGTGGTAGCGCGCCTGAAGAGCGAATACAACGTTGa                                                                             <  2:1057876/148‑1 (MQ=255)
                                                                       aGCTGCTCAAAGGGCTGGTACAGCTTTCCGAAGAGGGCGCGGTGCAGGTGTTCCGTCCAATCTCCAACAACGATCTGATCGTTGGTGCAGTTGGTGTGCTGCAGTTTGATGTGGTGGTAGCGCGCCTGAAGAGCGAATACAACGTTGa                                                                             <  4:649844/148‑1 (MQ=255)
                                                                         ctgctCAAAGGGCTGGTACAGCTTTCCGAAGAGGGCGCGGTGCAGGTGTTCCGTCCAATCTCCAACAACGATCTGATCGCTGGTGCAGTTGGTGTGCTGCAGTTTGATGTGGTGGTAGCGCGCCTGAAGAGCGAATACAACGTTGAAg                                                                           >  4:126875/1‑148 (MQ=255)
                                                                          tgctCAAAGGGCTGGTACAGCTTTCCGAAGAGGGCGCGGTGCAGGTGTTCCGTCCAATCTCCAACAACGATCTGATCGTTGGTGCAGTTGGTGTGCTGCAGTTTGATGTGGTGGTAGCGCGCCTGAAGAGCGAATACAACGTTGAAGc                                                                          <  1:47136/148‑1 (MQ=255)
                                                                              cAAAGGGCTGGTACAGCTTTCCGAAGAGGGCGCGGTGCAGGTGTTCCGTCCAATCTCCAACAACGATCTg                                                                                                                                                    >  1:614790/1‑70 (MQ=255)
                                                                              cAAAGGGCTGGTACAGCTTTCCGAAGAGGGCGCGGTGCAGGTGTTCCGTCCAATCTCCAACAACGATCTg                                                                                                                                                    >  1:1081648/1‑70 (MQ=255)
                                                                              cAAAGGGCTGGTACAGCTTTCCGAAGAGGGCGCGGTGCAGGTGTTCCGTCCAATCTCCAACAACGATCTGATCGTTGGTGCAGTTGGTGTGCTGCAGTTTGATGTGGTGGTAGCGCGCCTGAAGAGCGAATACAACGTTGAAGCAgtg                                                                      <  1:15509/148‑1 (MQ=255)
                                                                               aaaGGGCTGGTACAGCTTTCCGAAGAGGGCGCGGTGCAGGTGTTCCGTCCAATCTCCAACAACGATCTGATCGtt                                                                                                                                              >  2:84413/1‑75 (MQ=255)
                                                                               aaaGGGCTGGTACAGCTTTCCGAAGAGGGCGCGGTGCAGGTGTTCCGTCCAATCTCCAACAACGATCTGATCGtt                                                                                                                                              <  1:84413/75‑1 (MQ=255)
                                                                               aaaGGGCTGGTACAGCTTTCCGAAGAGGGCGCGGTGCAGGTGTTCCGTCCAATCTCCAACAACGATCTGATCGTTGGTGCAGTTGGTGTGCTGCAGTTTGATGTGGTGGTAGCGCGCCTGAAGAGCGAATACAACGTTGAAGCAgtgt                                                                     <  3:126875/148‑1 (MQ=255)
                                                                               aaaGGGCTGGTACAGCTTTCCGAAGAGGGCGCGGTGCAGGTGTTCCGTCCAATCTCCAACAACGATCTGATCGTTGGTGCAGTTGGTGTGCTGCAGTTTGATGTGGTGGTAGCGCGCCTGAAGAGCGAATACAACGTTGAAGCAgtgt                                                                     <  2:1108794/148‑1 (MQ=255)
                                                                                aaGGGCTGGTACAGCTTTCCGAAGAGGGCGCGGTGCAGGTGTTCCGTCCAATCTCCAACAACGATCTGATCGTTGGTGCAGTTGGTGTGCTGCAGTTTg                                                                                                                     >  2:618966/1‑99 (MQ=255)
                                                                                aaGGGCTGGTACAGCTTTCCGAAGAGGGCGCGGTGCAGGTGTTCCGTCCAATCTCCAACAACGATCTGATCGTTGGTGCAGTTGGTGTGCTGCAGTTTGATGTGGTGGTAGCGCGCCTGAAGAGCGAATACAACGTTGAAGCAGTGTa                                                                    >  1:33502/1‑148 (MQ=255)
                                                                                   ggCTGGTACAGCTTTCCGAAGAGGGCGCGGTGCAGGTGTTCCGTCCAATCTCCAACAACGATCTGATCGTTGGTGCAGTTGGTGTGCTGCAGTTTGATGTGGTGGTAGCGCGCCTGaa                                                                                               >  1:163059/1‑118 (MQ=255)
                                                                                   ggCTGGTACAGCTTTCCGAAGAGGGCGCGGTGCAGGTGTTCCGTCCAATCTCCAACAACGATCTGATCGTTGGTGCAGTTGGTGTGCTGCAGTTTGATGTGGTGGTAGCGCGCCTGAAGAGCGAATACAACGt                                                                                <  3:683210/133‑1 (MQ=255)
                                                                                   ggCTGGTACAGCTTTCCGAAGAGGGCGCGGTGCAGGTGTTCCGTCCAATCTCCAACAACGATCTGATCGTTGGTGCAGTTGGTGTGCTGCAGTTTGATGTGGTGGTAGCGCGCCTGAAGAGCGAATACAACGt                                                                                >  4:683210/1‑133 (MQ=255)
                                                                                    gCTGGTACAGCTTTCCGAAGAGGGCGCGGTGCAGGTGTTCCGTCCAATCTCCAACAACGATCTGATCGTTGGTGCAGTTGGTGTGCTGCAGttt                                                                                                                      >  2:564999/1‑94 (MQ=255)
                                                                                    gCTGGTACAGCTTTCCGAAGAGGGCGCGGTGCAGGTGTTCCGTCCAATCTCCAACAACGATCTGATCGTTGGTGCAGTTGGTGTGCTGCAGttt                                                                                                                      <  1:564999/94‑1 (MQ=255)
                                                                                          aCAGCTTTCCGAAGAGGGCGCGGTGCAGGTGTTCCGTCCAATCTCCAACAACGATCTGATCGTTGGTGCAGTTGGTGTGCTGCAGttt                                                                                                                      <  2:716393/88‑1 (MQ=255)
                                                                                          aCAGCTTTCCGAAGAGGGCGCGGTGCAGGTGTTCCGTCCAATCTCCAACAACGATCTGATCGTTGGTGCAGTTGGTGTGCTGCAGttt                                                                                                                      >  1:716393/1‑88 (MQ=255)
                                                                                             gCTTTCCGAAGAGGGCGCGGTGCAGGTGTTCCGTCCAATCTCCAACAACGATCTGATCGtt                                                                                                                                              <  4:770826/61‑1 (MQ=255)
                                                                                                ttCCGAAGAGGGCGCGGTGCAGGTGTTCCGTCCAATCTCCAACAACGATCTGATCGTTGGTGCAGTTGGTGTGCTGCAGTTTGATgtgg                                                                                                               >  3:557604/1‑89 (MQ=255)
                                                                                                   cGAAGAGGGCGCGGTGCAGGTGTTCCGTCCAATCTCCACCAACGATCAGATCGTTgtt                                                                                                                                           >  2:406648/1‑56 (MQ=37)
                                                                                                   cGAAGAGGGCGCGGTGCAGGTGTTCCGTCCAATCTCCACCAACGATCAGATCGTTgtt                                                                                                                                           <  1:406648/58‑3 (MQ=37)
                                                                                                     aaGAGGGCGCGGTGCAGGTGTTCCGTCCAATCTCCAACAACGATCTGATCGTTGGTGCAGTTGGTGTGCTGCAGTTTGATGTGGTGGTAGCGCGCCTGAAGAGCGAATACAACg                                                                                 >  2:382623/1‑114 (MQ=255)
                                                                                                     aaGAGGGCGCGGTGCAGGTGTTCCGTCCAATCTCCAACAACGATCTGATCGTTGGTGCAGTTGGTGTGCTGCAGTTTGATGTGGTGGTAGCGCGCCTGAAGAGCGAATACAACg                                                                                 <  1:382623/114‑1 (MQ=255)
                                                                                                          ggCGCGGTGCAGGTGTTCCGTCCAATCTCCAACAACGATCTGATCGTTGGTGCAGTTGGTGTGCTGCAGTTTGATGTGGTGGTAGCGCGCCTGAAGAGCGAATACAACGTTGAAGCAg                                                                        <  4:515077/118‑1 (MQ=255)
                                                                                                          ggCGCGGTGCAGGTGTTCCGTCCAATCTCCAACAACGATCTGATCGTTGGTGCAGTTGGTGTGCTGCAGTTTGATGTGGTGGTAGCGCGCCTGAAGAGCGAATACAACGTTGAAGCAg                                                                        >  3:515077/1‑118 (MQ=255)
                                                                                                                 tGCAGGTGTTCCGTCCAATCTCCAACAACGATCTGATCGTTGGTGCAGTTGGTGTGCTGCAGTTTGATGTGGTGGTAGCGCGCCTGAAGAGCGAATACAACGTTGAAGCAGTGTATGAGTCAGTCAACGTTGCCACTGCCCGCTGGGt                                   <  2:33502/148‑1 (MQ=255)
                                                                                                                     ggTGTTCCGTCCAATCTCCAACAACGATCTGATCGTTGGTGCAGTTGGTGTGCTGCAGTTTGATGTGGTGGTAGCGCGCCTGAAGAGCGAATACAACGTTGAAGCAGTGTATGAGTCAGTCAACGTTGCCACTGCCCGCTGGGTAGa                                <  1:823029/147‑1 (MQ=255)
                                                                                                                     ggTGTTCCGTCCAATCTCCAACAACGATCTGATCGTTGGTGCAGTTGGTGTGCTGCAGTTTGATGTGGTGGTAGCGCGCCTGAAGAGCGAATACAACGTTGAAGCAGTGTATGAGACAGTCAACGTTGCCACTGCCCGCTGGGTAGa                                >  2:823029/1‑147 (MQ=255)
                                                                                                                      gtgtTCCGTCCAATCTCCAACAACGATCTGATCGTTGGTGCAGTTGGTGTGCTGCAGTTTGATGTGGTGGTAGCGCGCCTGAAGAGCGAATACAACGTTGAAGCAGTGTATGAGTCAGTCAACGTTGCCACTGCCCGCt                                       <  3:347868/139‑1 (MQ=255)
                                                                                                                      gtgtTCCGTCCAATCTCCAACAACGATCTGATCGTTGGTGCAGTTGGTGTGCTGCAGTTTGATGTGGTGGTAGCGCGCCTGAAGAGCGAATACAACGTTGAAGCAGTGTATGAGTCAGTCAACGTTGCCACTGCCCGCt                                       >  4:347868/1‑139 (MQ=255)
                                                                                                                        gtTCCGTCCAATCTCCAACAACGATCTGATCGTTGGTGCAGTTGGTGTGCTGCAGTTTGATGTGGTGGTAGCGCGCCTGAAGAGCGAATACAACGTTGAAGCAGTGTATg                                                                  <  1:999444/110‑1 (MQ=255)
                                                                                                                          tCCGTCCAATCTCCAACAACGATCTGATCGTTGGTGCAGTTGGTGTGCTGCAGTTTGa                                                                                                                    >  2:162215/1‑58 (MQ=255)
                                                                                                                          tCCGTCCAATCTCCAACAACGATCTGATCGTTGGTGCAGTTGGTGTGCTGCAGTTTGATGTGGTGGTAGCGCGCCTGAAGAGCGAATACAACGTTGAAGCAGTGTATGAGTCAGTCAACGTTGCCACTGCCCGCTggg                                    <  1:763534/138‑1 (MQ=255)
                                                                                                                               ccAATCTCCAACGACGATCTGATCGTTGGTGCAGTTGGTGTGCTGCAGTTTGATGTGGTGGTAGCGCGCCTGAAGAGCGAATACAACGTTGa                                                                             >  4:11058/1‑92 (MQ=255)
                                                                                                                               ccAATCTCCAACGACGATCTGATCGTTGGTGCAGTTGGTGTGCTGCAGTTTGATGTGGTGGTAGCGCGCCTGAAGAGCGAATACAACGTTGa                                                                             <  3:11058/92‑1 (MQ=255)
                                                                                                                               ccAATCTCCAACAACGATCTGATCGTTGGTGCAGTTGGTGTGCTGCAGTTTGATGTGGTGGTAGCGCGCCTGAAGAGCGAATACAACGTTGa                                                                             <  1:756389/92‑1 (MQ=255)
                                                                                                                               ccAATCTCCAACAACGATCTGATCGTTGGTGCAGTTGGTGTGCTGCAGTTTGATGTGGTGGTAGCGCGCCTGAAGAGCGAATACAACGTTGa                                                                             >  2:756389/1‑92 (MQ=255)
                                                                                                                               ccAATCTCCAACAACGATCTGATCGTTGGTGCAGTTGGTGTGCTGCAGTTTGATGTGGTGGTAGCGCGCCTGAAGAGCGAATACAACGTTGAAGCAGTGTATGAGTCAGTCAACGTTGCCACTGCCCGCTg                                      >  2:1658/1‑131 (MQ=255)
                                                                                                                                cAATCTCCAACAACGATCTGATCGTTGGTGCAGTTGGTGTGCTGCAGTTTGATGTGGTGGTAGCGCGCCTGAAGAGCGAATACAACGTTGAAGCAGTGTATGAGTCAGTCAACGTTGCCACTGCCCGCTGGGTAGa                                >  1:926917/1‑136 (MQ=255)
                                                                                                                                cAATCTCCAACAACGATCTGATCGTTGGTGCAGTTGGTGTGCTGCAGTTTGATGTGGTGGTAGCGCGCCTGAAGAGCGAATACAACGTTGAAGCAGTGTATGAGTCAGTCAACGTTGCCACTGCCCGCTGGGTAGa                                <  2:926917/136‑1 (MQ=255)
                                                                                                                                 aaTCTCCAACAACGATCTGATCGTTGGTGCAGTTGGTGTGCTGCAGTTTGATGTGGTGGTAGCGCGCCTGAAGAGCGAATACAACGTTGAAGCAGTGTATGAGTCAGTCAACGTTGCCACTGCCCGCTGGGTAGAATGTGCAGAcgcg                   >  4:742898/1‑148 (MQ=255)
                                                                                                                                        aacaacGATCTGATCGTTGGTGCAGTTGGTGTGCTGCAGTTTGATGTGGTGGTAGCGCGCCTGAAGAGCGAATACAACGTTGAAGCAGTGTATGAGTCAGTCAACGTTGCCACTGCCCGCTGGGTAGAATGTGCAGACGCGAAGAAAt            >  3:1133060/1‑148 (MQ=255)
                                                                                                                                         ccaacGATCAGATCGTTGGTGCAGTTGGTGTGCTGCAGTTTGATGTGGTGGTAGCGCGCCTGAAGAGCGAATACAACGTTGAAGCAGTGTATGAGTCAGTCAACGTTGCCACTGcc                                           <  3:748335/115‑1 (MQ=255)
                                                                                                                                         acaacGATCTGATCGTTGGTGCAGTTGGTGTGCTGCAGTTTGATGTGGTGGTAGCGCGCCTGAAGAGCGAATACAACGTTGAAGCAGTGTATg                                                                  >  4:869841/1‑93 (MQ=255)
                                                                                                                                         acaacGATCTGATCGTTGGTGCAGTTGGTGTGCTGCAGTTTGATGTGGTGGTAGCGCGCCTGAAGAGCGAATACAACGTTGAAGCAGTGTATg                                                                  <  3:869841/93‑1 (MQ=255)
                                                                                                                                         acaacGATCTGATCGTTGGTGCAGTTGGTGTGCTGCAGTTTGATGTGGTGGTAGCGCGCCTGAAGAGCGAATACAACGTTGAAGCAGTGTATGAGTCAGTCAACGTTGCCACTGCCCGCTGGGTAGa                                <  3:851247/127‑1 (MQ=255)
                                                                                                                                         acaacGATCTGATCGTTGGTGCAGTTGGTGTGCTGCAGTTTGATGTGGTGGTAGCGCGCCTGAAGAGCGAATACAACGTTGAAGCAGTGTATGAGTCAGTCAACGTTGCCACTGCCCGCTGGGTAGAAtgtg                           <  1:179231/132‑1 (MQ=255)
                                                                                                                                         acaacGATCTGATCGTTGGTGCAGTTGGTGTGCTGCAGTTTGATGTGGTGGTAGCGCGCCTGAAGAGCGAATACAACGTTGAAGCAGTGTATGAGTCAGTCAACGTTGCCACTGCCCGCTGGGTAGAATGTGCAGACGCgaaga               <  4:179576/144‑1 (MQ=255)
                                                                                                                                          caacGATCTGATCGTTGGTGCAGTTGGTGTGCTGCAGTTTGATGTGGTGGTAGCGCGCCTGAAGAGCGAATACAACGTTGAAGCAGTGTATGAGTCAGTCAACGTTGcc                                                 <  2:242526/109‑1 (MQ=255)
                                                                                                                                          caacGATCTGATCGTTGGTGCAGTTGGTGTGCTGCAGTTTGATGTGGTGGTAGCGCGCCTGAAGAGCGAATACAACGTTGAAGCAGTGTATGAGTCAGTCAACGTTGcc                                                 >  1:242526/1‑109 (MQ=255)
                                                                                                                                          caacGATCTGATCGTTGGTGCAGTTGGTGTGCTGCAGTTTGATGTGGTGGTAGCGCGCCTGAAGAGCGAATACAACGTTGAAGCAGTGTATGAGTCAGTCAACGTTGCCACTGCCCGCTGGGTAGAATGTGCAGACGCGAAGAAATTc          >  3:178227/1‑148 (MQ=255)
                                                                                                                                          caacGATCAGATCGTTGGTGCAGTTGGTGTGCTGCAGTTTGATGTGGTGGTAGCGCGCCTGAAGAGCGaa                                                                                        >  1:501807/1‑70 (MQ=255)
                                                                                                                                          caacGATCAGATCGTTGGTGCAGTTGGTGTGCTGCAGTTTGATGTGGTGGTAGCGCGCCTGAAGAGCGaa                                                                                        <  2:501807/70‑1 (MQ=255)
                                                                                                                                           aacGATCTGATCGTTGGTGCAGTTGGTGTGCTGCAGTTTGATGTGGTGGTAGCGCGCCTGAAGAGCGa                                                                                         >  1:986965/1‑68 (MQ=255)
                                                                                                                                            acGATCTGATCGTTGGTGCAGTTGGTGTGCTGCAGTTTGATGTGGTGGTAGCGCGCCTGAAGAGCGAATACAACGTTGaa                                                                            <  3:1067483/80‑1 (MQ=255)
                                                                                                                                            acGATCTGATCGTTGGTGCAGTTGGTGTGCTGCAGTTTGATGTGGTGGTAGCGCGCCTGAAGAGCGAATACAACGTTGaa                                                                            >  4:1067483/1‑80 (MQ=255)
                                                                                                                                            acGATCTGATCGTTGGTGCAGTTGGTGTGCTGCAGTTTGATGTGGTGGTAGCGCGCCTGAAGAGCGAATACAACGTTGAAGCAGTGTATGAGTCAGTCAACGTTGCCACTGCCCGCTGGGTAGAATGTGCAGACGCga                  <  1:696235/138‑1 (MQ=255)
                                                                                                                                            acGATCTGATCGTTGGTGCAGTTGGTGTGCTGCAGTTTGATGTGGTGGTAGCGCGCCTGAAGAGCGAATACAACGTTGAAGCAGTGTATGAGTCAGTCAACGTTGCCACTGCCCGCTGGGTAGAATGTGCAGACGCga                  >  2:696235/1‑138 (MQ=255)
                                                                                                                                            acGATCTGATCGTTGGTGCAGTTGGTGTGCTGCAGTTTGATGTGGTGGTAGCGCGCCTGAAGAGCGAATACAACGTTGAAGCAGTGTATGAGTCAGTCAACGTTGCCACTGCCCGCTGGGTAGAATGTGCAGACGCGAAGAAATTCGa        <  4:1133060/148‑1 (MQ=255)
                                                                                                                                             cGATCTGATCGTTGGTGCAGTTGGTGTGCTGCAGTTTGATGTGGTGGTAGCGCGcc                                                                                                   <  2:666473/56‑1 (MQ=255)
                                                                                                                                             cGATCTGATCGTTGGTGCAGTTGGTGTGCTGCAGTTTGATGTGGTGGTAGCGCGcc                                                                                                   >  1:666473/1‑56 (MQ=255)
                                                                                                                                                tCTGATCGTTGGTGCAGTTGGTGTGCTGCAGTTTGATGTGGTGGTAGCGCGCCTGAAGAGCGa                                                                                         <  1:407655/63‑1 (MQ=255)
                                                                                                                                                tCTGATCGTTGGTGCAGTTGGTGTGCTGCAGTTTGATGTGGTGGTAGCGCGCCTGAAGAGCGa                                                                                         >  2:407655/1‑63 (MQ=255)
                                                                                                                                                tCTGATCGTTGGTGCAGTTGGTGTGCTGCAGTTTGATGTGGTGGTAGCGCGCCTGAAGAGCGAATACAACGTTGAAGc                                                                          >  3:656562/1‑78 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
CACCGGTATTCCGAACTTCGCACCAGAACTGTTCCGTCGTATCCGCCTGAAAGATCCGCTGAAGCAAAAACAGCTGCTCAAAGGGCTGGTACAGCTTTCCGAAGAGGGCGCGGTGCAGGTGTTCCGTCCAATCTCCAACAACGATCTGATCGTTGGTGCAGTTGGTGTGCTGCAGTTTGATGTGGTGGTAGCGCGCCTGAAGAGCGAATACAACGTTGAAGCAGTGTATGAGTCAGTCAACGTTGCCACTGCCCGCTGGGTAGAATGTGCAGACGCGAAGAAATTCGAAGAGTT  >  NC_000913/4610574‑4610867

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: