Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,259,286 T→C 5.1% K4R (AAG→AGG)  psuK ← pseudouridine kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,259,2860TC 5.1% 257.5 / 4.1 79K4R (AAG→AGG) psuKpseudouridine kinase
Reads supporting (aligned to +/- strand):  ref base T (31/44);  new base C (2/2);  total (33/46)
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.74e-01

TGTGCAATATTGCGCCCTACTCCACCAGGCGTAAATTTTATTTTACCTGGATTTGAATCCGCATAATTTAATGATTCATGTGAATATCCGGCGACATCAATATTCGCCGAACCTATAATTACGACATAATCCTTTTCGCGCATAGCACGTCCTCTGGTAAAAAATTACCAACCCTGTCAATCCACAAGAATATAGATACTGCAAAGACGGTTTAAACATGTGTTCATATTTAAACATGTGCTCATAGTATGCATGTTTTACGTAAAGTAAAT  >  NC_000913/2259154‑2259425
                                                                                                                                    |                                                                                                                                           
tgtgCAATATTGCGCCCTACTCCACCAGGCGTAAATTTTATTTTACCTGGATTTGAATCCGCATAATTTAATGATTCATGTGAATATCCGGCGACATCAATATTCGCCGAACCTATAATTACGACATAATCCTTTTcgcg                                                                                                                                      <  1:315434/140‑1 (MQ=255)
  tgCAATATTGCGCCCTACTCCACCAGGCGTAAATTTTATTTTACCTGGATTTGAATCCGCATAATTTAATGATTCATGTGAATATCCGGCGACATCAATATTCGCCGAACCTATAATTACGACATAATCCTTTTCGCGCa                                                                                                                                    >  1:1647090/1‑140 (MQ=255)
      atatTGCGCCCTACTCCACCAGGCGTAAATTTTATTTTACCTGGATTTGAATCCGCATAATTTAATGATTCATGTGAATATCCGGCGACATCAATATTCGCCGAACCTATAATTACGACATAATCCTTTTCGCGCATAGc                                                                                                                                >  1:1873146/1‑140 (MQ=255)
            cgcgCTACTCCACCAGGGGTAAATTTTATTTTACCTTGATTTGATTCCGCATATTTTAATGATTCATGTGAATATGCGTCGTAATAATTATTCGCAGTAGCTTTAATTTCGAGTTAATCGTTTTCGCTCATTGCTCGTcg                                                                                                                          >  2:197083/1‑139 (MQ=255)
            cgcCCTACTCCACCAGGCGTAAATTTTATTTTACCTGGATTTGAATCCTCATAATTTAATGATTCATGTGAATATCCGGCGACATCAATATTCGCCGAACCTAAAATTACGACATAATCCTTTTCGCGCATAGCACGTcc                                                                                                                          >  1:1169515/1‑140 (MQ=255)
                 tACTCCACCAGGCGTAAATTTTATTTTACCTGGATTTGAATCCGCATAATTTAATGATTCATGTGAATATCCGGCGACATCAATATTCGCCGAACCTATAATTACGACATAATCCCTTTCGCGCATAGCACGTCCTCTgg                                                                                                                     >  1:494885/1‑140 (MQ=255)
                 tACTCCACCAGGCGTAAATTTTATTTTACCTGGATTTGAATCCGCATAATTTAATGATTCATGTGAATATCCGGCGACATCAATATTCGCCGAACCTATAATTACGACATAATCCCTTTCGCGCATAGCACGTCCTCTgg                                                                                                                     >  1:494835/1‑140 (MQ=255)
                        ccaGGCGTAAATTTTATTTTACCTGGATTTGAATCCGCATAATTTAATGATTCATGTGAATATCCGGCGACATCAATATTCGCCGAACCTATAATTACGACATAATCCTTTTCGCGCATAGCACGTCCTCTGGTaaaaaa                                                                                                              >  1:1756354/1‑140 (MQ=255)
                         caGGCGTAAATTTTATTTTACCTGGATTTGAATCCGCATAATTTAATGATTCATGTGAATATCCGGCGACATCAATATTCGCCGAACCTATAATTACGACATAATCCTTTTCGCGCATAGCACGTCCTCTGGTAAAAAAt                                                                                                             >  1:1900446/1‑140 (MQ=255)
                           ggCGTAAATTTTATTTTACCTGGATTTGAATCCGCATAATTTAATGATTCATGTGAATATCCGGCGACATCAATATTCGCCGAACCTATAATTACGACATAATCCTTTTCGCGCATAGCACGTCCTCTGGTAAAAAATTa                                                                                                           >  1:429484/1‑140 (MQ=255)
                           ggCGTAAATTTTATTTTACCTGGATTTGAATCCGCATAATTTAATGATTCATGTGAATATCCGGCGACATCAATATTCGCCGAACCTATAATTACGACATAATCCTTTTCGCGCATAGCACGTCCTCTGGTAAAAAATTa                                                                                                           >  2:267590/1‑140 (MQ=255)
                           ggCGTAAATTTTATTTTACCTGGATTTGAATCCGCATAATTTAATGATTCATGTGAATATCCGGCGACATCAATATTCGCCGAACCTATAATTACGACATAATCCTTTTCGCGCATAGCACGTCCTCTGGTAAAAAATTa                                                                                                           >  1:860043/1‑140 (MQ=255)
                                aaaTTTTATTTTACCTGGATTTGAATCCGCATAATTTAATGATTCATGTGAATATCCGGCGACATCAATATTCGCCGAACCTATAATTACGACATAATCCTTTTCGCGCATAGCACGTCCTCTGGTAAAAAATTACCAAc                                                                                                      >  2:741016/1‑140 (MQ=255)
                                          ttaCCTGGATTTGAATCCGCATAATTTAATGATTCATGTGAATATCCGGCGACATCAATATTCGCCGAACCTATAATTACGACATAATCCTTTTCGCGCATAGCACGTCCTCTGGTAAAAAATTACCAACCCTGTCAATc                                                                                            <  2:584478/140‑1 (MQ=255)
                                             ccTGGATTTGAATCCGCATAATTTAATGATTCATGTGAATATCCGGCGACATCAATATTCGCCGAACCTATAATTACGACATAATCCTTTTCGCGCATAGCACGTCCTCTGGTAAAAAATTACCAACCCTGTCAATCcac                                                                                         <  1:1365082/140‑1 (MQ=255)
                                             ccTGGATTTGAATCCGCATAATTTAATGATTCATGTGAATATCCGGCGACATCAATATTCGCCGAACCTATAATTACGACATAATCCTTTTCGCGCATAGCACGTCCTCTGGTAAAAAATTACCAACCCTGTCAATCcac                                                                                         >  1:121366/1‑140 (MQ=255)
                                             ccTGGATTTGAATCCGCATAATTTAATGATTCATGTGAATATCCGGCGACATCAATATTCGCCGAACCTATAATTACGACATAATCCTTTTCGCGCATAGCACGTCCTCTGGTAAAAAATTACCAACCCTGTCAATCcac                                                                                         <  1:400780/140‑1 (MQ=255)
                                             ccTGGATTTGAATCCGCATAATTTAATGATTCATGTGAATATCCGGCGACATCAATATTCGCCGAACCTATAATTACGACATAATCCTTTTCGCGCATAGCACGTCCTCTGGTAAAAAATTACCAACCCTGTCAATCcac                                                                                         >  1:120848/1‑140 (MQ=255)
                                              cTGGATTTGAATCCGCATAATTTAATGATTCATGTGTATATCCGGCGACATCAATATTCGCCGAACCTATAATTACGACATAATCCTTTTCGCGCATAGCACGTCCTCTGGTAAAAAATTACCAACCCTGTCAATCcaca                                                                                        <  2:1070282/140‑1 (MQ=255)
                                                ggATTTGAATCCGCATAATTTAATGATTCATGTGAATATCCGGCGACATCAATATTCGCCGAACCTATAATTACGACATAATCCTTTTCGCGCATAGCACGTCCTCTGGTAAAAAATTACCAACCCTGTCAATCCACAAg                                                                                      <  1:120067/140‑1 (MQ=255)
                                                  aTTTGAATCCGCATAATTTAATGATTCATGTGAATATCCCGCGACATCAATATTCGCCGAACCTATAATTACGACATAAACCTATACGCGCATAGCACCTCCTCTGGTAAAAAATTACCAACCCTGTCAATCCACAAGaa                                                                                    >  2:255993/1‑140 (MQ=255)
                                                     tGAATCCGCATAATTTAATGATTCATGTGAATATCCGGCGACATCAATATTCGCCGAACCTATAATTACGACATAATCCTTTTCGCGCATAGCACGTCCTCTGGTAAAAAATTACCAACCCTGTCAATCCACAAGAatat                                                                                 <  1:539228/140‑1 (MQ=255)
                                                        aTCCGCATAATTTAATGATTCATGTGAATATCCGGCGACATCAATATTCGCCGAACCTATAATTACGACATAATCCTTTTCGCGCATAGCACGTCCTCTGGTAAAGAATTACCAACCCTGTCAATCCACAAGAATATAGa                                                                              <  1:1331277/140‑1 (MQ=255)
                                                        aTCCGCATAATTTAATGATTCATGTGAATATCCGGCGACATCAATATTCGCCGAACCTATAATTACGACATAATCCTTTTCGCGCATAGCACGTCCTCTGGTAAAAAATTACCAACCCTGTCAATCCACAAGAATATAGa                                                                              <  2:75373/140‑1 (MQ=255)
                                                        aTCCGCATAATTTAATGATTCATGTGAATATCCGGCGACATCAATATTCGCCGAACCTATAATTACGACATAATCCTTTTCGCGCATAGCACGTCCTCTGGTAAAAAATTACCAACCCTGTCAATCCACAAGAATATAGa                                                                              <  1:1664374/140‑1 (MQ=255)
                                                         tCCGCATAATTTAATGATTCATGTGAATATCCGGCGACATCAATATTCGCCGAACCTATAATTACGACATAATCCTTTTCGCGCATAGCACGTCCTCTGGTAAAAAATTACCAACCCTGTCAATCCACAAGAATATAGAt                                                                             <  1:1660833/140‑1 (MQ=255)
                                                           cGCATAATTTAATGATTCATGTGAATATCCGGCGACATCAATATTCGCCGAACCTATAATTACGACATAATCCTTTTCGCGCATAGCACGTCCTCTGGTAAAAAATTACCAACCCTGTCAATCCACAAGAATATAGATAc                                                                           <  1:1061396/140‑1 (MQ=255)
                                                           cGCATAATTTAATGATTCATGTGAATATCCGGCGACATCAATATTCGCCGAACCTATAATTACGACATAATCCTTTTCGCGCATAGCACGTCCTCTGGTAAAAAATTACCAACCCTGTCAATCCACAAGAATATAGATAc                                                                           >  1:1278662/1‑140 (MQ=255)
                                                             cATAATTTAATGATTCATGTGAATATCCGGCGACATCAATATTCGCCGAACCTATAATTACGACATAATCCTTTTCGCGCATAGCACGTCCTCTGGTAAAAAATTACCAACCCTGTCAATCCACAAGAATATAGATACTg                                                                         <  2:1162022/140‑1 (MQ=255)
                                                                  tttAATGATTCATGTGAATATCCGGCGACATCAATATTCGCCGAACCTATAATTACGACATAATCCTTTTCGCGCATAGCACGTCCTCTGGTAAAAAATTACCAACCCTGTCAATCCACAAGAATATAGATACTGCAAAg                                                                    <  2:1354662/140‑1 (MQ=255)
                                                                   ttAATGATTCATGTGAATATCCGGCGACATCAATATTCGCCGAACCTATAATTACGACATAATCCTTTTCGCGCATAGCACGTCCTCTGGTAAAAAATTACCAACCCTGTCAATCCACAAGAATATAGATACTGCAAAGa                                                                   <  1:1170175/140‑1 (MQ=255)
                                                                   ttAATGATTCATGTGAATATCCGGCGACATCAATATTCGCCGAACCTATAATTACGACATAATCCTTTTCGCGCATAGCACGTCCTCTGGTAAAAAATTACCAACCCTGTCAATCCACAAGAATATAGATACTGCAAAGa                                                                   <  1:579415/140‑1 (MQ=255)
                                                                     aaTGATTCATGTGAATATCCGGCGACATCAATATTCGCCGAACCTATAATTACGACATAATCCTTTTCGCGCATAGCACGTCCTCTGGTAAAAAATTACCAACCCTGTCAATCCACAAGAATATAGATACTGCAAAGACg                                                                 <  2:877142/140‑1 (MQ=255)
                                                                            cATGTGAATATCCGGCGACATCAATATTCGCCGAACCTATAATTACGACATAATCCTTTTCGCGCATAGCACGTCCTCTGGTAAAAAATTACCAACCCTGTCAATCCACAAGAATATAGATACTGCAAAGACGGTTTaaa                                                          <  2:1625324/140‑1 (MQ=255)
                                                                            cATGTGAATATCCGGCGACATCAATATTCGCCGAACCTATAATTACGACATAATCCTTTTCGCGCATAGCACGTCCTCTGGTAAAAAATTACCAACCCTGTCAATCCACAAGAATATAGATACTGCAAAGACGGTTTaaa                                                          >  2:540147/1‑140 (MQ=255)
                                                                            cATGTGAATATCCGGCGACATCAATATTCGCCGAACCTATAATTACGACATAATCCTTTTCGCGCATAGCACGTCCTCTGGTAAAAAATTACCAACCCTGTCAATCCACAAGAATATAGATACTGCAAAGACGGTTTaaa                                                          >  2:558781/1‑140 (MQ=255)
                                                                            cATGTGAATATCCGGCGACATCAATATTCGCCGAACCTATAATTACGACATAATCCTTTTCGCGCATAGCACGTCCTCTGGTAAAAAATTACCAACCCTGTCAATCCACAAGAATATAGATACTGCAAAGACGGTTTaaa                                                          >  1:1616890/1‑140 (MQ=255)
                                                                            cATGTGAATATCCGGCGACATCAATATTCGCCGAACCTATAATTACGACATAATCCTTTTCGCGCATAGCACGTCCTCTGGTAAAAAATTACCAACCCTGTCAATCCACAAGAATATAGATACTGCAAAGACGGTTTaaa                                                          <  2:1624713/140‑1 (MQ=255)
                                                                            cATGTGAATATCCGGCGACATCAATATTCGCCGAACCTATAATTACGACATAATCCTTTTCGCGCATAGCACGTCCTCTGGTAAAAAATTACCAACCCTGTCAATCCACAAGAATATAGATACTGCAAAGACGGTTTaaa                                                          >  1:1396698/1‑140 (MQ=255)
                                                                            cATGTGAATATCCGGCGACATCAATATTCGCCGAACCTATAATTACGACATAATCCTTTTCGCGCATAGCACGTCCTCTGGTAAAAAATTACCAACCCTGTCAATCCACAAGAATATAGATACTGCAAAGACGGTTTaaa                                                          >  2:1277187/1‑140 (MQ=255)
                                                                             aTGTGAATATCCGGCGACATCAATATTCGCCGAACCTATAATTACGACATAATCCTTTTCGCGCATAGCACGTCCTCTGGTAAAAAATTACCAACCCTGTCAATCCACAAGAATATAGATACTGCAAAGACGGTTTAAAc                                                         <  1:990832/140‑1 (MQ=255)
                                                                             aTGTGAATATCCGGCGACATCAATATTCGCCGAACCTATAATTACGACATAATCCTTTTCGCGCATAGCACGTCCTCTGGTAAAAAATTACCAACCCTGTCAATCCACAAGAATATAGATACTGCAAAGACGGTTTAAAc                                                         >  2:990832/1‑140 (MQ=255)
                                                                             aTGTGAATATCCGGCGACATCAATATTCGCCGAACCTATAATTACGACATAATCCTTTTCGCGCATAGCACGTCCTCTGGTAAAAAATTACCAACCCTGTCAATCCACAAGAATATAGATACTGCAAAGACGGTTTAAAc                                                         <  1:1740907/140‑1 (MQ=255)
                                                                               gtgAATATCCGGCGACATCAATATTCGCCGAACCTATAATTACGACATAATCCTTTTCGCGCATAGCACGTCCTCTGGTAAAAAATTACCAACCCTGTCAATCCACAAGAATATAGATACTGCAAAGACGGTTTAAACAt                                                       <  1:1848419/140‑1 (MQ=255)
                                                                               gtgAATATCCGGCGACATCAATATTCGCCGAACCTATAATTACGACATAATCCTTTTCGCGCATAGCACGTCCTCTGGTAAAAAATTACCAACCCTGTCAATCCACAAGAATATAGATACTGCAAAGACGGTTTAAACAt                                                       <  2:321488/140‑1 (MQ=255)
                                                                               gtgAATATCCGGCGACATCAATATTCGCCGAACCTATAATTACGACATAATCCTTTTCGCGCATAGCACGTCCTCTGGTAAAAAATTACCAACCCTGTCAATCCACAAGAATATAGATACTGCAAAGACGGTTTAAACAt                                                       <  2:321477/140‑1 (MQ=255)
                                                                                    tatCCGGCGACATCAATATTCGCCGAACCTATAATTACGACATAATCCTTTTCGCGCATAGCACGTCCTCTGGTAATAAATTACCAACCCTGTCAATCCACAAGAATATAGATACTGCAAAGACGGTTTAAACATGTGtt                                                  >  2:1092128/1‑140 (MQ=255)
                                                                                    tatCCGGCGACATCAATATTCGCCGAACCTATAATTACGACATAATCCTTTTCGCGCATAGCACGTCCTCTGGTAAAAAATTACCAACCCTGTCAATCCACAAGAATATAGATACTGCAAAGACGGTTTAAACATGTGtt                                                  >  2:1092151/1‑140 (MQ=255)
                                                                                    tatCCGGCGACATCAATATTCGCCGAACCTATAATTACGACATAATCCTTTTCGCGCATAGCACGTCCTCTGGTAAAAAATTACCAACCCTGTCAATCCACAAGAATATAGATACTGCAAAGACGGTTTAAACATGTGtt                                                  >  2:1327316/1‑140 (MQ=255)
                                                                                    tatCCGGCGACATCAATATTCGCCGAACCTATAATTACGACATAATCCTTTTCGCGCATAGCACGTCCTCTGGTAAAAAATTACCAACCCTGTCAATCCACAAGAATATAGATACTGCAAAGACGGTTTAAACATGTGtt                                                  >  2:1774545/1‑140 (MQ=255)
                                                                                    tatCCGGCGACATCAATATTCGCCGAACCTATAATTACGACATAATCCTTTTCGCGCATAGCACGTCCTCTGGTAAAAAATTACCAACCCTGTCAATCCACAAGAATATAGATACTGCAAAGACGGTTTAAACATGTGtt                                                  <  1:737957/140‑1 (MQ=255)
                                                                                     atCCGGCGACATCAATATTCGCCGAACCTATAATTACGACATAATCCTTTTCGCGCATAGCACGTCCTCTGGTAAAAAATTACCAACCCTGTCAATCCACAAGAATATAGATACTGCAAAGACGGTTTAAACATGTGTTc                                                 <  2:753280/140‑1 (MQ=255)
                                                                                            gACATCAATATTCGCCGAACCTATAATTACGACATAATCCTTTTCGCGCATAGCACGTCCTCTGGTAAAAAATTACCAACCCTGTCAATCCACAAGAATATAGATACTGCAAAGACGGTTTAAACATGTGTTCATATTTa                                          <  1:583386/140‑1 (MQ=255)
                                                                                            gACATCAATATTCGCCGAACCTATAATTACGACATAATCCTTTTCGCGCATAGCACGTCCTCTGGTAAAAAATTACCAACCCTGTCAATCCACAAGAATATAGATACTGCAAAGACGGTTTAAACATGTGTTCATATTTa                                          <  1:1829100/140‑1 (MQ=255)
                                                                                            gACATCAATATTCGCCGAACCTATAATTACGACATAATCCTTTTCGCGCATAGCACGTCCTCTGGTAAAAAATTACCAACCCTGTCAATCCACAAGAATATAGATACTGCAAAGACGGTTTAAACATGTGTTCATATTTa                                          <  1:535891/140‑1 (MQ=255)
                                                                                              cATCAATATTCGCCGAACCTATAATTACGACATAATCCTTTTCGCGCATAGCACGTCCTCTGGTAAAAAATTACCAACCCTGTCAATCCACAAGAATATAGATACTGCAAAGACGGTTTAAACATGTGTTCATATTTaaa                                        >  1:577884/1‑140 (MQ=255)
                                                                                              cATCAATATTCGCCGAACCTATAATTACGACATAATCCTTTTCGCGCATAGCACGTCCTCTGGTAAAAAATTACCAACCCTGTCAATCCACAAGAATATAGATACTGCAAAGACGGTTTAAACATGTGTTCATATTTaaa                                        <  1:291267/140‑1 (MQ=255)
                                                                                                tCAATATTCGCCGAACCTATAATTACGACATAATCCTTTTCGCGCATAGCACGTCCTCTGGTAAAAAATTACCAACCCTGTCAATCCACAAGAATATAGATACTGCAAAGACGGTTTAAACATGTGTTCATATTTAAACa                                      <  2:194076/140‑1 (MQ=255)
                                                                                                tCAATATTCGCCGAACCTATAATTACGACATAATCCTTTTCGCGCATAGCACGTCCTCTGGTAAAAAATTACCAACCCTGTCAATCCACAAGAATATAGATACTGCAAAGACGGTTTAAACATGTGTTCATATTTAAACa                                      <  2:194226/140‑1 (MQ=255)
                                                                                                tCAATATTCGCCGAACCTATAATTACGACATAATCCTTTTCGCGCATAGCACGTCCTCTGGTAAAAAATTACCAACCCTGTCAATCCACAAGAATATAGATACTGCAAAGACGGTTTAAACATGTGTTCATATTTAAACa                                      <  2:194232/140‑1 (MQ=255)
                                                                                                  aaTATTCGCCGAACCTATAATTACGACATAATCCTTTTCGCGCATAGCACGTCCTCTGGTAAAAAATTACCAACCCTGTCAATCCACAAGAATATAGATACTGCAAAGACGGTTTAAACATGTGTTCATATTTAAACAtg                                    <  1:1708392/140‑1 (MQ=255)
                                                                                                  aaTATTCGCCGAACCTATAATTACGACATAATCCTTTTCGCGCATAGCACGTCCTCTGGTAAAAAATTACCAACCCTGTCAATCCACAAGAATATAGATACTGCAAAGACGGTTTAAACATGTGTTCATATTTAAACAtg                                    <  1:1446436/140‑1 (MQ=255)
                                                                                                    tatTCGCCGAACCTATAATTACGACATAATCCTTTTCGCGCATAGCACGTCCTCTGGTAAAAAATTACCAACCCTGTCAATCCACAAGAATATAGATACTGCAAAGACGGTTTAAACATGTGTTCATATTTAAACAtgtg                                  <  1:252696/140‑1 (MQ=255)
                                                                                                      ttCGCCGAACCTATAATTACGACATAATCCTTTTCGCGCATAGCACGTCCTCTGGTAAAAAATTACCAACCCTGTCAATCCACAAGAATATAGATACTGCAAAGACGGTTTAAACATGTGTTCATATTTAAACATGTGCt                                <  1:525374/140‑1 (MQ=255)
                                                                                                      ttCGCCGAACCTATAATTACGACATAATCCTTTTCGCGCATAGCACGTCCTCTGGTAAAAAATTACCAACCCTGTCAATCCACAAGAATATAGATACTGCAAAGACGGTTTAAACATGTGTTCATATTTAAACATGTGCt                                <  1:106704/140‑1 (MQ=255)
                                                                                                        cGCCGAACCTATAATTACGACATAATCCTTTTCGCGCATAGCACGTCCTCTGGTAAAAAATTACCAACCCTGTCAATCCACAAGAATATAGATACTGCAAAGACGGTTTAAACATGTGTTCATATTTAAACATGTGCTCa                              <  2:554308/140‑1 (MQ=255)
                                                                                                         gCCGAACCTATAATTACGACATAATCCTTTTCGCGCATAGCACGTCCTCTGGTAAAAAATTACCAACCCTGTCAATCCACAAGAATATAGATACTGCAAAGACGGTTTAAACATGTGTTCATATTTAAACATGTGCTCAt                             <  1:1908875/140‑1 (MQ=255)
                                                                                                                cTATAATTACGACATAATCCTTTTCGCGCATAGCACGTCCTCTGGTAAAAAATTACCAACCCTGTCAATCCACAAGAATATAGATACTGCAAAGACGGTTTAAACATGTGTTCATATTTAAACATGTGCTCATAGTATGc                      <  1:635383/140‑1 (MQ=255)
                                                                                                                     aTTACGACATAATCCTTTTCGCGCATAGCACGTCCTCTGGTAAAAAATTACCAACCCTGTCAATCCACAAGAATATAGATACTGCACAGACGGTTTAAACATGTGTTCATATTTAAACATGTGCTCATAGTATGCATGtt                 <  1:264835/140‑1 (MQ=255)
                                                                                                                      ttACGACATAATCCTTTTCGCGCATAGCACGTCCTCTGGTAAAAAATTACCAACCCTGTCAATCCACAAGAATATAGATACTGCAAAGACGGTTTAAACATGTGTTCATATTTAAACATGTGCTCATAGTATGCATGttt                >  1:572134/1‑140 (MQ=255)
                                                                                                                      ttACGACATAATCCTTTTCGCGCATAGCACGTCCTCTGGTAAAAAATTACCAACCCTGTCAATCCACAAGAATATAGATACTGCAAAGACGGTTTAAACATGTGTTCATATTTAAACATGTGCTCATAGTATGCATGttt                >  1:572215/1‑140 (MQ=255)
                                                                                                                       tACGACATAATCCTTTTCGCGCATAGCACGTCCTCTGGTAAAAAATTACCAACCCTGTCAATCCACAAGAATATAGATACTGCAAAGACGGTTTAAACATGTGTTCATATTTAAACATGTGCTCATAGTATGCATGtttt               >  1:684938/1‑140 (MQ=255)
                                                                                                                        aCGACATAATCCTTTTCGCGCATAGCACGTCCTCTGGTAAAAAATTACCAACCCTGTCAATCCACAAGAATATAGATACTGCAAAGACGGTTTa                                                            >  2:580854/1‑94 (MQ=255)
                                                                                                                        aCGACATAATCCTTTTCGCGCATAGCACGTCCTCTGGTAAAAAATTACCAACCCTGTCAATCCACAAGAATATAGATACTGCAAAGACGGTTTa                                                            <  1:580854/94‑1 (MQ=255)
                                                                                                                          gACATAATCCTTTTCGCGCAAAGCACGTCCTCTGGTAAAAAATTACCAACCCTGTCAAGCCACATGTATAGAGATACTGCGAAGACGGGTGGAACATGTGTTCATAGTTGAACTTGGGCTCGTAGTGTGCATGATTTACg            >  2:493058/1‑140 (MQ=255)
                                                                                                                            cATAATCCTTTTCGCGCATAGCACGTCCTCTGGTAAAAAATTACCAACCCTGTCAATCCACAAGAATATAGATACTGCAAAGACGGTTTAAACATGTGTTCATATTTAAACATGTGCTCATAGTATGCATGTTTTACgta          <  1:1167287/140‑1 (MQ=255)
                                                                                                                            cATAATCCCTTTCGCGCATAGCACGTCCTCTGGTAAAAAATTACCAACCCTGTCAATCCACAAGAATATAGATACTGCAAAGACGGTTTAAACATGTGTTCATATTTAAACATGTGCTCATAGTATGCATGTTTTACgta          <  2:494885/140‑1 (MQ=255)
                                                                                                                            cATAATCCCTTTCGCGCATAGCACGTCCTCTGGTAAAAAATTACCAACCCTGTCAATCCACAAGAATATAGATACTGCAAAGACGGTTTAAACATGTGTTCATATTTAAACATGTGCTCATAGTATGCATGTTTTACgta          <  2:494835/140‑1 (MQ=255)
                                                                                                                                   cTTTTCGCGCATAGCACGTCCTCTGGTAAAAAATTACCAACCCTGTCAATCCACAAGAATATAGATACTGCAAAGACGGTTTAAACATGTGTTCATATTTAAACATGTGCTCATAGTATGCATGTTTTACgtaaagtaaa   >  1:4540/1‑140 (MQ=255)
                                                                                                                                    ttttCGCGCATAGCACGTCCTCTGGTAAAAAATTACCAACCCTGTCAATCCACAAGAATATAGATACTGCAAAGACGGTTTAAACATGTGTTCATATTTAAACATGTGCTCATAGTATGCATGTTTTACGTAAAGTAAAt  <  2:949218/140‑1 (MQ=255)
                                                                                                                                    |                                                                                                                                           
TGTGCAATATTGCGCCCTACTCCACCAGGCGTAAATTTTATTTTACCTGGATTTGAATCCGCATAATTTAATGATTCATGTGAATATCCGGCGACATCAATATTCGCCGAACCTATAATTACGACATAATCCTTTTCGCGCATAGCACGTCCTCTGGTAAAAAATTACCAACCCTGTCAATCCACAAGAATATAGATACTGCAAAGACGGTTTAAACATGTGTTCATATTTAAACATGTGCTCATAGTATGCATGTTTTACGTAAAGTAAAT  >  NC_000913/2259154‑2259425

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