Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 268,471 (C)5→6 coding (490/864 nt) ykfA ← CP4‑6 prophage; putative GTP‑binding protein YkfA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913268,4661.C98.9% 305.9 / ‑3.8 87coding (495/864 nt)ykfACP4‑6 prophage; putative GTP‑binding protein YkfA
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base C (43/43);  total (43/44)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TAATCATCCGCTCGGCCATCACCTTAAGCCCCCATTGTAGACGGACCGACACAGCACACACCGGATGGACGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGTCCCCCACCGCTGGTAGGTTCGGCCTTATCCGACTGGCTGATAACAAACAGCACCTTATGCCGGTATACCTCTCCAATCACCTGATGATAAAAATGTTCATCCACCGTCAGCGCCCGGTCATCGGCCTTAATCA  >  NC_000913/268327‑268599
                                                                                                                                            |                                                                                                                                      
tAATCATCCGCTCGGCCATCACCTTAAGCCCCCATTGTAGACGGACCGACACAGCACACACCGGATGGACGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGt                                                                                                                                         >  1:266310/1‑140 (MQ=255)
tAATCATCCGCTCGGCCATCACCTTAAGCCCCCATTGTAGACGGACCGACACAGCACACACCGGATGGACGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGt                                                                                                                                         <  1:1081576/140‑1 (MQ=255)
tAATCATCCGCTCGGCCATCACCTTAAGCCCCCATTGTAGACGGACCGACACAGCACACACCGGATGGACGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGt                                                                                                                                         <  1:607851/140‑1 (MQ=255)
       ccGCTCGGCCATCACCTTAAGCCCCCATTGTATACGGACCCACACAGCACACACCGGATGGACGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTGCGGCTCATATTCTGCTTCTGCGCCGTGGACAACTGTcccccca                                                                                                                                <  2:1090074/140‑3 (MQ=255)
       ccGCTCGGCCATCACCTTAAGCCCCCATTGTAGACGGACCGACACAGCACACACCGGATGGACGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGTcccccca                                                                                                                                <  1:11074/140‑3 (MQ=255)
       ccGCTCGGCCATCACCTTAAGCCCCCATTGTAGACGGACCGACACAGCACACACCGGATGGACGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGTcccccca                                                                                                                                >  1:832331/1‑138 (MQ=255)
         gCTCGGCCATCACCTTAAGCCCCCATTGTAGACGGACCGACACAGCACACACCGGATGGACGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGTccccccacc                                                                                                                              >  1:1015127/1‑136 (MQ=255)
         gCTCGGCCATCACCTTAAGCCCCCATTGTAGACGGACCGACACAGCACACACCGGATGGACGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGTccccccacc                                                                                                                              >  2:557256/1‑136 (MQ=255)
         gCTCGGCCATCACCTTAAGCCCCCATTGTAGACGGACCGACACAGCACACACCGGATGGACGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGTccccccacc                                                                                                                              >  1:78226/1‑136 (MQ=255)
          cTCGGCCATCACCTTAAGCCCCCATTGTAGACGGACCGACACAGCACACACCGGATGGACGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGTccccccaccg                                                                                                                             >  2:734796/1‑135 (MQ=255)
              gCCATCACCTTAAGCCCCCATTGTAGACGGACCGACACAGCACACACCGGATGGACGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGTCCGCCCACCGCTgg                                                                                                                          <  1:218346/140‑1 (MQ=255)
                         aaGCCCCCATTGTAGACGGACCGACACAGCACACACCGGATGGACGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGTCCCCCCACCGCTGGTAGGTTCGGcc                                                                                                               >  2:20435/1‑140 (MQ=255)
                              cccATTGTAGACGGACCGACACAGCACGCACCGGATGGACGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGTCCCCCCACCGCTgg                                                                                                                          <  2:206625/124‑1 (MQ=255)
                              cccATTGTAGACGGACCGACACAGCACGCACCGGATGGACGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGTCCCCCCACCGCTgg                                                                                                                          >  1:206625/1‑124 (MQ=255)
                               ccATTGTAGACGGACCGACACAGCACACACCGGATGGACGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGTCCCCCCACCGCTGGTAGGTTCGGCCTTATcc                                                                                                         >  1:155523/1‑140 (MQ=255)
                                cATTGTAGACGGACCGACACAGCACACACCGGATGGACGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGTCCCCCCACCGCTGGTAGGTTCGGCCTTATCCg                                                                                                        <  1:371233/140‑1 (MQ=255)
                                cATTGTAGACGGACCGACACAGCACACACCGGATGGACGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGTCCCCCCACCGCTGGTAGGTTCGGCCTTATCCg                                                                                                        >  1:180146/1‑140 (MQ=255)
                                cATTGTAGACGGACCGACACAGCACACACCGGATGGACGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGTCCCCCCACCGCTGGTAGGTTCGGCCTTATCCg                                                                                                        >  1:247745/1‑140 (MQ=255)
                                      aGACGGACCGACACAGCACACACCGGATGGACGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGTCCCCCCACCGCTGGTAGGTTCGGCCTTATCCGACTGGc                                                                                                  <  1:22494/140‑1 (MQ=255)
                                       gACGGACCGACACAGCACACACCGGATGGACGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGTCCCCCCACCGCTGGTAGGTTCGGCCTTATCCGAc                                                                                                      >  1:899930/1‑135 (MQ=255)
                                       gACGGACCGACACAGCACACACCGGATGGACGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGTCCCCCCACCGCTGGTAGGTTCGGCCTTATCCGAc                                                                                                      <  2:899930/135‑1 (MQ=255)
                                          ggACCGACACAGCACACACCGGATGGACGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGTCCCCCCACCGCTGGTAGGTTCGGCCTTATCCGa                                                                                                       <  1:1141083/131‑1 (MQ=255)
                                          ggACCGACACAGCACACACCGGATGGACGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGTCCCCCCACCGCTGGTAGGTTCGGCCTTATCCGa                                                                                                       >  2:1141083/1‑131 (MQ=255)
                                          ggACCGACACAGCACACACCGGATGGACGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGTCCCCCCACCGCTGGTAGGTTCGGCCTTATCCGACTGGCTGAt                                                                                              >  2:177816/1‑140 (MQ=255)
                                           gACCGACACAGCACACACCGGATGGACGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGTCCCCCCACCGCTGGTAGGTTCGGCCTTATCCGACTGGCTGATa                                                                                             <  1:348229/140‑1 (MQ=255)
                                           gACCGACACAGCACACACCGGATGGACGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGTCCCCCCACCGCTGGTAGGTTCGGCCTTATCCGACTGGCTGATa                                                                                             >  2:15795/1‑140 (MQ=255)
                                            aCCGACACAGCACACACCGGATGGACGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGTCCCCCCACCGCTGGTAGGTTCGGCCTTATCCGACTGGCTGATaa                                                                                            >  1:274582/1‑140 (MQ=255)
                                            aCCGACACAGCACACACCGGATGGACGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGTCCCCCCACCGCTGGTAGGTTCGGCCTTATCCGACTGGCTGATaa                                                                                            >  2:482168/1‑140 (MQ=255)
                                            aCCGACACAGCACACACCGGATGGACGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGTCCCCCCACCGCTGGTAGGTTCGGCCTTATCCGACTGGCTGATaa                                                                                            >  1:166774/1‑140 (MQ=255)
                                            aCCGACACAGCACACACCGGATGGACGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGTCCCCCCACCGCTGGTAGGTTCGGCCTTATCCGACTGGCTGATaa                                                                                            >  1:563697/1‑140 (MQ=255)
                                            aCCGACACAGCACACACCGGATGGACGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGTCCCCCCACCGCTGGTAGGTTCGGCCTTATCCGACTGGCTGATa                                                                                             >  1:1106413/1‑139 (MQ=255)
                                            aCCGACACAGCACACACCGGATGGACGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGTCCCCCCACCGCTGGTAGGTTCGGCCTTATCCGACTGGCTGATa                                                                                             <  2:1106413/139‑1 (MQ=255)
                                             ccGACACAGCACACACCGGATGGACGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGTCCCCCCACCGCTGGTAGGTTCGGCCTTATCCGACTGGCTGATaac                                                                                           >  2:570196/1‑140 (MQ=255)
                                             ccGACACAGCACACACCGGATGGACGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGTCCCCCCACCGCTGGTAGGTTCGGCCTTATCCGACTGGCTGATaac                                                                                           >  1:1168812/1‑140 (MQ=255)
                                             ccGACACAGCACACACCGGATGGACGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGTCCCCCCACCGCTGGTAGGTTCGGCCTTATCCGACTGGCTGATaac                                                                                           >  2:455040/1‑140 (MQ=255)
                                               gACACAGCACACACCGGATGGACGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGTCCCCCCACCGCTGGTAGGTTCGGCCTTATCCGACTGGCTGATaacaa                                                                                         >  1:111331/1‑140 (MQ=255)
                                                acacAGCACACACCGGATGGACGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGTCCCCCCACCGCTGGTAGGTTCGGCCTTATCCGACTGGCTGATaacaaa                                                                                        <  2:301278/140‑1 (MQ=255)
                                                   caGCACACACCGGATGGACGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGTCCCCCCACCGCTGGTAGGTTCGGCCTTATCCGACTGGCTGATAACAAACAg                                                                                     <  1:947935/140‑1 (MQ=255)
                                                   caGCACACACCGGATGGACGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGTCCCCCCACCGCTGGTAGGTTCGGCCTTATCCGACTGGCTGATAACAAACAg                                                                                     <  1:20435/140‑1 (MQ=255)
                                                      cacacaCCGGATGGACGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGTCCCCCCACCGCTGGTAGGTTCGGCCTTATCCGACTGGCTGATAACAAACAg                                                                                     <  2:295908/137‑1 (MQ=255)
                                                      cacacaCCGGATGGACGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGTCCCCCCACCGCTGGTAGGTTCGGCCTTATCCGACTGGCTGATAACAAACAg                                                                                     >  1:295908/1‑137 (MQ=255)
                                                       acacacCGGATGGACGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTGCGGCTGATATTCTGCTTCTGCTCCGTGGACAACTGTCCCCCCACCGCTGGTAGGTTCGGCCTTATCCGACTGGCTGATAACAAACAGCAcc                                                                                 <  2:273790/140‑1 (MQ=255)
                                                       acacacCGGATGGACGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGTCCCCCCACCGCTGGTAGGTTCGGCCTTATCCGACTGGCTGATAACAAACAGCAcc                                                                                 <  2:78226/140‑1 (MQ=255)
                                                          cacCGGATGGACGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGTCCCCCCACCGCTGGTAGGTTCGGCCTTATCCGACTGGCTGATAACAAACAGCACCTTa                                                                              <  2:247745/140‑1 (MQ=255)
                                                          cacCGGATGGACGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGTCCCCCCACCGCTGGTAGGTTCGGCCTTATCCGACTGGCTGATAACAAACAGCACCTTa                                                                              <  1:1093688/140‑1 (MQ=255)
                                                          cacCGGATGGACGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGTCCCCCCACCGCTGGTAGGTTCGGCCTTATCCGACTGGCTGATAACAAACAGCACCTTa                                                                              <  1:539795/140‑1 (MQ=255)
                                                          cacCGGATGGACGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGTCCCCCCACCGCTGGTAGGTTCGGCCTTATCCGACTGGCTGATAACAAACAGCACCTTa                                                                              <  1:482168/140‑1 (MQ=255)
                                                           acCGGATGGACGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGTCCCCCCACCGCTGGTAGGTTCGGCCTTATCCGACTGGCTGATAACAAACAGCACCTTAt                                                                             <  1:919521/140‑1 (MQ=255)
                                                           acCGGATGGACGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGTCCCCCCACCGCTGGTAGGTTCGGCCTTATCCGACTGGCTGATAACAAACAGCACCTTAt                                                                             >  2:436369/1‑140 (MQ=255)
                                                           acCGGATGGACGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGTCCCCCCACCGCTGGTAGGTTCGGCCTTATCCGACTGGCTGATAACAAACAGCACCTTAt                                                                             <  2:111331/140‑1 (MQ=255)
                                                           acCGGATGGACGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGTCCCCCCACCGCTGGTAGGTTCGGCCTTATCCGACTGGCTGATAACAAACAGCACCTTAt                                                                             <  2:1015127/140‑1 (MQ=255)
                                                           acCGGATGGACGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGTCCCCCCACCGCTGGTAGGTTCGGCCTTATCCGACTGGCTGATAACAAACAGCACCTTAt                                                                             >  2:627867/1‑140 (MQ=255)
                                                           acCGGATGGACGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGTCCCCCCACCGCTGGTAGGTTCGGCCTTATCCGACTGGCTGATAACAAACAGCACCTTAt                                                                             <  1:436369/140‑1 (MQ=255)
                                                           acCGGATGGACGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGTCCCCCCACCGCTGGTAGGTTCGGCCTTATCCGACTGGCTGATAACAAACAGCACCTTAt                                                                             <  1:734796/140‑1 (MQ=255)
                                                           acCGGATGGACGGGCTGGAACAGCTCATGCAGCAGGCAGAATTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGTCCCCCCACCGCTGGTAGGTTCGGCCTTATCCGACTGGCTGATaacaaaca                                                                                      >  2:1046492/1‑131 (MQ=255)
                                                           acCGGATGGACGGGCTGGAACAGCTCATGCAGCAGGCAGAATTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGTCCCCCCACCGCTGGTAGGTTCGGCCTTATCCGACTGGCTGATaacaaaca                                                                                      <  1:1046492/131‑1 (MQ=255)
                                                            ccGGATGGACGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGTCCCCCCACCGCTGGTAGGTTCGGCCTTATCCGACTGGCTGATAACAAACAGCACCTTATg                                                                            <  1:177816/140‑1 (MQ=255)
                                                                aTGGACGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGTCCCCCCACCGCTGGTAGGTTCGGCCTTATCCGACTGGCTGATAACAAACAGCACCTTATGCCgg                                                                        >  2:695746/1‑140 (MQ=255)
                                                                aTGGACGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGTCCCCCCACCGCTGGTAGGTTCGGCCTTATCCGACTGGCTGATAACAAACAGCACCTTATGCCgg                                                                        >  1:521456/1‑140 (MQ=255)
                                                                  ggACGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGTCCCCCCACCGCTGGTAGGTTCGGCCTTATCCGACTGGCTGATAACAAACAGCACCTTATGCCGGta                                                                      <  2:832331/140‑1 (MQ=255)
                                                                     cGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGTCCCCCCACCGCTGGTAGGTTCGGCCTTATCCGACTGGCTGATAACAAACAGCACCTTa                                                                              <  2:731488/129‑1 (MQ=255)
                                                                     cGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGTCCCCCCACCGCTGGTAGGTTCGGCCTTATCCGACTGGCTGATAACAAACAGCACCTTa                                                                              >  1:731488/1‑129 (MQ=255)
                                                                       ggCTGGAACAGCTCATGCAGCAGGCAGATTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGTCCCCCCACCGCTGGTAGGTTCGGCCTTATCCGACTGGCTGATAACAAACAGCACCTTATGCCGGTATACct                                                                 >  2:72765/1‑140 (MQ=255)
                                                                       ggCTGGAACAGCTCATGCAGCAGGCAGATTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGTCCCCCCACCGCTGGTAGGTTCGGCCTTATCCGACTGGCTGATAACAAACAGCACCTTATGCCGGTATACct                                                                 <  1:15795/140‑1 (MQ=255)
                                                                       ggCTGGAACAGCTCATGCAGCAGGCAGATTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGTCCCCCCACCGCTGGTAGGTTCGGCCTTATCCGACTGGCTGATAACAAACAGCACCTTATGCCGGTATACct                                                                 >  1:14905/1‑140 (MQ=255)
                                                                         cTGGAACAGCTCATGCAGCAGGCAGATTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGTCCCCCCACCGCTGGTAGGTTCGGCCTTATCCGACTGGCTGATAACAAACAGCACCTTAt                                                                             >  1:292558/1‑126 (MQ=255)
                                                                         cTGGAACAGCTCATGCAGCAGGCAGATTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGTCCCCCCACCGCTGGTAGGTTCGGCCTTATCCGACTGGCTGATAACAAACAGCACCTTAt                                                                             <  2:292558/126‑1 (MQ=255)
                                                                         cTGGAACAGCTCATGCAGCAGGCAGATTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGTCCCCCCACCGCTGGTAGGTTCGGCCTTATCCGACTGGCTGATAACAAACAGCACCTTATGCCgg                                                                        <  1:588004/131‑1 (MQ=255)
                                                                         cTGGAACAGCTCATGCAGCAGGCAGATTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGTCCCCCCACCGCTGGTAGGTTCGGCCTTATCCGACTGGCTGATAACAAACAGCACCTTATGCCgg                                                                        >  2:588004/1‑131 (MQ=255)
                                                                            gAACAGCTCATGCAGCAGGCAGATTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGTCCCCCCACCGCTGGTAGGTTCGGCCTTATCCGACTGGCTGATAACAAACAGCACCTTATGCCGGTATACCTCTCCa                                                            <  2:87041/140‑1 (MQ=255)
                                                                                  cTCATGCAGCAGGCAGATTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGTCCCCCCACCGCTGGTAGGTTCGGCCTTATCCGACTGGCTGATAACAAACAGCACCTTATGCCGGTATACCTCTCCAATCAcc                                                      >  1:1006629/1‑140 (MQ=255)
                                                                                  cTCATGCAGCAGGCAGATTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGTCCCCCCACCGCTGGTAGGTTCGGCCTTATCCGACTGGCTGATAACAAACAGCACCTTATGCCGGTATACCTCTCCAATCAcc                                                      >  1:842970/1‑140 (MQ=255)
                                                                                  cTCATGCAGCAGGCAGATTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACACTGTCCCCCCACCGCTGGTAGGTTCGGCCTTATCCGACTGGCTGATAACAAACAGCACCTTATGCCGGTATACCTCTCCAATCACCt                                                     <  1:263154/140‑1 (MQ=255)
                                                                                   tCATGCAGCAGGCAGATTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGTCCCCCCACCGCTGGTAGGTTCGGCCTTATCCGACTGGCTGATAACAAACAGCACCTTATGCCGGTATACCTCTCCAATCACCt                                                     <  1:695746/140‑1 (MQ=255)
                                                                                   tCATGCAGCAGGCAGATTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGTCCCCCCACCGCTGGTAGGTTCGGCCTTATCCGACTGGCTGATAACAAACAGCACCTTATGCCGGTATACCTCTCCAATCACCt                                                     <  1:476620/140‑1 (MQ=255)
                                                                                     aTGCAGCAGGCAGATTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGTCCCCCCACCGCTGGTAGGTTCGGCCTTATCCGACTGGCTGATAACAAACAGCACCTTATGCCGGTATACCTCTCCAATCACCtga                                                   <  2:563697/140‑1 (MQ=255)
                                                                                      tGCAGCAGGCAGATTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGTCCCCCCACCGCTGGTAGGTTCGGCCTTATCCGACTGGCTGATAACAAACAGCACCTTATGCCGGTATACCTCTCCAATCACCtgat                                                  <  1:570196/140‑1 (MQ=255)
                                                                                      tGCAGAAGGCAGATTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGTCCCCCCACCGCTGGTAGGTTCGGCCTTATCCGACTGGCTGATAACAAACAGCACCTTATGCCGGTATACCTCTCCAATCACCtgat                                                  >  1:1010170/1‑140 (MQ=255)
                                                                                       gcagcaGGCAGATTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGTCCCCCCACCGCTGGTAGGTTCGGCCTTATCCGACTGGCTGATAACAAACAGCACCTTATGCCGGTATACCTCTCCAATCACCtgatg                                                 >  2:426871/1‑140 (MQ=255)
                                                                                         agcagGCAGATTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGTCCCCCCACCGCTGGTAGGTTCGGCCTTATCCGACTGGCTGATAACAAACAGCACCTTATGCCGGTATACCTCTCCAATCACCTgatgat                                               <  2:166774/140‑1 (MQ=255)
                                                                                         agcagGCAGATTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGTCCCCCCACCGCTGGTAGGTTCGGCCTTATCCGACTGGCTGATAACAAACAGCACCTTATGCCGGTATACCTCTCCAATCACCTgatgat                                               <  2:316929/140‑1 (MQ=255)
                                                                                         agcagGCAGATTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGTCCCCCCACCGCTGGTAGGTTCGGCCTTATCCGACTGGCTGATAACAAACAGCACCTTATGCCGGTATACCTCTCCAATCACCTgatgat                                               <  1:426871/140‑1 (MQ=255)
                                                                                            aggcagATTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGTCCCCCCACCGCTGGTAGGTTCGGCCTTATCCGACTGGCTGATAACAAACAGCACCTTATGCCGGTATACCTCTCCAATCACCTGATGATaaa                                            >  1:577892/1‑140 (MQ=255)
                                                                                            aggcagATTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGTCCCCCCACCGCTGGTAGGTTCGGCCTTATCCGACTGGCTGATAACAAACAGCACCTTATGCCGGTATACCTCTCCAATCACCTGATGATaaa                                            >  2:293444/1‑140 (MQ=255)
                                                                                             ggcagATTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGTCCCCCCACCGCTGGTAGGTTCGGCCTTATCCGACTGGCTGATAACAAACAGCACCTTATGCCGGTATACCTCTCCAATCACCTGATGATaaaa                                           <  2:27025/140‑1 (MQ=255)
                                                                                                     ttGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGTCCCCCCACCGCTGGTAGGTTCGGCCTTATCCGACTGGCTGATAACAAACAGCACCTTATGCCGGTATACCTCTCCAATCACCTGATGATAAAAATGTTCAt                                   >  1:635521/1‑140 (MQ=255)
                                                                                                        cGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGTCCCCCCACCGCTGGTAGGTTCGGCCTTATCCGACTGGCTGATAACAAACAGCACCTTATGCCGGTATACCTCTCCAATCACCTGATGATAAAAATGTTCATCCa                                <  1:293444/140‑1 (MQ=255)
                                                                                                          gCTGATATTCTGCTTCTGCGCCGTGGACAACTGTCCCCCCACCGCTGGTAGGTTCGGCCTTATCCGACTGGCTGATAACAAACAGCACCTTATGCCGGTATACCTCTCCAATCACCTGATGATAAAAATGTTCATCCAcc                              >  2:226917/1‑140 (MQ=255)
                                                                                                               tatTCTGCTTCTGCGCCGTGGACAACTGTCCCCCCACCGCTGGTAGGTTCGGCCTTATCCGACTGGCTGATAACAAACAGCACCTTATGCCGGTATACCTCTCCAATCACCTGATGATAAAAATGTTCATCCACCGTCAg                         <  2:266310/140‑1 (MQ=255)
                                                                                                                         ctgcGCCGTGGACAACTGTCCCCCCACCGCTGGTAGGTTCGGCCTTATCCGACTGGCTGATAACAAACAGCACCTTATGCCGGTATACCTCTCCAATCACCTGATGATAAAAATGTTCATCCACCGTCAGCGCCCGGTCa               >  2:484155/1‑140 (MQ=255)
                                                                                                                         ctgcGCCGTGGACAACTGTCCCCCCACCGCTGGTAGGTTCGGCCTTATCCGACTGGCTGATAACAAACAGCACCTTATGCCGGTATACCTCTCCAATCACCTGATGATAAAAATGTTCATCCACCGTCAGCGCCCGGTCa               >  1:905831/1‑140 (MQ=255)
                                                                                                                           gcgcCGTGGACAACTGTCCCCCCACCGCTGGTAGGTTCGGCCTTATCCGACTGGCTGATAACAAACAGCACCTTATGCCGGTATACCTCTCCAATCACCTGATGATAAAAATGTTCATCCACCGTCAGCGCCCGGTCATc             >  1:591229/1‑140 (MQ=255)
                                                                                                                           gcgcCGTGGACAACTGTCCCCCCACCGCTGGTAGGTTCGGCCTTATCCGACTGGCTGATAACAAACAGCACCTTATGCCGGTATACCTCTCCAATCACCTGATGATAAAAATGTTCATCCACCGTCAGCGCCCGGTCATc             <  2:180146/140‑1 (MQ=255)
                                                                                                                            cgcCGTGGACAACTGTCCCCCCACCGCTGGTAGGTTCGGCCTTATCCGACTGGCTGATAACAAACAGCACCTTATGCCGGTATACCTCTCCAATCACCTGATGATAAAAATGTTCATCCACCGTCAGCGCCCGGTCATCg            >  2:583411/1‑140 (MQ=255)
                                                                                                                            cgcCGTGGACAACTGTCCCCCCACCGCTGGTAGGTTCGGCCTTATCCGACTGGCTGATAACAAACAGCACCTTATGCCGGTATACCTCTCCAATCACCTGATGATAAAAATGTTCATCCACCGTCAGCGCCCGGTCATCg            >  2:223860/1‑140 (MQ=255)
                                                                                                                                gTGGACAACTGTCCCCCCACCGCTGGTAGGTTCGGCCTTATCCGACTGGCTGATAACAAACAGCACCTTATGCCGGTATACCTCTCCAATCACCTGATGATAAAAATGTTCATCCACCGTCAGCGCCCGGTCATCGGCCt        >  1:390245/1‑140 (MQ=255)
                                                                                                                                     cAACTGTCCCCCCACCGCTGGTAGGTTCGGCCTTATCCGACTGGCTGATAACAAACAGCACCTTATGCCGGTATACCTCCCCAATCCCCTGATGATAAAAATGTTCCCCCACCGTCAGCGCCCGGTCATCGGCCTTAATc   <  2:14905/140‑1 (MQ=25)
                                                                                                                                      aaCTGTCCCCCCACCGCTGGTAGGTTCGGCCTTATCCGACTGGCTGATAACAAACAGCACCTTATGCCGGTATACCTCTCCAATCACCTGATGATAAAAATGTTCATCCACCGTCAGCGCCCGGTCATCGGCCTTAATCa  <  1:627867/140‑1 (MQ=255)
                                                                                                                                            |                                                                                                                                      
TAATCATCCGCTCGGCCATCACCTTAAGCCCCCATTGTAGACGGACCGACACAGCACACACCGGATGGACGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGTCCCCCACCGCTGGTAGGTTCGGCCTTATCCGACTGGCTGATAACAAACAGCACCTTATGCCGGTATACCTCTCCAATCACCTGATGATAAAAATGTTCATCCACCGTCAGCGCCCGGTCATCGGCCTTAATCA  >  NC_000913/268327‑268599

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