Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,138,902 C→A T59K (ACG→AAG)  yghU → disulfide reductase/organic hydroperoxide reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,138,9020CA100.0% 244.1 / NA 69T59K (ACG→AAG) yghUdisulfide reductase/organic hydroperoxide reductase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (28/41);  total (28/41)

GGGATAAATCCGCTGGCGGCGCGTTCGCCAATATCAATCGCCCGGTTTCTGGTCCGACGCATGAAAAAACGCTGCCCGTTGGCAAACACCCATTGCAACTTTATTCGCTGGGAACGCCGAACGGTCAGAAAGTAACGATTATGCTTGAGGAGCTGCTGGCGCTGGGCGTTACTGGTGCAGAGTACGACGCCTGGCTGATTCGTATTGGCGATGGCGATCAATTCTCCAGCGGCTTTGTCGAAGTGAACCCAAACTCGAAGATCCCGGCG  >  NC_000913/3138767‑3139035
                                                                                                                                       |                                                                                                                                     
gggATAAATCCGCTGGCGGCGCGTTCGCCAATATCAATCGCCCGGTTTCTGGTCCGACGCATGAAAAAACGCTGCCCGTTGGCAAACACCCATTGCAACTTTATTCGCTGGGAACGCCGAACGGTCAGAAAGTAAAGAtt                                                                                                                                   >  1:4331/1‑140 (MQ=255)
gggATAAATCCGCTGGCGGCGCGTTCGCCAATATCAATCGCCCGGTTTCTGGTCCGACGCATGAAAAAACGCTGCCCGTTGGCAAACACCCATTGCAACTTTATTCGCTGGGAACGCCGAACGGTCAGAAAGTAAAGAtt                                                                                                                                   >  2:549227/1‑140 (MQ=255)
gggATAAATCCGCTGGCGGCGCGTTCGCCAATATCAATCGCCCGGTTTCTGGTCCGACGCATGAAAAAACGCTGCCCGTTGGCAAACACCCATTGCAACTTTATTCGCTGGGAACGCCGAACGGTCAGAAAGTAAAGAtt                                                                                                                                   >  2:1628791/1‑140 (MQ=255)
      aaTCCGCTGGCGGCGCGTTCGCCAATATCAATCGCCCGGTTTCTGGTCCGACGCATGAAAAAACGCTGCCCGTTGGCAAACACCCATTGCAACTTTATTCGCTGGGAACGCCGAACGGTCAGAAAGTAAAGAt                                                                                                                                    >  2:392526/1‑133 (MQ=255)
      aaTCCGCTGGCGGCGCGTTCGCCAATATCAATCGCCCGGTTTCTGGTCCGACGCATGAAAAAACGCTGCCCGTTGGCAAACACCCATTGCAACTTTATTCGCTGGGAACGCCGAACGGTCAGAAAGTAAAGAt                                                                                                                                    <  1:392526/133‑1 (MQ=255)
        tCCGCTGGCGGCGCGTTCGCCAATATCAATCGCCCGGTTTCTGGTCCGACGCATGAAAAAACGCTGCCCGTTGGCAAACACCCATTGCAACTTTATTCGCTGGGAGCGCCGAACGGTCAGAAAGTAAAGATTATGCTTga                                                                                                                           <  1:257084/140‑1 (MQ=255)
            cTGGCGGCGCGTTCGCCAATATCAATCGCCCGGTTTCTGGTCCGACGCATGAAAAAACGCTGCCCGTTGGCAAACACCCATTGCAACTTTATTCGCTGGGAACGCCGAACGGTCAGAAAGTAAAGATTATGCTTgaggag                                                                                                                       >  2:1423445/1‑140 (MQ=255)
               gcggcgCGTTCGCCAATATCAATCGCCCGGTTTCTGGTCCGACGCATGAAAAAACGCTGCCCGTTGGCAAACACCCATTGCAACTTTATTCGCTGGGAACGCCGAACGGTCAGAAAGTAAAGATTATGCTTGAGGAgctg                                                                                                                    <  2:865790/140‑1 (MQ=255)
                 ggcgCGTTCGCCAATATCAATCGCCCGGTTTCTGGTCCGACGCATGAAAAAACGCTGCCCGTTGGCAAACACCCATTGCAACTTTATTCGCTGGGAACGCCGAACGGTCAGAAAGTAAAGATTATGCTTGAGGAgctgct                                                                                                                  <  1:907864/140‑1 (MQ=255)
                     cgTTCGCCAATATCAATCGCCCGGTTTCTGGTCCGACGCATGAAAAAACGCTGCCCGTTGGCAAACACCCATTGCAACTTTATTCGCTGGGAACGCCGAACGGTCAGAAAGTAAAGATTATGCTTGAGGAGCTGCTGgcg                                                                                                              <  2:1042331/140‑1 (MQ=255)
                     cgTTCGCCAATATCAATCGCCCGGTTTCTGGTCCGACGCATGAAAAAACGCTGCCCGTTGGCAAACACCCATTGCAACTTTATTCGCTGGGAACGCCGAACGGTCAGAAAGTAAAGATTATGCTTGAGGAGCTGCTGgcg                                                                                                              <  1:431730/140‑1 (MQ=255)
                           ccAATATCAATCGCCCGGTTTCTGGTCCGACGCATGAAAAAACGCTGCCCGTTGGCAAACACCCATTGCAACTTTATTCGCTGGGAACGCCGAACGGTCAGAAAGTAAAGATTATGCTTGAGGAGCTGCTGGCGCTGGGc                                                                                                        <  2:1722089/140‑1 (MQ=255)
                             aaTATCAATCGCCCGGTTTCTGGTCCGACGCATGAAAGAACGCTGCCCGTTGGCAAACACCCATTGCAACTTTATTCGCTGGGAACGCCGAACGGTCAGAAAGTAAAGATTATGCTTGAGGAGCTGCTGGCGCTGGGCGt                                                                                                      <  2:1450319/140‑1 (MQ=255)
                             aaTATCAATCGCCCGGTTTCTGGTCCGACGCATGAAAAAACGCTGCCCGTTGGCAAACACCCATTGCAACTTTATTCGCTGGGAACGCCGAACGGTCAGAAAGTAAAGATTATGCTTGAGGAGCTGCTGGCGCTGGGCGt                                                                                                      <  2:1478228/140‑1 (MQ=255)
                                atCAATCGCCCGGTTTCTGGTCCGACGCATGAAAAAACGCTGCCCGTTGGCAAACACCCATTGCAACTTTATTCGCTGGGAACGCCGAACGGTCAGAAAGTAAAGATTATGCTTGAGGAGCTGCTGGCGCTGGGCGTTAc                                                                                                   >  2:1006153/1‑140 (MQ=255)
                                    aTCGCCCGGTTTCTGGTCCGACGCATGAAAAAACGCTGCCCGTTGGCAAACACCCATTGCAACTTTATTCGCTGGGAACGCCGAACGGTCAGAAAGTAAAGATTATGCTTGAGGAGCTGCTGGCGCTGGGCGTTACTGGt                                                                                               >  2:1007967/1‑140 (MQ=255)
                                    aTCGCCCGGTTTCTGGTCCGACGCATGAAAAAACGCTGCCCGTTGGCAAACACCCATTGCAACTTTATTCGCTGGGAACGCCGAACGGTCAGAAAGTAAAGATTATGCTTGAGGAGCTGCTGGCGCTGGGCGTTACTGGt                                                                                               <  1:39981/140‑1 (MQ=255)
                                     tCGCCCGGTTTCTGGTCCGACGCATGAAAAAACGCTGCCCGTTGGCAAACACCCATTGCAACTTTATTCGCTGGGAACGCCGAACGGTCAGAAAGTAAAGATTATGCTTGAGGAGCTGCTGGCGCTGGGCGTTCCTGGTg                                                                                              <  2:246223/140‑1 (MQ=255)
                                     tCGCCCGGTTTCTGGTCCGACGCATGAAAAAACGCTGCCCGTTGGCAAACACCCATTGCAACTTTATTCGCTGGGAACGCCGAACGGTCAGAAAGTAAAGATTATGCTTGAGGAGCTGCTGGCGCTGGGCGTTACTGGTg                                                                                              <  2:1451508/140‑1 (MQ=255)
                                     tCGCCCGGTTTCTGGTCCGACGCATGAAAAAACGCTGCCCGTTGGCAAACACCCATTGCAACTTTATTCGCTGGGAACGCCGAACGGTCAGAAAGTAAAGATTATGCTTGAGGAGCTGCTGGCGCTGGGCGTTACTGGTg                                                                                              <  1:792490/140‑1 (MQ=255)
                                       gCCCGGTTTCTGGTCCGACGCATGAAAAAACGCTGCCCGTTGGCAAACACCCATTGCAACTTTATTCGCTGGGAACGCCGAACGGTCAGAAAGTAAAGATTATGCTTGAGGAGCTGCTGGCGCTGGGCGTTACTGGTGCa                                                                                            <  1:1220440/140‑1 (MQ=255)
                                        cccGGTTTCTGGTCCGACGCATGAAAAAACGCTGCCCGTTGGCAAACACCCATTGCAACTTTATTCGCTGGGAACGCCGAACGGTCAGAAAGTAAAGATTATGCTTGAGGAGCTGCTGGCGCTGGGCGTTACTGGTGCag                                                                                           <  2:1073444/140‑1 (MQ=255)
                                        cccGGTTTCTGGTCCGACGCATGAAAAAACGCTGCCCGTTGGCAAACACCCATTGCAACTTTATTCGCTGGGAACGCCGAACGGTCAGAAAGTAAAGATTATGCTTGAGGAGCTGCTGGCGCTGGGCGTTACTGGTGCag                                                                                           <  2:870947/140‑1 (MQ=255)
                                         ccGGTTTCTGGTCCGACGCATGAAAAAACGCTGCCCGTTGGCAAACACCCATTGCAACTTTATTCGCTGGGAACGCCGAACGGTCAGAAAGTAAAGATTATGCTTGAGGAGCTGCTGGCGCTGGGCGTTACTGGTGCaga                                                                                          >  2:538810/1‑140 (MQ=255)
                                         ccGGTTTCTGGTCCGACGCATGAAAAAACGCTGCCCGTTGGCAAACACCCATTGCAACTTTATTCGCTGGGAACGCCGAACGGTCAGAAAGTAAAGATTATGCTTGAGGAGCTGCTGGCGCTGGGCGTTACTGGTGCaga                                                                                          >  2:1512446/1‑140 (MQ=255)
                                         ccGGTTTCTGGTCCGACGCATGAAAAAACGCTGCCCGTTGGCAAACACCCATTGCAACTTTATTCGCTGGGAACGCCGAACGGTCAGAAAGTAAAGATTATGCTTGAGGAGCTGCTGGCGCTGGGCGTTACTGGTGCaga                                                                                          >  2:1346836/1‑140 (MQ=255)
                                          cGGTTTCTGGTCCGACGCATGAAAAAACGCTGCCCGTTGGCAAACACCCATTGCAACTTTATTCGCTGGGAACGCCGAACGGTCAGAAAGTAAAGATTATGCTTGAGGAGCTGCTGGCGCTGGGCGTTACTGGTGCagag                                                                                         >  1:1127436/1‑140 (MQ=255)
                                            gTTTCTGGTCCGACGCATGAAAAAACGCTGCCCGTTGGCAAACACCCATTGCAACTTTATTCGCTGGGAACGCCGAACGGTCAGAAAGTAAAGATTATGCTTGAGGAGCTGCTGGCGCTGGGCGTTACTGGTGCAGAGTa                                                                                       <  1:538810/140‑1 (MQ=255)
                                              ttCTGGTCCGACGCATGAAAAAACGCTGCCCGTTGGCAAACACCCATTGCAACTTTATTCGCTGGGAACGCCGAACGGTCAGAAAGTAAAGATTATGCTTGAGGAGCTGCTGGCGCTGGGCGTTACTGGTGCAGAGTacg                                                                                     <  2:1127436/140‑1 (MQ=255)
                                               tCTGGTCCGACGCATGAAAAAACGCTGCCCGTTGGCAAACACCCATTGCAACTTTATTCGCTGGGAACGCCGAACGGTCAGAAAGTAAAGATTATGCTTGAGGAGCTGCTGGCGCTGGGCGTTACTGGTGCAGAGTacga                                                                                    >  1:1285122/1‑140 (MQ=255)
                                                     ccGACGCATGAAAAAACGCTGCCCGTTGGCAAACACCCATTGCAACTTTATTCGCTGGGAACGCCGAACGGTCAGAAAGTAAAGATTATGCTTGAGGAGCTGCTGGCGCTGGGCGTTACTGGTGCAGAGTACGACGCCTg                                                                              <  2:1391653/140‑1 (MQ=255)
                                                      cGACGCATGAAAAAACGCTGCCCGTTGGCAAACACCCATTGCAACTTTATTCGCTGGGAACGCCGAACGGTCAGAAAGTAAAGATTATGCTTGAGGAGCTGCTGGCGCTGGGCGTTACTGGTGCAGAGTACGACGCCTgg                                                                             <  1:380837/140‑1 (MQ=255)
                                                      cGACGCATGAAAAAACGCTGCCCGTTGGCAAACACCCATTGCAACTTTATTCGCTGGGAACGCCGAACGGTCAGAAAGTAAAGATTATGCTTGAGGAGCTGCTGGCGCTGGGCGTTACTGGTGCAGAGTACGACGCCTgg                                                                             <  2:169451/140‑1 (MQ=255)
                                                        aCGCATGAAAAAACGCTGCCCGTTGGCAAACACCCATTGCAACTTTATTCGCTGGGAACGCCGAACGGTCAGAAAGTAAAGATTATGCTTGAGGAGCTGCTGGCGCTGGGCGTTACTGGTGCAGAGTACGACGCCTGGCt                                                                           >  2:336248/1‑140 (MQ=255)
                                                         cGCATGAAAAAACGCTGCCCGTTGGCAAACACCCATTGCAACTTTATTCGCTGGGAACGCCGAACGGTCAGAAAGTAAAGATTATGCTTGAGGAGCTGCTGGCGCTGGGCGTTACTGGTGCAGAGTACGACGCCTGGCTg                                                                          >  2:1395162/1‑140 (MQ=255)
                                                                    aCGCTGCCCGTTGGCAAACACCCATTGCAACTTTATTCGCTGGGAACGCCGAACGGTCAGAAAGTAAAGATTATGCTTGAGGAGCTGCTGGCGCTGGGCGTTACTGGTGCAGAGTACGACGCCTGGCTGATTCGTATtgg                                                               <  1:1247896/140‑1 (MQ=255)
                                                                      gCTGCCCGTTGGCAAACACCCATTGCAACTTTATTCGCTGGGAACGCCGAACGGTCAGAAAGTAAAGATTATGCTTGAGGAGCTGCTGGCGCTGGGCGTTACTGGTGCAGAGTACGACGCCTGGCTGATTCGTATtggcg                                                             >  2:484351/1‑140 (MQ=255)
                                                                         gCCCGTTGGCAAACACCCATTGCAACTTTATTCGCTGGGAACGCCGAACGGTCAGAAAGTAAAGATTATGCTTGAGGAGCTGCTGGCGCTGGGCGTTACTGGTGCAGAGTACGACGCCTGGCTGATTCGTATtggcgatg                                                          <  1:1395162/140‑1 (MQ=255)
                                                                          cccGTTGGCAAACACCCATTGCAACTTTATTCGCTGGGAACGCCGAACGGTCAGAAAGTAAAGATTATGCTTGAGGAGCTGCTGGCGCTGGGCGTTACTGGTGCAGAGTACGACGCCTGGCTGATTCGTATtggcgatgg                                                         <  1:446233/140‑1 (MQ=255)
                                                                          cccGTTGGCAAACACCCATTGCAACTTTATTCGCTGGGAACGCCGAACGGTCAGAAAGTAAAGATTATGCTTGAGGAGCTGCTGGCGCTGGGCGTTACTGGTGCAGAGTACGACGCCTGGCTGATTCGTATtggcgatgg                                                         >  1:1141997/1‑140 (MQ=255)
                                                                            cGTTGGCAAACACCCATTGCAACTTTATTCGCTGGGAACGCCGAACGGTCAGAAAGTAAAGATTATGCTTGAGGAGCTGCTGGCGCTGGGCGTTACTGGTGCAGAGTACGACGCCTGGCTGATTCGTATtggcgatggcg                                                       <  1:1006153/140‑1 (MQ=255)
                                                                            cGTTGGCAAACACCCATTGCAACTTTATTCGCTGGGAACGCCGAACGGTCAGAAAGTAAAGATTATGCTTGAGGAGCTGCTGGCGCTGGGCGTTACTGGTGCAGAGTACGACGCCTGGCTGATTCGTATtggcgatggcg                                                       <  1:1277311/140‑1 (MQ=255)
                                                                            cGTTGGCAAACACCCATTGCAACTTTATTCGCTGGGAACGCCGAACGGTCAGAAAGTAAAGATTATGCTTGAGGAGCTGCTGGCGCTGGGCGTTACTGGTGCAGAGTACGACGCCTGGCTGATTCGTATtggcgatggcg                                                       <  1:151125/140‑1 (MQ=255)
                                                                            cGTTGGCAAACACCCATTGCAACTTTATTCGCTGGGAACGCCGAACGGTCAGAAAGTAAAGATTATGCTTGAGGAGCTGCTGGCGCTGGGCGTTACTGGTGCAGAGTACGACGCCTGGCTGATTCGTATtggcgatggcg                                                       <  1:846350/140‑1 (MQ=255)
                                                                            cGTTGGCAAACACCCATTGCAACTTTATTCGCTGGGAACGCCGAACGGTCAGAAAGTAAAGATTATGCTTGAGGAGCTGCTGGCGCTGGGCGTTACTGGTGCAGAGTACGACGCCTGGCTGATTCGTATtggcgatggcg                                                       <  1:198218/140‑1 (MQ=255)
                                                                                   aaaCACCCATTGCAACTTTATTCGCTGGGAACGCCGAACGGTCAGAAAGTAAAGATTATGCTTGAGGAGCTGCTGGCGCTGGGCGTTACTGGTGCAGAGTACGACGCCTGGCTGATTCGTATTGGCGATGGCGATCAAtt                                                <  1:1628791/140‑1 (MQ=255)
                                                                                    aaCACCCATTGCAACTTTATTCGCTGGGAACGCCGAACGGTCAGAAAGTAAAGATTATGCTTGAGGAGCTGCTGGCGCTGGGCGTTACTGGTGCAGAGTACGACGCCTGGCTGATTCGTATTggcgatggcgat                                                     <  1:596438/134‑1 (MQ=255)
                                                                                    aaCACCCATTGCAACTTTATTCGCTGGGAACGCCGAACGGTCAGAAAGTAAAGATTATGCTTGAGGAGCTGCTGGCGCTGGGCGTTACTGGTGCAGAGTACGACGCCTGGCTGATTCGTATTggcgatggcgat                                                     >  2:596438/1‑134 (MQ=255)
                                                                                         ccATTGCAACTTTATTCGCTGGGAACGCCGAACGGTCAGAAAGTAAAGATTATGCTTGAGGAGCTGCTGGCGCTGGGCGTTACTGGTGCAGAGTACGACGCCTGGCTGATTCGTATTGGCGATGGCGATCAATTCTCCAg                                          >  2:735850/1‑140 (MQ=255)
                                                                                           aTTGCAACTTTATTCGCTGGGAACGCCGAACGGTCAGAAAGTAAAGATTATGCTTGAGGAGCTGCTGGCGCTGGGCGTTACTGGTGCAGAGTACGACGCCTGGCTGATTCGTATTGGCGATGGCGATCAATTCTCCAGCg                                        >  2:553489/1‑140 (MQ=255)
                                                                                            ttGCAACTTTATTCGCTGGGAACGCCGAACGGTCAGAAAGTAAAGATTATGCTTGAGGAGCTGCTGGCGCTGGGCGTTACTGGTGCAGAGTACGACGCCTGGCTGATTCGTATTGGCGATGGCGATCAATTCTCCAGTgg                                       <  2:74702/140‑1 (MQ=255)
                                                                                            ttGCAACTTTATTCGCTGGGAACGCCGAACGGTCAGAAAGTAAAGATTATGCTTGAGGAGCTGCTGGCGCTGGGCGTTACTGGTGCAGAGTACGACGCCTGGCTGATTCGTATTGGCGATGGCGATCAATTCTCCAGCgg                                       <  2:177899/140‑1 (MQ=255)
                                                                                            ttGCAACTTTATTCGCTGGGAACGCCGAACGGTCAGAAAGTAAAGATTATGCTTGAGGAGCTGCTGGCGCTGGGCGTTACTGGTGCAGAGTACGACGCCTGGCTGATTCGTATTGGCGATGGCGATCAATTCTCCAGCgg                                       <  2:693482/140‑1 (MQ=255)
                                                                                              gCAACTTTATTCGCTGGGAACGCCGAACGGTCAGAAAGTAAAGATTATGCTTGAGGAGCTGCTGGCGCTGGGCGTTACTGGTGCAGAGTACGACGCCTGGCTGATTCGTATTGGCGATGGCGATCAATTCTCCAGCGGCt                                     <  1:1706748/140‑1 (MQ=255)
                                                                                                 aCTTTATTCGCTGGGAACGCCGAACGGTCAGAAAGTAAAGATTATGCTTGAGGAGCTGCTGGCGCTGGGCGTTACTGGTGCAGAGTACGACGCCTGGCTGATTCGTATTGGCGATGGCGATCAATTCTCCAGCGGCTTTg                                  >  1:568905/1‑140 (MQ=255)
                                                                                                   tttATTCGCTGGGAACGCCGAACGGTCAGAAAGTAAAGATTATGCTTGAGGAGCTGCTGGCGCTGGGCGTTACTGGTGCAGAGTACGACGCCTGGCTGATTCGTATTGGCGATGGCGATCAATTCTCCAGCGGCTTTGTc                                <  1:549227/140‑1 (MQ=255)
                                                                                                   tttATTCGCTGGGAACGCCGAACGGTCAGAAAGTAAAGATTATGCTTGAGGAGCTGCTGGCGCTGGGCGTTACTGGTGCAGAGTACGACGCCTGGCTGATTCGTATTGGCGATGGCGATCAATTCTCCAGCGGCTTTGTc                                >  1:1617414/1‑140 (MQ=255)
                                                                                                     tATTCGCTGGGAACGCCGAACGGTCAGAAAGTAAAGATTATGCTTGAGGAGCTGCTGGCGCTGGGCGTTACTGGTGCAGAGTACGACGCCTGGCTGATTCGTATTGGCGATGGCGATCAATTCTCCAGCGGCTTTGTCGa                              >  1:1504271/1‑140 (MQ=255)
                                                                                                        tCGCTGGGAACGCCGAACGGTCAGAAAGTAAAGATTATGCTTGAGGAGCTGCTGGCGCTGGGCGTTACTGGTGCAGAGTACGACGCCTGGCTGATTCGTATTGGCGATGGCGATCAATTCTCCAGCGGCTTTGTCGAAGt                           <  2:1427889/140‑1 (MQ=255)
                                                                                                            tGGGAACGCCGAACGGTCAGAAAGTAAAGATTATGCTTGAGGAGCTGCTGGCGCTGGGCGTTACTGGTGCAGAGTACGACGCCTGGCTGATTCGTATTGGCGATGGCGATCAATTCTCCAGCGGCTTTGTCGAAGTGAAc                       >  2:534111/1‑140 (MQ=255)
                                                                                                            tGGGAACGCCGAACGGTCAGAAAGTAAAGATTATGCTTGAGGAGCTGCTGGCGCTGGGCGTTACTGGTGCAGAGTACGACGCCTGGCTGATTCGTATTGGCGATGGCGATCAATTCTCCAGCGGCTTTGTCGAAGTGAAc                       >  1:1218095/1‑140 (MQ=255)
                                                                                                            tGGGAACGCCGAACGGTCAGAAAGTAAAGATTATGCTTGAGGAGCTGCTGGCGCTGGGCGTTACTGGTGCAGAGTACGACGCCTGGCTGATTCGTATTGGCGATGGCGATCAATTCTCCAGCGGCTTTGTCGAAGTGAAc                       >  2:1620541/1‑140 (MQ=255)
                                                                                                            tGGGAACGCCGAACGGTCAGAAAGTAAAGATTATGCTTGAGGAGCTGCTGGCGCTGGGCGTTACTGGTGCAGAGTACGACGCCTGGCTGATTCGTATTGGCGATGGCGATCAATTCTACAGCGGCTTTGTCGAAGTGAAc                       >  2:658414/1‑140 (MQ=255)
                                                                                                                 aCGCCGAACGGTCAGAAAGTAAAGATTATGCTTGAGGAGCTGCTGGCGCTGGGCGTTACTGGTGCAGAGTACAACGCCTGGCTGATTCGTATTGGCGATGGCGATCAATTCTCCAGCGGCTTTGTCGAAGTGAACCCaaa                  >  2:86416/1‑140 (MQ=255)
                                                                                                                      gAACGGTCAGAAAGTAAAGATTATGCTTGAGGAGCTGCTGGCGCTGGGCGTTACTGGTGCAGAGTACGACGCCTGGCTGATTCGTATTGGCGATGGCGATCAATTCTCCAGCGGCTTTGTCGAAGTGAACCCa                    <  2:808076/133‑1 (MQ=255)
                                                                                                                      gAACGGTCAGAAAGTAAAGATTATGCTTGAGGAGCTGCTGGCGCTGGGCGTTACTGGTGCAGAGTACGACGCCTGGCTGATTCGTATTGGCGATGGCGATCAATTCTCCAGCGGCTTTGTCGAAGTGAACCCa                    >  1:808076/1‑133 (MQ=255)
                                                                                                                           gTCAGAAAGTAAAGATTATGCTTGAGGAGCTGCTGGCGCTGGGCGTTACTGGTGCAGAGTACGACGCCTGGCTGATTCGTATTGGCGATGGCGATCAATTCTCCAGCGGCTTTGTCGAAGTGAACCCAAACTCGAAGATc        <  1:180600/140‑1 (MQ=255)
                                                                                                                                 aaGTAAAGATTATGCTTGAGGAGCTGCTGGCGCTGGGCGTTACTGGTGCAGAGTACGACGCCTGGCTGATTCGTATTGGCGATGGCGATCAATTCTCCAGCGGCTTTGTCGAAGTGAACCCAAACTCGAAGATCCCGgcg  <  2:1170119/140‑1 (MQ=255)
                                                                                                                                 aaGTAAAGATTATGCTTGAGGAGCTGCTGGCGCTGGGCGTTACTGGTGCAGAGTACGACGCCTGGCTGATTCGTATTGGCGATGGCGATCAATTCTCCAGCGGCTTTGTCGAAGTGAACCCAAACTCGAAGATCCCGgcg  <  1:484351/140‑1 (MQ=255)
                                                                                                                                       |                                                                                                                                     
GGGATAAATCCGCTGGCGGCGCGTTCGCCAATATCAATCGCCCGGTTTCTGGTCCGACGCATGAAAAAACGCTGCCCGTTGGCAAACACCCATTGCAACTTTATTCGCTGGGAACGCCGAACGGTCAGAAAGTAACGATTATGCTTGAGGAGCTGCTGGCGCTGGGCGTTACTGGTGCAGAGTACGACGCCTGGCTGATTCGTATTGGCGATGGCGATCAATTCTCCAGCGGCTTTGTCGAAGTGAACCCAAACTCGAAGATCCCGGCG  >  NC_000913/3138767‑3139035

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