Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,574,101 T→C pseudogene (183/2937 nt) yjiV → pseudogene, conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,574,1010TC100.0% 345.7 / NA 95pseudogene (183/2937 nt)yjiVpseudogene, conserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (42/53);  total (42/53)

GATCGATGTCTGTTCAGGGACTCAACCTAATGCAAAGGATCCCGCATTTCTTAAAGTCAGGGCGCATATCTTCCAGCGTAATACCTAGGGAATATGGGCTTGTGTTGATAAAGATTGCAGATTAAAGCACGGTACACCGCTCGACAAAGGCTGGCCCTTTGGCTATGTGTATGTGAACCAGCGACAAAATTGTGACTGTGGAAGCCCTGTATACGAAGTTGCATTCTGTAATGATTGTAATGAGCCTC  >  NC_000913/4573969‑4574216
                                                                                                                                    |                                                                                                                   
gATCGATGTCTGTTCAGGGACTCAACCTAATGCAAAGGATCCCGCATTTCTTAAAGTCAGGGCGCATATCTTCCAGCGTAATACCTAGGGAATATGGGCTCGTGTTGATAAAGATTGCAGATTAAAGCACGGCACACCg                                                                                                               >  1:151900/1‑139 (MQ=255)
 aTCGATGTCTGTTCAGGGACTCAACCTAATGCAAAGGATCCCGCATTTCTTAAAGTCAGGGCGCATATCTTCCAGCGTAATACCTAGGGAATATGGGCTTGTGTTGATAAAGATTGCAGATTAAAGCACGGCACACCGc                                                                                                              <  1:1332076/139‑1 (MQ=255)
     aTGTCTGTTCAGGGACTCAACCTAATGCAAAGGATCCCGCATTTCTTAAAGTCAGGGCGCATATCTTCCAGCGTAATACCTAGGGAATATGGGCTTGTGTTGATAAAGATTGCAGATTAAAGCACGGCACACCGCTCGa                                                                                                          >  1:828228/1‑139 (MQ=255)
     aTGTCTGTTCAGGGACTCAACCTAATGCAAAGGATCCCGCATTTCTTAAAGTCAGGGCGCATATCTTCCAGCGTAATACCTAGGGAATATGGGCTTGTGTTGATAAAGATTGCAGATTAAAACACGGCACACCGCTCGa                                                                                                          >  2:1323016/1‑139 (MQ=255)
     aTGTCTGTTCAGGGACTCAACCTAATGCAAAGGACCCCGCATTTCTTAAAGTCAGGGCGCATATCTTCCAGCGTAATACCTAGGGAATATGGGCTTGTGTTGATAAAGATTGCAGATTAAAGCACGGCACACCGCTCGa                                                                                                          >  2:976655/1‑139 (MQ=255)
             tCAGGTACTCAACCTAATGCAAAGGATCCCGCATTTCTTAAAGTCAGGGCGCATATCTTCCAGCGTAATACCTAGGGAATATGGGCTTGTGTTGATAAAGATTGCAGATTAAAGCACGGCACACCGCTCGACAAAGGCt                                                                                                  <  2:900586/139‑1 (MQ=255)
                    cTCAACCTAATGCAAAGGATCCCGCATTTCTTAAAGTCAGGGCGCATATCTTCCAGCGTAATACCTAGGGAATATGGGCTTGTGTTGATAAAGATTGCAGATTAAAGCACGGCACACCGCTCGACAAAGGCTGGCCCtt                                                                                           >  2:464375/1‑139 (MQ=255)
                      cAACCTAATGCAAAGGATCCCGCATTTCTTAAAGTCAGGGCGCATATCTTCCAGCGTAATACCTAGGGAATATGGGCTTGTGTTGATAAAGATTGCAGATTAAAGCACGGCACACCGCTCGACAAAGGCTGGCCCTTTg                                                                                         >  2:1131266/1‑139 (MQ=255)
                      cAACCTAATGCAAAGGATCCCGCATTTCTTAAAGTCAGGGCGCATATCTTCCAGCGTAATACCTAGGGAATATGGGCTTGTGTTGATAAAGATTGCAGATTAAAGCACGGCACACCGCTCGACAAAGGCTGGCCCTTTg                                                                                         >  2:853638/1‑139 (MQ=255)
                       aaCCTAATGCAAAGGATCCCGCATTTCTTAAAGTCAGGGCGCATATCTTCCAGCGTAATACCTAGGGAATATGGGCTTGTGTTGATAAAGATTGCAGATTAAAGCACGGCACACCGCTCGACAAAGGCTGGCCCTTTgg                                                                                        <  1:1231307/139‑1 (MQ=255)
                       aaCCTAATGCAAAGGATCCCGCATTTCTTAAAGTCAGGGCGCATATCTTCCAGCGTAATACCTAGGGAATATGGGCTTGTGTTGATAAAGATTGCAGATTAAAGCACGGCACACCGCTCGACAAAGGCTGGCCCTTTgg                                                                                        <  2:477738/139‑1 (MQ=255)
                        aCCTAATGCAAAGGATCCCGCATTTCTTAAAGTCAGGGCGCATATCTTCCAGCGTAATACCTAGGGAATATGGGCTTGTGTTGATAAAGATTGCAGATTAAAGCACGGCACACCGCTCGACAAAGGCTGGCCCTTTGGc                                                                                       <  2:50975/139‑1 (MQ=255)
                         ccTAATGCGAAGGATCCCGCATTTCTTAAAGTCAGGGCGCATATCTTCCAGCGTAATACCTAGGGAATATGGGCTTGTGTTGATAAAGATTGCAGATTAAAGCACGGCACACCGCTCGACAAAGGCTGGCCCTTTGGCt                                                                                      <  2:917481/139‑1 (MQ=255)
                         ccTAATGCAAAGGATCCCGCATTTCTTAAAGTCAGGGCGCATATCTTCCAGCGTAATACCTAGGGAATATGGGCTTGTGTTGATAAAGATTGCAGATTAAAGCACGGCACACCGCTCGACAAAGGCTGGCCCTTTGGCt                                                                                      <  1:386672/139‑1 (MQ=255)
                         ccTAATGCAAAGGATCCCGCATTTCTTAAAGTCAGGGCGCATATCTTCCAGCGTAATACCTAGGGAATATGGGCTTGTGTTGATAAAGATTGCAGATTAAAGCACGGCACACCGCTCGACAAAGGCTGGCCCTTTGGCt                                                                                      <  1:1492259/139‑1 (MQ=255)
                            aaTGCAAAGGATCCCGCATTTCTTAAAGTCAGGGCGCATATCTTCCAGCGTAATACCTAGGGAATATGGGCTTGTGTTGATAAAGATTGCAGATTAAAGCACGGCACACCGCTCGACAAAGGCTGGCCCTTTGGCtatg                                                                                   <  1:1080390/139‑1 (MQ=255)
                             aTGCAAAGGATCCCGCATTTCTTAAAGTCAGGGCGCATATCTTCCAGCGTAATACCTAGGGAATATGGGCTTGTGTTGATAAAGATTGCAGATTAAAGCACGGCACACCGCTCGAc                                                                                                         >  2:58990/1‑116 (MQ=255)
                             aTGCAAAGGATCCCGCATTTCTTAAAGTCAGGGCGCATATCTTCCAGCGTAATACCTAGGGAATATGGGCTTGTGTTGATAAAGATTGCAGATTAAAGCACGGCACACCGCTCGAc                                                                                                         <  1:58990/116‑1 (MQ=255)
                             aTGCAAAGGATCCCGCATTTCTTAAAGTCAGGGCGCATATCTTCCAGCGTAATACCTAGGGAATATGGGCTTGTGTTGATAAAGATTGCAGATTAAAGCACGGCACACCGCTCGACAAAGGCTGGCCCTTTGGCt                                                                                      <  2:1058318/135‑1 (MQ=255)
                             aTGCAAAGGATCCCGCATTTCTTAAAGTCAGGGCGCATATCTTCCAGCGTAATACCTAGGGAATATGGGCTTGTGTTGATAAAGATTGCAGATTAAAGCACGGCACACCGCTCGACAAAGGCTGGCCCTTTGGCt                                                                                      >  1:1058318/1‑135 (MQ=255)
                                 aaaGGATCCCGCATTTCTTAAAGTCAGGGCGCATATCTTCCAGCGTAATACCTAGGGAATATGGGCTTGTGTTGATAAAGATTGCAGATTAAAGCACGGCACACCGCTCGACAAAGGCTGGCCCTTTGGCtatgtgtat                                                                              >  2:1167696/1‑139 (MQ=255)
                                    ggATCCCGCATTTCTTAAAGTCAGGGCGCATATCTTCCAGCGTAATACCTAGGGAATATGGGCTTGTGTTGATAAAGATTGCAGATTAAAGCACGGCACACCGCTCGACAAAGGCTGGCCCt                                                                                            <  2:1417374/122‑1 (MQ=255)
                                    ggATCCCGCATTTCTTAAAGTCAGGGCGCATATCTTCCAGCGTAATACCTAGGGAATATGGGCTTGTGTTGATAAAGATTGCAGATTAAAGCACGGCACACCGCTCGACAAAGGCTGGCCCt                                                                                            >  1:1417374/1‑122 (MQ=255)
                                          cGCATTTCTTAACGTCAGGGCGCATATCTTCCAGCGTAATACCTAGGGAATAAGGGTGTGTGTTGATAAAGATTGCAGATTAAAGCACGGCACACCGCTCGACAAAGGCTGGCCCTTGGGCTATGCGTATGTGAACCAg                                                                     >  2:1297666/1‑139 (MQ=255)
                                               ttCTTAAAGTCAGGGCGCATATCTTCCAGCGTAATACCTAGGGAATATGGGCTTGTGTTGATAAAGATTGCAGATTAAAGCACGGCACACCGCTCGACAAAGGCTGGCCCTTTGGCTATGTGTATGTGAACCAGCGACa                                                                >  2:1439923/1‑139 (MQ=255)
                                                   tAAAGTCAGGGCGCATATCTTCCAGCGTAATACCTAGGGAATATGGGCTTGTGTTGATAAAGATTGCAGATTAAAGCACGGCACACCGCTCGACAAAGGCTGGCCCTTTGGCTATGTGTATGTGAACCAGCGACAAAAt                                                            <  1:1330200/139‑1 (MQ=255)
                                                    aaaGTCAGGGCGCATATCTTCCAGCGTAATACCTAGGGAATATGGGCTTGTGTTGATAAAGATTGCAGATTAAAGCACGGCACACCGCTCGACAAAGGCTGGCCCTTTGGCTATGTGTATGTGAACCAGCGACAAAAtt                                                           >  1:1162062/1‑139 (MQ=255)
                                                    aaaGTCAGGGCGCATATCTCCCAGCGTAATACCTAGGGAATATGGGCTTGTGTTGATAAAGATTGCAGATTAAAGCACGGCACACCGCTCGACAAAGGCTGGCCCTTTGGCTATGTGTATGTGAACCAGCGACAAAAtt                                                           <  2:1033136/139‑1 (MQ=255)
                                                           gggCGCATATCTTCCAGCGTAATACCTAGGGAATATGGGCTTGTGTTGATAAAGATTGCAGATTAAAGCACGGCACACCGCTCGACAAAGGCTGGCCCTTTGGCTATGTGTATGTGAACCAGCGACAAAATTGTGACtg                                                    >  2:554486/1‑139 (MQ=255)
                                                            ggCGCATATCTTCCAGCGTAATACCTAGGGAATATGGGCTTGTGTTGATAAAGATTGCAGATTAAAGCACGGCACAcc                                                                                                                >  1:133388/1‑78 (MQ=255)
                                                            ggCGCATATCTTCCAGCGTAATACCTAGGGAATATGGGCTTGTGTTGATAAAGATTGCAGATTAAAGCACGGCACAcc                                                                                                                <  2:133388/78‑1 (MQ=255)
                                                            ggCGCATATCTTCCAGCGTAATACCTAGGGAATATGGGCTTGTGTTGATAAAGATTGCAGATTAAAGCACGGCACACCGCTCGACAAAGGCTGGCCCTTTGGCt                                                                                      >  2:774509/1‑104 (MQ=255)
                                                            ggCGCATATCTTCCAGCGTAATACCTAGGGAATATGGGCTTGTGTTGATAAAGATTGCAGATTAAAGCACGGCACACCGCTCGACAAAGGCTGGCCCTTTGGCt                                                                                      <  1:774509/104‑1 (MQ=255)
                                                            ggCGCATATCTTCCAGCGTAATACCTAGGGAATATGGGCTTGTGTTGATAAAGATTGCAGATTAAAGCACGGCACACCGCTCGACAAAGGCTGGCCCTTTGGCTATGTGTATGTGAACCAGCGACAAAATTGTGACtgt                                                   >  2:751514/1‑139 (MQ=255)
                                                                cATATCTTCCAGCGTAATACCTAGGGAATATGGGCTTGTGTTGATAAAGATTGCAGATTAAAGCACGGCACACCGCTCGACAAAGGCTGGCCCTTTGGCTATGTGTATGTGAACCAGCGACAAAATTGTGACTGTGGaa                                               <  2:959187/139‑1 (MQ=255)
                                                                   atCTTCCAGCGTAATACCTAGGGAATATGGGCTTGTGTTGATAAAGATTGCAGATTAAAGCACGGCACACCGCTCGACAAAGGCTGGCCCTTTGGCTATGTGTATGTGAACCAGCGACAAAATTGTGACTGTGGAAGcc                                            <  2:978652/139‑1 (MQ=255)
                                                                   atCTTCCAGCGTAATACCTAGGGAATATGGGCTTGTGTTGATAAAGATTGCAGATTAAAGCACGGCACACCGCTCGACAAAGGCTGGCCCTTTGGCTATGTGTATGTGAACCAGCGACAAAATTGTGACTGTGGAAGcc                                            <  2:986246/139‑1 (MQ=255)
                                                                   atCTTCCAGCGTAATACCTAGGGAATATGGGCTTGTGTTGATAAAGATTGCAGATTAAAGCACGGCACACCGCTCGACAAAGGCTGGCCCTTTGGCTATGTGTATGTGAACCAGCGACAAAATTGTGACTGTGGAAGcc                                            <  1:853638/139‑1 (MQ=255)
                                                                   atCTTCCAGCGTAATACCTAGGGAATATGGGCTTGTGTTGATAAAGATTGCAGATTAAAGCACGGCACACCGCTCGACAAAGGCTGGCCCTTTGGCTATGTGTATGTGAACCAGCGACAAAATTGTGACTGTGGAAGcc                                            <  2:848129/139‑1 (MQ=255)
                                                                   atCTTCCAGCGTAATACCTAGGGAATATGGGCTTGTGTTGATAAAGATTGCAGATTAAAGCACGGCACACCGCTCGACAAAGGCTGGCCCTTTGGCTATGTGTATGTGAACCAGCGACAAAATTGTGACTGTGGAAGcc                                            <  1:841674/139‑1 (MQ=255)
                                                                   atCTTCCAGCGTAATACCTAGGGAATATGGGCTTGTGTTGATAAAGATTGCAGATTAAAGCACGGCACACCGCTCGACAAAGGCTGGCCCTTTGGCTATGTGTATGTGAACCAGCGACAAAATTGTGACTGTGGAAGcc                                            <  1:1131266/139‑1 (MQ=255)
                                                                   atCTTCCAGCGTAATACCTAGGGAATATGGGCTTGTGTTGATAAAGATTGCAGATTAAAGCACGGCACACCGCTCGACAAAGGCTGGCCCTTTGGCTATGTGTATGTGAACCAGCGACAAAATTGTGACTGTGGAAGcc                                            <  1:751514/139‑1 (MQ=255)
                                                                   atCTTCCAGCGTAATACCTAGGGAATATGGGCTTGTGTTGATAAAGATTGCAGATTAAAGCACGGCACACCGCTCGACAAAGGCTGGCCCTTTGGCTATGTGTATGTGAACCAGCGACAAAATTGTGACTGTGGAAGcc                                            <  1:1167696/139‑1 (MQ=255)
                                                                   atCTTCCAGCGTAATACCTAGGGAATATGGGCTTGTGTTGATAAAGATTGCAGATTAAAACACGGCACACCGCTCGACAAAGGCTGGCCCTTTGGCTATGTGTATGTGAACCAGCGACAAAATTGTGACTGTGGAAGcc                                            <  1:1323016/139‑1 (MQ=255)
                                                                       tCCAGCGTAATACCTAGGGAATATGGGCTTGTGTTGATAAAGATTGCAGATTAAAGCACGGCACACCGCTCGACAAAGGCTGGCCCttt                                                                                          >  2:562191/1‑89 (MQ=255)
                                                                       tCCAGCGTAATACCTAGGGAATATGGGCTTGTGTTGATAAAGATTGCAGATTAAAGCACGGCACACCGCTCGACAAAGGCTGGCCCttt                                                                                          <  1:562191/89‑1 (MQ=255)
                                                                       tCCAGCGTAATACCTAGGGAATATGGGCTTGTGTTGATAAAGATTGCAGATTAAAGCACGGCACACCGCTCGACAAAGGCTGGCCCTTTGGCtatgt                                                                                  >  2:135605/1‑97 (MQ=255)
                                                                       tCCAGCGTAATACCTAGGGAATATGGGCTTGTGTTGATAAAGATTGCAGATTAAAGCACGGCACACCGCTCGACAAAGGCTGGCCCTTTGGCtatgt                                                                                  <  1:135605/97‑1 (MQ=255)
                                                                           gCGTAATACCTAGGGAATATGGGCTTGTGTTGATAAAGATTGCAGATTAAAGCACGGCACACCGCTCGACAAAGGCTGGCCCTTTGGCTATGTGTATGTGAACCAGCGACAAAATTGTGACTGTGGAAGCCCTGTATAc                                    >  2:1301566/1‑139 (MQ=255)
                                                                              tAATACCTAGGGAATATGGGCTTGTGTTGATAAAGATTGCAGATTAAAGCACGGCACACCGCTCGACAAAGGCTGGCCCTTTGGCTATGTGTATGTGAACCAGCGACAAAATTGTGACTGTGGAAGCCCTGTATACGaa                                 >  2:271908/1‑139 (MQ=255)
                                                                               aaTACCTAGGGAATATGGGCTTGTGTTGATAAAGATTGCAGATTAAAGCACGGCACACCGCTCGACAAAGGCTGGCCCTTTGGCTATGTGTATGTGAACCAGCGACAAAATTGTGACTGTGGAAGCCCTGTATACGAAg                                <  1:554486/139‑1 (MQ=255)
                                                                                aTACCTAGGGAATATGGGCTTGTGTTGATAAAGATTGCAGATTAAAGCACGGCACACCGCTCGACAAAGGCTGGCCCTTTGGCTATGTGTATGTGAACCAGCGACAAAATTGTGa                                                       <  1:1258080/115‑1 (MQ=255)
                                                                                aTACCTAGGGAATATGGGCTTGTGTTGATAAAGATTGCAGATTAAAGCACGGCACACCGCTCGACAAAGGCTGGCCCTTTGGCTATGTGTATGTGAACCAGCGACAAAATTGTGa                                                       >  2:1258080/1‑115 (MQ=255)
                                                                                aTACCTAGGGAATATGGGCTTGTGTTGATAAAGATTGCAGATTAAAGCACGGCACACCGCTCGACAAAGGCTGGCCCTTTGGCTATGTGTATGTGAACCAGCGACAAAATTGTGACTGTGGAAGCCCTGTATACGAAGt                               >  2:112896/1‑139 (MQ=255)
                                                                                 tACCTAGGGAATATGGGCTTGTGTTGATAAAGATTGCAGATTAAAGCACGGCACACCGCTCGACAAAGGCTGGCCCTTTgg                                                                                        >  2:1471419/1‑81 (MQ=255)
                                                                                 tACCTAGGGAATATGGGCTTGTGTTGATAAAGATTGCAGATTAAAGCACGGCACACCGCTCGACAAAGGCTGGCCCTTTgg                                                                                        <  1:1471419/81‑1 (MQ=255)
                                                                                 tACCTAGGGAATATGGGCTTGTGTTGATAAAGATTGCAGATTAAAGCACGGCACACCGCTCGACAAAGGCTGGCCCTTTGGCTATGTGTATGTGAACCAGCGACAAAATTGTGACTGTGGAAGCCCTGTATACGAAGtt                              >  2:1361851/1‑139 (MQ=255)
                                                                                    cTAGGGAATATGGGCTTGTGTTGATAAAGATTGCAGATTAAAGCACGGCACACCGCTCGACAAAGGCTGGCCCTTTGGCTATGTGTATGTGAACCAGCGACAAAATTGTGACTGTGGAAGCCCTGTATACGAAGTTGCa                           <  1:1301566/139‑1 (MQ=255)
                                                                                    cTAGGGAATATGGGCTTGTGTTGATAAAGATTGCAGATTAAAGCACGGCACACCGCTCGACAAAGGCTGGCCCTTTGGCTATGTGTATGTGAACCAGCGACAAAATTGTGACTGTGGAAGCCCTGTATACGAAGTTGCa                           <  2:1162062/139‑1 (MQ=255)
                                                                                            tatGGGCTTGTGTTGATAAAGATTGCAGATTAAAGCACGGCACACCGCTCGACAAAGGCTGGCCCTTTGGCTATGTGTATGTGAACCAGCGACAAAATTGTGACTGTGGAAGCCCTGTATACGAAGTTGCATTCTGTaa                   <  1:1361851/139‑1 (MQ=255)
                                                                                            tatGGGCTTGTGTTGATAAAGATTGCAGATTAAAGCACGGCACACCGCTCGACAAAGGCTGGCCCTTTGGCTATGTGTATGTGAACCAGCGACAAAATTGTGACTGTGGAAGCCCTGTATACGAAGTTGCATTCTGTaa                   <  2:405363/139‑1 (MQ=255)
                                                                                            tatGGGCTTGTGTTGATAAAGATTGCAGATTAAAGCACGGCACACCGCTCGACAAAGGCTGGCCCTTTGGCTATGTGTATGTGAACCAGCGACAAAATTGTGACTGTGGAAGCCCTGTATACGAAGTTGCATTCTGTaa                   <  2:828228/139‑1 (MQ=255)
                                                                                             atGGGCTTGTGTTGATAAAGATTGCAGATTAAAGCACGGCACACCGCTCGACAAAGGCTGGCCCTTTGGCTATGTGTATGTGAACCAGCGACAAAATTGTGACTGTGGAAGCCCTGTATACGAAGTTGCATTCTGTAAt                  >  2:403304/1‑139 (MQ=255)
                                                                                                ggCTTGTGTTGATAAAGATTGCAGATTAAAGCACGGCACACCGCTCGACAAAGGCTGGCCCTTTGGCTATGTGTATGTGAACCAGCGACAAAATTGTGACTGTGGAAGCCCTGTATACGAAGTTGCATTCTGTAATGAt               <  2:433970/139‑1 (MQ=255)
                                                                                                  cTTGTGTTGATAAAGATTGCAGATTAAAGCACGGCACACCGCTCGACAAAGGCTGGCCCTTTGGCTATGTGTATGTGAACCAGCGACAAAATTGTGACTGTGGAAGCCCTGTATACGAAGTTGCATTCTGTAATg                 <  1:601071/135‑1 (MQ=255)
                                                                                                  cTTGTGTTGATAAAGATTGCAGATTAAAGCACGGCACACCGCTCGACAAAGGCTGGCCCTTTGGCTATGTGTATGTGAACCAGCGACAAAATTGTGACTGTGGAAGCCCTGTATACGAAGTTGCATTCTGTAATg                 >  2:601071/1‑135 (MQ=255)
                                                                                                  cTTGTGTTGATAAAGATTGCAGATTAAAGCACGGCACACCGCTCGACAAAGGCTGGCCCTTTGGCTATGTGTATGTGAACCAGCGACAAAATTGTGACTGTGGAAGCCCTGTATACGAAGTTGCATTCTGTAATGa                >  2:601072/1‑136 (MQ=255)
                                                                                                  cTTGTGTTGATAAAGATTGCAGATTAAAGCACGGCACACCGCTCGACAAAGGCTGGCCCTTTGGCTATGTGTATGTGAACCAGCGACAAAATTGTGACTGTGGAAGCCCTGTATACGAAGTTGCATTCTGTAATGa                <  1:601072/136‑1 (MQ=255)
                                                                                                  cTTGTGTTGATAAAGATTGCAGATTAAAGCACGGCACACCGCTCGACAAAGGCTGGCCCTTTGGCTATGTGTATGTGAACCAGCGACAAAATTGTGACTGTGGAAGCCCTGTATACGAAGTTGCATTCTGTAATGATTg             >  2:1455088/1‑139 (MQ=255)
                                                                                                    tgtgTTGATAAAGATTGCAGATTAAAGCACGGCACACCGCTCGACa                                                                                                        <  2:1240890/46‑1 (MQ=255)
                                                                                                    tgtgTTGATAAAGATTGCAGATTAAAGCACGGCACACCGCTCGACa                                                                                                        >  1:1240890/1‑46 (MQ=255)
                                                                                                    tgtgTTGATAAAGATTGCAGATTAAAGCACGGCACACCGCTCGACAAAGGCt                                                                                                  >  2:252502/1‑52 (MQ=255)
                                                                                                    tgtgTTGATAAAGATTGCAGATTAAAGCACGGCACACCGCTCGACAAAGGCt                                                                                                  <  1:252502/52‑1 (MQ=255)
                                                                                                    tgtgTTGATAAAGATTGCAGATTAAAGCACGGCACACCGCTCGACAAAGGCTGGCCCTTTGGc                                                                                       >  1:614527/1‑63 (MQ=255)
                                                                                                    tgtgTTGATAAAGATTGCAGATTAAAGCACGGCACACCGCTCGACAAAGGCTGGCCCTTTGGc                                                                                       <  2:614527/63‑1 (MQ=255)
                                                                                                    tgtgTTGATAAAGATTGCAGATTAAAGCACGGCACACCGCTCGACAAAGGCTGGCCCTTTGGCTATGTGTATGTGAACCAGCGACAAAGTTGTGa                                                       >  1:956742/1‑95 (MQ=255)
                                                                                                    tgtgTTGATAAAGATTGCAGATTAAAGCACGGCACACCGCTCGACAAAGGCTGGCCCTTTGGCTATGTGTATGTGAACCAGCGACAAAATTGTGa                                                       <  2:956742/95‑1 (MQ=255)
                                                                                                    tgtgTTGATAAAGATTGCAGATTAAAGCACGGCACACCGCTCGACAAAGGCTGGCCCTTTGGCTATGTGTATGTGAACCAGCGACAAAATTGTGACTGTGGAAGCCCTGTATACGAAGTTGCATTCTGTa                    >  2:844721/1‑130 (MQ=255)
                                                                                                    tgtgTTGATAAAGATTGCAGATTAAAGCACGGCACACCGCTCGACAAAGGCTGGCCCTTTGGCTATGTGTATGTGAACCAGCGACAAAATTGTGACTGTGGAAGCCCTGTATACGAAGTTGCATTCTGTa                    <  1:844721/130‑1 (MQ=255)
                                                                                                    tgtgTTGATAAAGATTGCAGATTAAAGCACGGCACACCGCTCGACAAAGGCTGGCCCTTTGGCTATGTGTATGTGAACCAGCGACAAAATTGTGACTGTGGAAGCCCTGTATACGAAGTTGCATTCTGTAATGATTGTa           >  1:761386/1‑139 (MQ=255)
                                                                                                     gggtttgaaaaGATTGCAGATTAAAGCACGGCACACCGCTCGACAAAGGCTGGCCCTTTGGCTATGTGTATGTGAACCAGCGACAAAATTGTGACTGTGGAAGCCCTGTATACGAAGTTGCATTCTGTAATGATTGTaa          <  1:679092/131‑1 (MQ=255)
                                                                                                      tgtTGATAAAGATTGCAGATTAAAGCACGGCACACCGCTCGACAAAGGCTGGCCCTTTGGCTATGTGTATGTGAACCAGCGACAAAATTGTGACTGTGGAAGCCCTGTATACGAAGTTGCATTCTGTAATGATTGTAAt         <  1:112896/139‑1 (MQ=255)
                                                                                                             aaaGATTGCAGCTTAAAGCACGGCACACCGCTCGACAAAGGCTGGCCCTTTGGCTATGTGTATGTGAACCAGCGACAAAATTGTGACTGTGGAAGCCCTGTATACGAAGTTTCATTCTGTAATGATTGTAATGAGCCTc  <  1:976655/139‑1 (MQ=255)
                                                                                                                gATTGCAGATTAAAGCACGGCACACCGCTCGACAAAGGCTGGCCCTTTGGCtatgtgt                                                                                <  2:483960/58‑1 (MQ=255)
                                                                                                                gATTGCAGATTAAAGCACGGCACACCGCTCGACAAAGGCTGGCCCTTTGGCtatgtgt                                                                                >  1:483960/1‑58 (MQ=255)
                                                                                                                  ttGCAGATTAAAGCACGGCACACCGCTCGACAAAGGCTGGCCCTTTGGCtatg                                                                                   >  2:711163/1‑53 (MQ=255)
                                                                                                                  ttGCAGATTAAAGCACGGCACACCGCTCGACAAAGGCTGGCCCTTTGGCtatg                                                                                   <  1:711163/53‑1 (MQ=255)
                                                                                                                  ttGCAGATTAAAGCACGGCACACCGCTCGACAAAGGCTGGCCCTTTGGCTATGTGTATGTGAACCAGCGACAAAATTGTGACTGTGGAAGCCCTGTATACGAAGTTGCATTCTGTaa                   >  1:1389000/1‑117 (MQ=255)
                                                                                                                  ttGCAGATTAAAGCACGGCACACCGCTCGACAAAGGCTGGCCCTTTGGCTATGTGTATGTGAACCAGCGACAAAATTGTGACTGTGGAAGCCCTGTATACGAAGTTGCATTCTGTaa                   <  2:1389000/117‑1 (MQ=255)
                                                                                                                   tGCAGATTAAAGCACGGCACACCGCTCGACAAAGGCTGGCCCTTTGGCTATGTGTATGTGAACCAGCGACAAAATTGTGACTGTGGaa                                               >  1:300852/1‑88 (MQ=255)
                                                                                                                   tGCAGATTAAAGCACGGCACACCGCTCGACAAAGGCTGGCCCTTTGGCTATGTGTATGTGAACCAGCGACAAAATTGTGACTGTGGaa                                               <  2:300852/88‑1 (MQ=255)
                                                                                                                     cAGATTAAAGCACGGCACACCGCTCGACAAAGGCTGGCCCTTTGGCTATGTGTATGTGAACCAGCGACAAAATTGTGACTGTGGAAGcc                                            >  2:815584/1‑89 (MQ=255)
                                                                                                                     cAGATTAAAGCACGGCACACCGCTCGACAAAGGCTGGCCCTTTGGCTATGTGTATGTGAACCAGCGACAAAATTGTGACTGTGGAAGcc                                            <  1:815584/89‑1 (MQ=255)
                                                                                                                     cAGATTAAAGCACGGCACACCGCTCGACAAAGGCTGGCCCTTTGGCTATGTGTATGTGAACCAGCGACAAAATTGTGACTGTGGAAGCCCTGTATACGAAGTTGCa                           >  2:1104357/1‑106 (MQ=255)
                                                                                                                     cAGATTAAAGCACGGCACACCGCTCGACAAAGGCTGGCCCTTTGGCTATGTGTATGTGAACCAGCGACAAAATTGTGACTGTGGAAGCCCTGTATACGAAGTTGCa                           <  1:1104357/106‑1 (MQ=255)
                                                                                                                                    |                                                                                                                   
GATCGATGTCTGTTCAGGGACTCAACCTAATGCAAAGGATCCCGCATTTCTTAAAGTCAGGGCGCATATCTTCCAGCGTAATACCTAGGGAATATGGGCTTGTGTTGATAAAGATTGCAGATTAAAGCACGGTACACCGCTCGACAAAGGCTGGCCCTTTGGCTATGTGTATGTGAACCAGCGACAAAATTGTGACTGTGGAAGCCCTGTATACGAAGTTGCATTCTGTAATGATTGTAATGAGCCTC  >  NC_000913/4573969‑4574216

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑