Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,439,845 T→A 9.1% 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 9.1% 338.6 / 16.3 128intergenic (+52/+27)ytfK/ytfLDUF1107 family protein/UPF0053 family inner membrane protein
Reads supporting (aligned to +/- strand):  ref base T (59/57);  new base A (7/5);  total (66/62)
Fisher's exact test for biased strand distribution p-value = 7.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.82e-01

Only 100 of 130 total aligned reads displayed.

GTAAGATTCTTATCCCAAAAGTGAAGGATAAGCTGCATTTGTCAGTGATGTCCGAAGTTAACCGTCAGGTTATGCGTCTGCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGTTTTGGCGAAAGGGCGGTGGCCTTGCTGTCGATCCGGGTCACCAGCAGCTGGTCGATGCGGTAGTTATCGATATCC  >  NC_000913/4439696‑4439986
                                                                                                                                                     |                                                                                                                                             
gTAAGATTCTTATCCCAAAAGTGAAGGATAAGCTGCATTTGTCAGTGATGTCCGAAGTTAACCGTCAGGTTATGCGTCTGCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGAt                                                                                                                                               <  1:1137319/150‑1 (MQ=255)
 tAAGATTCTTATCCCAAAAGTGAAGGATAAGCTGCATTTGTCAGTGATGTCCGAAGTTAACCGTCAGGTTATGCGTCTGCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGAtt                                                                                                                                              <  2:1204378/150‑1 (MQ=255)
  aaGATTCTTATCCCAAAAGTGAAGGATAAGCTGCATTTGTCAGTGATGTCCGAAGTTAACCGTCAGGTTAAGCGTCTGCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTc                                                                                                                                             <  2:399096/150‑1 (MQ=255)
  aaGATTCATATCCCAAAAGTGAAGGATAAGCTGCATTTGTCAGTGATGTCCGAAGTTAACCGTCAGGTTATGCGTCTGCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTc                                                                                                                                             <  1:947318/150‑1 (MQ=255)
    gATTCTTATCCCAAAAGTGAAGGATAAGCTTCATTTGTCAGTGAAGTCCGAAGTTAACAGTCAGGTTATGCGTCTGCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAg                                                                                                                                           >  1:67612/1‑150 (MQ=255)
    gATTCTTATCCCAAAAGTGAAGGATAAGCTGCATTTGTCAGTGATGTCCGAAGTTAACCGTCAGGTTATGCGTCTGCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAg                                                                                                                                           >  2:185559/1‑150 (MQ=255)
    gATTCTTATCCCAAAAGTGAAGGATAAGCTGCATTTGTCAGTGATGTCCGAAGTTAACCGTCAGGTTATGCGTCTGCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAg                                                                                                                                           >  2:68858/1‑150 (MQ=255)
        cTTATCCCAAAAGTGAAGGATAAGCTGCATTTGTCAGTGATGTCCGAAGTTAACCGTCAGGTTATGCGTCTGCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGc                                                                                                                                       >  2:964592/1‑150 (MQ=255)
        cTTATCCCAAAAGTGAAGGATAAGCTGCATTTGTCAGTGATGTCCGAAGTTAACCGTCAGGTTATGCGTCTGCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGAATCAGGAGc                                                                                                                                       <  1:1236287/150‑1 (MQ=255)
        cTTATCCCAAAAGTGAAGGATAAGCTGCATTTGTCAGTGATGTCCGAAGTTAACCGGCAGGTTATGCGTCTGCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGc                                                                                                                                       <  1:872604/150‑1 (MQ=255)
         ttATCCCAAAAGTGAAGGATAAGCTGCATTTGTCAGTGATGTCCGAAGTTAACCGTCAGGTTATGCGTCTGCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGcc                                                                                                                                      >  1:30596/1‑150 (MQ=255)
         ttATCCCAAAAGTGAAGGAAAAGCTGCATTTGTCAGTGATGTCCGAAGTTAACCGTCAGGTTATGCGTCTGCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGcc                                                                                                                                      <  2:1490289/150‑1 (MQ=255)
           aTCCCAAAAGTGAAGGATAAGCTGCATTTGTCAGGGATGTCCGAAGTTAACCGTCAGGTTATGCGTCTGCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGAATCAGGAGCCGt                                                                                                                                    <  2:325604/150‑1 (MQ=255)
             cccAAAAGTGAAGGATAAGCTGCATTTGTCAGTGATGTCCGAAGTTAACCGTCAGGTTATGCGTCTGCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTATAAAGTGGTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGTTg                                                                                                                                  >  2:102026/1‑150 (MQ=255)
             cccAAAAGTGAAGGATAAGCTGCATTTGTCAGTGATGTCCGAAGTTAACCGTCAGGTTATGCGTCTGCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGAATCAGGAGCCGtct                                                                                                                                  <  1:1140203/150‑3 (MQ=255)
              ccAAAAGTGAAGGATAAGCTGCATTTGTCAGTGATGTCCGAAGTTAACCGTCAGGTTATGCGTCTGCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGTTGa                                                                                                                                 >  2:548203/1‑150 (MQ=255)
              ccAAAAGTGAAGGATAAGCTGCATTTGTCAGTGATGTCCGAAGTTAACCGTCAGGTTATGCGTCTGCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGTTGa                                                                                                                                 >  2:548202/1‑150 (MQ=255)
              ccAAAAGTGAAGGATAAGCTGCATTTGTCAGTGATGTCCGAAGTTAACCGTCAGGTTATGCGTCTGCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGAATCAGGAGCCGtctt                                                                                                                                 <  2:1705275/150‑4 (MQ=255)
               cAAAAGTGAAGGATAAGCTGCATTTGTCAGTGATGTCCGAAGTTAACCGTCAGGTTATGCGTCTGCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGTTGAt                                                                                                                                >  2:743778/1‑150 (MQ=255)
                aaaaGTGAAGGATAAGCTGCATTTGTCAGTGATGTCCGAAGTTACCCGTCAGGTTATGCGTCTGCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTATTAAGACGGCTCCTGATTCAGGAGCCGTTGATg                                                                                                                               >  2:255428/1‑150 (MQ=255)
                  aaGTGAAGGATAAGCTGCATTTGTCAGTGATGTCCGAAGTTAACCGTCAGGTTATGCGTCTGCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGTTGATGtt                                                                                                                             <  1:68858/150‑1 (MQ=255)
                     tGAAGGATAAGCTGCATTTGTCAGTGATGTCCGAAGTTAACCGTCAGGTTATGCGTCTGCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCt                                                                                                                          >  2:1675189/1‑150 (MQ=255)
                      gAAGGATAAGCTGCATTTGTCAGTGATGTCCGAAGTTAACCGTCAGGTTATGCGTCTGCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCGTGATTCAGGAGCCGTTGATGTTTCTg                                                                                                                         >  2:255107/1‑150 (MQ=255)
                         ggATAAGCTGCATTTGTCAGTGATGTCCGAAGTTAACCGTCATGTTATGCGTCTGCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTgggg                                                                                                                      <  1:1052957/150‑1 (MQ=255)
                         ggATAAGCTGCATTTGTCAGTGATGTCCGAAGTTAACCGTCAGGTTATGCGTCTGCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTgggg                                                                                                                      >  1:899386/1‑150 (MQ=255)
                          gATAAGCTGCATTTGTCAGTGATGTCCGAAGTTAACCGTCAGGTTATGCGTCTGCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTggggg                                                                                                                     <  1:1000156/150‑1 (MQ=255)
                           aTAAGCTGCATTTGTCAGTGATGTCCGAAGTTAACCGTCAGGTTATGCGTCTGCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGt                                                                                                                    <  2:814207/150‑1 (MQ=255)
                                cTGCATTTGTCAGTGATGTCCGAAGTTTACCGTCAGGTTATGCGTCTGCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTAcgc                                                                                                               >  2:1013347/1‑150 (MQ=255)
                                 tGCATTTGTCAGTGATGTCCGAAGTTAACCGTCATGTTATGCGTCTGCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTAcgcg                                                                                                              <  1:1294042/150‑1 (MQ=255)
                                  gCATTTGTCAGTGATGTCCGAAGTTAACCGTCAGGTTATGCGTCTGCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGa                                                                                                             >  1:1334525/1‑150 (MQ=255)
                                       tGTCAGTGATGTCCGAAGTTAACCGTCAGGTTATGCGTCTGCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCtt                                                                                                        >  2:1271914/1‑150 (MQ=255)
                                          cAGTGATGTCCGAAGTTAACCGTCAGGTTATGCGTCTGCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTttct                                                                                                     >  1:1338580/1‑150 (MQ=255)
                                          cAGTGATGTCCGAAGTTAACCGTCAGGTTATGCGTCTGCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTttct                                                                                                     >  1:1224966/1‑150 (MQ=255)
                                           cgggATGTCCGAAGTTAACCGTCAGGTTATGCGTCGGCAAACAGAGATGGCTAAACCAAAGGGCAATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTttctt                                                                                                    <  1:306349/147‑1 (MQ=255)
                                            gTGATGTCCGAAGTTAACCGTCAGGTTATGCGTCTGCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTttcttc                                                                                                   <  1:21020/150‑1 (MQ=255)
                                            gTGATGTCCGAAGTTAACCGTCAGGTTATGCGTCTGCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTttcttc                                                                                                   <  1:210867/150‑1 (MQ=255)
                                             tGATGTCCGAAGTTAACCGTCAGGTTATGCGTCTGCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTtcttct                                                                                                  >  1:1206521/1‑150 (MQ=255)
                                              gATGTCCGAAGTTAACCGTCAGGTTATGCGTCTGCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTcttctt                                                                                                 <  1:1361835/150‑1 (MQ=255)
                                                tgaccGAAGTTAACCGTCAGGTTATGCGTCTGCAAACTGAGATGGCTTCACCAAAGTGATATGCTGTAATAAGTGCTATCCAGTAATAAGACGGCTCCTCAATCAGGAGACGTTGATGTTTCTGGGGGTTAAGAAAAGCTTTcttcttaa                                                                                               >  1:97577/4‑148 (MQ=255)
                                                tGTCCGAAGTTAACCGTCAGGTTATGCGTCTGCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTtcttta                                                                                               <  2:481851/150‑1 (MQ=255)
                                                 gTCCGAAGTTAACCGTCAGGTTATGCGTCTGCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTtctttat                                                                                              >  1:580000/1‑150 (MQ=255)
                                                  tCCGAAGTTAACCGTCAGGTTATGCGTCTGCAAACAGAGATTGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTtctttatc                                                                                             <  2:1084204/150‑1 (MQ=255)
                                                   ccGAAGTTAACCGTCAGGTTATGCGTCTGCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTtctttatct                                                                                            >  1:1618118/1‑150 (MQ=255)
                                                     gAAGTTAACCGTCAGTTTATGCGTCTGCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTtctttatcttt                                                                                          <  2:98749/150‑1 (MQ=255)
                                                           aaCCGTCAGGTTATGCGTCTGCAAACACAGAAGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTTGGGGTTACGCGACGGTATCTACTTTATCTTAAGCGta                                                                                    >  1:166256/1‑149 (MQ=255)
                                                              cGTCAGGTTATGCGTCTGCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAgg                                                                                 <  2:1462112/150‑1 (MQ=255)
                                                              cGTCAGGTTATGCGTCTGCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAgg                                                                                 <  2:899386/150‑1 (MQ=255)
                                                              cGTCAGGTTATGCGTCTGCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAgg                                                                                 <  1:252805/150‑1 (MQ=255)
                                                                 cAGGTTATGCGTCTGCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAg                                                                              <  2:1191529/150‑1 (MQ=255)
                                                                   ggTTATGCGTCTGCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTTGCGTCAGGAAGtt                                                                            >  2:1265106/1‑150 (MQ=255)
                                                                   ggTTATGCGTCTGCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGtt                                                                            <  1:668885/150‑1 (MQ=255)
                                                                   ggTTATGCGTCTGCAAACAGAGATGGCATAACCAAAGTGCTATGCAGTAAAAAGTGCTATTCAGTAATAAGACGGCTCCTGATTCAGGAGCCGTTCATGTATCTGGGGGTTAACCAACTCTTTCTTCGTTATCTGTAGAGTCAGGCatta                                                                            >  2:849743/1‑147 (MQ=255)
                                                                      tATGCGTCTGCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTGCGTTTTCTTTAGCGTCAGGCAGTTTTg                                                                         >  2:697955/1‑150 (MQ=255)
                                                                        tGCGTCTGCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGTTTTGGc                                                                       >  1:949317/1‑150 (MQ=255)
                                                                         gCGTCTGCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGTTTTGGCg                                                                      <  2:580000/150‑1 (MQ=255)
                                                                         gCGTCTGCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGTTTTGGCg                                                                      >  2:527364/1‑150 (MQ=255)
                                                                              tGCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCa                                                                                   >  1:144857/1‑132 (MQ=255)
                                                                              tGCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGTTTTGGCGAAAgg                                                                 >  1:1457760/1‑150 (MQ=255)
                                                                              tGCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGTTTTGGCGAAAgg                                                                 >  1:19888/1‑150 (MQ=255)
                                                                              tGCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGTTTCGGCGAAAgg                                                                 >  1:490998/1‑150 (MQ=255)
                                                                               gCAAACAGAGATGGCTTAACCAAAGTGCTCTGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGTTTTGGCGAAAggg                                                                <  1:1199050/150‑1 (MQ=255)
                                                                               gCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGTTTTGGCGAAAggg                                                                <  2:1224966/150‑1 (MQ=255)
                                                                               gCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGTTTTGGCGAAAggg                                                                >  1:1639347/1‑150 (MQ=255)
                                                                               gCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGTTTTGGCGAAAggg                                                                >  1:713342/1‑150 (MQ=255)
                                                                                 aatcAGAGATCTCTTATCCAAAGTGCTATGCTGTAAAAAGGTCTATGCAGTAATCAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGTGTTTCGCGTCTCTTTCTTCTTTATCTTTAGCGTCAGGCAGATTTGGCGAAAGGGCg                                                              <  2:97577/147‑1 (MQ=255)
                                                                                        gaTGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGTTTTGGCGAAAGGGCGGTGGCCt                                                       <  2:1052219/150‑1 (MQ=255)
                                                                                          tGGCTTAATCAAAGGGCTATGGCGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGTTTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGTTTTGGCTAAAGGGCGGTGGCCTTg                                                     <  1:961141/150‑1 (MQ=255)
                                                                                           ggCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGGTTTGGCGAAAGGGCGGTGGCCTTGc                                                    >  2:114757/1‑150 (MQ=255)
                                                                                            gCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGTTTTGGCGAAAGGGCGGTGGCCTTGCt                                                   <  2:949317/150‑1 (MQ=255)
                                                                                            gCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGTTTTGGCGAAAGGGCGGTGGCCTTGCt                                                   >  2:709902/1‑150 (MQ=255)
                                                                                            gCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGTTTTGGCGAAAGGGCGGTGGCCTTGCt                                                   <  1:365574/150‑1 (MQ=255)
                                                                                            gCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGTTTTGGCGAAAGGGCGGTGGCCTTGCt                                                   <  1:359283/150‑1 (MQ=255)
                                                                                               tAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGTTTTGGCGAAAGGGCGGTGGCCTTGCTGTc                                                <  1:1675189/150‑1 (MQ=255)
                                                                                                aaCCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGTTTTGGCGAAAGGGCGGTGGCCTTGCTGTCg                                               >  2:744351/1‑150 (MQ=255)
                                                                                                aaCCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTATATCTTTAGCGTCATGCAGTTTTGGCGAAAGGGCGGTGGCCTTGCTGTCg                                               >  1:10037/1‑150 (MQ=255)
                                                                                                   cAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGTCAGTTTTGGCGAAAGGGCGGTGGCCTTGCTGTCGATc                                            >  2:1221527/1‑150 (MQ=255)
                                                                                                   cAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGTTTTGGCGAAAGGGCGGTGGCCTTGCTGTCGATc                                            >  2:1642/1‑150 (MQ=255)
                                                                                                                 gtaaaaagtgctatgcagtaaTAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGTTTTGGCGAAAGGGCGGTGGCCTTGCTGTCGATCCGGGTCACcagcag                              >  1:1461155/1‑150 (MQ=255)
                                                                                                                    aaaagtgctatgcagtaaTAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGTTTTGGCGAAAGGGCGGTGGCCTTGCTGTCGATCCGGGTCACCAGCAGCtt                           >  2:1343714/1‑149 (MQ=255)
                                                                                                                     aaagtgctatgcagtaaTAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGTTTTGGCGAAAGGGCGGTGGCCTTGCTGTCGATCCGGGTCACCAGCAGCTgg                          >  1:1190788/1‑150 (MQ=255)
                                                                                                                      aagtgctatgcagtaaTAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACTCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGTTTTGGCGAAAGGGCGGTGGCCTTGCTGTCGATCCGGGTCACCAGCAGCTGGt                         <  1:1509096/150‑1 (MQ=255)
                                                                                                                       agtgctatgcagtaaTAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGTTTTGGCGAAAGGGCGGTGGCCTTGCTGTCGATCCGGGTCACCAGCAGCTGGTc                        <  2:1339263/150‑1 (MQ=255)
                                                                                                                       actgctatgcagaaaGAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGTTTTGGCGAAAGGGCGGTGGCCTTGCTGTCGATCCGGGTCACCAGCAGCTGGTc                        <  1:102026/148‑1 (MQ=255)
                                                                                                                        gtgctatgcagtaaTAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGTTTTGGCGAAAGGGCGGTGGCCTTGCTGTCGATCCGGGTCACCAGCAGCTGGTCg                       <  2:1136988/150‑1 (MQ=255)
                                                                                                                          gctatgcagtaaTAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGTTTTGGCGAAAGGGCGGTGGCCTTGCTGTCGATCAGGGTCACCAGCAGCTGGTCGAt                     >  1:225951/1‑150 (MQ=255)
                                                                                                                           ctatgcagtaaTAAGACGGCTCCTGATTCGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGTTTTGGCGAAAGGGCGGTGGCCTTGCTGTCGATCCGGGTCACCAGCAGCTGGTCGATGc                   >  2:1712393/1‑150 (MQ=255)
                                                                                                                            tatgcagtaaTAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGTTTTGGCGAAAGGGCGGTGGCCTTGCTGTCGATCCGGGTCACCAGCTGCTGGtgacccg                   <  2:1457760/150‑7 (MQ=255)
                                                                                                                              tgcagtaaTAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGTTTTGGCGAAAGGGCGGTGGCCTTGCTGTCGATCCGGGTCACCAGCAGCTGGTCGATGCgg                 >  2:216430/1‑150 (MQ=255)
                                                                                                                               gcagtaaTAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGTTTTGGCGAAAGGGCGGTGGCCTTGCTGTCGATCCGGGTCACCAGCAGCTGGTCGATGCGGt                <  1:743778/150‑1 (MQ=255)
                                                                                                                                cagtaaTAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGTTTTGGCGAAAGGGCGGTGGCCTTGCTGTCGATCCGGGTCACCAGCAGCTGGTCGATGCGGTa               >  1:1243300/1‑150 (MQ=255)
                                                                                                                                cagtaaTAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTATTTATCTTTAGCGTCAGGCAGTTTTGGCGAAAGGGCGGTGGCCTTGCTGTCGATCCGGGTCACCAGCAGCAGGTCGATGCGGTa               >  1:346448/1‑150 (MQ=255)
                                                                                                                                 ataaatacGCCGGCGCCTTATTCAGGAGCGGTTTATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGGCAGGCAGTTTGGGCGAACGGGCGGTGGCCTTGCTGGCGAACCGGGTCACCAGCAGCTGGTCGAGGCGGTAg              <  2:1190788/147‑1 (MQ=255)
                                                                                                                                 agtaaTAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGTTTTGGCGAAAGGGCGGTGGCCTTGCTGTCGATCCGGGTCACCAGCAGCTGGTCGATGCGGTAg              <  1:183483/150‑1 (MQ=255)
                                                                                                                                      atcaaCGGCTCCTGAATCAGGAGCCGTTGATGTATCTGGGGGTTACGCGACGCTTTCTACTTTAACATTAGCGACAGGCAGTTTTGGCGAAAGCGCGGTGGCCTTGCTGTCGATCCGGGTCACAAGCAGCTGGTCGATGCGGTATTTATc         >  2:1405086/5‑150 (MQ=255)
                                                                                                                                         gACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGTGTTACGCTACGCATTCATCTATATCATTAGCGTCAGGCAGTTTTGGCGAAAGGGCGGTGGCCTTGCTGTAGATCAGGGTCACCAGCAGCTGGTCGATGCGCTAGTTATCGat      >  1:858691/1‑150 (MQ=255)
                                                                                                                                         gACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTGTATCTTTAGCGTCAGGCAGTTTTGGCGAAAGGGCGGTGGCCTTGCTGTCGATCCGGGTCACCAGCAGCTGGTCGATGCGGTAGTTATCGat      >  1:1168649/1‑150 (MQ=255)
                                                                                                                                          aCGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGTTTTGGCGAAAGGGCGGTGGCCTTGCTGTCGATCCGGGTCACCAGCAGCTGGTCGATGCGGTAGTTATCGata     >  2:1647668/1‑150 (MQ=255)
                                                                                                                                           cGGCTCCTGAATCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGTTTTGGCGAAAGGGCGGTGGCCTTGCTGTCGATCCGGGTCACCAGCAGCTGGTCGATGCGGTAGTTATCGatat    >  1:1523776/1‑150 (MQ=255)
                                                                                                                                             gCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGTTTTGGCGAAAGGGCGGTGGCCTTGCTGTCGATCCGGGTCACCAGCAGCTGGTCGATGCGGTAGTTATCGATATcc  <  2:713342/150‑1 (MQ=255)
                                                                                                                                             gCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGTTTTGGCGAAAGGGCGGTGGCCTTGCTGTCGATCCGGGTCACCAGCAGCTGGTCGATGCGGTAGTTATCGATATcc  <  2:1639347/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                             
GTAAGATTCTTATCCCAAAAGTGAAGGATAAGCTGCATTTGTCAGTGATGTCCGAAGTTAACCGTCAGGTTATGCGTCTGCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGTTTTGGCGAAAGGGCGGTGGCCTTGCTGTCGATCCGGGTCACCAGCAGCTGGTCGATGCGGTAGTTATCGATATCC  >  NC_000913/4439696‑4439986

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