Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,439,845 T→A 6.5% intergenic (+52/+27) ytfK → / ← ytfL DUF1107 family protein/UPF0053 family inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,439,8450TA 6.5% 390.0 / 11.3 142intergenic (+52/+27)ytfK/ytfLDUF1107 family protein/UPF0053 family inner membrane protein
Reads supporting (aligned to +/- strand):  ref base T (71/61);  new base A (7/3);  total (78/64)
Fisher's exact test for biased strand distribution p-value = 5.12e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

Only 100 of 143 total aligned reads displayed.

TAAGATTCTTATCCCAAAAGTGAAGGATAAGCTGCATTTGTCAGTGATGTCCGAAGTTAACCGTCAGGTTATGCGTCTGCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGTTTTGGCGAAAGGGCGGTGGCCTTGCTGTCGATCCGGGTCACCAGCAGCTGGTCGATGCGGTAGTTATCGATATCCACC  >  NC_000913/4439697‑4439989
                                                                                                                                                    |                                                                                                                                                
      tCTTATCCCAAAAGTGAAGGATAAGCTGCATTTGTCAGTGATGTCCGAAGTTAACCGTCAGGTTATGCGTCTGCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAg                                                                                                                                           <  1:1213318/150‑1 (MQ=255)
         tATCCCAAAAGTGAAGGATAAGCTGCATTTGTCAGTGATGTCCGAAGTTAACCGTCAGGTTATGCGTCTGCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCg                                                                                                                                        <  1:378811/150‑1 (MQ=255)
          aTCACAAAAGTGAAGTATAAGCTGCATTTGTCAGTGATGTCCGAAGTTAACCGTCAGGTTATGCGTCTGCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGAATCAGGAGCCGt                                                                                                                                       <  2:201386/150‑1 (MQ=255)
            cccAAAAGTGAAGGATAAGCTGCATTTGTCAGTGATGTCCGAAGTTAACCGTCAGGTTATGCGTCTGCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGAATCAGGAGCCGtct                                                                                                                                     <  1:1855104/150‑3 (MQ=255)
             ccAAAAGTGAAGGATAAGCTGCATTTGTCAGTGATGTCCGAAGTTAACCGTCAGGTTATGCGTCTGCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGTTGa                                                                                                                                    <  2:858952/150‑1 (MQ=255)
             ccAAAAGTGAAGGATAAGCTGCATTTGTCAGTGATGTACGAAGTTAACCGTCAGGTTATGCGTCTGCAAACAGAGATGGCTTAACAAAAGTACTATGCAGTAAAAAGTGATATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGTTGa                                                                                                                                    >  1:305827/1‑150 (MQ=255)
              cAAAAGTGAAGGATAAGCTGCATTTGTCAGTGATGTCCGAAGTTAACCGTCAGGTTATGCGTCTGCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGTTGAt                                                                                                                                   >  2:1361896/1‑150 (MQ=255)
               aaaaGTGAAGGATAAGCTGCATTTGTCAGTGATGTCCGAAGTTAACCGTCAGGTTATGCGTCTGCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGTTGATg                                                                                                                                  <  1:1463906/150‑1 (MQ=255)
               aaaaGTGAAGGATAAGCTGCATTTGTCAGTGATGTCCGAAGTTAACCGTCAGGTTATGCGTCTGCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGAATCAGGAGCCGtcttat                                                                                                                                  <  2:1291538/150‑6 (MQ=255)
                   gTGAAGGATAAGCTGCATTTGTCAGTGATGTCCGAAGTTAACCGTCAGGTTATGCGTCTGCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTc                                                                                                                              >  2:1843919/1‑150 (MQ=255)
                   gTGAAGGATAAGCTGCATTTGTCAGTGATGTCCGAAGTTAACCGTCAGGTTATGCGTCTGCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTc                                                                                                                              >  2:1216824/1‑150 (MQ=255)
                   ctGAAGGATAATCTGCATTTGTCAGGGATGTCCGAAGTTAACCGACAGGTTATGCGTCTGCACACAGAGATGGCTAAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTc                                                                                                                              <  2:85576/149‑1 (MQ=255)
                     gAAGGATAAGCTGCATTTGTCAGTGATGTCCGAAGTTAACCGTCAGGTTATGCGTCTGCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTg                                                                                                                            >  1:1452115/1‑150 (MQ=255)
                       aGGATAAGCTGCATTTGTCAGTGATGTCCGAAGTTAACCGTCAGGTTATGCGTCTGCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTggg                                                                                                                          >  2:410832/1‑150 (MQ=255)
                       aGGATAAGCTGCATTTGTCAGTGATGTCCGAAGTTAACCGTCAGGTTATGCGTCTGCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTggg                                                                                                                          <  1:1727528/150‑1 (MQ=255)
                        tgttAAGCTGCATTTGTCAGTGATGTCAGAATTTAACCGTCAGGTTATGCGTCTGCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAACAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTgggg                                                                                                                         <  2:1258917/147‑1 (MQ=255)
                        ggATAAGCTGCATTTGTCAGTGATGTCCGAAGTTAACCGTCAGGTTATGCGTCTGCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTgggg                                                                                                                         >  2:244625/1‑150 (MQ=255)
                          aTAAGCTGCATTTGTCAGTGATGTCCGAAGTTAACCGTCAGGTTATGCGTCTGCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGt                                                                                                                       <  2:465949/150‑1 (MQ=255)
                           tAAGCTGCATTTGTCAGTGATGTCCGAAGTTAACCGTCAGGTTATGCGTCTGCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGtt                                                                                                                      <  1:1485857/150‑1 (MQ=255)
                            aaGCTGCATTTGTCAGTGATGTCCGAAGTTAACCGTCAGGTTATGCGTCTGCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTa                                                                                                                     <  1:66967/150‑1 (MQ=255)
                            aaGCTGCATTTGTCAGTGATGTCCGAAGTTAACCGTCAGGTTATGCGTCTGCAAACAGAGATGGCTTAACAAAAGTGCTATGCAGTAAAAAGTCATATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGTTGAAGTTTCTGGGGGTTa                                                                                                                     >  1:98242/1‑150 (MQ=255)
                             aGCTGCATTTGTCAGTGATGTCCGAAGTTAACCGTCAGGTTATGCGTCTGCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTAc                                                                                                                    >  2:1472759/1‑150 (MQ=255)
                             aGCTGCATTTGTCAGTGATGTCCGAAGTTAACCGTCAGGTTATGCGTCTGCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTAc                                                                                                                    >  1:1754924/1‑150 (MQ=255)
                                  cATTTGTCAGTGATGTCCGAAGTTAACCGTCAGGTTATGCGTCTGCAAACAGAGATGGGTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCATGATTCAGGAGCCGTTGATGTTTCTGGGGGTTActccac                                                                                                               >  1:1198514/1‑145 (MQ=255)
                                  cATTTGTCAGTGATGTCCGAAGTTAACCGTCAGGTTATGCGTCTGCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGAc                                                                                                               >  2:471686/1‑150 (MQ=255)
                                      tGTCAGTGATGTCCGAAGTTAACCGTCAGGTTATGCGTCTGCAAACAGAGATGGCTTACCCAAAGTGCTATGCAGTAAAAAGTCCTATGTAGTAATAAGACGGCTCATGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCtt                                                                                                           >  1:71524/1‑150 (MQ=255)
                                          aGTGATGTCCGAAGTTAACCGTCAGGTTATGCGTCTGCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGTTGATCTTTCTGGGGGTTACGCGCCGCTttctt                                                                                                       >  2:900408/1‑150 (MQ=255)
                                           gggATGTCCGAAGTTAACCGTCAGGTTATGCGTCTGCAAACAGAGATGGCGTAACCAAAGTGCTATGCAGTAAAAATTGCTATTAATTATTAATACGGCTCCTGATTCAGGAGACGTTTATTTTTATTGTGGTTACTCGACGCTttcttc                                                                                                      <  2:1062812/148‑1 (MQ=255)
                                           gTGATGTCCGAAGTTAACCGTCAGGTTATGCGTCTGCAAACAGAGATGGCGTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTttcttc                                                                                                      <  1:90108/150‑1 (MQ=255)
                                             gATGTCCGAAGTTAACCGTCAGGTTATGCGTCTGCAAACAGAGATGGCTTAACCAAAGTGCTAAGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTcttctt                                                                                                    <  1:1073723/150‑1 (MQ=255)
                                               tGTCCGAAGTTAACCGTCAGGTTATGCGTCTGCAAACAGAGATGGCTTACCCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTtcttta                                                                                                  <  1:397943/150‑1 (MQ=255)
                                               tGTCCGAAGTTAACCGTCAGGTTATGCGGCTGCAAACAGAGATCGCTTAACCACAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTGAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTtcttta                                                                                                  <  2:144135/150‑1 (MQ=255)
                                                 tCCGAAGTTAACCGTCAGGTTATGCGTCTGCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTtctttatc                                                                                                >  2:691371/1‑150 (MQ=255)
                                                  ccGAAGTTAACCGTCAGGTTATGCGTCTGCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTtctttatct                                                                                               >  1:1079279/1‑150 (MQ=255)
                                                  ccGAAGTTAACCGTCAGGTTATGCGTCTGCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTtctttatct                                                                                               <  2:510299/150‑1 (MQ=255)
                                                   cGAAGTTAACCGTCAGGTTATGCGTCTGCAAACAGAGTTGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTtctttatctt                                                                                              >  2:1879600/1‑150 (MQ=255)
                                                    gAAGTTAACCGTCAGGTTATGCGTCTGCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTtctttatcttt                                                                                             <  2:1022153/150‑1 (MQ=255)
                                                         tAACCGTCAGGTTATGCGTCTGCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGt                                                                                        >  1:335407/1‑150 (MQ=255)
                                                           aCCGTCAGGTTATGCGTCTGCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCa                                                                                      <  2:1412482/150‑1 (MQ=255)
                                                             cGTCAGGTTATGCGTCTGCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAgg                                                                                    >  2:1304431/1‑150 (MQ=255)
                                                             cGTCAGGTTATGCGTCTGCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAgg                                                                                    >  1:968382/1‑150 (MQ=255)
                                                               tCAGGTTATGCGTCTGCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCa                                                                                  <  1:649793/150‑1 (MQ=255)
                                                                 aGGTTATGCGTCTGCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATTAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGt                                                                                >  1:1385865/1‑150 (MQ=255)
                                                                  ggTTATGCGTCTGCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGtt                                                                               >  2:1029804/1‑150 (MQ=255)
                                                                     tATGCTTCTGCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAATAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGTGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGTTTTg                                                                            <  1:1584569/150‑1 (MQ=255)
                                                                       tGCGTCTGCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGTTTTGGc                                                                          >  2:319865/1‑150 (MQ=255)
                                                                        gCGTCTGCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGTTTTGGCg                                                                         <  1:588154/150‑1 (MQ=255)
                                                                            cTGCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGTTTTGGCGAAAg                                                                     <  2:1087093/150‑1 (MQ=255)
                                                                             tGCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGTTTTGGCGAAAgg                                                                    <  1:151229/150‑1 (MQ=255)
                                                                              gCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGTTTTGGCGAAAggg                                                                   >  2:540152/1‑150 (MQ=255)
                                                                              gCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTATTTATCTTTAGCGTCAGGCAGTTTTGGCGAAAggg                                                                   >  2:1746449/1‑150 (MQ=255)
                                                                               cAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGTTTTGGCGAAAGGGc                                                                  >  1:1578814/1‑150 (MQ=255)
                                                                               cAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGTTTTGGCGAAAGGGc                                                                  <  1:1469183/150‑1 (MQ=255)
                                                                                aaaCAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGTTTTGGCGAAAGGGCg                                                                 >  2:632551/1‑150 (MQ=255)
                                                                                aaaCAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAAAAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGTTTTGGCGAAAGGGCg                                                                 <  2:1665587/150‑1 (MQ=255)
                                                                                   cAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAATTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGTTTTGGCGAAAGGGCGGTg                                                              <  2:1173697/150‑1 (MQ=255)
                                                                                   cAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGTTTTGGCGAAAGGGCGGTg                                                              <  2:1170796/150‑1 (MQ=255)
                                                                                   cAGAGATGGCTTAACCAAAGTGCTATGATGTAAAAAGTGCTATGCAGTATTAAGACGGCTCCTGATTCAGGAGCAGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCTGTTTTGGCGAAAGGGCGGTg                                                              >  1:1015922/1‑150 (MQ=255)
                                                                                     gagaTGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGTTTTGGCGAAAGGGCGGTGGc                                                            >  2:1058264/1‑150 (MQ=255)
                                                                                     gagaTGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGTTTTGGCGAAAGGGCGGTGGc                                                            >  1:17025/1‑150 (MQ=255)
                                                                                     gagaTGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGTTTTGGCGAAAGGGCGGTGGc                                                            >  2:906862/1‑150 (MQ=255)
                                                                                           gCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGTTTTGGCGAAAGGGCGGTGGCCTTGCt                                                      >  1:509680/1‑150 (MQ=255)
                                                                                              tAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGTTTTGGCGAAAGGGCGGTGGCCTTGCTGTc                                                   <  1:207948/150‑1 (MQ=255)
                                                                                                 ccAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGTTTTGGCGAAAGGGCGGTGGCCTTGCTGTCGAt                                                >  2:1125588/1‑150 (MQ=255)
                                                                                                 ccAAAGTGCTATGAAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGTTTTGGCGAAAGGGCGGTGGCCTTGCTGTCGAt                                                >  2:539292/1‑150 (MQ=255)
                                                                                                      gtgctatgcagtaaaaagtgctatgcagtaaTAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGTTTTGGCGAAAGGGCGGTGGCCTTGCTGTCGATCCggg                                           <  2:35483/150‑1 (MQ=255)
                                                                                                        gctatgcagtaaaaagtgctatgcagtaaTAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCGTTATCTTTAGCGTCAGGCAGGTTTGGCGAAAGGGCGGTGGCCTTGCTGTCGATCCGGGTc                                         >  1:906872/1‑150 (MQ=255)
                                                                                                         ctatgcagtaaaaagtgctatgcagtaaTAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGTTTTGGCGAAAGGGCGGTGGCCTTGCTGTCGATCCGGGTCa                                        <  2:1557320/150‑1 (MQ=255)
                                                                                                          tatgcagtaaaaagtgctatgcagtaaTAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGTTTTGGCGAAAGGGCGGTGGCCTTGCTGTCGATCCGGGTCAc                                       <  1:1668812/150‑1 (MQ=255)
                                                                                                          tatgcagtaaaaagtgctatgcagtaaTAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGTTTTGGCGAAAGGGAGGTGGCCTTGCTGTCGATCCGGGTCAc                                       >  1:1108428/1‑150 (MQ=255)
                                                                                                          tatgcagtaaaaagtgctaagcagtaaTAAGACCGCTCCTTATTCAGGAGCCGTAGAAGTTTCTGGGGGTTACGCGACGCTATCTTATTTATATTAAGCGACACGCAGGTTAGGCGAAAGGGCGGTGCCCATCCAGTCGATCCCGGTCAc                                       >  2:795204/1‑150 (MQ=255)
                                                                                                               agtaaaaagtgctatgcagtaaTAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGTTTTGGCGAAAGGGCGGTGGCCTTGCTGTCGATCCGGGTCACcagca                                  >  1:1868808/1‑150 (MQ=255)
                                                                                                                gtaaaaagtgctatgcagtaaTAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGTTTTGGCGAAAGGGCGGTGGCCTTGCTGTCGATCCGGGTCACcagcag                                 <  1:1029804/150‑1 (MQ=255)
                                                                                                                        tgctatgcagtaaTAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGTTTTGGCGAAAGGGCGGTGGCCTTGCTGTCGATCCGGGTCACCAGCAGCTGGTCGa                         >  1:528269/1‑150 (MQ=255)
                                                                                                                          ctatgcagtaaTAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCATGCAGTTTTGGCGAAAGGGCGGTGGCCTTGCTGTCGATCCGGGTCACCAGCAGCTGGTCGATg                       <  1:999096/150‑1 (MQ=255)
                                                                                                                            atgcagtaaTAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGTTTTGGCGAAAGGGCGGTGGCCTTGCTGTCGATCCGGGTCACCAGCAGCTGGTCGATGCg                     >  1:1125338/1‑150 (MQ=255)
                                                                                                                             tgcagtaaTAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGTTTTGGCGAAAGGGCGGTGGCCTTGCTGTCGATCCGGGTCACCAGCAGCTGGTCGATGCgg                    >  1:807062/1‑150 (MQ=255)
                                                                                                                              gcagtaaTAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGTTTTGGCGAAAGGGCGGTGGCCTTGCTGTCGATCCGGGTCACCAGCAGCTGGTCGATGCGGt                   <  2:1528949/150‑1 (MQ=255)
                                                                                                                              gcagtaaTAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGTTTTGGCGAAAGGGCGGTGGCCTTGCTGTCGATCCGGGTCACCAGCAGCTGGTCGATGCGGt                   >  1:1150039/1‑150 (MQ=255)
                                                                                                                              gcagtaaTAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGTTTTGGCGAAAGGGCGGTGGCCTTGCTGTCGATCCGGGTCACCAGCAGCTGGTCGATGCGGt                   <  1:410832/150‑1 (MQ=255)
                                                                                                                               cagtaaTAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGTTTTGGCGAAAGGGCGGTGGCCTTGCTGTCGATCCGGGTCACCAGCAGCTGGTCGATGCGGTa                  >  2:39671/1‑150 (MQ=255)
                                                                                                                               cagtaaTAAGACGGCGCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCCTCAGGCAGTTTTGGCGAAAGGGCGGTGGCCTTGCTGTCGATCCGGGTCACCAGCAGCTGGTCGATGCGGTa                  <  1:1050297/150‑1 (MQ=255)
                                                                                                                                agtaaTAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGTTTTGGCGAAAGGGCGGTGGCCTTGCTGTCGATCCGGGTCACCAGCAGCTGGTCGATGCGGTAg                 >  1:1899366/1‑150 (MQ=255)
                                                                                                                                 gtaaTAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGTTTTGGCGAAAGGGCGGTGGCCTTGCTGTCGATCCGGGTCACCAGCAGCTGGTCGATGCGGTAGt                <  1:1361896/150‑1 (MQ=255)
                                                                                                                                    catcaaCGGCTCCTGAATCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGTTTTGGCGAAAGGGCGGTGGCCTTGCTGTCGATCCGGGTCACCAGCAGCTGGTCGATGCGGTAGTTAt             >  1:840869/6‑150 (MQ=255)
                                                                                                                                    ataaGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGTTTTGGCGAAAGTGCGGTGGCCTTGCTGTCGATCCGGGTCACCAGCAGCTGGTCGATGCGGTAGTTAt             >  2:980929/1‑150 (MQ=255)
                                                                                                                                     acaGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGGTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGGTTTGGCGAACGGGCGGTGGCCTTGCTGTCGATCCGGGTCACCAGCAGCTGGTCGATGAGGTAGTTATc            <  1:94782/148‑1 (MQ=255)
                                                                                                                                      aaGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGTTTTGGCGAAAGGGCGGTGGCCTTGCTGTCGATCCGGGTCACCAGCAGCTGGTCGATGCGGTAGTTATCg           <  2:199232/150‑1 (MQ=255)
                                                                                                                                       cgccGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGTTTTGGCGAAAGGGCGGTGGCCTTGCTGTCGATCCGGGTCACCAGCAGCTGGTCGATGCGGTAGTTATCGa          <  1:69861/147‑1 (MQ=255)
                                                                                                                                       cgACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGTTTTGGCGAAAGGGCGGTGGCCTTGCTGGCGATCCGGGTCACCAGCAGCTGGTCGATGCGGTAGTTATCGa          <  2:1617658/149‑1 (MQ=255)
                                                                                                                                       caaCGGCTCCTGAATCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGAGTCAGGCAGTTTTGGCGAAAGGGCGGTGGCCTTGCTGTCGATCCGGGTCACCTGCAGCTGGTCGATGCGGTAGTTATCGa          >  2:715737/3‑150 (MQ=255)
                                                                                                                                       aGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTTGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGTTTTGGCGAAAGGGCGGTGGCCTTGCTGTCGAGCCGGGTCACCAGCAGCTGGTCGATGCGGTAGTTATCGa          >  2:1148298/1‑150 (MQ=255)
                                                                                                                                       aGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTCGCGTCAGGCAGTTTTGGCGAAAGGGCGGTGGCCTTGCTGTCGATCCGGGTCACCAGCAGCTGGTCGATGCGGTAGTTATCGa          <  1:999834/150‑1 (MQ=255)
                                                                                                                                       aGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGTTTTGGCGAAAGGGCGGTGGCCTTGCTGTCGATCCGGGTCACCAGCAGCTGGTCGATGCGGTAGTTATCGa          <  1:1094872/150‑1 (MQ=255)
                                                                                                                                       aGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGTTATGGCGAAAGGGCGGTGGCCTTGCAGTCGATCCTGGTCACCAGCAGCTGGTCGATGCGGTAGTTATCGa          >  1:780376/1‑150 (MQ=255)
                                                                                                                                         aCGGCTCCTGAATCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGTTTTGGCGAAAGGGCGGTGGCCTTGCTGTCGATCCGGGTCACCAGCAGCTGGTCGATGCGGTAGTTATCGata        >  1:1000208/1‑150 (MQ=255)
                                                                                                                                           ggCTCCTGAATCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGTTTTGGCGAAAGGGCGGTGGCCTTGCTGTCGATCCGGGTCACCAGCAGCTGGTCGATGCGGTAGTTATCGATATc      >  2:626910/1‑150 (MQ=255)
                                                                                                                                            ccTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATATTTAGCGTCAGGCAGTATTGGCGAAAGGGCGTTCGCCTTGCTCTCGCTCCGGGTCACCAGCAGCTGGTCGCTGCGGTAGTTATCGATTTcc     >  2:983308/2‑150 (MQ=255)
                                                                                                                                               ccTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGTTTTGGCGAAAGGGCGGTGGCCTTGCTGTCGATCCGGGTCACCAGCAGCTGGTCGATGCGGTAGTTATCGATATccacc  <  1:1860689/150‑1 (MQ=255)
                                                                                                                                               ccTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGTTTTGGCGAAAGGGCGGTGGCCTTGCTGTCGATCCGGGTCACCAGCAGCTGGTCGATGCGGTAGTTATCGATATccacc  >  2:1349010/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                
TAAGATTCTTATCCCAAAAGTGAAGGATAAGCTGCATTTGTCAGTGATGTCCGAAGTTAACCGTCAGGTTATGCGTCTGCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGTTTTGGCGAAAGGGCGGTGGCCTTGCTGTCGATCCGGGTCACCAGCAGCTGGTCGATGCGGTAGTTATCGATATCCACC  >  NC_000913/4439697‑4439989

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: