Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,323,128 A→T 7.3% Y42N (TAC→AAC)  mrcB ← transglycosylase/transpeptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,323,1280AT 7.3% 190.2 / 3.4 64Y42N (TAC→AAC) mrcBtransglycosylase/transpeptidase
Reads supporting (aligned to +/- strand):  ref base A (28/31);  new base T (3/2);  total (31/33)
Fisher's exact test for biased strand distribution p-value = 6.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCGAGCGCATCGAGCTCGGTGAGGAAGTCGTTCTTGCTCAACTTCTGGCCGGTGAACAGCTCCAGCGGTCGGGCATACACTTTGGCCGGGATGGTCCAGCGCTTGCCGGAGAACTTTTCCTGGACGACGGCATCGAGGTAAACCGCGAAGCCGGCAACGATCACCAGGCCAACCAGGCTGAGCTTCAAAGCCCAGCCCAGCCAGGCGCGGGAGCGGCCGGTCGGGCGTTTCTGAGGGGTGCGGGGATTTCGGGTTCGGGTCATGGCGGCGGATTA  >  NC_002947/5322989‑5323263
                                                                                                                                           |                                                                                                                                       
ccGAGCGCATCGAGCTCGGTGAGGAAGTCGTTCTTGCTCAACTTCCGGCCGGTGAACAGCTCCAGCGGTCGGGCATACACTTTGGCCGGGATGGTCCAGCGCTTGCCGGAGAACTTTTCCTGGACGACGGCATCGAGGTAAACCGCGaa                                                                                                                                <  2:274466/149‑1 (MQ=255)
      gcATCGAGCTCGGTGAGGAAGTCGTTCTTGCTCAACTTCTGGCCGGTGAACAGCTCCAGCGGTCGGGCATACACTTTGGCCGGGATGGTCCAGCGCTTGCCGGAGAACTTTTCCTGGACGACGGCATCGAGGTAAACCGCGAAGCCGGc                                                                                                                          <  3:88226/149‑1 (MQ=255)
         tCGAGCTCGGTGAGGAAGTCGTTCTTGCTCAACTTCTGGCCGGTGAACAGCTCCAGCGGTCGGGCATACACTTTGGCCGGGATGGTCCAGCGCTTGCCGGAGAACTTTTCCTGGACGACGGCATCGAGGTAAACCGCGAAGCCGGCAAc                                                                                                                       >  2:295230/1‑149 (MQ=255)
            aGCTCGGTGAGGAAGTCGTTCTTGCTCAACTTCTGGCCGGTGAACAGCTCCAGCGGTCGGGCATACTCTTTTGCATGGATGGTCCAGCGCTTGCCGGAGAACTTTTCCTGGACGAGGGCATCGAGGTTAACCGCGAAGCCGGCAACGAt                                                                                                                    <  8:277935/149‑1 (MQ=255)
            aGCTCGGTGAGGAAGTCGTTCTTGCTCAACTTCTGGCCGGTGAACAGCTCCAGCGGTCGGGCATACACTTTGGCCGGGATGGTCCAGCGCTTGCCGGAGAACTTTTCCTGGACGACGGCTTCGCGGTTTACCTCGATGCCgtcgtcc                                                                                                                      <  3:202120/147‑7 (MQ=255)
            aGCTCGGTGAGGAAGTCGTTCTTGCTCAACTTCTGGCCGGTGAACAGCTCCAGCGGTCGGGCATACACTTTGGCCGGGATGGTCCAGCGCTTGCCGGAGAACTTTTCCTGGACGACGGCTTCGCGGTTTACCTCGATGCCgtcgtcc                                                                                                                      >  4:202120/1‑141 (MQ=255)
            aGCTCGGTGAGGAAGTCGTTCTTGCTCAACTTCTGGCCGGTGAACAGCTCCAGCGGTCGGGCATACACTTTGGCCGGGATGGTCCAGCGCTTGCCGGAGAACTTTTCCTGGACGACGGCATCGAGGTAAACCGCGAAGCCGGCAACGAt                                                                                                                    >  1:164507/1‑149 (MQ=255)
                 ggTGAGGAAGTCGTTCTTGCTCAACTTCTGGCCGGTGAACAGCTCCAGCGGTCGGGCATACACTTTGGCCGGGATGGTCCAGCGCTTGCCGGAGAACTTTTCCTGGACGACGGCATCGAGGTAAACCGCGAAGCCGGCAACGATCAc                                                                                                                 <  6:89010/147‑1 (MQ=255)
                       gAAGTCGTTCTTCCTCAACTTCTGGCCGGTGAACAGCTCCAGCGGTCGGGCATACACTTTGGCCGGGATGGTCCAGCGCTTGCCGGAGAACTTTTCCTGGACGACGGCATCGAGGTAAACCGCGAAGCCGGCAACGATCACCAGGCCaa                                                                                                         >  3:179058/1‑149 (MQ=255)
                           tCGTTCTTTCTAAACTTCTGGCAGGTGAACAGCTCCAGCGGTCGGGCATACACGTTGGACGGGATGGTCCAGAGCTTGCCGGAGAACTTTTCCTGGACTACGGCATCGACGTAAACCGCGAAGCCGGCAACTATCACCAGGCCAACCAg                                                                                                     >  6:351369/1‑149 (MQ=255)
                             gTTCTTGCTCAACTTCTGGCCGGTGAACAGCTCCAGCGGTCGGGCATACACTTTGGCCGGGATGGTCCAGCGCTTGCCGGAGAACTTTTCCTGGACGACGGCATCGAGGTAAACCGCGAAGCCGGCAACGATCACCAGGCCAACCAgg                                                                                                    <  5:274325/148‑1 (MQ=255)
                                cTTGCTCAACTTCTGGCCGGTGAACAGCTCCAGCGGTCGGGCATACACTTTGGCCGGGATGGTCCAGCGCTTGCCGGAGAACTTTTCCTGGACGACGGCATCGAGGTAAACCGCGAAGCCGGCAACGATCACCAGGCCAACCAGGCTg                                                                                                 <  2:164507/148‑1 (MQ=255)
                                    cTCAACTTCTGGCCGGTGAACAGCTCCAGCGGTCGGGCATACACTTTGGCCGGGATGGTCCAGCGCTTGCCGGAGAACTTTTCCTGGACGACGGCATCGAGGTAAACCGCGAAGCCGGCAACGATCACCAGGCCAACAAGGCTGAGCt                                                                                             >  1:324306/1‑148 (MQ=255)
                                     tCAACTTCTGGCCGGTGAACAGCTCCAGCGGTCGGGCATACACTTTGGCCGGGATGGTCCAGCGCTTGCCGGAGAACTTTTCCTGGACGACGGCATCGAGGTAAACCGCGAAGCCGGCAACGATCACCAGGCCAACCAGGCTGAGCtt                                                                                            >  3:90783/1‑148 (MQ=255)
                                     tCAACTTCTGGCCGGTGAACAGCTCCAGCGGTCGGGCATACACTTTGGCCGGGATGGTCCAGCGCTTGCCGGAGAACTTTTCCTGGACGACGGCATCGAGGTAAACCGCGAAGCCGGCAACGATCACCAGGCCAACCAGGCTGAGCtt                                                                                            >  2:6643/1‑148 (MQ=255)
                                       aaCTTCTGGCCGGTGAACAGCTGCAGCGGTCGGGCATACACTTTGGCCGGGATGGTCCAGCGCTTGCCGGAGAACTTTTCCTGGACGACGTCATCGAGGTAAAACGCGAAGCCGGCAACGATCACCAGGCCAACCa                                                                                                      <  5:31811/136‑1 (MQ=255)
                                       aaCTTCTGGCCGGTGAACAGCTCCAGCGGTCGGGCATACACTTTGGCCGGGATGGTCCAGCGCTTGCCGGAGAACTTTTCCTGGACGACGGCATCGAGGTAAACCGCGAAGCCGGCAACGATCACCAGGCCAACCAGGCTGAGCTTCa                                                                                          >  4:22962/1‑148 (MQ=255)
                                       aaCTTCTGGCCGGTGAACAGCTCCAGCGGTCGGGCATACACTTTGGCCGGGATGGTCCAGCGCTTGCCGGAGAACTTTTCCTGGACGACGGCATCGAGGTAAACCGCGAAGCCGGCAACGATCACCAGGCCAACCAGGCTGAGCTTCa                                                                                          >  8:288391/1‑148 (MQ=255)
                                       aaCTTCTGGCCGGTGAACAGCTCCAGCGGTCGGGCATACACTATGGCCGGGATGGTCCAGCGCTTGCCGGAGAACTTTTCCTGGACGACGGCATCGAGGTAAACCGAGAAGCCGGCAACGATCACCAGGCCAACCa                                                                                                      >  6:31811/1‑136 (MQ=255)
                                        aCTTCTGGCCGGTGAACAGCTCCAGCGGTCGGGCATACACTTTGGCCGGGATGGTCCAGCGCTTGCCGGAGAACTTTTCCTGGACGACGGCATCGAGGTAAACCGCGAAGCCGGCAACGATCACCAGGCCAACCAGGCTGAGCTTCaaa                                                                                        <  4:90783/149‑1 (MQ=255)
                                         cTTCTGGCCGGTGAACAGCTCCAGCGGTCGGGCATACACTTTGGCCGGGATGGTCCAGCGCTTGCCGGAGAACTTTTCCTGGACGACGGCATCGAGGTAAACCGCGAAGCCGGCAACGATCACCAg                                                                                                              >  4:142982/1‑126 (MQ=255)
                                           tCTGGCCGGTGAACAGCTCCAGCGGTCGGGCATACACTTTGGCCGGGATGGTCCAGCGCTTGCCGGAGAACTTTTCCTGGACGACGGCATCGAGGTAAACCGCGAAGCCGGCAACGATCACCAGGCCAACCAGGCTGAGCTTCAAagcc                                                                                     >  2:281907/1‑149 (MQ=255)
                                               gCCGGTGAACAGCTCCAGCGGTCGGGCATACACTTTGGCCGGGATGGTCCAGCGCTTGCCGGAGAACTTTTCCTGGACGACGGCATCGAGGTAAACCGCGAAGCCGGCAACGATCACCAGGCCAACCAGGCTGAGCTTCaaa                                                                                        <  1:13269/142‑1 (MQ=255)
                                               gCCGGTGAACAGCTCCAGCGGTCGGGCATACACTTTGGCCGGGATGGTCCAGCGCTTGCCGGAGAACTTTTCCTGGACGACGGCATCGAGGTAAACCGCGAAGCCGGCAACGATCACCAGGCCAACCAGGCTGAGCTTCaaa                                                                                        >  2:13269/1‑142 (MQ=255)
                                                ccGGTGAACAGCTCCAGCGGTCGGGCATACACTTTGGCCGGGATGGTCCAGCGCTTGCCGGAGAACTTTTCCTGGACGACGGCATCGAGGTAAACCGCGAAGCCGGCAACGATCACCAGGCCAACCAGGCTGAGCTTCAAagcccagcc                                                                                <  3:22962/149‑1 (MQ=255)
                                                 cGGTGAACAGCTCCAGCGGTCGGGCATACACTTTGGCCGGGATGGTCCAGCGCTTGCCGGAGAACTTTTCCTGGACGACGGCATCGAGGTAAACCGCGAAGCCGGCAACGATCACCAGGCCAACCAGGCTGAGCTTCAAagcccagcc                                                                                >  5:130064/1‑148 (MQ=255)
                                                     gAACAGCTCCAGCGGTCGGGCATACACTTTGGCCGGGATGGTCCAGCGCTTGCCGGAGAACTTTTCCTGGACGACGGCATCGAGGTAAACCGCGAAGCCGGCAACGATCACCAGGCCAACCAGGCTGAGCTTCAAAGCccagcccagc                                                                            <  7:245148/148‑1 (MQ=255)
                                                      aaCAGCTCCAGCGGTCGGGCATACACTTTGGCCGGGATGGTCCAGCGCTTGCCGGAGAACTTTTCCTGGACGACGGCATCGAGGTAAACCGCGAAGCCGGCAACGATCACCAGGCCAACCAGGCTGGGCTTTGAagctcagcctggtt                                                                           <  6:213146/148‑6 (MQ=255)
                                                      aaCAGCTCCAGCGGTCGGGCATACACTTTGGCCGGGATGGTCCAGCGCTTGCCGGAGAACTTTTCCTGGACGACGGCATCGAGGTAAACCGCGAAGCCGGCAACGATCACCAGGCCAACCAGGCTGAGCTTCaaa                                                                                        >  2:233650/1‑135 (MQ=255)
                                                      aaCAGCTCCAGCGGTCGGGCATACACTTTGGCCGGGATGGTCCAGCGCTTGCCGGAGAACTTTTCCTGGACGACGGCATCGAGGTAAACCGCGAAGCCGGCAACGATCACCAGGCCAACCAGGCTGAGCTTCaaa                                                                                        <  1:233650/135‑1 (MQ=255)
                                                           cTCCAGCGGTCGGGCATACACTTTGGCCGGGATGGTCCAGCGCTTGCCGGAGAACTTTTCCTGGACGACGGCATCGAGGTAAACCGCGAAGCCGGCAACGATCACCAGGCCAACCAGGCTGAGCTTCAAAGCCCAGCCCAGCCAGgcgc                                                                     >  4:190365/1‑149 (MQ=255)
                                                              cAGCGGTCGGGCATACACTTTGGCCGGGATGGTCCAGCGCTTGCCGGAGAACTTTTCCTGGACGACGGCATCGAGGTAAACCGCGAAGCCGGCAACGATCACCAGGCCAACCAGGCTGAGCTTCAAAGCCCAGCccagccag                                                                         <  5:54907/142‑1 (MQ=255)
                                                              cAGCGGTCGGGCATACACTTTGGCCGGGATGGTCCAGCGCTTGCCGGAGAACTTTTCCTGGACGACGGCATCGAGGTAAACCGCGAAGCCGGCAACGATCACCAGGCCAACCAGGCTGAGCTTCAAAGCCCAGCccagccag                                                                         >  6:54907/1‑142 (MQ=255)
                                                              cAGCGGTCGGGCATACACTTTGGCCGGGATGGTCCAGCGCTTGCCGGAGAACTTTTCCTGGACGACGGCATCGAGGTAAACCGCGAAGCCGGCAACGATCACCAGGCCAACCAGGCTGAGCTTCAAAGCCCAGCCCAGCCAGGCGCgg                                                                   <  1:6643/148‑1 (MQ=255)
                                                              cAGCGGTCGGGCATACACTTTGGCCGGGATGGTCCAGCGCTTGCCGGAGAACTTTTCCTGGACGACGGCATCGAGGTAAACCGCGAAGCCGGCAACGATCACCAGGCCAACCAGGCTGAGCTTCAAAGCCCAGCCCAGCCAGGCGCgg                                                                   >  4:343363/1‑148 (MQ=255)
                                                              cAGCGGTCGGGCATACACTTTGGCCGGGATGGTCCAGCGCTTGCCGGAGAACTTTTCCTGGACGACGGCATCGAGGTAAACCGCGAAGCCGGCAACGATCACCAGGCCAACCAGGCTGAGCTTCAAAGCCCAGCCCAGCCAGGCGCgg                                                                   >  3:339658/1‑148 (MQ=255)
                                                                gCGGTCGGGCATACACTTTGGCCGGGATAGTCCAGCGCTTGCCGGAGAACTTTTCCTGGACGACGGCATCGAGGTAAACCGCGAAGCCGGCAACGATCACCAGGCCAACCAGGCTGAGCTTCAAAGCCCAGCCCAGCCAGGCGCGGGAg                                                                >  1:196461/1‑149 (MQ=255)
                                                                 cGGTCGGGCATACACTTTGGCCGGGATGGTCCAGCGCTTGCCGGAGAACTTTTCCTGGACGACGGCATCGAGGTAAACCGCGAAGCCGGCAACGATCACCAGGCCAACCAGGCTGAGCTTCAAAGCCCAGCCCAGCCAGGCGCGGGAGc                                                               <  8:215688/149‑1 (MQ=255)
                                                                   gTCGGGCATACACTTTGGCCGGGATGGTCCAGCGCTTGCCGGAGAACTTTTCCTGGACGACGGCATCGAGGTAAACCGCGAAGCCGGCAACGATCACCAGGCCAACCAGGCTGAGCTTCAAAGCCCAGCCCAGCCAGGCGCGGGAGCgg                                                             <  7:288391/149‑1 (MQ=255)
                                                                          aTACACTTTGGCCGGGATGGTCCAGCGCTTGCCGGAGAACTTTTCCTGGACGACGGCATCGAGGTAAACCGCGAAGCCGGCAACGATCACCAGGCCAACCAGGCt                                                                                                  <  8:161723/105‑1 (MQ=255)
                                                                          aTACACTTTGGCCGGGATGGTCCAGCGCTTGCCGGAGAACTTTTCCTGGACGACGGCATCGAGGTAAACCGCGAAGCCGGCAACGATCACCAGGCCAACCAGGCt                                                                                                  >  7:161723/1‑105 (MQ=255)
                                                                              acTTTGGCCGGGATGGTCCAGCGCTTGCCGGAGAACTTTTCCTGGACGACGGCATCGAGGTAAACCGCGAAGCCGGCAACGATCACCAGGCCAACCAGGCTGAGCTTCAAAGCCCAGCCCAGCCAGGCGCGGGAGCGGCCGGTCGGGCg                                                  >  7:126953/1‑149 (MQ=255)
                                                                              acTTTGGCCGGGATGGTCCAGCGCTTGCCGGAGAACTTTTCCTGGACGACGGCATCGAGGTAAACCGCGAAGCCGGCAACGATCACCAGGCCAACCAGGCTGAGCTTCAAAGCCCAGCCCAGCCAGGCGCGGGAGCGGCCGGTCGGGCg                                                  >  2:290493/1‑149 (MQ=255)
                                                                                tttGGCCGGGATGGTCCAGCGCTTGCCGGAGAACTTTTCCTGGACGACGGCATCGAGGTAAACCGCGAAGCCGGCAACGATCACCAGGCCAACCAGGCTGAGCTTCAAAGCCCAGCCCAGCCAGGCGCGGGAGCGGCCGGTCGGGCGtt                                                >  3:301141/1‑149 (MQ=255)
                                                                                   ggCCGGGATGGTCCAGCGCTTGCCGGAGAACTTTTCCTGGACGACGGCATCGAGGTAAACCGCGAAGCCGGCAACGATCACCAGGCCAACCAGGCTGAGCTTCAAAGCCCAGCCCAGCCAGGCGCGGGAGCGGCCGGTCGGGCGTTTc                                              <  6:130064/148‑1 (MQ=255)
                                                                                      cGGGATGGTCCAGCGCTTGCCGGAGAACTTTTCCTGGACGACGGCATCGAGGTAAACCGCGAAGCCGGCAACGATCACCAGGCCAACCAGGCTGAGCTTCAAAGCCCAGCCCAGCCAGGCGCGGGAGCGGCCGGTCGGGCGTTTCTGAg                                          <  2:324306/149‑1 (MQ=255)
                                                                                        ggATGGTCCAGCGCTTGCCGGAGAACTTTTCCTGGACGACGGCATCGAGGTAAACCGCGAAGCCGGCAACGATCACCAGGCCAACCAGGCTGAGCTTCAAAGCCCAGCCCAGCCAGGCGCGGGAGCGGCCGGTCGGGCGTTTCTGAggg                                        <  3:190365/149‑1 (MQ=255)
                                                                                              tCCAGCGCTTGCCGGAGAACTTTTCCTGGACGACGGCATCGAGGTAAACCGCGAAGCCGGCAACGATCACCAGGCCAACCAGGCTGAGCTTCAAAGCCCAGCCCAGCCAGGCGCGGGAGCGGCCGGTCGGGCGTTTCTGAGGGGTGCgg                                  <  8:126953/149‑1 (MQ=255)
                                                                                                cAGCGCTTGCCGGAGAACTTTTCCTGGACGACGGCATCGAGGTAAACCGCGAAGCCGGCAACGATCACCAGGCCAACCAGGCTGAGCTTCAAAGCCCAGCCCAGCCAGGCGCGGGAGCGGCCGGTCGGGCGTTTCTGAGGGGTGCggg                                 <  1:290493/148‑1 (MQ=255)
                                                                                                  gcgcTTGCCGGAGAACTTTTCCTGGACGACGGCATCGAGGTAAACCGCGAAGCCGGCAACGATCACCAGGCCAACCAGGCTGAGCTTCAAAGCCCAGCCCAGCCAGGCGCGGGAGCGGc                                                            <  3:49964/119‑1 (MQ=255)
                                                                                                  gcgcTTGCCGGAGAACTTTTCCTGGACGACGGCATCGAGGTAAACCGCGAAGCCGGCAACGATCACCAGGCCAACCAGGCTGAGCTTCAAAGCCCAGCCCAGCCAGGCGCGGGAGCGGc                                                            >  4:49964/1‑119 (MQ=255)
                                                                                                   cgcTTGCCGGAGAACTTTTCCTGGACGACGGCATCGAGGTAAACCGCGAAGCCGGCAACGATCACCAGGCCAACCAGGCTGAGCTTCAAAGCCCAGCCCAGCCAGGCGCGGGAGCGGCCGGTCGGGCGTTTCTGAGGGGTGCGGGGAtt                             >  4:160828/1‑149 (MQ=255)
                                                                                                     cTTGCCGGAGAACTTTTCCTGGACGACGGCATCGAGGTAAACCGCGAAGCCGGCAACGATCACCAGGCCAACCAGGCTGAGCTTCAAAGCCCAGCCCAGCCAGGCGCGGGAGCGGCCGGTCGGGCGTTTCTGAGGGGTGCGGGGATttc                           <  4:339658/149‑1 (MQ=255)
                                                                                                         ccGGAGAAATTTTCCTGGACGACGGCATCGTGGTAAACCGCTAAGCCGGAAACGTTCATCCGGCCAACCTGGCTGAGCTTCAAAGCCCAGCCCAGCCAGGCGCGGGAGCGGCCGGTCGGGCGTTTCTGAGGGGTGCGGGGATttcggg                        <  1:281907/148‑1 (MQ=255)
                                                                                                           ggAGAACTTTTCCTGGACGACGGCATCGAGGTAAACCGCGAAGCCGGCAACGATCACCAGGCCAACCAGGCTGAGCTTCAAAGCCCAGCCCAGCCAGGCGCGGGAGCGGCCGGTCGGGCGTTTCTGAGGGGTGCGGGGATttcgggttc                     >  6:58835/1‑149 (MQ=255)
                                                                                                             agaACTTTTCCTGGACGACGGCATCGAGGTAAACCGCGAAGCCGGCAACGATCACCAGGCCAACCAGGCTTAGCTTCAAAGCCCAGCCCAGCCAGTCGCTGGATCGGCCTGTCGGTCGTTTCTGAGGGGTGCGGGGATttctggttcgg                   <  5:351369/149‑1 (MQ=255)
                                                                                                             agaACTTTTCCTGGACGACGGCATCGAGGTAAACCGCGAAGCCGGCAACGATCACCAGGCCAACCAGGCTGAGCTTCAAAGCCCAGCCCAGCCAGGCGCGGGAGCGGCCGGTCGGGCGTTTCTGAGGGGTGCGGGGATttcgggttcgg                   <  1:295230/149‑1 (MQ=255)
                                                                                                              gaACTTTTCCTGGACGACGGCATCGAGGTAAACCGCGAAGCCGGCAACGATCACCAGGCCAACCAGGCTGAGCTTCAAAGCCCAGCCCAGCCAGGCGCGGGAGCGGCCGGTCGGGCGTTTCTGAGGGGTGCGGGGATttcgggttcggg                  <  5:58835/149‑1 (MQ=255)
                                                                                                                        tGGACGACGGCATCGAGGTAAACCGCGAAGCCGGCAACGATCACCAGGCCAACCAGGCTGAGCTTCAAAGCCCAGCCCAGCCAGGCGCGGGAGCGGCCGGTCGGGCGTTTCTGAGGGGTGCGGGGATttcgggttcgg                   >  6:62385/1‑138 (MQ=255)
                                                                                                                        tGGACGACGGCATCGAGGTAAACCGCGAAGCCGGCAACGATCACCAGGCCAACCAGGCTGAGCTTCAAAGCCCAGCCCAGCCAGGCGCGGGAGCGGCCGGTCGGGCGTTTCTGAGGGGTGCGGGGATttcgggttcgg                   <  5:62385/138‑1 (MQ=255)
                                                                                                                           acgacgGTATCGTGTTATACCGCGAAGCCGTCAACGATCACCAGGCCAACCAGGCTGAGCTTCAAAGCCCAGGCCAGCCAGTCGCGGTAGCGGCCGGTCGGGCGTTTCTGAGGGGTGCGGGGATTTCGGGTTCGGGTCATGGCGGCGGa     <  4:179058/149‑1 (MQ=255)
                                                                                                                            cgacgGCATCGAGGTAAACCGCGAAGCCGGCAACGATCACCAGGCCAACCAGGCTGAGCTTCAAAGCCCAGCCCAGCCAGGCGCGGGAGCGGCCGGTCGGGCGTTTCTGAGGGGTGCGGGGATTTCGGGTTCGGGTCATGGCGGCGGAt    <  4:301141/149‑1 (MQ=255)
                                                                                                                              ccgGCTTCGCGGTTTACCTCGATGCCGGCAACGATCACCAGGCCAACCAGGCTGAGCTTCAAAGCCCAGCCCAGCCAGGCGCGGGAGCGGCCGGTCGGGCGTTTCTGAGGGGTGCGGGGATTTCGGGTTCGGGTCATGGCGGCGGATta  >  7:317602/2‑149 (MQ=255)
                                                                                                                               cgGCTTCGCGGTTTACCTCGATGCCGGCAACGATCACCAGGCCAACCAGGCTGAGCTTCAAAGCCCAGCCCAGCCAGGCGCGGGAGCGGCCGGTCGGGCGTTTCTGAGGGGTGCGGGGATTTCGGGTTCGGGTCATGGCGGCGGATta  >  7:122186/1‑148 (MQ=255)
                                                                                                                                           |                                                                                                                                       
CCGAGCGCATCGAGCTCGGTGAGGAAGTCGTTCTTGCTCAACTTCTGGCCGGTGAACAGCTCCAGCGGTCGGGCATACACTTTGGCCGGGATGGTCCAGCGCTTGCCGGAGAACTTTTCCTGGACGACGGCATCGAGGTAAACCGCGAAGCCGGCAACGATCACCAGGCCAACCAGGCTGAGCTTCAAAGCCCAGCCCAGCCAGGCGCGGGAGCGGCCGGTCGGGCGTTTCTGAGGGGTGCGGGGATTTCGGGTTCGGGTCATGGCGGCGGATTA  >  NC_002947/5322989‑5323263

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: