Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 548,908 G→T 7.0% L586L (CTG→CTT)  fusA → elongation factor G 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947548,9080GT 7.0% 211.1 / 6.6 72L586L (CTG→CTT) fusAelongation factor G 1
Reads supporting (aligned to +/- strand):  ref base G (34/33);  new base T (3/2);  total (37/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.89e-02

GACAAGGGCAACATCACAGAAGGTCTGGTGTTCTCCAACGAAGTCGTTGGTGGTGTGATTCCTAAGGAATACATCCCTGCGATCCAGAAGGGTATCGAAGAGCAGATGAAAAACGGCGTTGTTGC‑‑‑CGGTTATCCTCTGATC‑GGCCTGAAGGCCGCGGTATTCGACGGTTCGTACCACGACGTCGACTCCAACGAAATGGCGTTCAAGATCGCTGCTTCGA‑TGGCGACCAAGCAACTGGCCCAGAAGGGCGGTGGCGTGGTTCTCGAGCCGATCATGAAAGTCGAAGTTGT  >  NC_002947/548762‑549051
                                                                                                                                                      |                                                                                                                                                
gacAAGGGCAACATCACAGAAGGTCTGGTGTTCTCCAACGAAGTCGTTGGTGGTGTGATTCCTAAGGAATACATCCCTGCGATCCAGAAGGGTATCGAAGAGCAGATGAAAAACGGCGTTGTTGC‑‑‑CGGTTATCCTCTGATC‑GGCCTg                                                                                                                                                  >  7:130688/1‑147 (MQ=255)
gacAAGGGCAACATCACAGAAGGTCTGGTGTTCTCCAACGAAGTCGTTGGTGGTGTGATTCCTAAGGAATACATCCCTGCGATCCAGAAGGGTATCGAAGAGCAGATGAAAAACGGCGTTGTTGC‑‑‑CGGTTATCCTCTGATC‑GGCCTGaa                                                                                                                                                <  1:184476/149‑1 (MQ=255)
    aGGGCAACATCACAGAAGGTCTGGTGTTCTCCAACGAAGTCGTTGGTGGTGTGATTCCTAAGGAATACATCCCTGCGATCCAGAAGGGTATCGAAGAGCAGATGAAAAACGGCGTTGTTGC‑‑‑CGGTTATCCTCTGATC‑GGCCTGAAGGcc                                                                                                                                            >  4:127614/1‑149 (MQ=255)
       gCAACATCACAGAAGGTCTGGTGTTCTCCAACGAAGTCGTTTGTGGTGTGATTCCTAAGGAATACATCCCTGCGATCCAGAAGGGTATCGAAGAGCAGATGAAAAACGGCGTTGTTGC‑‑‑CGGTTATCCTCTGATC‑GGCCTGAAGGCcgcg                                                                                                                                         <  5:244995/149‑1 (MQ=255)
        cAACATCACAGAAGGTCTGGTGTTCTCCAACGAAGTCGTTGGTGGTGTGATTCCTAAGGAATACATCCCTGCGATCCAGAAGGGTATCGAAGAGCAGATGAAAAACGGCGTTGTTGC‑‑‑CGGTTATCCTCTGATC‑GGCCTGAAGGCCGCgg                                                                                                                                        <  7:188626/149‑1 (MQ=255)
         aaCATCACAGCAGGTCTGGTGTTCTCCAACGAAGTCGTTGGTGGTGTGATTCCTAAGGAATACATCCCTGCGATCCAGAAGGGTATCGAAGAGCAGATGAAAAACGGCGTTGTTGC‑‑‑CGGTTATCCTCTGATC‑GGCCTGAAg                                                                                                                                               <  2:187507/141‑1 (MQ=255)
         aaCATCACAGAAGGTCTGGTGTTCTCCAACGAAGTCGTTGGTGGTGTGATTCCTAAGGAATACATCCCTGCGATCCAGAAGGGTATCGAAGAGCAGATGAAAAACGGCGTTGTTGC‑‑‑CGGTTATCCTCTGATC‑GGCCTGAAggagaacg                                                                                                                                        >  1:187507/1‑142 (MQ=255)
            aTCACAGAAGGTCTGGTGTTCTCCAACGAAGTCGTTGGTGGTGTGATTCCTAAGGAATACATCCCTGCGATCCAGAAGGGTATCGAAGAGCAGATGAAAAACGGCGTTGTTGC‑‑‑CGGTTATCCTCTGATC‑GGCCTGAAGGCCGCGGTAtt                                                                                                                                    >  2:128202/1‑149 (MQ=255)
                        cTGGTGTTCTCCAACGAAGTCGTTGGTGGTGTGATTCCTAAGGAATACATCCCTGCGATCCAGAAGGGTATCGAAGAGCAGATGAAAAACGGCGTTGTTGC‑‑‑CGGTTATCCTCTGATC‑GGCCTGAAGGCCGCGGTATTCGACGGTTCGTa                                                                                                                        >  5:138292/1‑149 (MQ=255)
                          ggTGTTCTCCAACGAAGTCGTTGGTGGTGTGATTCCTAAGGAATACATCCCTGCGATCCAGAAGGGTATCGAAGAGCAGATGAAAAACGGCGTTGTTGC‑‑‑CGGTTATCCTCTGATC‑GGCCTGAAGGCCGCGGTATTCGACGGTTCGTAc                                                                                                                       >  1:15744/1‑148 (MQ=255)
                            tgtTCTCCAACGAAGTCGTTGGTGGTGTGATTCCTAAGGAATACATCCCTGCGATCCAGAAGGGTATCGAAGAGCAGATGAAAAACGGCGTTGTTGC‑‑‑CGGTTATCCTCTGATC‑GGCCTGAAGGCcgcg                                                                                                                                         <  6:23093/128‑1 (MQ=255)
                            tgtTCTCCAACGAAGTCGTTGGTGGTGTGATTCCTAAGGAATACATCCCTGCGATCCAGAAGGGTATCGAAGAGCAGATGAAAAACGGCGTTGTTGC‑‑‑CGGTTATCCTCTGATC‑GGCCTGAAGGCcgcg                                                                                                                                         >  5:23093/1‑128 (MQ=255)
                            tgtTCTCCAACGAAGTCGTTGGTGGTGTGATTCCTAAGGAATACATCCCTGCGATCCAGAAGGGTATCGAAGAGCAGATGAAAAACGGCGTTGTTGC‑‑‑CGGTTATCCTCTGATC‑GGCCTGAAGGCCGCGGTATTCGACGGTTCGTACCac                                                                                                                    <  6:138292/149‑1 (MQ=255)
                             gtTCTCCAACGAAGTCGTTGGTGGTGTGATTCCTAAGGAATACATCCCTGCGATCCAGAAGGGTATCGAAGAGCAGATGAAAAACGGCGTTGTTGC‑‑‑CGGTTATCCTCTGATC‑GGCCTGAAGGCCGCGGTATTc                                                                                                                                   <  2:194174/133‑1 (MQ=255)
                             gtTCTCCAACGAAGTCGTTGGTGGTGTGATTCCTAAGGAATACATCCCTGCGATCCAGAAGGGTATCGAAGAGCAGATGAAAAACGGCGTTGTTGC‑‑‑CGGTTATCCTCTGATC‑GGCCTGAAGGCCGCGGTATTc                                                                                                                                   >  1:194174/1‑133 (MQ=255)
                                caccAACGACTTCGTTGGTGGTGTGATTCCTAAGGAATACATCCCTGCGATCCAGAAGGGTATCGAAGAGCAGATGAAAAACGGCGTTGTTGC‑‑‑CGGTTATCCTCTGATC‑GGCCTGAAGGCCGCGGTATTCGACGGTTCGTACCacgacg                                                                                                                >  2:78173/3‑149 (MQ=255)
                                 tcCAACGAAGTCGTTGGTGGTGTGATTCCTAAGTAATACATCCCTGCGATCCAGAAGGGTATCGAAGAGCAGATGAAAAACGGCGTTGTTGC‑‑‑CGGTTATCCTCTGATC‑GGCCTGAAGGCcgcg                                                                                                                                         <  4:15180/123‑1 (MQ=255)
                                 tcCAACGAAGTCGTTGGTGGTGTGATTCCTAAGGAATACATCCCTGCGATCCAGAAGGGTATCGAAGAGCAGATGAAAAACGGCGTTGTTGC‑‑‑CGGTTATCCTCTGATC‑GGCCTGAAGGCCGCGGTATTCGACGGTTCGTACCACGACGt                                                                                                               >  3:149546/1‑149 (MQ=255)
                                  ccAACGAAGTCGTTGGTGGTGTGATTCCTAAGGAATACATCCCTGCGATCCAGAAGGGTATCGAAGAGCAGATGAAAAACGGCGTTGTTGC‑‑‑CGGTTATCCTCTGATC‑GGCCTGAAGGCCGCGGTATTCGACGGTTCGt                                                                                                                         >  5:212996/1‑138 (MQ=255)
                                  ccAACGAAGTCGTTGGTGGTGTGATTCCTAAGGAATACATCCCTGCGATCCAGAAGGGTATCGAAGAGCAGATGAAAAACGGCGTTGTTGC‑‑‑CGGTTATCCTCTGATC‑GGCCTGAAGGCCGCGGTATTCGACGGTTCGt                                                                                                                         <  6:212996/138‑1 (MQ=255)
                                  ccAACGAAGTCGTTAGTGGTGTGATTCCTAAGGAATACATCCCTGCGATCCAGAAGGGTATCGAAGAGCAGATGAAAAACGGCGTTGTTGC‑‑‑CGGTTATCCTCTGATC‑GGCCTGAAGGCCGCGGTATTCGACGGTTCGTACCACGACGTc                                                                                                              <  8:272046/149‑1 (MQ=255)
                                   cAACGAAGTCGTTGGTGGTGTGATTCCTAAGGAATACATCCCTGTGATCCAGAAGGGTATCGAAGAGCAGATGAAAAACGGCGTTGTTGC‑‑‑CGGTTATCCTCTGATC‑GGCCTTCAGGCcgat                                                                                                                                         >  4:62389/1‑119 (MQ=255)
                                   cAACGAAGTCGTTGGTGGTGTGATTCCTAAGGAATACATCCCTGTGATCCAGAAGGGTATCGAAGAGCAGATGAAAAACGGCGTTGTTGC‑‑‑CGGTTATCCTCTGATC‑GGCCTTCAGGCcgat                                                                                                                                         <  3:62389/121‑3 (MQ=255)
                                   cAACGAAGTCGTTGGTGGTGTGATTCCTAAGGAATACATCCCTGCGATACAGAAGGGTATCGAAGAGCAGATGAAAAACGGCGTTGTTGC‑‑‑CGGTTATCCTCTGATC‑GGCCTGAAGGCCGCGGTATTAGACGGTTCGTACCACGACGTc                                                                                                              >  3:294237/1‑148 (MQ=255)
                                    aaCGAAGTCGTTGGTGGTGTGATTCCTAAGGAATACATCCCTGCGATCCAGAAGGGTATCGAAGAGCAGATGAAAAACGGCGTTGTTGC‑‑‑CGGTTATCCTCTGATC‑GGCCTGAAGGCCGCGGTATTCGACGGTTCGTACCACGACGTCg                                                                                                             >  8:104736/1‑148 (MQ=255)
                                    aaCGAAGTCGTTGGTGGTGTGATTCCTAAGGAATACATCCCTGCGATCCAGAAGGGTATCGAAGAGCAGATGAAAAACGGCGTTGTTGC‑‑‑CGGTTATCCTCTGATC‑GGCCTGAAGGCCGCGGTATTCGACGGTTCGTACCACGACGTCGa                                                                                                            >  4:301014/1‑149 (MQ=255)
                                    aaCGAAGTCGTTGGTGGTGTGATTCCTAAGGAATACATCCCTGCGATCCAGAAGGGTATCGAAGAGCAGATGAAAAACGGCGTTGTTGC‑‑‑CGGTTATCCTCTGATA‑GGCCTGAAGGCCGCGGTATTCGACGGTTCGTACCACGACGTCGa                                                                                                            >  1:241619/1‑149 (MQ=255)
                                    aaCGAAGTCGTTGGTGGTGTGATTCCTAAGGAATACATCCCTGCGATCCAGAAGGGTATCGAAGAGCAGATGAAAAACGGCGTTGTTGC‑‑‑CGGTTATCCTATGATC‑GGCCTGAAGGCCGCGGTATTCGACGGTTCGTACCACGACGTCGa                                                                                                            >  1:245963/1‑149 (MQ=255)
                                                  tggtgTGATTCCTAAGGAATACATCCCTGCGATCCAGAAGGGTATCGAAGAGCAGATGAAAAACGGCGTTGTTGC‑‑‑CGGTTATCCTCTGATC‑GGCCTGAAGGCCGCGGTATTCGACGGTTCGTACCACGACGTCGACTCCAACGAAATgg                                                                                              >  4:8151/1‑149 (MQ=255)
                                                       tgATTCCTAAGGAATACATCCCTGCGATCCAGAAGGGTATCGAAGAGCAGATGAAAAACGGCGTTGTTGCATTAAGTTATCCTCTGATCAGGCCTGAAGGCCGCGGTATTCGACGGTTCGTACCACGACGTCGACTCCAACGAAATGGc                                                                                             >  1:267144/1‑149 (MQ=255)
                                                       tgATTCCTAAGGAATACATCCCTGCGATCCAGAAGGGTATCGAAGAGCAGATGAAAAACGGCGTTGTTGC‑‑‑CGGTTATCCTCTGATC‑GGCCTGAAGGCCGCGGTATTCGACGGTTCGTACCACGACGTCGACTCCAACGAAATGGCGtt                                                                                          <  3:8151/148‑1 (MQ=255)
                                                            ccTAAGCAATACATCCCTGCGATCCAGAAGGGTTTCGAATAGCAGATGAAACACGGCGTTTTTGC‑‑‑GGGTTATCCTCTGAAC‑GGCCTGATGGCCGCGGTATTCGACGGGTCGTACCACGACGTCGACTCCAACGAAATGGCGTGCAAGa                                                                                     <  6:323825/148‑1 (MQ=255)
                                                              tAAGGAATACATCCCTGCGATCCAGAAGGGTATCGAAGAGCAGATGAAAAACGGCGTTGTTGC‑‑‑CGGTTATCCTCTGATC‑GGCCTTCAGGCcgatca                                                                                                                                       <  3:1071/96‑5 (MQ=255)
                                                              tAAGGAATACATCCCTGCGATCCAGAAGGGTATCGAAGAGCAGATGAAAAACGGCGTTGTTGC‑‑‑CGGTTATCCTCTGATC‑GGCCTTCAGGCcgatca                                                                                                                                       >  4:1071/1‑92 (MQ=255)
                                                              tAAGGAATACATCCCTGCGATCCAGAAGGGTATCGAAGAGCAGATGAAAAACGGCGTTGTTGC‑‑‑CGGTTATCCTCTGATC‑GGCCTGAAGGCCGCGGTATTCGACGGTTCGTACCACGACGTCGACTCCAACGAAATGGCGTTCAAGATCg                                                                                  >  3:200686/1‑149 (MQ=255)
                                                                  gAATACATCCCTGCGATCCAGAAGGGTATCGAAGAGCAGATGAAAAACGGCGTTGTTGC‑‑‑CGGTTATCCTCTGATC‑GGCCTGAAGGCCGCGGTATTCGACGGTTc                                                                                                                           >  2:105833/1‑104 (MQ=255)
                                                                  gAATACATCCCTGCGATCCAGAAGGGTATCGAAGAGCAGATGAAAAACGGCGTTGTTGC‑‑‑CGGTTATCCTCTGATC‑GGCCTGAAGGCCGCGGTATTCGACGGTTc                                                                                                                           <  1:105833/104‑1 (MQ=255)
                                                                        aTCCCTGCGATCCAGAAGGGTATCGAAGAGCAGATGAAAAACGGCGTTGTTGC‑‑‑CGGTTATCCTCTGATC‑GGCCTGAAGGCCGCGGTATTCGACGGTTCGTACCACGACGTCGACTCCAACGAAATGGCGTTCAAGATCGCTGCTTCGA‑t                                                                       <  2:15744/149‑1 (MQ=255)
                                                                        aTCCCTGCGATCCAGAAGGGTATCGAAGAGCAGATGAAAAACGGCGTTGTTGC‑‑‑CGGTTATCCTCTGATC‑GGCCTGAAGGCCGCGGTATTCGACGGTTCGTACCACGACGTCGACTCCAACGAAATGGCGTTCAAGATCGCTGCTTCGA‑t                                                                       <  8:130688/149‑1 (MQ=255)
                                                                         tCCCTGCGATCCAGAAGGGTATCGAAGAGCAGATGAAAAACGGCGTTGTTGC‑‑‑CGGTTATCCTCTGATC‑GGCCTGAAGGCCGCGGTATTCGACGGTTCGTACCACGACGTCGACTCCAACGAAATGGCGTTCAAGATCGCTGCTTCGA‑Tg                                                                      <  1:128202/149‑1 (MQ=255)
                                                                         tCCCTGCGATCCAGAAGGGTATCGAAGAGCAGATGAAAAACGGCGTTGTTGC‑‑‑CGGTTATCCTCTGATC‑GGCCTGAAGGCCGCGGTATTCGACGGTTCGTACCACGACGTCGACTCCAACGAAATGGCGTTCAAGATCGCTGCTTCGA‑Tg                                                                      >  7:54870/1‑149 (MQ=255)
                                                                         tCCCTGCGATCCAGAAGGGTATCGAAGAGCAGATGAAAAACGGCGTTGTTGC‑‑‑CGGTTATCCTCTGATC‑GGCCTGAAGGCCGCGGTATTCGACGGTTCGTACCACGACGTCGACTCCAACGAAATGGCGTTCAAGATCGCTGCTTCGA‑Tg                                                                      >  2:139156/1‑149 (MQ=255)
                                                                          cccTGCGATCCAGAAGGGTATCGAAGAGCAGATGAAAAACGGAGTTGTTGC‑‑‑CGCTTATCCTCTGATC‑GGCCTGAAGGCCGCGGTATTCGACGGTTCGTACCACGACGTCGACTCCATCGATATGGCGTTCAAGATCGCTGCTTCGA‑Tgg                                                                     >  5:234760/1‑149 (MQ=255)
                                                                            cTGCGATCCAGAAGGGTATCGAAGAGCAGATGAAAAACGGCGTTGTTGC‑‑‑CGGTTATCCTCTGATC‑GGCCTGAAGGCCGCGGTATTCGACGGTTCGTACCACGACGTCGACTCCAACGAAATGGCGTTCAAGATCGCTGCTTCGA‑TGGc                                                                    >  5:328770/1‑148 (MQ=255)
                                                                            cTGCGATCCAGAAGGGTATCGAAGAGCAGATGAAAAACGGCGTTGTTGC‑‑‑CGGTTATCCTCTGATC‑GGCCTGAAGGCCGCGGTATTCGAAGGTTCGTACCACGACGTCGACTCCAACGAAATGGCGTTCAAGATCGCTGCTTCGA‑TGGc                                                                    >  8:208526/1‑148 (MQ=255)
                                                                               cGATCCAGAAGGGTATCGAAGAGCAGATGAAAACCGGCGTGGTTGC‑‑‑CGGTTATCCTCTGATC‑GGCCTGAAGGCCGCGGTATTCGACGGTTCGTACCACGACGTCGACTCCAACGAAATGGCGTTCAAGATCGCTGCTTCGA‑TGGCGAcc                                                                <  6:78348/149‑1 (MQ=255)
                                                                                  tCCAGAAGGGTATCGAAGAGCAGATGAAAAACGGCGTTGTTGC‑‑‑CGGTTATCCTCTGATC‑GGCCTGAAGGCCGCGGTATTCGACGGTTCGTACCACGACGTCGACTCCAACGAAATGGCGTTCAAGATCGCTGCTTCGA‑TGGCGACCAAg                                                             <  1:139156/149‑1 (MQ=255)
                                                                                      gAAGGGTATCGAAGAGCAGATGAAAAACGGCGTTGTTGC‑‑‑CGGTTATCCTCTGATC‑GGCCTGAAGGCCGCGGTATTCGACGGTTCGTACCACGACGTCGACTCCAACGAAATGGCGTTCAAGATCGCTGCTTCGA‑TGGCGACCAAGCAAc                                                         <  5:229567/149‑1 (MQ=255)
                                                                                           gTATCGAAGAGCAGATGAAAAACGGCGTTGTTGC‑‑‑CGGTTATCCTCTGATC‑GGCCTGAAGGCCGCGGTATTCGACGGTTCGTACCACGACGTCGACTCCAACGAAATGGCGTTCAAGATCGCTGCTTCGA‑TGGCGACCAAGCAACTGGcc                                                    >  8:230527/1‑149 (MQ=255)
                                                                                                gAAGAGCAGATGAAAAACGGCGTTGTTGC‑‑‑CGGTTATCCTCTGATC‑GGCCTGAAGGCCGCGGTATTCGACGGTTCGTACCACGACGTCGACTCCAACGAAATGGCGTTCAAGATCGCTGCTTCGA‑TGGCGACCAAGCAACTGGCCCAGaa                                               <  7:208526/149‑1 (MQ=255)
                                                                                                     gCAGATGAAAAACGGCGTTGTTGC‑‑‑CGGTTATCCTCTGATC‑GGCCTGAAGGCCGCGGTATTCGACGGTTCGTACCACGACGTCGACTCCAACGAAATGGCGTTCAAGATCGCTGCTTCGA‑TGGCGACCAAGCAACTGGCCCAGAAGGGCg                                          >  4:316686/1‑149 (MQ=255)
                                                                                                              aaaCGGCGTTGTTGC‑‑‑CGGTTATCCTCTGATC‑GGCCTGAAGGCCGCGGTATTCGACGGTTCGTACCACGACGTCGACTCCAACGAAATGGCGTTCAAGATCGCTGCTTCGA‑TGGCGACCAAGCAACTGGCCCAGAAGGGCGGTGGCGTgg                                 <  3:316686/149‑1 (MQ=255)
                                                                                                                     gttgttGC‑‑‑CGGTTATCCTCTGATC‑GGCCTGAAGGCCGCGGTATTCGACGGTTCGTACCACGACGTCGACTCCAACGAAATGGCGTTCAAGATCGCTGCTTCGA‑TGGCGACCAAGCAACTGGCCCAGAAGGGCGGTGGCGTGGTTCTCGa                          <  2:241619/149‑1 (MQ=255)
                                                                                                                     gttgttGC‑‑‑CGGTTATCCTCTGATC‑GGCCTGAAGGCCGCGGTATTCGACGGTTCGTACCACGACGTCGACTCCAACGAAATGGCGTTCAAGATCGCTGCTTCGA‑TGGCGACCAAGCAACTGGCCCAGAAGGGCGGTGGCGTGGTTCTCGa                          <  2:245963/149‑1 (MQ=255)
                                                                                                                        gttgC‑‑‑CGGTTATCCTCTGATC‑GGCCTGAAGGCCGCGGTATTCGACGGTTCGTACCACGACGTCGACTCCAACGAAATGGCGTTCAAGATCGCTGCTTCGA‑TGGCGACCAAGCAACTGGCCCAGAAGGGCGGTGGCGTGGTTCTCGAGc                        <  4:149546/148‑1 (MQ=255)
                                                                                                                        gttgC‑‑‑CGGTTATCCTCTGATC‑GGCCTGAAGGCCGCGGTATTCGACGGTTCGTACCACGACGTCGACTCCAACGAAATGGCGTTCAAGATCGCTGCTTCGA‑TGGCGACCAAGCAACTGGCCCAGAAGGGCGGTGGCGTGGTTCTCGAGc                        <  7:230527/148‑1 (MQ=255)
                                                                                                                         ttgC‑‑‑CGGTTATCCTCTGATC‑GGCCTGAAGGCCGCGGTATTCGACGGTTTGTACCACGACGTCGACTCCAACGAAATGGCGTTCAAGATCGCTGCTTCGA‑TGGCGACCAAGCAACTGGCCCAGAAGGGCGGTGGCGTg                                  <  4:287222/137‑1 (MQ=255)
                                                                                                                         ttgC‑‑‑CGGTTATCCTCTGATC‑GGCCTGAAGGCCGCGGTATTCGACGGTTCGTACCACGACGTCGACTCCAACGAAATGGCGTTCAAGATCGCTGCTTCGATTGGCGACCAAGCAACTGGCCCAGAAGGGCGGTGGCGTGGTTCTCGAGc                        <  1:78173/148‑1 (MQ=255)
                                                                                                                         ttgC‑‑‑CGGTTATCCTCTGATC‑GGCCTGAAGGCCGCGGTATTCGACGGTTAGTACCACGACGTCGACTCCAACGAAATGGCGTTCAAGATCGCTGCTTCGA‑TGGCGACCAAGCAACTGGCCCAGAAGGGCGGTGGCGTg                                  >  3:287222/1‑137 (MQ=255)
                                                                                                                           gC‑‑‑CGGTTATCCTCTGATC‑GGCCTGAAGGCCGCGGTATTCGACGGTTCGTACCACGACGTCGACTCCAACGAAATGGCGTTCAAGATCGCTGCTTCGA‑TGGCGACCAAGCAACTGGCCCAGAAGGGCGGTGGCGTGGTTCTCGAGCCGa                     <  3:127614/148‑1 (MQ=255)
                                                                                                                            c‑‑‑cGGTTATCCTCTGATC‑GGCCTGAAGGCCGCGGTATTCGACGGTTCGTACCACGACGTCGACTCCAACGAAATGGCGTTCAAGATCGCTGCTTCGA‑TGGCGACCAAGCAACTGGCCCAGAAGGGCGGTGGCGTGGTTCTCGAGCCGATc                   <  3:301014/149‑1 (MQ=255)
                                                                                                                            c‑‑‑cGGTAATCCTCTGATC‑GGCCAGAACGCCGCGGTATTCGACGGTTCGTACCACGACGTCGACTCCAACGAAATGCCGTTCAAGATCGCTGCTTCGA‑TGGCGACCAAGCACCTGGCCCTGAAGGGCGGTGGCGAGGTGCTCGAGCCgct                    >  8:271710/1‑146 (MQ=255)
                                                                                                                                      tCCTCTGATC‑GGCCTGAAGGCCGCGGTATTCGACGGTTCGTACCACGACGTCGACTCCAACGAAATGGCGTTCAAGATCGCTGCTTCGA‑TGGCGACCAAGCAACTGGCCCAGAAGGGCGGTGGCGTGGTTCTCGAGCCGATCATGaaa             <  7:104736/148‑1 (MQ=255)
                                                                                                                                        ctatGATC‑GGCCTGAAGGCCGCGGTATTCGACGGTTCGTACCACGACGTCGACTCCAACGAAATTGCGTTCAAGATCGCTGCTTCGA‑TGGCGACCAAGCAACTGGCCCAGAAGGGCGGTGGCGTGGTTCTCTAGCCGATCATGAAAGTc          <  6:328770/146‑1 (MQ=255)
                                                                                                                                         tctGATC‑GGCCTGAAGGCCGCGGTATTCGACGGTTCGTACCACGACGTCGACTCCAACGAAATGGCGTTCAAGATCGCTGCTTCGA‑TGGCGACCAAGCAACTGGCCCAGAAGGGCGGTGGCGTGGTTCTCGAGCCGATCATGAAAGTc          >  1:80167/1‑148 (MQ=255)
                                                                                                                                         tctGATC‑GGCCTGAAGGCCGCGGTATTCGACGGTTCGTACCACGACGTCGACTCCAACGAAATGGCGTTCAAGATCGCTGCTTCGA‑TGGCGACCAAGCAACTGGCCCAGAAGGGCGGTGGCGTGGTTCTCGAGCCGATCATGAAAGTc          >  2:26743/1‑148 (MQ=255)
                                                                                                                                            ccgc‑GGCCTTCAGGCCGCGGTATTCGACGGTTCGTACCACGACGTCGACTCCAACGAAATGGCGTTCAAGATCGCTGCTTCGA‑TGGCGCCCAAGCAACTGGCCCAGAAGGGCGGTGGCGTGGTTCTCGAGCAGATACTGAAAGTCGAAg      >  8:209948/4‑149 (MQ=255)
                                                                                                                                             aTC‑GGCCTGACGGCCGCGGTATTCGACGGTTCGTACCACGACGTCGACTCCAACGAAATGGCGTTCAAGATCGCTGCTTCGA‑TGGCGACCAAGCAACTGGCGCAGAAGGGCTGTGGCGTGGTTCTCGAGCCGATCATGAAAGTCgta       <  4:294237/147‑3 (MQ=255)
                                                                                                                                             aTC‑GGCCTGAAGGCCGCGGTATTCGACGGTTCGTACCACGACGTCGACTCCAACGAAATGGCGTTCAAGATCGCTGCTTCGA‑TGGCGACCAAGCAACTGGCCCAGAAGGGCGGTGGCGTGGTTCTCGAGCCGATCATGAAAGTCGAAGt     >  7:11284/1‑149 (MQ=255)
                                                                                                                                             aTC‑GGCCTGAAGGCCGCGGTATTCGACGGTTCGTACCACGACGTCGACTCCAACGAAATGGCGTTCAAGATCGCTGCTTCGA‑TGGCGACCAAGCAACTGGCCCAGAAGGGCGGTGGCGTGGTTCTCGAGCCGATCATGAAAGTCGAAGt     >  3:280693/1‑149 (MQ=255)
                                                                                                                                                  gCCTGAAGGCCGCGGTATTCGACGGTTCGTACCACGACGTCGACTCCAACGAAATGGCGTTCAAGATCGCTGCTTCGA‑TGGCGACCAAGCAACTGGCCCAGAAGGGCGGTGGCGTGGTTCTCGAGCCGATCATGAAAGTCGAAGTTGt  >  7:267840/1‑148 (MQ=255)
                                                                                                                                                     tGAAGGCCGCGGTATTCGACGGTTCGTACCACGACGTCGACTCCAACGAAATGGCGTTCAAGATCGCTGCTTCGA‑TGGCGACCAAGCAACTGGCCCAGAAGGGCGGTGGCGTGGTTCTCGAGc                        <  4:227318/123‑1 (MQ=255)
                                                                                                                                                     tGAAGGCCGCGGTATTCGACGGTTCGTACCACGACGTCGACTCCAACGAAATGGCGTTCAAGATCGCTGCTTCGA‑TGGCGACCAAGCAACTGGCCCAGAAGGGCGGTGGCGTGGTTCTCGAGc                        >  3:227318/1‑123 (MQ=255)
                                                                                                                                                      |                                                                                                                                                
GACAAGGGCAACATCACAGAAGGTCTGGTGTTCTCCAACGAAGTCGTTGGTGGTGTGATTCCTAAGGAATACATCCCTGCGATCCAGAAGGGTATCGAAGAGCAGATGAAAAACGGCGTTGTTGC‑‑‑CGGTTATCCTCTGATC‑GGCCTGAAGGCCGCGGTATTCGACGGTTCGTACCACGACGTCGACTCCAACGAAATGGCGTTCAAGATCGCTGCTTCGA‑TGGCGACCAAGCAACTGGCCCAGAAGGGCGGTGGCGTGGTTCTCGAGCCGATCATGAAAGTCGAAGTTGT  >  NC_002947/548762‑549051

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