Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 273,814 T→C I323T (ATC→ACC)  ykfC → CP4‑6 prophage; PF00078 domain‑containing protein YkfC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913273,8140TC98.5% 225.7 / ‑3.4 67coding (968/1146 nt)ykfCCP4‑6 prophage; PF00078 domain‑containing protein YkfC
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (33/33);  minor base . (1/0);  total (34/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.77e-01

GTACTGCCGCTATGCCGATGATTTTGTCCTCATCGTCAAAGGCACCAAAGCACAGGTGGAAGCCATCAGGGAGGAGTGTCGGGGTGTGCTCGAAGGCAGTCTGAAACTCAGGCTGAACATGGATAAGACTAAAATCCCCCATGTTAATGACGGCTTTATCTTTCTGGGGCACAGGCTCATTCGCAAACGCAGTCGTTATGGCGAGATGCGAGTGGTCTCAACGATCCCGCAGGAGAAAGCCAGAAACTTCGCCGCATCGCTGACAGCAC  >  NC_000913/273680‑273948
                                                                                                                                      |                                                                                                                                      
gTACTGCCGCTATGCCGATGATTTTGTCCTCATCGTCAAAGGCACCAAAGCACAGGTGGAAGCCATCAGGGAGGAGTGTCGGGGTGTGCTCGAAGGCAGTCTGAAACTCAGGCTGAACATGGATAAGACTAAAAcccccc                                                                                                                                   <  1:1056133/140‑1 (MQ=255)
gTACTGCCGCTATGCCGATGATTTTGTCCTCATCGTCAAAGGCACCAAAGCACAGGTGGAAGCCATCAGGGAGGAGTGTCGGGGTGTGCTCGAAGGCAGTCTGAAACTCAGGCTGAACATGGATAAGACTAAAAcccccc                                                                                                                                   <  2:3506/140‑1 (MQ=255)
gTACTGCCGCTATGCCGATGATTTTGTCCTCATCGTCAAAGGCACCAAAGCACAGGTGGAAGCCATCAGGGAGGAGTGTCGGGGTGTGCTCGAAGGCAGTCTGAAACTCAGGCTGAACATGGATAAGACTAAAAcccccc                                                                                                                                   <  1:882153/140‑1 (MQ=255)
   cTGCCGCTATGCCGATGATTTTGTCCTCATCGTCAAAGGCACCAAAGCACAGGTGGAAGCCATCAGGGAGGAGTGTCGGGGTGTGCTCGAAGGCAGTCTGAAACTCAGGCTGAACATGGATAAGACTAAAACCCCCCATg                                                                                                                                <  2:710405/140‑1 (MQ=255)
        gCTATGCCGATGATTTTGTCCTCATCGTCAAAGGCACCAAAGCACAGGTGGAAGCCATCAGGGAGGAGTGTCGGGGTGTGCTCGAAGGCAGTCTGAAACTCAGGCTGAACATGGATAAGACTAAAACCCCCCATGTCAAt                                                                                                                           >  1:259405/1‑140 (MQ=255)
          tATGCCGATGATTTTGTCCTCATCGTCAAAGGCACCAAAGCACAGGTGGAAGCCATCAGGGAGGAGTGTCGGGGTGTGCTCGAAGGCAGTCTGAAACTCAGGCTGAACATGGATAAGACTAAAACCCCCCATGTTAATGa                                                                                                                         <  1:1351254/140‑1 (MQ=255)
           aTGCCGATGATTTTGTCCTCATCGTCAAAGGCACCAAAGCACAGGTTGAAGCCATCAGGGAGGAGTGTCGGGGTGTGCTCGAAGGCAGTCTGAAACTCAGGCTGAACATGGATAAGACTAAAACCCCCCATGTTAATGAc                                                                                                                        >  2:609417/1‑140 (MQ=255)
           aTGCCGATGATTTTGTCCTCATCGTCAAAGGCACCAAAGCACAGGTGGAAGCCATCAGGGAGGAGTGTCGGGGTGTGCTCGAAGGCAGTCTGAAACTCAGGCTGAACATGGATAAGACTAAAACCCCCCATGTTAATGAc                                                                                                                        <  1:482014/140‑1 (MQ=255)
                gatgatTTTGTCCTCATCGTCAAAGGCACCAAAGCACAGGTGGAAGCCATCAGGGAGGAGTGTCGGGGTGTGCTCGAAGGCAGTCTGAAACTCAGGCTGAACATGGATAAGACTAAAACCCCCCATGTTAATGACGGCtt                                                                                                                   >  2:1284505/1‑140 (MQ=255)
                  tgatTTTGTCCTCATCGTCAAAGGCACCAAAGCACAGGTGGAAGCCATCAGGGAGGAGTGTCGGGGTGTGCTCGAAGGCAGTCTGAAACTCAGGCTGAACATGGATAAGACTAAAACCCCCCATGTTAATGACGGCTTTa                                                                                                                 <  2:414513/140‑1 (MQ=255)
                     ttttGTCCTCATCGTCAAAGGCACCAAAGCACAGGTGGAAGCCATCAGGGAGGAGTGTCGGGGTGTGCTCGAAGGCAGTCTGAAACTCAGGCTGAACATGGATAAGACTAAAACCCCCCATGTTAATGACGGCTTTATCt                                                                                                              <  1:541670/140‑1 (MQ=255)
                           ccTCATCGTCAAAGGCACCAAAGCACAGGTGGAAGCCATCAGGGAGGAGTGTCGGGGTGTGCTCGAAGGCAGTCTGAAACTCAGGCTGAACATGGATAAGACTAAAACCCCCCATGTTAATGACGGCTTTATCTTTCTgg                                                                                                        >  2:349267/1‑140 (MQ=255)
                               aTCGTCAAAGGCACCAAAGCACAGGTGGAAGCCATCAGGGAGGAGTGTCGGGGTGTGCTCGAAGGCAGTCTGAAACTCAGGCTGAACATGGATAAGACTAAAACCCCCCATGTTAATGACGGCTTTATCTTTCTGGGGca                                                                                                    <  2:1372067/140‑1 (MQ=255)
                                      aaGGCACCAAAGCACAGGTGGAAGCCATCAGGGAGGAGTGTCGGGGTGTGCTCGAAGGCAGTCTGAAACTCAGGCTGAACATGGATAAGACTAAAACCCCCCATGTTAATGACGGCTTTATCTTTCTGGGGCACAGGCTc                                                                                             >  1:455995/1‑140 (MQ=255)
                                         gCACCAAAGCACAGGTGGAAGCCATCAGGGAGGAGTGTCGGGGTGTGCTCGAAGGCAGTCTGAAACTCAGGCTGAACATGGATAAGACTAAAACCCCCCATGTTAATGACGGCTTTATCTTTCTGGGGCACAGGCTCAtt                                                                                          >  1:1052354/1‑140 (MQ=255)
                                           aCCAAAGCACAGGTGGAAGCCATCAGGGAGGAGTGTCGGGGTGTGCTCGAAGGCAGTCTGAAACTCAGGCTGAACATGGATAAGACTAAAACCCCCCATGTTAATGACGGCTTTATCTTTCTGGGGca                                                                                                    >  1:1167966/1‑128 (MQ=255)
                                           aCCAAAGCACAGGTGGAAGCCATCAGGGAGGAGTGTCGGGGTGTGCTCGAAGGCAGTCTGAAACTCAGGCTGAACATGGATAAGACTAAAACCCCCCATGTTAATGACGGCTTTATCTTTCTGGGGca                                                                                                    <  2:1167966/128‑1 (MQ=255)
                                           aCCAAAGCACAGGTGGAAGCCATCAGGGAGGAGTGTCGGGGTGTGCTCGAAGGCAGTCTGAAACTCAGGCTGAACATGGATAAGACTAAAACCCCCCATGTTAATGACGGCTTTATCTTTCTGGGGCACAGGCTCATTCg                                                                                        >  1:641355/1‑140 (MQ=255)
                                                aGCACAGGTGGAAGCCATCAGGGAGGAGTGTCGGGGTGTGCTCGAAGGCAGTCTGAAACTCAGGCTGAACATGGATAAGACTAAAACCCCCCATGTTAATGACGGCTTTATCTTTCTGGGGCACAGGCTCATTCGCAAAc                                                                                   >  1:825268/1‑140 (MQ=255)
                                                     aGGTGGAAGCCATCAGGGAGGAGTGTCGGGGTGTGCTCGAAGGCAGTCTGAAACTCAGGCTGAACATGGATAAGACTAAAACCCCCCATGTTAATGACGGCTTTATCTTTCTGGGGCACAGGCTCATTCGCAAACGCAGt                                                                              >  2:11432/1‑140 (MQ=255)
                                                     aGGTGGAAGCCATCAGGGAGGAGTGTCGGGGTGTGCTCGAAGGCAGTCTGAAACTCAGGCTGAACATGGATAAGACTAAAACCCCCCATGTTAATGACGGCTTTATCTTTCTGGGGCACAGGCTCATTCGCAAACGCAGt                                                                              >  2:1080130/1‑140 (MQ=255)
                                                     aGGTGGAAGCCATCAGGGAGGAGTGTCGGGGTGTGCTCGAAGGCAGTCTGAAACTCAGGCTGAACATGGATAAGACTAAAACCCCCCATGTTAATGACGGCTTTATCTTTCTGGGGCACAGGCTCATTCGCAAACGCAGt                                                                              >  1:495780/1‑140 (MQ=255)
                                                         ggAAGCCATCAGGGAGGAGTGTCGGGGTGTGCTCGAAGGCAGTCTGAAACTCAGGCTGAACATGGATAAGACTAAAACCCCCCATGTTAATGACGGCTTTATCTTTCTGGGGCACAGGCTCATTCGCAAACGCAGTCGtt                                                                          <  1:1284505/140‑1 (MQ=255)
                                                                 tCAGGGAGGAGTGTCGGGGTGTGCTCGAAGGCAGTCTGAAACTCAGGCTGAACATGGATAAGACTAAAACCCCCCATGTTAATGACGGCTTTATCTTTCTGGGGCACAGGCTCATTCGCAAACGCAGTCGTTATGGCgag                                                                  >  1:1007311/1‑140 (MQ=255)
                                                                 tCAGGGAGGAGTGTCGGGGTGTGCTCGAAGGCAGTCTGAAACTCAGGCTGAACATGGATAAGACTAAAACCCCCCATGTTAATGACGGCTTTATCTTTCTGGGGCACAGGCTCATTCGCAAACGCAGTCGTTATGGCgag                                                                  >  2:1295557/1‑140 (MQ=255)
                                                                 tCAGGGAGGAGTGTCGGGGTGTGCTCGAAGGCAGTCTGAAACTCAGGCTGAACATGGATAAGACTAAAACCCCCCATGTTAATGACGGCTTTATCTTTCTGGGGCACAGGCTCATTCGCAAACGCAGTCGTTATGGCgag                                                                  >  1:746710/1‑140 (MQ=255)
                                                                 tCAGGGAGGAGTGTCGGGGTGTGCTCGAAGGCAGTCTGAAACTCAGGCTGAACATGGATAAGACTAAAACCCCCATGTTAATGACGGCTTTATCTTTCTGGGGCACAGGCTCATTCGCAAACGCAGTCGTTATGGCgaga                                                                 >  2:221805/1‑140 (MQ=255)
                                                                       aggagTGTCGGGGTGTGCTCGAAGGCAGTCTGAAACTCAGGCTGAACATGGATAAGACTAAAACCCCCCATGTTAATGACGGCTTTATCTTTCTGGGGCACAGGCTCATTCGCAAACGCAGTCGTTATGGCGAGATGCGa                                                            <  2:154601/140‑1 (MQ=255)
                                                                       aggagTGTCGGGGTGTGCTCGAAGGCAGTCTGAAACTCAGGCTGAACATGGATAAGACTAAAACCCCCCATGTTAATGACGGCTTTATCTTTCTGGGGCACAGGCTCATTCGCAAACGCAGTCGTTATGGCGAGATGCGa                                                            <  2:455995/140‑1 (MQ=255)
                                                                       aggagTGTCGGGGTGTGCTCGAAGGCAGTCTGAAACTCAGGCTGAACATGGATAAGACTAAAACCCCCCATGTTAATGACGGCTTTATCTTTCTGGGGCACAGGCTCATTCGCAAACGCAGTCGTTATGGCGAGATGCGa                                                            <  2:746710/140‑1 (MQ=255)
                                                                        ggagTGTCGGGGTGTGCTCGAAGGCAGTCTGAAACTCAGGCTGAACATGGATAAGACTAAAACCCCCCATGTTAATGACGGCTTTATCTTTCTGGGGCACAGGCTCATTCGCAAACGCAGTCGTTATGGCGAGATGCGAg                                                           <  2:884005/140‑1 (MQ=255)
                                                                        ggagTGTCGGGGTGTGCTCGAAGGCAGTCTGAAACTCAGGCTGAACATGGATAAGACTAAAACCCCCCATGTTAATGACGGCTTTATCTTTCTGGGGCACAGGCTCATTCGCAAACGCAGTCGTTATGGCGAGATGCGAg                                                           <  1:1295557/140‑1 (MQ=255)
                                                                            tgtCGGGGTGTGCTCGAAGGCAGTCTGAAACTCAGGCTGAACATGGATAAGACTAAAACCCCCCATGTTAATGACGGCTTTATCTTTCTGGGGCACAGGCTCATTCGCAAACGCAGTCGTTATGGCGAGATGCGAGTGGt                                                       >  1:699104/1‑140 (MQ=255)
                                                                            tgtCGGGGTGTGCTCGAAGGCAGTCTGAAACTCAGGCTGAACATGGATAAGACTAAAACCCCCCATGTTAATGACGGCTTTATCTTTCTGGGGCACAGGCTCATTCGCAAACGCAGTCGTTATGGCGAGATGCGAGTGGt                                                       <  2:360242/140‑1 (MQ=255)
                                                                            tgtCGGGGTGTACTCGAAGGCAGTCTGAAACTCAGGCTGAACATGGATAAGACTAAAACCCCCCATGTTAATGACGGCTTTATCTTTCTGGGGCACAGGCTCATTCGCAAACGCAGTCGTTATGGCGAGATGCGAGTGGt                                                       >  2:1162529/1‑140 (MQ=255)
                                                                             gtCGGGGTGTGCTCGAAGGCAGTCTGAAACTCAGGCTGAACATGGATAAGACTAAAACCCCCCATGTTAATGACGGCTTTATCTTTCTGGGGCACAGGCTCATTCGCAAACGCAGTCGTTATGGCGAGATGCGAGTGGtc                                                      >  1:620699/1‑140 (MQ=255)
                                                                               cGGGGTGTGCTCGAAGGCAGTCTGAAACTCAGGCTGAACATGGATAAGACTAAACCCCCCCATGTTAATGACGGTTTTATCTTTCTGGGGCACAGGCTCATTCGCAAACGCAGTCGTTATGGCGAGATGCGAGTGGtctc                                                    <  2:699104/140‑1 (MQ=255)
                                                                               cGGGGTGTGCTCGAAGGCAGTCTGAAACTCAGGCTGAACATGGATAAGACTAAAACCCCCCATGTTAATGACGGCTTTATCTTTCTGGGGCACAGGCTCATTCGCAAACGCAGTCGTTATGGCGAGATGCGa                                                            <  1:1413025/132‑1 (MQ=255)
                                                                               cGGGGTGTGCTCGAAGGCAGTCTGAAACTCAGGCTGAACATGGATAAGACTAAAACCCCCCATGTTAATGACGGCTTTATCTTTCTGGGGCACAGGCTCATTCGCAAACGCAGTCGTTATGGCGAGATGCGa                                                            >  2:1413025/1‑132 (MQ=255)
                                                                               cGGGGTGTGCTCGAAGGCAGTCTGAAACTCAGGCTGAACATGGATAAGACTAAAACCCCCCATGTTAATGACGGCTTTATCTTTCTGGGGCACAGGCTCATTCGCAAACGCAGTCGTTATGGCGAGATGCGAg                                                           <  2:1124056/133‑1 (MQ=255)
                                                                               cGGGGTGTGCTCGAAGGCAGTCTGAAACTCAGGCTGAACATGGATAAGACTAAAACCCCCCATGTTAATGACGGCTTTATCTTTCTGGGGCACAGGCTCATTCGCAAACGCAGTCGTTATGGCGAGATGCGAg                                                           >  1:1124056/1‑133 (MQ=255)
                                                                               cGGGGTGTGCTCGAAGGCAGTCTGAAACTCAGGCTGAACATGGATAAGACTAAAACCCCCCATGTTAATGACGGCTTTATCTTTCTGGGGCACAGGCTCATTCGCAAACGCAGTCGTTATGGCGAGATGCGAGTGGtctc                                                    >  2:787946/1‑140 (MQ=255)
                                                                                      tgCTCGAAGGCAGTCTGGAACTCAGGCTGAACATGGATAAGACTAAAACCCCCCATGTTAATGACGGCTTTATCTTTCTGGGGCACAGGCTCATTCGCAAACGCAGTCGTTATGGCGAGATGCGAGTGGTCTCAACGATc                                             >  1:479299/1‑140 (MQ=255)
                                                                                      tgCTCGAAGGCAGTCTGAAACTCAGGCTGAACATGGATAAGACTAAAACCCCCCATGTTAATGACGGCTTTATCTTTCTGGGGCACAGGCTCATTCGCAAACGCAGTCGTTATGGCGAGATGCGAGTGGTCTCAACGATc                                             >  1:959068/1‑140 (MQ=255)
                                                                                      tgCTCGAAGGCAGTCTGAAACTCAGGCTGAACATGGATAAGACTAAAACCCCCCATGTTAATGACGGCTTTATCTTTCTGGGGCACAGGCTCATTCGCAAACGCAGTCGTTATGGCGAGATGCGAGTGGTCTCAACGATc                                             >  1:178094/1‑140 (MQ=255)
                                                                                        cTCGAAGGCAGTCTGAAACTCAGGCTGAACATGGATAAGACTAAAACCCCCCATGTTAATGACGGCTTTATCTTTCTGGGGCACAGGCTCATTCGCAAACGCAGTCGTTATGGCGAGATGCGAg                                                           <  2:1379841/124‑1 (MQ=255)
                                                                                        cTCGAAGGCAGTCTGAAACTCAGGCTGAACATGGATAAGACTAAAACCCCCCATGTTAATGACGGCTTTATCTTTCTGGGGCACAGGCTCATTCGCAAACGCAGTCGTTATGGCGAGATGCGAg                                                           >  1:1379841/1‑124 (MQ=255)
                                                                                              ggCAGTCTGAAACTCAGGCTGAACATGGATAAGACTAAAACCCCCCATGTTAATGACGGCTTTATCTTTCTGGGGCACAGGCTCATTCGCAAACGCAGTCGTTATGGCGAGATGCGAGTGGTCTCAACGATCCCGCAGga                                     <  1:609417/140‑1 (MQ=255)
                                                                                                cAGTCTGGAACTCAGGCTGAACATGGATAAGACTAAAACCCCCCATGTTAATGACGGCTTTATCTTTCTGGGGCACAGGCTCATTCGCAAACGCAGTCGTTATGGCGAGATGCGAGTGGTCTCAACGATCCCGCAGgaga                                   <  2:479299/140‑1 (MQ=255)
                                                                                                cAGTCTGAAACTCAGGCTGAACATGGATAAGACTAAAACCCCCCATGTTAATGACGGCTTTATCTTTCTGGGGCACAGGCTCATTCGCAAACGCAGTCGTTATGGCGAGATGCGAGTGGTCTCAACGATCCCGCAGgaga                                   <  2:199536/140‑1 (MQ=255)
                                                                                                   tCTGAAACTCAGGCTGAACATGGATAAGACTAAAACCCCCCATGTTAATGACGGCTTTATCTTTCTGGGGCACAGGCTCATTCGCAAACGCAGTCGTTATGGCGAGATGCGAGTGGTCTCAACGATCCCGCAGGAGAAAg                                >  2:353813/1‑140 (MQ=255)
                                                                                                         aCTCAGGCTGAACATGGATAAGACTAAAACCCCCCATGTTAATGACGGCTTTATCTTTCTGGGGCACAGGCTCATTCGCAAACGCAGTCGTTATGGCGAGATGCGAGTGGTCTCAACGATCCCGCAGGAGAAAGCCAGaa                          >  1:342372/1‑140 (MQ=255)
                                                                                                             aGGCTGAACATGGATAAGACTAAAACCCCCCATGTTAATGACGGCTTTATCTTTCTGGGGCACAGGCTCATTCGCAAACGCAGTCGTTATGGCGAGATGCGAGTGGTCTCAACGATCCCGCAGGAGAAAGCCAGAAACtt                      >  1:393546/1‑140 (MQ=255)
                                                                                                                 tGAACATGGATAAGACTAAAACCCCCCATGTTAATGACGGCTTTATCTTTCTGGGGCACAGGCTCATTCGCAAACGCAGTCGTTATGGCGAGATGCGAGTGGTCTCAACGATCCCGCAGGAGAAAGCCAGAAACTTcgcc                  <  1:353813/140‑1 (MQ=255)
                                                                                                                    aCATGGATAAGACTAAAACCCCCCATGTTAATGACGGCTTTATCTTTCTGGGGCACAGGCTCATTCGCAAACGCAGTCGTTATGGCGAGATGCGAGTGGTCTCAACGATCCCGCAGGAGAAAGCCAGAAACTTCGCCGCa               >  1:170109/1‑140 (MQ=255)
                                                                                                                    aCATGGATAAGACTAAAACCCCCCATGTTAATGACGGCTTTATCTTTCTGGGGCACAGGCTCATTCGCAAACGCAGTCGTTATGGCGAGATGCGAGTGGTCTCAACGATCCCGCAGGAGAAAGCCAGAAACTTCGCCGCa               <  1:349267/140‑1 (MQ=255)
                                                                                                                     cATGGATAAGACTAAAACCCCCCATGTTAATGACGGCTTTATCTTTCTGGGGGACAGGCTCATTCGCAAACGCAGTCGTTATGGCGAGATGCGAGTGGTCTCAACGATCCCGCAGGAGAAAGCCAGAAACTTCGCCGCAt              >  1:1268214/1‑140 (MQ=255)
                                                                                                                       tGGATAAGACTAAAACCCCCCATGTTAATGACGGCTTTATCTTTCTGGGGCACAGGCTCATTCGCAAACGCAGTCGTTATGGCGAGATGCGAGTGGTCTCAACGATCCCGCAGGAGAAAGCCAGAAACTTCGCCGCATCg            >  2:189863/1‑140 (MQ=255)
                                                                                                                         gATAAGACTAAAACCCCCCATGTTAATGACGGCTTTATCTTTCTGGGGCACAGGCTCATTCGCAAACGCAGTCGTTATGGCGAGATGCGAGTGGTCTCAACGATCCCGCAGGAGAAAGCCAGAAACTTCGCCGCATCGCt          <  2:170109/140‑1 (MQ=255)
                                                                                                                         gATAAGACTAAAACCCCCCATGTTAATGACGGCTTTATCTTTCTGGGGCACAGGCTCATTCGCAAACGCAGTCGTTATGGCGAGATGCGAGTGGTCTCAACGATCCCGCAGGAGAAAGCCAGAAACTTCGCCGCATCGCt          <  2:620699/140‑1 (MQ=255)
                                                                                                                             aGACTAAAACCCCCCATGTTAATGACGGCTTTATCTTTCTGGGGCACAGGCTCATTCGCAAACGCAGTCGTTATGGCGAGTTGCGAGTGGTCTCAACGATCCCGCAGGAGAAAGCCAGAAACTTCGCCGCATCGCTGa        >  2:1097309/1‑138 (MQ=255)
                                                                                                                             aGACTAAAACCCCCCATGTTAATGACGGCTTTATCTTTCTGGGGCACAGGCTCATTCGCAAACGCAGTCGTTATGGCGAGTTGCGAGTGGTCTCAACGATCCCGCAGGAGAAAGCCAGAAACTTCGCCGCATCGCTGa        <  1:1097309/138‑1 (MQ=255)
                                                                                                                              gACTAAAACCCCCCATGTTAATGACGGCTTTATCTTTCTGGGGCACAGGCTCATTCGCAAACGCAGTCGTTATGGCGAGATGCGAGTGGTCTCAACGATCCCGCAGGAGAAAGCCAGAAACTTCGCCGCATCGCTGACAg     <  2:825268/140‑1 (MQ=255)
                                                                                                                               aCTAAAACCCCCCATGTTAATGACGGCTTTATCTTTCTGGGGCACAGGCTCATTCGCAAACGCAGTCGTTATGGCGAGATGCGAGTGGTCTCAACGATCCCGCAGGAGAAAGCCAGAAACTTCGCCGCATCGCTGACAGc    <  1:899362/140‑1 (MQ=255)
                                                                                                                                ccaaaaCCCCCCATGTTAATGACGGCTTTATCTTTCTGGGGCACAGGCTCATTCGCAAACGCAGTCGTTATGGCGAGATGCGAGTGGTCTCAACGATCCCGCAGGAGAAAGCCAGAAACTTCGCCGCATCGCTGACAGCa   <  1:1463457/138‑1 (MQ=255)
                                                                                                                                cTAAAACCCCCCATGTTAATGACGGCTTTATCTTTCTGGGGCACAGGCTCATTCGCAAACGCAGTCGTTATGGCGAGATGCGAGTGGTCTCAACGATCCCGCAGGAGAAAGCCAGAAACTTCGCCGCATCGCTGACAGCa   <  2:1052354/140‑1 (MQ=255)
                                                                                                                                cTAAAACCCCCCATGTTAATGACGGCTTTATCTTTCTGGGGCACAGGCTCATTCGCAAACGCAGTCGTTATGGCGAGATGCGAGTGGTCTCAACGATCCCGCAGGAGAAAGCCAGAAACTTCGCCGCATCGCTGACAGCa   <  2:641355/140‑1 (MQ=255)
                                                                                                                                cTAAAACCCCCCATGTCAATGACGGCTTTATCTTTCTGGGGCACAGGCTCATTCGCAAACGCAGTCGTTATGGCGAGATGCGAGTGGTCTCAACGATCCCGCAGGAGAAAGCCAGAAACTTCGCCGCATCGCTGACAGCa   <  2:259405/140‑1 (MQ=255)
                                                                                                                                 tAAAACCCCCCATGTTAATGACGGCTTTATCTTTCTGGGGCACAGGCTCATTCGCAAACGCAGTCGTTATGGCGAGATGCGAGTGGTCTCAACGATCCCGCAGGAGAAAGCCAGAAACTTCGCCGCATCGCTGACAGCAc  <  2:959068/140‑1 (MQ=255)
                                                                                                                                  aaaaCCCCCCATGTTAATGACGGCTTTATCTTTCTGGGGCACAGGCTCATTCGCAAACGCAGTCGTTATGGCGAGATGCGAGTGGTCTCAACGATCCCGCAGGAGAAAGCCAGAAACTTCGCCGCATCGCTGa        <  2:1397820/133‑1 (MQ=255)
                                                                                                                                  aaaaCCCCCCATGTTAATGACGGCTTTATCTTTCTGGGGCACAGGCTCATTCGCAAACGCAGTCGTTATGGCGAGATGCGAGTGGTCTCAACGATCCCGCAGGAGAAAGCCAGAAACTTCGCCGCATCGCTGa        >  1:1397820/1‑133 (MQ=255)
                                                                                                                                      |                                                                                                                                      
GTACTGCCGCTATGCCGATGATTTTGTCCTCATCGTCAAAGGCACCAAAGCACAGGTGGAAGCCATCAGGGAGGAGTGTCGGGGTGTGCTCGAAGGCAGTCTGAAACTCAGGCTGAACATGGATAAGACTAAAATCCCCCATGTTAATGACGGCTTTATCTTTCTGGGGCACAGGCTCATTCGCAAACGCAGTCGTTATGGCGAGATGCGAGTGGTCTCAACGATCCCGCAGGAGAAAGCCAGAAACTTCGCCGCATCGCTGACAGCAC  >  NC_000913/273680‑273948

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: