Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 548,909 A→C 7.0% K587Q (AAG→CAG)  fusA → elongation factor G 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947548,9090AC 7.0% 216.5 / 6.6 71K587Q (AAG→CAG) fusAelongation factor G 1
Reads supporting (aligned to +/- strand):  ref base A (34/32);  new base C (3/2);  total (37/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.08e-01

GACAAGGGCAACATCACAGAAGGTCTGGTGTTCTCCAACGAAGTCGTTGGTGGTGTGATTCCTAAGGAATACATCCCTGCGATCCAGAAGGGTATCGAAGAGCAGATGAAAAACGGCGTTGTTGC‑‑‑CGGTTATCCTCTGATC‑GGCCTGAAGGCCGCGGTATTCGACGGTTCGTACCACGACGTCGACTCCAACGAAATGGCGTTCAAGATCGCTGCTTCGA‑TGGCGACCAAGCAACTGGCCCAGAAGGGCGGTGGCGTGGTTCTCGAGCCGATCATGAAAGTCGAAGTTGTAACCCC  >  NC_002947/548762‑549057
                                                                                                                                                       |                                                                                                                                                     
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                                                                                                                                               >  5:23093/1‑128 (MQ=255)
                            tgtTCTCCAACGAAGTCGTTGGTGGTGTGATTCCTAAGGAATACATCCCTGCGATCCAGAAGGGTATCGAAGAGCAGATGAAAAACGGCGTTGTTGC‑‑‑CGGTTATCCTCTGATC‑GGCCTGAAGGCcgcg                                                                                                                                               <  6:23093/128‑1 (MQ=255)
                            tgtTCTCCAACGAAGTCGTTGGTGGTGTGATTCCTAAGGAATACATCCCTGCGATCCAGAAGGGTATCGAAGAGCAGATGAAAAACGGCGTTGTTGC‑‑‑CGGTTATCCTCTGATC‑GGCCTGAAGGCCGCGGTATTCGACGGTTCGTACCac                                                                                                                          <  6:138292/149‑1 (MQ=255)
                             gtTCTCCAACGAAGTCGTTGGTGGTGTGATTCCTAAGGAATACATCCCTGCGATCCAGAAGGGTATCGAAGAGCAGATGAAAAACGGCGTTGTTGC‑‑‑CGGTTATCCTCTGATC‑GGCCTGAAGGCCGCGGTATTc                                                                                                                                         >  1:194174/1‑133 (MQ=255)
                             gtTCTCCAACGAAGTCGTTGGTGGTGTGATTCCTAAGGAATACATCCCTGCGATCCAGAAGGGTATCGAAGAGCAGATGAAAAACGGCGTTGTTGC‑‑‑CGGTTATCCTCTGATC‑GGCCTGAAGGCCGCGGTATTc                                                                                                                                         <  2:194174/133‑1 (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                                                                            >  2:139156/1‑149 (MQ=255)
                                                                         tCCCTGCGATCCAGAAGGGTATCGAAGAGCAGATGAAAAACGGCGTTGTTGC‑‑‑CGGTTATCCTCTGATC‑GGCCTGAAGGCCGCGGTATTCGACGGTTCGTACCACGACGTCGACTCCAACGAAATGGCGTTCAAGATCGCTGCTTCGA‑Tg                                                                            >  7:54870/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:245963/149‑1 (MQ=255)
                                                                                                                     gttgttGC‑‑‑CGGTTATCCTCTGATC‑GGCCTGAAGGCCGCGGTATTCGACGGTTCGTACCACGACGTCGACTCCAACGAAATGGCGTTCAAGATCGCTGCTTCGA‑TGGCGACCAAGCAACTGGCCCAGAAGGGCGGTGGCGTGGTTCTCGa                                <  2:241619/149‑1 (MQ=255)
                                                                                                                        gttgC‑‑‑CGGTTATCCTCTGATC‑GGCCTGAAGGCCGCGGTATTCGACGGTTCGTACCACGACGTCGACTCCAACGAAATGGCGTTCAAGATCGCTGCTTCGA‑TGGCGACCAAGCAACTGGCCCAGAAGGGCGGTGGCGTGGTTCTCGAGc                              <  7:230527/148‑1 (MQ=255)
                                                                                                                        gttgC‑‑‑CGGTTATCCTCTGATC‑GGCCTGAAGGCCGCGGTATTCGACGGTTCGTACCACGACGTCGACTCCAACGAAATGGCGTTCAAGATCGCTGCTTCGA‑TGGCGACCAAGCAACTGGCCCAGAAGGGCGGTGGCGTGGTTCTCGAGc                              <  4:149546/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                              >  3:227318/1‑123 (MQ=255)
                                                                                                                                                     tGAAGGCCGCGGTATTCGACGGTTCGTACCACGACGTCGACTCCAACGAAATGGCGTTCAAGATCGCTGCTTCGA‑TGGCGACCAAGCAACTGGCCCAGAAGGGCGGTGGCGTGGTTCTCGAGc                              <  4:227318/123‑1 (MQ=255)
                                                                                                                                                       aaGGCCGCGGTATTCGACGGTTCGTACCACGACGTCGACTCCAACGAAATGGCGTTCAAGATCGCTGCTTCGA‑TGGCGACCAAGCAACTGGccc                                                         <  2:213094/94‑1 (MQ=255)
                                                                                                                                                       aaGGCCGCGGTATTCGACGGTTCGTACCACGACGTCGACTCCAACGAAATGGCGTTCAAGATCGCTGCTTCGA‑TGGCGACCAAGCAACTGGCCCAGAAGGGCGGTGGCGTGGTTCTCGAGCCGATCATGAAAGTCGAAGTTGTAAcccc  >  2:78033/1‑149 (MQ=255)
                                                                                                                                                       |                                                                                                                                                     
GACAAGGGCAACATCACAGAAGGTCTGGTGTTCTCCAACGAAGTCGTTGGTGGTGTGATTCCTAAGGAATACATCCCTGCGATCCAGAAGGGTATCGAAGAGCAGATGAAAAACGGCGTTGTTGC‑‑‑CGGTTATCCTCTGATC‑GGCCTGAAGGCCGCGGTATTCGACGGTTCGTACCACGACGTCGACTCCAACGAAATGGCGTTCAAGATCGCTGCTTCGA‑TGGCGACCAAGCAACTGGCCCAGAAGGGCGGTGGCGTGGTTCTCGAGCCGATCATGAAAGTCGAAGTTGTAACCCC  >  NC_002947/548762‑549057

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: ‑