Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,008,391 T→C 6.2% W39R (TGG→CGG)  yedE → UPF0394 family sulphur transport domain‑containing inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,008,3910TC 6.2% 195.5 / 2.6 64W39R (TGG→CGG) yedEUPF0394 family sulphur transport domain‑containing inner membrane protein
Reads supporting (aligned to +/- strand):  ref base T (26/34);  new base C (2/2);  total (28/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.81e-01

TGTTTCAAACCTTCTGAGTTCAGAGGCTACATTCATGTCATGGCAGCAATTCAAACACGCCTGGTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGGTATTCTCTCTACTTACTATTTTGGCATTACTGGCACCTTTTGGGCTGTCACGGGTGAATTTACCCGTTGGGGCGGCCAGCTCCTGCAACTGTTCGGCGTCCATGCTGAAGAGTGGGGTTACTTTAAAATTATCCATCTGGAAGGATCGCCATTAACCCGCATCGACGGGATGATGATCCTCGGTA  >  NC_000913/2008243‑2008535
                                                                                                                                                    |                                                                                                                                                
tGTTTCAAACATTCTGAGTTCAGAGGCTACATTCATGTCATGGCAGCAATTCAAACACGCCTGGTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGGTATTCTCTCTACTTACTATTTTGGCATTACTGGCACCtttt                                                                                                                                                  >  4:313329/1‑149 (MQ=255)
 gattCAAACCTTCTGAGTTCAGAGGCTACATTCTTGTCATGGCAGCAATTCAAACACGCCTGGTTGATTAAATTCTGTGACCCCATCCCTGCGGTCATCGCTGCTGGTATTCTCTCTACTTACTATTTTGGCATTACTGGCACCtttt                                                                                                                                                  <  3:207747/146‑1 (MQ=255)
  tttCAAACCTTCTGAGTTCAGAGGCTACATTCATGTCATGGCAGCAATTCAAACACGCCTGGTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGGTATTCTCTCTACTTACTATTTTGGCATTACTGGCACCTTTTgg                                                                                                                                                <  6:189357/149‑1 (MQ=255)
  tttCAAACCTTCTGAGTTCAGAGGCTACATTCATGTCATGGCAGCAATTCAAACACGCCTGGTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGGTATTCTCTCTACTTACTATTTTGGCATTACTGGCACCTTTTgg                                                                                                                                                <  7:32721/149‑1 (MQ=255)
  tttCAAACCTTCTGAGTTCAGAGGCTACATTCATGCCATGGCAGCAATTCAAACACTCCTGGTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGGTAATCTCTCTACTTACTATTTTGGCATTACAGGCACCTTTTg                                                                                                                                                 <  4:85639/148‑1 (MQ=255)
      aaaCCTTCTGAGTTCAGAGGCTACATTCATGTCATGGCAGCAATTCAAACACGCCTGGTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGGTATTCTCTCTACTTACTATTTTGGCATTACTGGCACCTTTTGGGCTg                                                                                                                                            >  6:232767/1‑149 (MQ=255)
       aaCCTTCTGAGTTCAGAGGCTACATTCATGTCATGGCAGCAATTCAAACACGCCTGGTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGGTATTCTCTCTACTTACTATTTTGGCATTACTGGCACCTTTTGGGCTGt                                                                                                                                           <  5:232767/149‑1 (MQ=255)
        aCCTTCTGAGTTCAGAGGCTACATTCATGTCATGGCAGCAATTCAAACACGCCTGGTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGGTATTCTCTCTACTTACTATTTTGGCATTACTGGCACCTTTTGGGCTGTc                                                                                                                                          <  3:232268/149‑1 (MQ=255)
           ttCTGAGTTCAGAGGCTACATTCATGTCATGGCAGCAATTCAAACACGCCTGGTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGGTATTCTCTCTACTTACTATTTTGGCATTACTGGCACCTTTTGGGCTGTCACg                                                                                                                                       >  7:379140/1‑149 (MQ=255)
             cTGAGTTCAGAGGCTACATTCATGTCATGGCAGCAATTCAAACACGCCTGGTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGGTATTCTCTCTACTTACTATTTTGGCATTACTGGCACCTTTTGGGCTGTCACgg                                                                                                                                      <  6:210852/148‑1 (MQ=255)
                  ttCAGAGGCTACATTCATGTCATGGCAGCAATTCAAACACGCCTGGTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGGTATTCTCTCTACTTACTATTTTGGCATTACTGGCACCTTTTGGGCTGTCACGGGTGAAt                                                                                                                                >  8:245785/1‑149 (MQ=255)
                    cAGAGGCTACATTCATGTCATGGCAGCAATTCAAACACGCCTGGTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGTTATTCTCTCTACTTACTATTTTGGCATTACTGGCACCTTTTGGGCTGTCACGGGTGAAtt                                                                                                                               <  2:118282/148‑1 (MQ=255)
                    aagagGCTACATTCATGTCATGGCAGCAATTCAAACACGCCTGGTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGGTATTCTCTCTACTTACTATTTTGGCATTACTGGCACCTTTTGGGCTGTCACGGGTGAAtt                                                                                                                               >  5:256518/2‑148 (MQ=255)
                      gagCCTACATTCATGTCATGGCAGCAATTCAAACACGCCTGGTTGATTAAATTCTGGGCCCCCATCCCTGCGGGCATCGCGGCGGGTATTCTCTCTACTTACTATTTTGGCATTACTGGCACCTTTTGGGCTGTCACGGGTGAATTTAc                                                                                                                            <  8:330698/149‑1 (MQ=255)
                             cATTCATGTCATGGCAGCAATTCAAACACGCCTGGTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGGTATTCTCTCTACTTACTATTTTGGCATTACTGGCACCTTTTGGGCTGTCACGGGTGAATTTAc                                                                                                                            >  6:13407/1‑142 (MQ=255)
                             cATTCATGTCATGGCAGCAATTCAAACACGCCTGGTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGGTATTCTCTCTACTTACTATTTTGGCATTACTGGCACCTTTTGGGCTGTCACGGGTGAATTTAc                                                                                                                            <  5:13407/142‑1 (MQ=255)
                             cATTCATGTCATGGCAGCAATTCAAACACGCCTGGTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGGTATTCTCTCTACTTACTATTTTGGCATTACTGGCACCTTTTGGGCTGTCACGGGTGAATTTACCCGTTgg                                                                                                                     >  4:12877/1‑149 (MQ=255)
                             cATTCATGTCATGGCAGCAATTCAAACACGCCTGGTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGGTATTCTCTCTACTTACTATTTTGGCATTACTGGCACCTTTTGGGCTGTCACGGGTGAATTTACCCGTTg                                                                                                                      <  7:374525/148‑1 (MQ=255)
                              aTTCATGTCATGGCAGCAATTCAAACACGCCTGGTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGGTATTCTCTCTACTTACTATTTTGGCATTACTGGCACCTTTTGGGCTGTCACGGGTGAATTTACCCGTTggg                                                                                                                    >  2:253837/1‑149 (MQ=255)
                                        tgGCAGCAATTCAAACACGCCTGGTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGGTATTCTCTCTACTTACTATTTTGGCATTACTGGCACCTTTTGGGCTGTCACGGGTGAATTTACCCGTTGGGGCGGCCAGCt                                                                                                          >  7:3821/1‑149 (MQ=255)
                                                     aaCACGCCTGGTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGGTATTCTCTCTACTTACTATTTTGGCATTACTGGCACCTTTTGGGCTGTCACGGGTGAATTTACCCGTTGGGGCGGCCAGCTCCTGCAACTGTTc                                                                                             >  1:177199/1‑149 (MQ=255)
                                                     aaCACGCCTGGTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGGTATTCTCTCTACTTACTATTTTGGCATTACTGGCACCTTTTGGGCTGTCACGGGTGAATTTACCCGTTGGGGCGGCCAGCTCCTGCAACTGTTc                                                                                             <  7:245785/149‑1 (MQ=255)
                                                      acacGCCTGGTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGGTATTCTCTCTACTTACTATTTTGGCATTACTGGCACCTTTTGGGCTGTCACGGGTGAATTTACCCGTTGGGGCGGCCAGCTCCTGCAACTGTTCg                                                                                            <  8:379140/149‑1 (MQ=255)
                                                         cGCCTGGTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGGTATTCTCTCTACTTACTATTTTGGCATTACTGGCACCTTTTGGGCTGTCACGGGTGAATTTACCCGTTGGGGCGGCCAGCTCCTGCAACTGTTCGGCg                                                                                         <  3:313329/149‑1 (MQ=255)
                                                            cTGGTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGGTATTCTCTCTACTTACTATTTTGGCATTACTGGCACCTTTTGGGCTGTCACGGGTGAATTTACCCGTTGGGGCGGCCAGCTCCTGCAACTGTTCGGCGTc                                                                                       >  3:316351/1‑148 (MQ=255)
                                                            cTGGTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGCGTATTCTCTCTACTTACTATTTTGGCATTACTGGCACCTTTTGGGCTGTCACGGGTGAATTTACCCGTTGGGGCGGCCAGCTCCTGCAACTGTTCGGCGTcc                                                                                      <  5:364054/149‑1 (MQ=255)
                                                                  gATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGGTATTCTCTCTACTTACTATTTTGGCATTACTGGCACCTTTTGGGCTGTCACGGGTGAATTTACCCGTTg                                                                                                                      <  5:20700/111‑1 (MQ=255)
                                                                   aTTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGGTATTCTCTCTACTTACTATTTTGGCATTACTGGCACCTTTTGGGCTGTCACGGGTGAATTTACCCGTTGGGGCGGCCAGCTCCTGCAACTGTTCGGCGTCCAt                                                                                    >  1:155113/1‑144 (MQ=255)
                                                                   aTTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGGTATTCTCTCTACTTACTATTTTGGCATTACTGGCACCTTTTGGGCTGTCACGGGTGAATTTACCCGTTGGGGCGGCCAGCTCCTGCAACTGTTCGGCGTCCAt                                                                                    <  2:155113/144‑1 (MQ=255)
                                                                     tAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGGTATTCTCTCTACTTACTATTTTGGCATTACTGGCACCTTTTGGGCTGTCACGGGTGAATTTACCCGTTGGGGCGGCCAGCTCCTGCAACTGtt                                                                                              >  1:167683/1‑132 (MQ=255)
                                                                     tAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGGTATTCTCTCTACTTACTATTTTGGCATTACTGGCACCTTTTGGGCTGTCACGGGTGAATTTACCCGTTGGGGCGGCCAGCTCCTGCAACTGtt                                                                                              <  2:167683/132‑1 (MQ=255)
                                                                         ttCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGGTATTCTCTCTACTTACTATTTTGGCATTACTGGCACCTTTTGGGCTGTCACGGGTGAATTTACCCGTTGGGGCGGCCAGCTCCTGCAACTGTTCGGCGTCCATGCTGAAGAGTg                                                                         <  3:315235/149‑1 (MQ=255)
                                                                              ggCCCCCATCCCTGCGGTCATCGCGGCGGGTATTCTCTCTACTTACTATTTTGGCATTACTGGCACCTTTTGGGCTGTCACGGGTGAATTTACCCGTTGGGGCGGACAGCTCCTGCAACTGTTCGGCGTCCATGCTGAAGAGTggg                                                                       <  2:244339/146‑1 (MQ=255)
                                                                                  cccATCCCTGCGGTCATCGCGGCGGGTATTCTCTCTACTTACTATTTTGGCATTACTGGCACCTTTTGGGCTGTCACGGGTGAATTTACCCGTTGGGGCGGCCAGCTCCTGCAACTGTTCGGCGTCCATGCt                                                                                 >  2:268116/1‑132 (MQ=255)
                                                                                  cccATCCCTGCGGTCATCGCGGCGGGTATTCTCTCTACTTACTATTTTGGCATTACTGGCACCTTTTGGGCTGTCACGGGTGAATTTACCCGTTGGGGCGGCCAGCTCCTGCAACTGTTCGGCGTCCATGCt                                                                                 <  1:268116/132‑1 (MQ=255)
                                                                                   ccATCCCTGCGGTCATCGCGGCGGGTATTCTCTCTACTTACTATTTTGGCATTACTGGCACCTTTTGGGCTGTCACGGGTGAATTTACCCGTTGGGGCGGCCAGCTCCTGCAACTGTTCGGCGTCCATGCTGAAGAGTGGGGTTACttt                                                               <  1:253837/149‑1 (MQ=255)
                                                                                        ccTGCGGTCATCGCGGCGGGTATTCTCTCTACTTACTATTTTGGCATTACTGGCACCTTTTGGGCTGTCACGGGTGAATTTACCCGTTGGGGCGGCCAGCTCCTGCAACTGTTCGGCGTCCATGCTGAAGAGTggg                                                                       <  4:216804/136‑1 (MQ=255)
                                                                                        ccTGCGGTCATCGCGGCGGGTATTCTCTCTACTTACTATTTTGGCATTACTGGCACCTTTTGGGCTGTCACGGGTGAATTTACCCGTTGGGGCGGCCAGCTCCTGCAACTGTTCGGCGTCCATGCTGAAGAGTggg                                                                       >  3:216804/1‑136 (MQ=255)
                                                                                          tGCGGTCATCGCGGCGGGTATTCTCTCTACTTACTATTTTGGCATTACTGGCACCTTTTGGGCTGTCACGGGTGAATTTACCCGTTGGGGCGGCCAGCTCCTGCAACTGTTCGGCGTCCATg                                                                                   >  5:176354/1‑122 (MQ=255)
                                                                                          tGCGGTCATCGCGGCGGGTATTCTCTCTACTTACTATTTTGGCATTACTGGCACCTTTTGGGCTGTCACGGGTGAATTTACCCGTTGGGGCGGCCAGCTCCTGCAACTGTTCGGCGTCCATg                                                                                   <  6:176354/122‑1 (MQ=255)
                                                                                                 aTCGCGGCGGGTATTCTCTCTACTTACTATTTTGGCATTACTGGCACCTTTTGGGCTGTCACGGGTGAATTTACCCGTTGGGGCGGCCAGCTCCTGCAACTGTTCGGCGTCCATGCTGAAGAGTGGGGTTACTTTAAAATTATCCATCt                                                 <  7:15514/149‑1 (MQ=255)
                                                                                                 aTCGCGGCGGGTATTCTCTCTACTTACTATTTTGGCATTACTGGCACCTTTTGGGCTGTCACGGGTGAATTTACCCGTTGGGGCGGCCAGCTCCTGCAACTGTTCGGCGGCCATGCTGAAGAGTGGGGTTACTTTAAAATTATCCATc                                                  <  3:12877/148‑1 (MQ=255)
                                                                                                    gcggcgGGTATTCTCTCTACTTACTATTTTGGCATTACTGGCACCTTTTGGGCTGTCACGGGTGAATTTACCCGTTGGGGCGGCCAGCTCCTGCAACTGTTCGGCGTCCATGCTGAAGAGTGGGGTTACTTTAAAATTATCCATCTgg                                               >  1:49761/1‑148 (MQ=255)
                                                                                                     cggcggGTATTCTCTCTACTTACTATTTTGGCAGTACTGGCACCTTTTGGGCTGTCACGGGTGAATTTACCCGTTGGGGCGGCCAGCTCCTGCAACTGTTCGGCGTCCATGCTGAAGAGTGGGGTTACTTTAAAATTATCCATCTGGaa                                             <  6:201620/149‑1 (MQ=255)
                                                                                                        ggggTATTCTCTCTACTTACTATTTTGGCATTACTGGCACCTTTTGGGCTGTCACGGGTGAATTTACCCGTTGGGGCGGCCAGCTCCTGCAACTGTTCGGCGTCCATGCTGAAGAGTGGGGTTACTTTAAAATTATCCATCTGGAAGGa                                          >  6:272143/2‑149 (MQ=255)
                                                                                                          ggTATTCTCTCTACTTACTATTTTGGCATTACTGGCACCTTTTGGGCTGTCACGGGTGAATTTACCCGTTGGGGCGGCCAGCTCCTGCAACTGTTCGGCGTCCATTCTGAAGAGTGGTGTTACTTTAAAATTATCCATCTGGAAGGATc                                        <  2:49761/149‑1 (MQ=255)
                                                                                                                 tctctACTTACTATTTTGGCATTACTGGCACCTTTTGGGCTGTCACGGGTGAATTTACCCGTTGGGGCGGCCAGCTCCTGCAACTGTTCGGCGTCCATGCTGAAGAGTGGGGTTACTTTAAAATTATCCATCTGGAAGGATCGCCAtt                                  >  7:300890/1‑148 (MQ=255)
                                                                                                                      acttactATTTTGGCATTACTGGCACCTTTTGGGCTGTCACGGGTGAATTTACCCGTTGGGGCGGCCAGCTCCTGCAACTGTTCGGCGTCCATGCTGAAGAGTGGGGTTACTTTAAAATTATCCATCTGGaa                                             >  5:311930/1‑132 (MQ=255)
                                                                                                                      acttactATTTTGGCATTACTGGCACCTTTTGGGCTGTCACGGGTGAATTTACCCGTTGGGGCGGCCAGCTCCTGCAACTGTTCGGCGTCCATGCTGAAGAGTGGGGTTACTTTAAAATTATCCATCTGGaa                                             <  6:311930/132‑1 (MQ=255)
                                                                                                                       cttactATTTTGGCATTACTGGCACCTTTTGGGCTGTCACGGGTGAATTTACCCGTTGGGGCGGCCAGCTCCTGCAACTGTTCGGCGTCCATGCTGAAGAGTGGGGTTACTTTAAAATTATCCATCTGGAAGGATCGCCATTAACCCg                            >  7:41479/1‑148 (MQ=255)
                                                                                                                        ttactATTTTGGCATTACTGGCACCTTTTGGGCTGTCACGGGTGAATTTACCCGTTGGGGCGGCCAGCTCCTGCAACTGTTCGGCGTCCATGCTGAAGAGTGGGGTTACTTTAAAATTATCCATCTGGAAGGATCGCCATTAACCCGCa                          >  8:155803/1‑149 (MQ=255)
                                                                                                                        ttactATTTTGGCATTACTGGCACCTTTTGGGCTGTCACGGGTGAATTTACCCGTTGGGGCGGCCAGCTCCTGCAACTGTTCGGCGTCCATGCTGAAGAGTGGGGTTACTTTAAAATTATCCATCTGGAAGGATCGCCATTAACCCGCa                          >  1:102347/1‑149 (MQ=255)
                                                                                                                                 tGGCATTACTGGCACCTTTCGGGCTGTCACGGGTGAATTTACCCGTTGGGGCGGCCAGCTCCTGCAACTGTTCGGCGTCCATGCTGAAGAGTGGGGTTACTTTAAAATTATCCATCTGGAAGGATCGCCATTAcccc                             >  7:208012/1‑137 (MQ=255)
                                                                                                                                 tGGCATTACTGGCACCTTTCGGGCTGTCACGGGTGAATTTACCCGTTGGGGCGGCCAGCTCCTGCAACTGTTCGGCGTCCATGCTGAAGAGTGGGGTTACTTTAAAATTATCCATCTGGAAGGATCGCCATTAAccc                             >  7:213522/1‑137 (MQ=255)
                                                                                                                                 tGGCATTACTGGCACCTTTCGGGCTGTCACGGGTGAATTTACCCGTTGGGGCGGCCAGCTCCTGCAACTGTTCGGCGTCCATGCTGAAGAGTGGGGTTACTTTAAAATTATCCATCTGGAAGGATCGCCATTAAccc                             <  8:213522/137‑1 (MQ=255)
                                                                                                                                 tGGCATTACTGGCACCTTTCGGGCTGTCACGGGTGAATTTACCCGTTGGGGCGGCCAGCTCCTGCAACTGTTCGGCGTCCATGCTGAAGAGTGGGGTTACTTTAAAATTATCCATCTGGAAGGATCGCCATTAAccc                             <  8:208012/137‑1 (MQ=255)
                                                                                                                                  ggCATTACTGGCACCTTTTGGGCTGTCACGGGTGAATTTACCCGTTGGGGCGGCCAGCTCCTGCAACTGTTCGGCGTCCATGCTGAAGAGTGGGGTTACTTTAAAATTATCCATCTGGAAGGATCGCCATTAACCCGCATCGACGGga                 <  8:300890/148‑1 (MQ=255)
                                                                                                                                  ggCATTACTGGCACCTTTTGGGCTGTCACGGGTGAATTTACCCGTTGGGGCGGCCAGCTCCTGCAACTGTTCGGCGTCCATGCTGAAGAGTGGGGTTACTTTAAAATTATCCATCTGGAAGGATCGCCATTAACCCGCATCGACGGga                 <  7:155803/148‑1 (MQ=255)
                                                                                                                                         cTGGCACCTTTTGGGCTGTCACGGGTGAATTTACCCGTTGGGGCGGCCAGCTCCTGCAACTGTTCGGCGTCCATGCTGAAGAGTGGGGTTACTTTAAAATTATCCATCTGGAAGGATCGCCATTAACCCGCa                          >  5:169126/1‑132 (MQ=255)
                                                                                                                                         cTGGCACCTTTTGGGCTGTCACGGGTGAATTTACCCGTTGGGGCGGCCAGCTCCTGCAACTGTTCGGCGTCCATGCTGAAGAGTGGGGTTACTTTAAAATTATCCATCTGGAAGGATCGCCATTAACCCGCa                          <  6:169126/132‑1 (MQ=255)
                                                                                                                                         cTGGCACCTTTTGGGCTGTCACGGGTGAATTTACCCGTTGGGGCGGCCAGCTCCTGCAACTGTTCGGCGTCCATGCTGAAGAGTGGGGTTACTTTAAAATTATCCATCTGGAAGGATCGCCATTAACCCGCATCGACGGgatgatgat          <  8:41479/148‑1 (MQ=255)
                                                                                                                                           ggCACCTTTTGGGCTGTCACGGGGGAATTTACCCGTTGGGGCGGCCAGCTCCTGCAACTGTTCGGCGTCCATGCTGAAGAGTGGGGTTACTTTAAAATTATCCATCTGGAAGGATCGCCATTAACCCGCATCGCCGGGATGATGATCCt       >  3:211421/1‑149 (MQ=255)
                                                                                                                                            gCACCTTTTGGGCTGTCACGGGTGAATTTACCCGTTGGGGCGGCCAGCTCCTGCAACTGTTCGGCGTCCATGCTGAAGAGTGGGGTTACTTTAAAATTATCCATCTGGAAGGATCGCCATTAACCCGCATCGACGGGATGATGATCCt       <  5:272143/148‑1 (MQ=255)
                                                                                                                                              aCCTTTTGGGCTGTCACGGGTGAATTTACCCGTTGGGGCGGCCAGCTCCTGCAACTGTTCGGCGTCCATGCTGAAGAGTGGGGTTACTTTAAAATTATCCATCTGGa                                              <  6:85874/107‑1 (MQ=255)
                                                                                                                                              aCCTTTTGGGCTGTCACGGGTGAATTTACCCGTTGGGGCGGCCAGCTCCTGCAACTGTTCGGCGTCCATGCTGAAGAGTGGGGTTACTTTAAAATTATCCATCTGGa                                              >  5:85874/1‑107 (MQ=255)
                                                                                                                                                cTTTTGGGCTGTCACGGGTGAATTTACCCGTTGGGGCGGCCAGCTCCTGCAACTGTTCGGCGTCCATGCTGAAGAGTGGGGTTACTTTAAAATTATCCATCTGGAAGGATCGCCATTAACCCGCATCGACGGGATGATGATCCTCGGTa  >  5:27293/1‑149 (MQ=255)
                                                                                                                                                 ttttGGGCTGTCACGGGTGAATTTACCCGTTGGGGCGGCCAGCTCCTGCAACTGTTCGGCGTCCATGCTGAAGAGTGGGGTTACTTTAAAATTATCCATCTGGAAGGATCGCCATTAACCCGCATCGACGGgatgatga           >  5:359527/1‑139 (MQ=255)
                                                                                                                                                 ttttGGGCTGTCACGGGTGAATTTACCCGTTGGGGCGGCCAGCTCCTGCAACTGTTCGGCGTCCATGCTGAAGAGTGGGGTTACTTTAAAATTATCCATCTGGAAGGATCGCCATTAACCCGCATCGACGGgatgatga           <  6:359527/139‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                
TGTTTCAAACCTTCTGAGTTCAGAGGCTACATTCATGTCATGGCAGCAATTCAAACACGCCTGGTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGGTATTCTCTCTACTTACTATTTTGGCATTACTGGCACCTTTTGGGCTGTCACGGGTGAATTTACCCGTTGGGGCGGCCAGCTCCTGCAACTGTTCGGCGTCCATGCTGAAGAGTGGGGTTACTTTAAAATTATCCATCTGGAAGGATCGCCATTAACCCGCATCGACGGGATGATGATCCTCGGTA  >  NC_000913/2008243‑2008535

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