Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,709,415 (C)7→9 coding (274/1059 nt) rsxD → SoxR [2Fe‑2S] reducing system protein RsxD

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,709,4081.C88.3% 452.0 / 46.2 145coding (267/1059 nt)rsxDSoxR [2Fe‑2S] reducing system protein RsxD
Reads supporting (aligned to +/- strand):  ref base . (5/12);  new base C (57/71);  total (62/83)
Fisher's exact test for biased strand distribution p-value = 3.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
*NC_0009131,709,4082.C84.8% 439.9 / 64.3 145coding (267/1059 nt)rsxDSoxR [2Fe‑2S] reducing system protein RsxD
Reads supporting (aligned to +/- strand):  ref base . (8/14);  new base C (54/69);  total (62/83)
Fisher's exact test for biased strand distribution p-value = 6.41e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 168 total aligned reads displayed.

CTCGTTCAGATCCTGTTGGCATCGGTTAGTGCTCTGTTAGCCGAAGCTCTCGTACTCAAACTACGCAAGCAGTCGGTAGCCGCAACGTTGAAAGATAACTCAGCATTGCTGACAGGCTTATTGCTGGCGGTAAGTATT‑‑‑CCCC‑‑C‑‑CC‑‑TCGCGCCATGGTGGATGGTCGTGCTGGGTACGGTGTTTGCGGTGATCATCGCTAAACAGTTGTATGGCGGTCTGGGACAAAACCCGTTTAATCCGGCAATGATTGGTTATGTGGTCTTACTGA  >  NC_000913/1709271‑1709538
                                                                                                                                          ||                                                                                                                                          
ctcGTTCAGATCCTGTTGGCATCGGTTAGTGCTCTGTTAGCCGAAGCTCTCGTACTCAAACTACGCAAGCAGTCGGTAGCCGCAACGTTGAAAGATAACTCATCATTGCTAACAAGCTTATTGCAGGCGGTAAGTATT‑‑‑cc                                                                                                                                         >  1:1182851/1‑140 (MQ=255)
   gTTCAGATCCTGTTGGCATCGGTTAGTGCTCTGTTAGCCGAAGCTCTCGTACTCAAACTACGCAAGCAGTCGGTAGCCGCAACGTTGAAAGATAACTCAGCATTGCTGACAGGCTTATTGCTGGCGGTAAGTATT‑‑‑cccc‑‑c                                                                                                                                   <  2:57708/140‑1 (MQ=255)
      cAGATCCTGTTGGCATCGGTTAGTGCTCTGTTAGCCGAAGCTCTCGTACTCAAACTACGCAAGCAGTCGGTAGCCGCAACGTTGAAAGATAACTCAGCATTGCTGACAGGCTTATTGCTGGCGGTAAGTATT‑‑‑cccc‑‑c‑‑ccc                                                                                                                              >  2:280575/1‑139 (MQ=255)
       aGATCCTGTTGGCATCGGTTAGTGCTCTGTTAGCCGAAGCTCTCGTACTCAAACTACGCAAGCAGTCGGTAGCCGCAACGTTGAAAGATAACTCAGCATTGCTGACAGGCTTATTGCTGGCGGTAAGTATT‑‑‑cccc‑‑c‑‑cccc                                                                                                                             <  1:2031488/140‑3 (MQ=255)
        gATCCTGTTGGCATCGGTTAGTGCTCTGTTAGCCGAAGCTCTCGTACTCAAACTACGCAAGCAGTCGGTAGCCGCAACGTTGAAAGATAACTCAGCATTGCTGACAGGCTTATTGCTGGCGGTAAGTATT‑‑‑cccc‑‑c‑‑cccct                                                                                                                            <  2:41191/140‑4 (MQ=255)
         aTCCTGTTGGCATCGGTTAGTGCTCTGTTAGCCGAAGCTCTCGTACTCAAACTACGCAAGCAGTCGGTAGCCGCAACGTTGAAAGATAACTAAGCATTGCTGACAGGCCTATTGCTGGCGGTCAGTATT‑‑‑cccc‑‑c‑‑cccctc                                                                                                                           >  2:990703/1‑136 (MQ=255)
           ccTGTTGGCATCGGTTAGTGCTCTGTTAGCCGAAGCTCTCGTACTCAAACTACGCAAGCAGTCGGTAGCCGCAACGTTGAAAGATAACTCAGCATTGCTGACAGGCTTATTGCTGGCGGTAAGTATT‑‑‑CCCC‑‑C‑‑CC‑‑CCtcgc                                                                                                                       >  1:210642/1‑140 (MQ=255)
            cTGTTGGCATCGGTTAGTGCTCTGTTAGCGGAAGCTCTGTTTCTCAAATTGCGCAAGCGGTCGGTAGTCGCAACGTTGAAAGATAACTCAGCATTGCTGACAGGCTTATTGCTGGCGGTAAGTATT‑‑‑CCCC‑‑C‑‑CC‑‑CCtcgcg                                                                                                                      <  2:901169/140‑2 (MQ=255)
            cTGTTGGCATCGGTTAGTGCTCTGTTAGCCGAAGCTCTCGTACTCAAACTACGCAAGCAGTCGGTAGCCGCAACGTTGAAAGATAACTCAGCATTGCTGACAGGCTTATTGCTGGCGGTAAGTATT‑‑‑CCCC‑‑C‑‑CC‑‑CCtcgcg                                                                                                                      <  2:1954008/140‑2 (MQ=255)
             tGTTGGCATCGGTTAGTGCTCTGTTAGCCGAAGCTCTCGTACTCAAACTACGCAAGCAGTCGGTAGCCGCAACGTTGAAAGATAACTCAGCATTGCTGACAGGCTTATTGCTGGCGGTAAGTATT‑‑‑CCCC‑‑C‑‑CC‑‑CCtcgcgc                                                                                                                     >  2:1602877/1‑138 (MQ=255)
              gTTGGCATCGGTTAGTGCTCTGTTAGCCGAAGCTCTCGTACTCAAACTACGCAAGCAGTCGGTAGCCGCAACGTTGAAAGATAACTCAGCATTGCTGACAGGCTTATTGCTGGCGGTAAGTATT‑‑‑CCCC‑‑C‑‑CC‑‑CCtcgcgcc                                                                                                                    <  1:2230774/140‑4 (MQ=255)
               ttGGCATCGGTTAGTGCTCTGTTAGGCGATGCTCTCGTACTCAAACTACGCAAGCAGTCGGTAGCCGCAACGTTGAAAGATAACTCAGCATTGCTGACAGGCTTATTGCTGGCGGTAAGTATTCCCCCC‑‑C‑‑CC‑‑TCGCGCCa                                                                                                                     <  1:996711/140‑1 (MQ=255)
               ttGGCATCGGTTAGTGCTCTGTTAGCCGAAGCTCTCGTACTCAAACTACGCAAGCAGTCGGTAGCCGCAACGTTGAAAGATAACTCAGCATTGCTGACAGGCTTATTGCTGGCGGTAAGTATTCCCCCC‑‑C‑‑CC‑‑TCGCGCCa                                                                                                                     >  2:2151425/1‑140 (MQ=255)
               ttGGCATCGGTTAGTCCTCTGTTAGCCGAGGCTCTCGTACTCAAACTACGCAAGCAGTCGGTAGCCGCAACGTTGAAAGATAACTCAGCATTGCTGACAGGCTTATTGCTGGCGGTAAGTATTCCCCCC‑‑C‑‑CC‑‑TCGCGCCa                                                                                                                     <  1:728774/140‑1 (MQ=255)
                   cATCGGTTAGTGCTCTGTTAGCCGAAGCTCTCGTACTCAAACTACGCAAGCAGTCGGTAGCCGCAACGTTGAAAGATAACTCAGCATTGCTGACAGGCTTATTGCTGGCGGTAAGTATT‑‑‑cccc‑‑c‑‑ccccc                                                                                                                            <  1:2320083/129‑4 (MQ=255)
                   cATCGGTTAGTGCTCTGTTAGCCGAAGCTCTCGTACTCAAACTACGCAAGCAGTCGGTAGCCGCAACGTTGAAAGATAACTCAGCATTGCTGACAGGCTTATTGCTGGCGGTAAGTATT‑‑‑cccc‑‑c‑‑ccccc                                                                                                                            >  2:2320083/1‑126 (MQ=255)
                   cATCGGTTAGTGCTCTGTTAGCCGAAGCTCTCGTACTCAAACTACGCAAGCAGTCGGTAGCCGCAACGTTGAAAGATAACTCAGCATTGCTGACAGGCTTATTGCTGGCGCTAACTAGT‑‑‑CCCCCTC‑‑CC‑‑TCGCGCCAtggg                                                                                                                 >  1:1145736/1‑139 (MQ=255)
                    aTCGGTTAGTGCTCTGTTAGCCGAAGCTCTCGTACTCAAACTACGCAAGCAGTCGGTAGCCGCAACGTTGAAAGATAACTCAGCATTGCTGACAGGCTTATTGCTGGCGGTAAGTATTCCCCCC‑‑C‑‑CC‑‑TCGCGCCa                                                                                                                     <  1:421311/135‑1 (MQ=255)
                         ttAGTGCTGTGTTAGCATAAGCTCTCGTACCAAAACTACGCAAGCAGTTGGTAGCCGCAACGTTGAAAGATAACTCAGCATTGCTGACAGGCTTATTGCTGGCGGTAAGTATTCCCCCC‑‑C‑‑CC‑‑TCGCGCCATGGTGGATgg                                                                                                           <  2:1182851/140‑1 (MQ=255)
                         ttAGTGCTCTGTTAGCCGAAGCTCTCGTACTCAAACTACGCAAGCAGTCGGTAGCCGCAACGTTGAAAGATAACTCAGCATTGCTGACAGGCTTATTGCTGGCGGTAAGTATTCCCCCC‑‑C‑‑CC‑‑TCGCGCCATGGTGGATgg                                                                                                           >  1:1864042/1‑140 (MQ=255)
                           aGTGCTCTGTTAGCCGAAGCTCTCGTACTCAAACTACGCAAGCAGTCGGTAGCCGCAACGTTGAAAGATAACTCAGCATTGCTGACAGGCTTATTGCTGGCGGTAAGTATTCCCCCC‑‑C‑‑CC‑‑TCGCGCCATGGTGGATGGTc                                                                                                         <  1:1697218/140‑1 (MQ=255)
                             tGCTCTGTTAGCCGAAGCTCTCGTACTCAAACTACGCAAGCAGTCGGTAGCCGCAACGTTGAAAGATAACTCAGCATTGCTGACAGGCTTATTGCTGGCGGTAAGTATTCCCCCC‑‑C‑‑CC‑‑TCGCGCCATGGTGGATGGTCGt                                                                                                       <  1:1816640/140‑1 (MQ=255)
                               ctctGTTAGCCGAAGCTCTCGTACTCAAACTACGCAAGCAGTCGGTAGCCGCAACGTTGAAAGATAACTCAGCATTGCTGACAGGCTTATTGCTGGCGGTAAGTATTCCCCCC‑‑C‑‑CC‑‑TCGCGCCATGGTGGATGGTCGTGc                                                                                                     <  1:1624008/140‑1 (MQ=255)
                               ctctGTTAGCCGAAGCTCTCGTACTCAAACTACGCAAGCAGTCGGTAGCCGCAACGTTGAAAGATAACTCAGCATTGCTGACAGGCTTATTGCTGGCGGTAAGTATTCCCCCC‑‑C‑‑CC‑‑TCGCGCCATGGTGGATGGTCGTGc                                                                                                     <  1:1700207/140‑1 (MQ=255)
                                  tGTTAGCCGAAGCTCTCGTACTCAAACTACGCAAGCAGTCGGTAGCCGCAACGTTGAAAGATAACTCAGCATTGCTGACAGGCTTATTGCTTGCGGTAAGTATTGCCCCC‑‑C‑‑CC‑‑CCGCGCCATAAGGGATGGTCGTGCTgg                                                                                                  >  2:1060889/1‑140 (MQ=255)
                                  tGTTAGCCGAAGCTCTCGTACTCAAACTACGCAAGCAGTCGGTAGCCGCAACGTTGAAAGATAACTCAGCATTGCTGACAGGCTTATTGCTGGCGGTAAGTATTCCCCCC‑‑C‑‑CC‑‑TGGCACCATGGGGGATGGTCGTGCTgg                                                                                                  >  1:919302/1‑140 (MQ=255)
                                  tGTTAGCCGAAGCTCTCGTACTCAAACTACGCAAGCAGTCGGTAGCCGCAACGTTGAAAGATAACTCAGCATTGCTGACAGGCTTATTGCTGGCGGTAAGTATTCCCCCC‑‑C‑‑CC‑‑TCGCGCCATGGTGGATGGTCGTGCTgg                                                                                                  >  2:1610579/1‑140 (MQ=255)
                                       gCCGAAGCTCTCGTACTCAAACTACGCAAGCAGTCGGTAGCCGCAACGTTGAAAGATAACTCAGCATTGCTGACAGGCTTATTGCTGGCGGTAAGTATTCCCCCC‑‑C‑‑CC‑‑TCGCGCCATGTTGGATGGTCGTGCTGGGTACg                                                                                             >  2:1665375/1‑140 (MQ=255)
                                       gCCGAAGCTCTCGTACTCAAACTACGCAAGCAGTCGGTAGCCGCAACGTTGAAAGATAACTCAGCATTGCTGACAGGCTTATTGCTGGCGGTAAGTATTCCCCCC‑‑C‑‑CC‑‑TCGCGCCATGGTGGATGGTCGTGCTGGGTACg                                                                                             >  1:1398403/1‑140 (MQ=255)
                                            aGCTCTCGTACTCAAACTACGCAAGCAGTCGGTAGCCGCAACGTTGAAAGATAACTCAGCATTGCTGACAGGCTTATTGCTGGCGGTAAGTATTCCCCCC‑‑C‑‑CC‑‑TCGCGCCATGGTGGATGGTCGTGCTGGGTACGGTGtt                                                                                        <  1:151067/140‑1 (MQ=255)
                                            aGCTCTCGTACTCAAACTACGCAAGCAGTCGGTAGCCGCAACGTTGAAAGATAACTCAGCATTGCTGACAGGCTTATTGCTGGCGGTAAGTATT‑‑‑CCCC‑‑CCTCC‑‑TTGCGCCATGGTGGGTGGTAGTGCTGGGTACGGTGtt                                                                                        >  2:854649/1‑140 (MQ=255)
                                              ctctCGTACTCAAACTACGCAAGCAGTCGGTAGCCGCAACGTTGAAAGATAACTCAGCATTGCTGACAGGCTTATTGCTGGCGGTAAGTATTCCCCCC‑‑C‑‑CC‑‑TCGCGCCATGGTGGATGGTCGTGCTGGGTACGGTGTTTg                                                                                      >  2:117488/1‑140 (MQ=255)
                                                 tcGTACTCAAACTACGCAAGCAGTCGGTAGCCGCAACGTTGAAAGATAACTCAGCATTGCTGACAGGCTTATTGCTGGCGGTAAGTATTCCCCCC‑‑C‑‑CC‑‑TCGCGCCATGGTTGATAGTAGTGCTGCGTACGGTGTTTGCgg                                                                                   >  1:690998/1‑140 (MQ=255)
                                                 tcGTACTCAAACTACGCAAGCAGTCGGTAGCCGCAACGTTGAAAGATAACTCAGCATTGCTGACAGGCTTATTGCTGGCGGTAAGTATTCCCCCC‑‑C‑‑CC‑‑TCGCGCCATGGTGGATGGTCGTGCTGGGTACGGTGTTTGCgg                                                                                   >  2:2381512/1‑140 (MQ=255)
                                                     aCTCAAACTACGCAATCAGTCGCTATCCGCAACGTTGAAAAATAACTCAGCATTGCTGACAGGCTTATTGCTGGCGGTAAGTATTCCCCCC‑‑C‑‑CC‑‑TCGCGCCATGGTGGATGGTCGTGCTGGGTACGGTGTTTGCGGTGat                                                                               <  1:1060889/140‑1 (MQ=255)
                                                       tCAAACTACGCAAGCAGTCGGTTGCCGCAACTTTGCAAGATAACTCAGCATTGCTGACAGGCTTATTGCTGGCGGTAAGTATTCCCCCC‑‑C‑‑CC‑‑TCGCGCCATGGTGGATGGTCGTGCTGGGTACGGTGTTTGCGGTGatca                                                                             <  2:780512/140‑1 (MQ=255)
                                                        cAAACTACGCAAGCAGTCGGTATCCGCAACGCTGATAGATAACTCAGCCTTGCTGACAGGCTTATTGCTGGCGGTAAGTATTCCCCCC‑‑C‑‑CC‑‑TCACGCCATGGTGGATGGTCGTGCTGGGTACGGTGTTTGCGGTGatcat                                                                            <  2:937693/140‑1 (MQ=255)
                                                          aaCTACGCAAGCAGTCGGTAGCCGCAACGTTGAAAGATAACTCAGCATTGCTGACAGGCTTATTGCTGGCGGTAAGTATTCCCCCC‑‑C‑‑CC‑‑TCGCGCCATGGTGGATGGTCGTGCTGGGTACGGTGTTTGCGGTGATCATCg                                                                          <  1:2433829/140‑1 (MQ=255)
                                                          aaCTACGCAAGCAGTCGGTAGCCGCAACGTTGAAAGATAACTCAGCATTGCTGACAGGCTTATTGCTGGCGGTAAGTATTCCCCCC‑‑C‑‑CC‑‑TCGCGCCATGGTGGATGGTCGTGCTGGGTACGGTGTTTGCGGTGATCATCg                                                                          <  1:2311342/140‑1 (MQ=255)
                                                             tACGCAAGCAGTCGGTAGCCGCAACGTTGAAAGATAACTCAGCATTGCTGACAGGCTTATTGCTGGCGGTAAGTATTCCCCCC‑‑C‑‑CC‑‑TCGCGCCATGGTGGATGGTCGTGCTGGGTACGGTGttt                                                                                       <  2:2412012/124‑1 (MQ=255)
                                                             tACGCAAGCAGTCGGTAGCCGCAACGTTGAAAGATAACTCAGCATTGCTGACAGGCTTATTGCTGGCGGTAAGTATTCCCCCC‑‑C‑‑CC‑‑TCGCGCCATGGTGGATGGTCGTGCTGGGTACGGTGTTTGCGGTGATCATCGCTa                                                                       <  1:1602877/140‑1 (MQ=255)
                                                             tACGCAAGCAGTCGGTAGCCGCAACGTTGAAAGATAACTCAGCATTGCTGACAGGCTTATTGCTGGCGGTAAGTATTCCCCCC‑‑C‑‑CC‑‑TCGCGCCATGGGGGATGGTCGTGCTGGGTACGGTGttt                                                                                       >  1:2412012/1‑124 (MQ=255)
                                                               cGCAAGCAGTCGGTAGCCGCAACGTTGAAAGATAACTCAGCATTGCTGACAGGCTTATTGCTGGCGGTAAGTATTCCCCCC‑‑C‑‑CC‑‑TCGCGCCATGGTGGATGGTCGTGCTGGGTACGGTGTTTGCGGTGATCATCGCTaaa                                                                     >  2:406894/1‑140 (MQ=255)
                                                               cGCAAGCAGTCGGTAGCCGCAACGTTGAAAGATAACTCAGCATTGCTGACAGGCTTATTGCTGGCGGTAAGTATTCCCCCC‑‑C‑‑CC‑‑TCGCGCCATGGGGGATGGTCGTGCTGGGTACGGTGTTTGCGGTGATCATCGCTaaa                                                                     >  2:281619/1‑140 (MQ=255)
                                                               cGCAAGCAGTCGGTAGCCGCAACGTTGAAAGATAACTCAGCATTGCGGACAGGCTTATTGCTGGCGGTAAGTATTCCCCCC‑‑C‑‑CC‑‑TCGCGCCATGGTGGATGGTCGTGCTGGGTACGGTGTTTGCGGTGATCATCGCTaaa                                                                     >  2:1485248/1‑140 (MQ=255)
                                                                gCAAGCAGTCGGTAGCCGCAACGTTGAAAGATAACTCAGCATTGCTGACAGGCTTATTGCTGGCGGTAAGTATTCCCCCC‑‑C‑‑CC‑‑TCGCGCCATGGTGGATGGTCGTGCTGGGTACGGTGTTTGTGGTGATCATCGCTAAAc                                                                    <  2:1320219/140‑1 (MQ=255)
                                                                 cAAGCAGTCGGTAGCCGCAACGTTGAAAGATAACTCAGCATTGTTGACAGGCTTATTGCTGGCGGTAAGTATTCCCCCC‑‑C‑‑CC‑‑TCGCGCCATGGTGGATGGTCGTGATGGGTACGGTGTTTGCGGTGATCATCGCTaaa                                                                     <  2:1705735/138‑1 (MQ=255)
                                                                   aGCAGTCGGTAGCCGCAACGTTGAAAGATAACTCAGCATTGCTGACAGGCTTATTGCTGGCGGTAAGTATTCCCCCC‑‑C‑‑CC‑‑TCGCGCCATGGTGGATGGTCGTGCTGGGTACGGTGTTTGCGGTGATCATCg                                                                          >  2:1713054/1‑131 (MQ=255)
                                                                   aGCAGTCGGTAGCCGCAACGTTGAAAGATAACTCAGCATTGCTGACAGGCTTATTGCTGGCGGTAAGTATTCCCCCC‑‑C‑‑CC‑‑TCGCGCCATGGTGGATGGTCGTGCTGGGTACGGTGTTTGCGGTGATCATCg                                                                          <  1:1713054/131‑1 (MQ=255)
                                                                     cAGTCGGTAGCCGCAACGTTGAAAGATAACTCAGCATTGCTGACAGGCTTATTGCTGGCGGTAAGTATTCCCCCC‑‑C‑‑CC‑‑TCGCGCCATGGTGGATGGTCGTGCTGGGTACGGTGTTTGCGGTGATCATCGCTAAACAGTTg                                                               >  2:2329930/1‑140 (MQ=255)
                                                                          ggTAGCCGCAACGTTGAAAGATAACTCAGCATTGCTGACAGGCTTATTGCTGGCGGTAAGTATTCCCCCC‑‑C‑‑CC‑‑TCGCGCCATGGTGGATGGTCGTGCTGGGTACGGTGTTTGCGGTGATCATCGCTAAACAGTTGTATgg                                                          <  2:1883147/140‑1 (MQ=255)
                                                                          ggTAGCCGCAACGTTGAAAGATAACTCAGCATTGCTGACAGGCTTATTGCTGGCGGTAAGTATTCCCCCC‑‑C‑‑CC‑‑TCGCGCCATGGTGGATGGTCGTGCTGGGTACGGTGTTTGCGGTGATCATCGCTAAACAGTTGTATgg                                                          <  2:210642/140‑1 (MQ=255)
                                                                          ggTAGCCGCAACGTTGAAAGATAACTCAGCATTGCTGACAGGCTTATTGCTGGCGGTAAGTATTCCCCCC‑‑C‑‑CC‑‑TCGCGCCATGGTGGATGGTCGTGCTGGGTACGGTGTTTGCGGTGATCATCGCTAAACAGTTGTATgg                                                          <  1:1656285/140‑1 (MQ=255)
                                                                          ggTAGCCGCAACGTTGAAAGATAACTCAGCATTGCTGACAGGCTTATTGCTGGCGGTAAGTATTCCCCCC‑‑C‑‑CC‑‑TCGCGCCATGGTGGATGGTCGTGCTGGGTACGGTGTTTGCGGTGATCATCGCTAAACAGTTGTATgg                                                          <  1:1500374/140‑1 (MQ=255)
                                                                          ggTAGCCGCAACGTTGAAAGATAACTCAGCATTGCTGACAGGCTTATTGCTGGCGGTAAGTATTCCCCCC‑‑C‑‑CC‑‑TCGCGCCATGGTGGATGGTCGTGCTGGGTACGGTGTTTGCGGTGATCATCGCTAAACAGTTGTATgg                                                          <  2:1294549/140‑1 (MQ=255)
                                                                           gTAGCCGCAACGTTGAAAGATAACTCAGCATTGCTGACAGGCTTATTGCTGGCGGTAAGTATTCCCCCC‑‑C‑‑CC‑‑TCGCGCCATGGTGGATGGTCGTGCTGGGTACGGTGTTTGCGGTGATCATCGCTAAACAGTTGTATGGc                                                         <  2:1517661/140‑1 (MQ=255)
                                                                             aGCCGCAACGTTGAAAGATAACTCAGCATTGCTGACAGGCTTATTGCTGGCGGTAAGTATTCCCCCC‑‑C‑‑CC‑‑CCGCGCCCTGGGGGGTGGGCGGGGTTGGGACGGGGTTTTGGGGGGTCCTCCCTTAACAGGTGTATGgggg                                                       >  2:424136/1‑140 (MQ=255)
                                                                                cGCAACGTTGAAAGATAACTCAGCATTGCTGACAGGCTTATTGCTGGCGGTAAGTATTCCCCCC‑‑C‑‑CC‑‑TCGCGCCATGGTGGATGGTCGTGCTGGGTACGGTGTTTGCGGTGATCATCGCTAAACAGTTGTATGGCGGTCt                                                    <  1:117488/140‑1 (MQ=255)
                                                                                  cAACGTTGAAAGATAACTCAGCATTGCTGACAGGCTTATTGCTGGCGGTAAGTATTCCCCCC‑‑C‑‑CC‑‑TCGCGCCATGGTGGATGGTCGTGCTGGGTACGGTGTTTGCGGTGATCATCGCTAAACAGTTGTATGGCGGTCTgg                                                  >  2:39556/1‑140 (MQ=255)
                                                                                    aCGTTGAAAGATAACTCAGCATTGCTGACAGGCTTATTGCTGGCGGTAAGTATTCCCCCC‑‑C‑‑CC‑‑TCGCGCCATGGTGGATGGTCGTGCTGGGTACGGTGTTTGCGGTGATCATCGCTAAACAGTTGTATGGCGGTCTGGGa                                                <  1:1775977/140‑1 (MQ=255)
                                                                                    aCGTTGAAAGATAACTCAGCATTGCTGACAGGCTTATTGCTGGCGGTAAGTATTCCCCCC‑‑C‑‑CC‑‑TCGCGCCATGGGGGATGGTCGTGCTGGGTACGGTGTTTGCGGTGATCATCGCTAAACAGTTGTATGGCGGTCTGGGa                                                >  1:1598540/1‑140 (MQ=255)
                                                                                         gAAAGATAACTCAGCATTGCTGACAGGCTTATTGCTGGCGGTAAGTATTC‑‑CCCC‑‑C‑‑CC‑‑TCGCGCCATGGGGGATGGTCGTGCTGGGTACGGTGTTTGCGGTGATCATCGCTAAACAGTTGTATGGCGGTCTGGGACTAAAc                                          >  2:2450285/1‑140 (MQ=255)
                                                                                             gATAACTCAGCATCGCTGACAGGCTTATTGCTGGCGGTAAGTATTCCCCCC‑‑C‑‑CC‑‑TCGCGCCATGGTGGATGGTCGTGCTGGGTACGGTGTTTGCGGTGATCATCGCTAAACAGTTGTATGGCGGTCTGGGACAAAACCCg                                       <  2:2285148/140‑1 (MQ=255)
                                                                                               tAACTCAGCATTGCTGTCAGGCTTATTGCTGGCGGTAAGTATTCCCCCC‑‑C‑‑CC‑‑TCGCGCCATGGTGGATGGTCGTGCTGGGTACGGTGTTTGCGGTGATCATCGCTAAACAGTTGTATGGCGGTCTGGGACAAAACCCGtt                                     <  2:1504109/140‑1 (MQ=255)
                                                                                               tAACTCAGCATTGCTGACAGGCTTATTGCTGGCGGTAAGTATTCCCCCC‑‑C‑‑CC‑‑TCGCGCCATGGTGGATGGTCGTGCTGGGTACGGTGTTTGCGGTGATCATCGCTAAACAGTTGTATGGCGGTCTGGGACAAAACCCGtt                                     <  1:428435/140‑1 (MQ=255)
                                                                                               tAACTCAGCATTGCTGACAGGCTTATTGCTGGCGGTAAGTATTCCCCCC‑‑C‑‑CC‑‑TCGCGCCATGGTGGATGGTCGGGCTGGGGACGGGGTTTGCGGCGATCATCCCTAAACAGTTGTATGGCGGTCTGGGCCAAAAccccgt                                     >  2:1530344/1‑137 (MQ=255)
                                                                                               tAACTCAGCATTGCTGACAGGCTTATTGCTGGCGGTAAGTATTCCCCCC‑‑C‑‑CC‑‑TCGAGCCATGGTGGATGGTCGTGCTGGGTACGGTGTTTGCGGTGATCATCGCTAAACAGTTGTATGGCGGTCTGGGACAAAACCCGtt                                     >  1:58138/1‑140 (MQ=255)
                                                                                                  cTCAGCATTGCTGACCGGCTTACTGCTGGTGGTAAGTATTCCCCCC‑‑C‑‑CC‑‑TCGCGCCATGGTGGATGGTCGTGCTGGGTACGGTGTTTGCGGTGATCATCGCTAAACAGTTGTATGGCGGTCTGGGACAAAACCCGTTTaa                                  <  1:687885/140‑1 (MQ=255)
                                                                                                   tCAGCATTGCTGACAGGCTTATTGCTGGCGGTAAGTATTCCCCCC‑‑C‑‑CC‑‑TCGCGCCATGGTGGATGGTCGTGCTGGGTACGGTGTTTGCGGTGATCATCGCTAAACAGTTGTATGGCGGTCTGGGACAAAACCCGTTTAAt                                 <  1:1530344/140‑1 (MQ=255)
                                                                                                       cATTGCTGACAGGCTTATTGCTGGCGGTAAGTATTCCCCCC‑‑C‑‑CC‑‑TCGCGCCATGGTGGATGGTCGTGCTGGGTACGGTGTTTGCGGTGATCATCGCTAAACAGTTGTATGGCGGTCTGGGACAAAACCCGTTTAATCCgg                             <  2:1738961/140‑1 (MQ=255)
                                                                                                       cATTGCTGACAGGCTTATTGCTGGCGGTAAGTATTCCCCCC‑‑C‑‑CC‑‑TCGCGCCATGGTGGATGGTCGTGCTGGGTACGGTGTTTGCGGTGATCATCGCTAAACAGTTGTATGGCGGTCTGGGACAAAACCCGTTTAATCCgg                             >  2:2032831/1‑140 (MQ=255)
                                                                                                       cATTGCTGACAGGCTTATTGCTGGCGGTAAGTATTCCCCCC‑‑C‑‑CC‑‑TCGCGCCATGGTGGATGGTCGTGCTGGGTACGGTGTTTGCGGTGATCATCGCTAAACAGTTGTATGGCGGTCTGGGACAAAACCCGTTTAATCCgg                             <  1:259481/140‑1 (MQ=255)
                                                                                                        aTTGCTGACAGGCTTTTTGCTGTCGCTAATTATT‑‑‑CCCTC‑‑C‑‑CC‑‑TCGCGCCATGGTGGATGGTCGTGCTGGGTACGGTGTTTGCGGTGATCATCGCTAAACAGTTGTATGGCGGTCTGGGACAAAACCCGTTTAATCCGGCa                           <  1:823962/140‑1 (MQ=255)
                                                                                                         ttGCTGACAGGCTTATTGCTGGCGGTAAGTATTCCCCCC‑‑C‑‑CC‑‑TCGCGCCATGGTGGATGGTCGTGCTGGGTACGGTGTTTGCGGTGATCATCGCTAAACAGTTGTATGGCGGTCTGGGACAAAACCCGTTTAATCCGGCa                           <  2:2346387/140‑1 (MQ=255)
                                                                                                          tGCTGACAGGCTTATTGCTGGCGGTAAGTATTCCCCCCC‑‑C‑‑CC‑‑TCGCGCCATGGGGGATGGTCGTGCTGGGTACGGGGTTTGCGGGGATCATCGCTAAACAGTCGTATGGCAGTCAGGAACAAAATCCGTTTAATCCGGCa                           >  1:1028838/1‑140 (MQ=255)
                                                                                                          tGCTGACAGGCTTATTGCTGGCGGTAAGTATTCCCCCC‑‑C‑‑CC‑‑TCGCGCCATGGTGGATGGTCGTGCTGGGTACGGTGTTTGCGGTGATCATCGCTAAACAGTTGTATGGCGGTCTGGGACAAAACCCGTTTAATCCGGCaa                          >  2:1883414/1‑140 (MQ=255)
                                                                                                            cTGACAGGCTTATTGCTCGCGGAAAGTATT‑‑‑CCCCCGC‑‑CC‑‑ACGCGCCATCGTGGATGGTCGTGCTGGGTACGGTGTTTGCGGTGATCATCGCTAAACAGTTGTATGGCGGTCTGGGACAAAACCCGTTTAATCCGGCAATg                        <  2:1037986/140‑1 (MQ=255)
                                                                                                            cTGACAGGATTATTGCGTGCGGTAAGTATTAGCTCC‑‑C‑‑CC‑‑ACGCGCCATAGTGGATGGTCGTGCTGGGTACGGTGTTTGCGGTGATCATCGCTAAACAGTTGTATGGCGGTCTGGGACAAAACCCGTTTAATCCGGCAATg                        <  2:1204720/140‑1 (MQ=255)
                                                                                                                   gCTTATTGCTGGCGGTAAGTATTCCCCCC‑‑C‑‑CC‑‑TCGCGCCATGGTGGATGGTCGTGCTGGGTACGGTGTTTGCGGTGATCATCGCTAAACAGTTGTATGGCGGTCTGGGACAAAACCCGTTTAATCCGGCAATGAt                      <  1:1373042/135‑1 (MQ=255)
                                                                                                                   gCTTATTGCTGGCGGTAAGTATTCCCCCC‑‑C‑‑CC‑‑TCGCGCCATGGTGGATGGTCGTGCTGGGTACGGTGTTTGCGGTGATCATCGCTAAACAGTTGTATGGCGGTCTGGGACAAAACCCGTTTAATCCGGCAATGAt                      >  2:1373042/1‑135 (MQ=255)
                                                                                                                    cTTATTGCTGGCGGTAAGTATTCCCCCC‑‑C‑‑CC‑‑TCGCGCCATGGTGGATGGTCGTGCTGGGTACGGTGTTTGCGGTGATCATCGCTAAACAGTTGTATGGCGGTCTGGGACAAAACCCGTTTAAt                                 >  2:1918555/1‑123 (MQ=255)
                                                                                                                    cTTATTGCTGGCGGTAAGTATTCCCCCC‑‑C‑‑CC‑‑TCGCGCCATGGTGGATGGTCGTGCTGGGTACGGTGTTTGCGGTGATCATCGCTAAACAGTTGTATGGCGGTCTGGGACAAAACCCGTTTAAt                                 <  1:1918555/123‑1 (MQ=255)
                                                                                                                      tATTGCTGGCGGTAAGTATTCCCCCC‑‑C‑‑CC‑‑TCGCGCCATGGTGGATGGTCGTGCTGGGTACGGTGTTTGCGGTGATCATCGCTAAACAGTTGTATGGCGGTCTGGGACAAAACCCGTTTAATCCGGCAATGATTGGTTAtg              <  1:39556/140‑1 (MQ=255)
                                                                                                                      tATTGCTGGCGGTAAGTATTCCCCCC‑‑C‑‑CC‑‑TCGCGCCATGGTGGATGGTCGTGCTGGGTACGGTGTTTGCGGTGATCATCGCTAAACAGTTGTATGGCGGTCTGGGACAAAACCCGTTTAATCCGGCAATGATTGGTTAtg              <  1:367318/140‑1 (MQ=255)
                                                                                                                      tATTGCTGGCGGTAAGTATTCCCCCC‑‑C‑‑CC‑‑TCGCGCCATGGTGGATGGTCGTGCTGGGTACGGTGTTTGCGGTGATCATCGCTAAACAGTTGTATGGCGGTCTGGGACAAAACCCGTTTAATCCGGCAATGATTGGTTAtg              <  1:2111474/140‑1 (MQ=255)
                                                                                                                      tATTGCTGGCGGTAAGTATTCCCCCC‑‑C‑‑CC‑‑TCGCGCCATGGTGGATGGTCGTGCTGGGTACGGTGTTTGCGGTGATCATCGCTAAACAGTTGTATGGCGGTCTGGGACAAAACCCGTTTAATCCGGCAATGATTGGTTAtg              <  1:1883414/140‑1 (MQ=255)
                                                                                                                         tGCTGGCGGTAAGTATTCCCCCC‑‑C‑‑CC‑‑TCGCGCCATGGTGGATGGTCGTGCTGGGTACGGTGTTTGCGGTGATCATCGCTAAACAGTTGTATGGCGGTCTGGGACAAAACCCGTTTAATCCGGCAATGATTGGTTATGTgg           <  1:1485248/140‑1 (MQ=255)
                                                                                                                            tGGCGGTAAGTATT‑‑‑CCCC‑‑C‑‑CCCTTCGCGCAATGGTGGTTGGTCGTGCTGGGTACGGTGTTTGCGGTGATCATCGCTAAACAGTTGTATGGCGGTCTGGGACAAAACCCGTTTAATCCGGCAATGATTGGTTATGTGGTCt        <  1:851924/140‑1 (MQ=255)
                                                                                                                             ggCGGTAAGTATTCCCCCC‑‑C‑‑CC‑‑TCGCGCCATGGTGGATGGTCGTGCTGGGTACGGTGTTTGCGGTGATCATCGCTAAACAGTTGTATGGCGGTCTGGGACAAAACCCGTTTAATCCGGCAATGATTGGTTAtg              >  1:449280/1‑133 (MQ=255)
                                                                                                                             ggCGGTAAGTATTCCCCCC‑‑C‑‑CC‑‑TCGCGCCATGGTGGATGGTCGTGCTGGGTACGGTGTTTGCGGTGATCATCGCTAAACAGTTGTATGGCGGTCTGGGACAAAACCCGTTTAATCCGGCAATGATTGGTTAtg              <  1:1253140/133‑1 (MQ=255)
                                                                                                                             ggCGGTAAGTATTCCCCCC‑‑C‑‑CC‑‑TCGCGCCATGGTGGATGGTCGTGCTGGGTACGGTGTTTGCGGTGATCATCGCTAAACAGTTGTATGGCGGTCTGGGACAAAACCCGTTTAATCCGGCAATGATTGGTTAtg              <  2:449280/133‑1 (MQ=255)
                                                                                                                             ggCGGTAAGTATTCCCCCC‑‑C‑‑CC‑‑TCGCGCCATGGTGGATGGTCGTGCTGGGTACGGTGTTTGCGGTGATCATCGCTAAACAGTTGTATGGCGGTCTGGGACAAAACCCGTTTAATCCGGCAATGATTGGTTATGTGGTCtt       >  2:2251105/1‑140 (MQ=255)
                                                                                                                             ggCGGTAAGTATTCCCCCC‑‑C‑‑CC‑‑TCGCGCCATGGTGGATGGTCGTGCTGGGTACGGTGTTTGCGGTGATCATCGCTAAACAGTTGTATGGCGGTCTGGGACAAAACCCGTTTAATCCGGCAATGATTGGTTATGTGGTCtt       >  1:1266093/1‑140 (MQ=255)
                                                                                                                             ggCGGTAAGTATTCCCCCC‑‑C‑‑CC‑‑TCGCGCCACGGTGGATGGTCGTGCTGGGTACGGTGTTTGCGGTGATCATCGCTAAACAGTTGTATGGCGAGCTGGGACCAAACCCGGTTAATCCGGCAATGATTGGTTATGGGGTCtt       >  2:883278/1‑140 (MQ=255)
                                                                                                                                 gTAAGTATTCCCCCC‑‑C‑‑CC‑‑TCGCGCCATGGTGGATGGTCGTGCTGGGTACGGTGTTTGCGGTGATCATCGCTAAACAGTTGTATGGCGGTCTGGGACAAAACCCGTTTAATCCGGCAATGATTGGTTATGTGGTCTTACTg   <  2:2118493/140‑1 (MQ=255)
                                                                                                                                 gTAAGTATTCCCCCC‑‑C‑‑CC‑‑TCGCGCCATGGTGGATGGTCGTGCTGGGTACGGTGTTTGCGGTGATCATCGCTAAACAGTTGTATGGCGGTCTGGGACAAAACCCGTTTAATCCGGCAATGATTGGTTATGTGGTCTTACTg   <  2:58138/140‑1 (MQ=255)
                                                                                                                                 gTAAGTATTCCCCCC‑‑C‑‑CC‑‑TCGCGCCATGGTGGATGGTCGTGCTGGGTACGGTGTTTGCGGTGATCATCGCTAAACAGTTGTATGGCGGTCTGGGACAAAACCCGTTTAATCCGGCAATGATTGGTTATGTGGTCTTACTg   <  1:406894/140‑1 (MQ=255)
                                                                                                                                 gTAAGTATT‑‑‑CCCC‑‑CCGCC‑‑TCGCTCCATGGTGGATGGTTGTGCTGGGTACGGTGTTTGCGGTGATCATCGCTAAACAGTTGTATGGCGGTCTGGGACAAAACCCGTTTAATCCGGCAATGATTGGTTATGTGGTCTTACTg   <  1:1031798/140‑1 (MQ=255)
                                                                                                                                 gTAAGTATT‑‑‑CCCC‑‑C‑‑CCCGTCGCTCCATGGTGCATGTTCGTGCTGGGTATGGTGTTTGCGGTGATCATCGCTAAACAGTTGTATGGCGGTCTGGGACAAAACCCGTTTAATCCGGCAATGATTGGTTATGTGGTCTTACTg   <  2:816693/140‑1 (MQ=255)
                                                                                                                                  tAAGTATTCCCCCC‑‑C‑‑CC‑‑TCGCGCCATGGTGGATGGTCGTGCTGGGTACGGTGTTTGCGGTGATCATCGCTAAACAGTTGTATGGCGGTCTGGGACAAAACCCGTTTAATCCGGCAATGATTGGTTATGTGGTCTTACTGa  <  1:1294996/140‑1 (MQ=255)
                                                                                                                                          ||                                                                                                                                          
CTCGTTCAGATCCTGTTGGCATCGGTTAGTGCTCTGTTAGCCGAAGCTCTCGTACTCAAACTACGCAAGCAGTCGGTAGCCGCAACGTTGAAAGATAACTCAGCATTGCTGACAGGCTTATTGCTGGCGGTAAGTATT‑‑‑CCCC‑‑C‑‑CC‑‑TCGCGCCATGGTGGATGGTCGTGCTGGGTACGGTGTTTGCGGTGATCATCGCTAAACAGTTGTATGGCGGTCTGGGACAAAACCCGTTTAATCCGGCAATGATTGGTTATGTGGTCTTACTGA  >  NC_000913/1709271‑1709538

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