Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 708,760 G→A 7.9% V388I (GTC→ATC)  murJ → lipid II flippase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947708,7600GA 7.9% 179.0 / 4.2 63V388I (GTC→ATC) murJlipid II flippase
Reads supporting (aligned to +/- strand):  ref base G (31/27);  new base A (2/3);  total (33/30)
Fisher's exact test for biased strand distribution p-value = 6.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.99e-01

TGGCAAGTTCAGCGCGTTTGATGCGGCGATGACCCAGCGGGCGTTGATCGCCTATTCCGTGGGCTTGCTGGCGATCATTCTGATCAAGGTACTGGCACCTGGCTTCTATGCGCAGCAGAATATTCGCACGCCGGTGAAGATCGCAGTCTTCACCCTGGTCTGCACCCAGCTGTTCAACCTTGCCCTGGTCGGCCCGCTTGCACATGCCGGCCTGGCATTGGCGATCAGCCTCGGTGCCTGCCTGAATGCCGGTTTGCTGTTCTGGAAGCTGCGCAGTCAGCAGCTGTTC  >  NC_002947/708615‑708903
                                                                                                                                                 |                                                                                                                                               
tGGCAAGTTCAGCGCGTTTGATGCGGCGATGACCCAGCGGGCGTTGATCGCCTATTCCGTGGGCTTGCTGGCGATCATTCTGATCAAGGTACTGGCACCTGGCTTCTATGCGCAGCAGAATATTCGCACGCCGGTGAAGATCGCAGTc                                                                                                                                               <  1:104926/148‑1 (MQ=255)
    aaGTTCAGCGCGTTTGATGCGGCGATGACCCAGCGGGCGTTGATCGCCTATTCCGTGGGCTTGCTGGCGATCATTCTGATCAAGGTACTGGCACCTGGCTTCTATGCGCAGCAGAATATTCGCACGCCGGTGAAGATCGCAGTCTTCag                                                                                                                                          >  1:123317/1‑148 (MQ=255)
     aGTTCAGCGCGTTTGATGCGGCGATGACCCAGCGGGCGTTGATCGCCTATTCCGTGGGCTTGCTGGCGATCATTCTGATCAAGGTACTGGCACCTGGCTTCTATGCGCAGCAGAATATTCGCACGCCGGTGAAGACTGCGATCTTCAcc                                                                                                                                         <  6:53276/149‑1 (MQ=255)
      gTTCAGCGCGTTTGATGCGGCGATGACCCAGCGGGCGTTGATCGCCTATTCCGTGGGCTTGCTGGCGATCATTCTGATCAAGGTACTGGCACCTGGCTTCTATGCGCAGCAGAATATTCGCACGCCGGTGAAGACTGCGATCTTCAccg                                                                                                                                        <  1:346142/149‑2 (MQ=255)
          aGCGCGTTTGATGCGGCGATGACCCAGCGGGCGTTGATCGCCTATTCCGTGGGCTTGCTGGCGATCATTCTGATCAAGGTACTGGCACCTGGCTTCTATGCGCAGCAGAATATTCGCACGCCGGTGAAGATCGCAGTCTTCACCCTGGt                                                                                                                                    <  2:123317/149‑1 (MQ=255)
           gcgcGTTTGATGCGGCGATGACCCAGCGGGCGTTGATCGCCTATTCCGTGGGCTTGCTGGCGATCATTCTGATCAAGGTACTGGCACCTGGCTTCTATGCGCAGCAGAATATTCGCACGCCGGTGAAGATCGCAGTCTTCACCCTGGTc                                                                                                                                   >  1:180012/1‑149 (MQ=255)
           gcgcGTTTGATGCGGCGATGACCCAGCGGGCGTTGATCGCCTATTCCGTGGGCTTGCTGGCGATCATTCTGATCAAGGTACTGGCACCTGGCTTCTATGCGCAGCAGAATATTCGCACGCCGGTGAAGATCGCAGTCTTCACCCTGGTc                                                                                                                                   <  6:259891/149‑1 (MQ=255)
            cgcgTTTGATGCGGCGATGACCCAGCGGGCGTTGATCGCCTATTCCGTGGGCTTGCTGGCGATCATTCTGATCAAGGTACTGGCACCTGGCTTCTATGCGCAGCAGAATATTCGCACGCCGGTGAAGATCGCAGTCTTCACCCTGGTCt                                                                                                                                  >  2:172160/1‑149 (MQ=255)
               gTTTGATGCGGCGATGACCCAGCGGGCGTTGATCGCCTATTCCGTGGGCTTGCTGGCGATCATTCTGATCAAGGTACTGGCACCTGGCTTCTATGCGCAGCAGAATATTCGCACGCCGGTGAAGATCGCAGTCTTCACCCTGGTCTGCa                                                                                                                               <  8:250997/149‑1 (MQ=255)
                                cccAGCGGGCGTTGATCGCCTATTCCGTGGGCTTGCTGGCGATCATTCTGATCAAGGTACTGGCACCTGGCTTCTATGCGCAGCAGAATATTCGCACGCCGGTGAAGACTGCGATCTTCAccg                                                                                                                                        >  4:221639/1‑122 (MQ=255)
                                cccAGCGGGCGTTGATCGCCTATTCCGTGGGCTTGCTGGCGATCATTCTGATCAAGGTACTGGCACCTGGCTTCTATGCGCAGCAGAATATTCGCACGCCGGTGAAGACTGCGATCTTCAccg                                                                                                                                        <  3:221639/123‑2 (MQ=255)
                                 ccAGCGGGCGTTGATCGCCTATTCCGTGGGCTTGCTGGCGATCATTCTGATCAAGGTACTGGCACCTGGCTTCTATGCGCAGCAGAATATTCGCACGCCGGTGAAGATCGCAGTCTTCACCCTGGTCTGCACCCAGCTGTTCAACCTTg                                                                                                             >  7:262413/1‑149 (MQ=255)
                                 ccAGCGGGCGTTGATCGCCTATTCCGTGGGCTTGCTGGCGATCATTCTGATCAAGGTACTGGCACCTGGCTTCTATGCGCAGCAGAATATTCGCACGCCGGTGAAGATCGCAGTCTTCACCCTGGTCTGCACCCAGCTGTTCAACCTTg                                                                                                             >  5:248032/1‑149 (MQ=255)
                                 ccAGCGGGCGTTGATCGCCTATTCCGTGGGCTTGCTGGCGATCATTCTGATCAAGGTACTGGCACCTGGCTTCTATGCGCAGCAGAATATTCGCACGCCGGTGAAGATCGCAGTCTTCACCCTGGTCTGCACCCAGCTGTTCAACCTTg                                                                                                             >  3:46912/1‑149 (MQ=255)
                                    gCGGGCGTTGATCGCCTATTCCGTGGGCTTGCTGGCGATCATTCTGATCAAGGTACTGGCACCTGGCTTCTATGCGCAGCAGAATATTCGCACGCCGGTGAAGATCGCAGTCTTCACCCTGGTCTGCACCCAGCTGTTCAACCTTGccc                                                                                                          <  6:248032/149‑1 (MQ=255)
                                             gATCGCCTATTCCGTGGGCTTGCTGGCGATCATTCTGATCAAGGTACTGGCACCTGTCTTCTATGCGCAGCAGAATATTCGCACGCCCGTGAAGATCGCAGTCTACACACTGGTCTGCACCCTGCTGTTCAACATTGCCCTGGTCGGcc                                                                                                 >  8:324462/1‑149 (MQ=255)
                                              aTCGCCTATTCCGTGGGCTTGCTGGCGATCATTCTGATCAAGGTACTGGCACCTGGCTTCTATGCGCAGCAGAATATTCGCACGCCGGTGAAGATCGCAGTCTTCACCCTGGTCTGCACCCAGCTGTTCAACCTTGCCCTGGTCGGccc                                                                                                <  6:338560/149‑1 (MQ=255)
                                               tCGCCTATTCCGTGGGCTTGCTGGCGATCATTCTGATCAAGGTACTGGCACCTGGCTTCTATGCGCAGCAGAATATTCGCACGCCGGTGAAGATCGCAGTCTTCACCCTGGTCTGCACCCAGc                                                                                                                         >  3:27641/1‑123 (MQ=255)
                                               tCGCCTATTCCGTGGGCTTGCTGGCGATCATTCTGATCAAGGTACTGGCACCTGGCTTCTATGCGCAGCAGAATATTCGCACGCCGGTGAAGATCGCAGTCTTCACCCTGGTCTGCACCCAGc                                                                                                                         <  4:27641/123‑1 (MQ=255)
                                                  ccTATTCCGTGGGCTTGCTGGCGATCATTCTGATCAAGGTACTGGCACCTGGCTTCTATGCGCAGCAGAATATTCGCACGCCGGTGAAGATCGCAGTCTTCACCCTGGTCTGCACCCAGCTGTTCAACCTTGCCCTGGTCGGCCCGCtt                                                                                            >  1:106485/1‑149 (MQ=255)
                                                  ccTATTCCGTGGGCTTGCTGGCGATCATTCTGATCAAGGTACTGGCACCTGGCTTCTATGCGCAGCAGAATATTCGCACGCCGGTGAAGATCGCAGTCTTCACCCTGGTCTGCACCCAGCTGTTCAACCTTGCCCTGGTCGGCCCGCtt                                                                                            >  4:156471/1‑149 (MQ=255)
                                                      ttCCGTGGGCTTGCTGGCGATCATTCTGATCAAGGTACTGGCACCTGGCTTCTATGCGCAGCAGAATATTCGCACGCCGGTGAAGATCGCAGTCTTCACCCTGGTCTGCACCCAGCTGTTCAACCTTGCCCTGGTCGGCCCGCTTGc                                                                                          <  8:262413/147‑1 (MQ=255)
                                                               cTTGCTGGCGATCATTCTGATCAAGGTACTGGCACCTGGCTTCTATGCGCAGCAGAATATTCGCACGCCGGTGAAGATCGCAGTCTTCACCCTGGTCTGCACCCAGCTGTTCAACCTTGCCCTGGTCGGCCCGCTTGCACATGCCGGcc                                                                               >  7:105980/1‑149 (MQ=255)
                                                                  gCTGGCGATCATTCTGATCAAGGTACTGGCACCTGGCTTCTATGCGCAGCAGAATATTCGCACGCCGGTGAAGATCGCAGTCTTCACCCTGGTCTGCACCCAGCTGTTCAACCTTGCCCTGGTCGGCCCGCTTGCACATGCCGGCCTgg                                                                            >  7:27745/1‑149 (MQ=255)
                                                                     ggCGATCATTCTGATCAAGGTACTGGCACCTGGCTTCTATGCGCAGCAGAATATTCGCACGCCGGTGAAGATCGCAGTCTTCACCCTGGTCTGCACCCAGCTGTTCAACCTTGCCCTGGTCGGCCCGCTTGCACATGCCGGCCTGGCAt                                                                         <  2:169134/149‑1 (MQ=255)
                                                                        gATCATTCTGATCAAGGTACTGGCACCTGGCTTCTATGCGCAGCAGAATATTCGCACGCCGGTGAAGATCGCAGTCTTCACCCTGGTCTGCACCCAGCTGTTCAACCtt                                                                                                              <  2:194504/109‑1 (MQ=255)
                                                                        gATCATTCTGATCAAGGTACTGGCACCTGGCTTCTATGCGCAGCAGAATATTCGCACGCCGGTGAAGATCGCAGTCTTCACCCTGGTCTGCACCCAGCTGTTCAACCTTGCCCTGGTCGGCCCGCTTGCACATGCCGGCCTGGCATTgg                                                                      <  2:106485/149‑1 (MQ=255)
                                                                        gATCATTCTGATCAAGGTACTGGCACCTCGCTTCTATGCGCAGCAGAATATTCGCACGCCGGTGAAGATCGCAGTCTTCACCCTGGTCTGCACCCAGCTGTTCAACCtt                                                                                                              >  1:194504/1‑109 (MQ=255)
                                                                         aTCATTCTGATCAAGGTACTGGCACCTGGCTTCTATGCGCAGCAGAATATTCGCACGCCGGTGAAGATCGCAGTCTTCACCCTGGTCTGCACCCAGCTGTTCAACCTTGCCCTGGTCGGCCCGCTTGCACATGCCGGCCTGGCATTGGc                                                                     <  8:27745/149‑1 (MQ=255)
                                                                            aTTCTGATCAAGGTACTGGCACCTGGCTTCTATGCGCAGCAGAATATTCGCACGCCGGTGAAGATCGCAGTCTTCACCCTGGTCTGCACCCAGCTGTTCAACCTTGCCCTGGTCGGCCCGCTTGCACATGCCGGCCTGGCATTGGCGAt                                                                  >  5:197059/1‑149 (MQ=255)
                                                                             ttCTGATCAAGGTACTGGCACCTGGCTTCTATGCGCAGCAGAATATTCGCACGCCGGTGAAGATCGCAGTCTTCACCCTGGTCTGCACCCAGCTGTTCAACCTTGCCCTGGTCGGCCCGCTTGCACATGCCAGGCCGGCAtgt                                                                       <  3:211865/143‑3 (MQ=255)
                                                                             ttCTGATCAAGGTACTGGCACCTGGCTTCTATGCGCAGCAGAATATTCGCACGCCGGTGAAGATCGCAGTCTTCACCCTGGTCTGCACCCAGCTGTTCAACCTTGCCCTGGTCGGCCCGCTTGCACATGCCAGGCCGGCAtgt                                                                       >  4:211865/1‑141 (MQ=255)
                                                                                            tGGCACCTGGCTTCTATGCGCAGCAGAATTTTTGCACGCCGGTGAAGATCGCAGTCTTCACCCTGGTCTGCACCCAGCTGTTCAACCTTGCCCTGGTCGGCCCGCTTGCACATGCCGGCCTGGCATTGGCGATCAGCCTCGGtgcctgc                                                  >  3:153559/1‑149 (MQ=255)
                                                                                            tGGCACCTGGCTTCTATGCGCAGCAGAATATTCGCACGCCGGTGAAGATCGCAGTCTTCACCCTGGTCTGCACCCAGCTGTTCAACCTTGCCCTg                                                                                                        <  7:134464/95‑1 (MQ=255)
                                                                                            tGGCACCTGGCTTCTATGCGCAGCAGAATATTCGCACGCCGGTGAAGATCGCAGTCTTCACCCTGGTCTGCACCCAGCTGTTCAACCTTGCCCTg                                                                                                        >  8:134464/1‑95 (MQ=255)
                                                                                            tGGCACCTGGCTTCTATGCGCAGCAGAATATTCGCACGCCGGTGAAGATCGCAGTCTTCACCCTGGTCTGCACCCAGCTGTTCAACCTTGCCCTGGTCGGCCCGCTTGCACATGCCGGCCTGGCATTGGCGATCAGCCTCGGtgcctgc                                                  >  8:294857/1‑149 (MQ=255)
                                                                                             ggCACCTGGCTTCTATGTGCAGCAGAATATTCGCACGCCGGTGAAGATCGCAGTCTTCACCCTGGTCTGCACCCAGCTGTTCAACCTTGCCCTGGTCGGCCCGCTTGCACATGCCGGCCTGGCATTGGCGATCAGCCTCGGtgcctgcc                                                 >  2:81781/1‑149 (MQ=255)
                                                                                                     gCTTCTATGCGCAGCAGAATATTCGCACGCCGGTGAAGATCGCAGTCTTCACCCTGGTCTGCACCCAGCTGTTCAACCTTGCCCTGGTCGGCCCGCTTGCACATGCCg                                                                                  <  7:264259/108‑1 (MQ=255)
                                                                                                     gCTTCTATGCGCAGCAGAATATTCGCACGCCGGTGAAGATCGCAGTCTTCACCCTGGTCTGCACCCAGCTGTTCAACCTTGCCCTGGTCGGCCCGCTTGCACATGCCg                                                                                  >  8:264259/1‑108 (MQ=255)
                                                                                                     gCTTCTATGCGCAGCAGAATATTCGCACGCCGGTGAAGATCGCAGTCTTCACCCTGGTCTGCACCCAGCTGTTCAACCTTGCCCTGGTCGGCCCGCTTGCACATGCCGGCCTGGCATTGGCGATCAGCCTCGGTGCCTGCCTGAATGcc                                         <  1:172160/149‑1 (MQ=255)
                                                                                                      cTTCTATGCGCAGCAGAATATTCGCACGCCGGTGAAGATCGCAGTCTTCACCCTGGTCTGCACCCAGCTGTTCAACCTTGCCCTGGTCGGCCCGCTTGCACATGCCGGCCTGGCATTGGCGATCAGCCTCGGTGCCTGCCTGAATGCCg                                        <  3:156471/149‑1 (MQ=255)
                                                                                                        tCTATGCGCAGCAGAATATTCGCACGCCGGTGAAGATCGCAGTCTTCACCCTGGTCTGCACCCAGCTGTTCAACCTTGCCCTGGTCGGCCCGCTTGCACATGCCGGCCTGGCATTGGCGATCAGCCTCGGTGCCTGCCTGAATGCCGGt                                      >  7:232842/1‑149 (MQ=255)
                                                                                                         cTATGCGCAGCAGAATATTCGCACGCCGGTGAAGATCGCAGTCTTCACCCTGGTCTGCACCCAGCTGTTCAACCTTGCCCTGGTCGGCCCGCTTGCACATGCCGGCCTGGCATTGGCGATCAGCCTCGGTGCCTGCCTGAAt                                            <  3:46237/142‑1 (MQ=255)
                                                                                                         cTATGCGCAGCAGAATATTCGCACGCCGGTGAAGATCGCAGTCTTCACCCTGGTCTGCACCCAGCTGTTCAACCTTGCCCTGGTCGGCCCGCTTGCACATGCCGGCCTGGCATTGGCGATCAGCCTCGGTGCCTGCCTGAAt                                            >  4:46237/1‑142 (MQ=255)
                                                                                                         cTATGCGCAGCAGAATATTCGCACGCCGGTGAAGATCGCAGTCTTCACCCTGGTCTGCACCCAGCTGTTCAACCTTGCCCTGGTCGGCCCGCTTGCACATGCCGGCCTGGCATTGGCGATCAGCCTCGGTGCCTGCCTGAATGCCGGtt                                     >  7:11955/1‑149 (MQ=255)
                                                                                                               gcagcaGAATATTCGCACGCCGGTGAAGATCGCAGTCTTCACCCTGGTCTGCACCCAGCTGTTCAACCTTGCCCTGGTCGGCCCGCTTGCACATGCCGGCCTGGCATTGGCGATCAGCCTCGGTGCCTGCCTGAATGCCGGTTTGCTGt                               >  6:84191/1‑149 (MQ=255)
                                                                                                                 agcagAATATTCGCACGCCGGTGAAGATCGCAGTCTTCACCCTGGTCTGCACCCAGCTGTTCAACCTTGCCCTGGTCGGCCCGCTTGCACATGCCGGCCTGGCATTGGCGATCAGCCTCGGTGCCTGCCTGAATGCCGGTTTGCTGTTc                             >  5:124070/1‑149 (MQ=255)
                                                                                                                  gcagAATATTCGCACGCCGGTGAAGATCGCAGTCTTCACCCTGGTCTGCACCCAGCTGTTCAACCTTGCCCTGGTCGGCCCGCTTGCACATGCCGGCCTGGCATTGGCGATCAGCCTCGGTGCCTGCCTGAATGCCGGTTTGCTGTTCt                            <  2:180012/149‑1 (MQ=255)
                                                                                                                  gcagAATATTCGCACGCCGGTGAAGATCGCAGTCTTCACCCTGGTCTGCACCCAGCTGTTCAACCTTGCCCTGGTCGGCCCGCTTGCACATGCCGGCCTGGCATTGGCGATCAGCCTCGGTGCCTGCCTGAATGCCGGTTTGCTGTTCt                            <  8:11955/149‑1 (MQ=255)
                                                                                                                   cagAATATTCGCACGCCGGTGAAGATCGCAGTCTTCACCCTGGTCTGCACCCAGCTGTTCAACCTTGCCCCGGTCGGCCCGCTTGCACATGCCGGCCTGGCATTGGCTATGAGCCTCGGTGCCTGCCTGAATGCCGGTTTGCTGTTCTg                           <  4:139311/149‑1 (MQ=255)
                                                                                                                      aaTATTCGCACGCCGGTGAAGATCGCAGTCTTCACCCTGGTCTGCACCCAGCTGTTCATCCTTGCCCTGGTCGGCCCGCTTGCACATGCCGGCCTGGCATTGGCGATCAGCCTCGGTGCCTGCCTGAATGCCGGTTTTCTGTTCTGGaa                        <  7:324462/149‑1 (MQ=255)
                                                                                                                           tcacACGCCGGTGAAGATCGCAGTCTTCACCCTGGTCTGCACCCAGCTGTTCAACCTTGCCCTGGTCGGCCCGCTTGCACATGCCGGCCTGGCATTGGCGATCAGCCTCGGTGCCTGCCTGAATGCCGGTTTGCTGTTCTGGAAGCTgc                   >  3:250714/4‑149 (MQ=255)
                                                                                                                           tCGCACGCCGGTGAAGATCGCAGTCTTCACCCTGGTCTGCACCCAGCTGTTCAACCTTGCCCTGGTCGGCCCGCTTGCACATGCCGGCCTGGCATTGGCGATCAGCCTCGGTGCCTGCCTGAATGCCGGTTTGCTGTTCTGGAAGCTgc                   >  7:288628/1‑149 (MQ=255)
                                                                                                                             gCACGCCGGTGAAGATCGCAGTCTTCACCCTGGTCTGCACCCAGCTGTTCAACCTTGCCCTGGTCGGCCCGCTTGCACATGCCGGCCTGGCATTGGCGATCAGCCTCGGTGCCTGCCTGAATGCCGGTTTGCTGTTCTGGAAGCTgcgc                 <  1:81781/149‑1 (MQ=255)
                                                                                                                             gCACGCCGGTGAAGATCGCAGTCTTCACCCTGGTCTGCACCCAGCTGTTCAACCTTGCCCTGGTCGGCCCGCTTGCACATGCCGGCCTGGCATTGGCGATCAGCCTCGGTGCCTGCCTGAATGCCGGTTTGCTGTTCTGGAAGCTgcg                  >  2:270525/1‑148 (MQ=255)
                                                                                                                             acACGCCGGTGAAGATCGCAGTCTTCACCCTGGTCTGCACCCAGCTGTTCAACCTTGCCCTGGTCGGCCCGCTTGCACATGCCGGCCTGGCATTGGCGATCAGCCTCGGTGCCTGCCTGAATGCCGGTTTGCTGTTCTGGAAGCTgcg                  <  4:250714/147‑1 (MQ=255)
                                                                                                                                cGCCGGTGAAGATCGCAGTCTTCACCCTGGTCTGCACCCAGCTGTTCAACCTTGCCCTGGTCGGCCCGCTTGCACATGCCGGCCTGGCATTGGCGATCAGCCTCGGTGCCTGCCTGAATGCCGGTTTGCTGTTCTGGAAGCTGCGCAGt              <  4:46912/149‑1 (MQ=255)
                                                                                                                                 gCCGGTGAAGATCGCAGTCTTCACCCTGGTCTGCACCCAGCTGTTCAACCTTGCCCTGGTCGGCCCGCTTGCACATGCCGGCCTGGCATTGGCGATCAGCCTCGGTGCCTGCCTGAATGCCGGTTTGCTGTTCTGGAAGCTGCGCAGTc             <  4:153559/149‑1 (MQ=255)
                                                                                                                                    ggTGAAGACTGCGATCTTCACCCTGGTCTGCACCCAGCTGTTCAACCTTGCCCTGGTCGGCCCGCTTGCACATGCCGGCCTGGCATTGGCGATCAGCCTCGGTGCCTGCCTGAATGCCGGTTTGCTGTTCTGGAAGCTGCGCAGTcag           >  3:298927/1‑148 (MQ=255)
                                                                                                                                       gAAGATCGCAGTCTTCACCCTGGTCTGCACCCAGCTGTTCAACCTTGCCCTGGTCGGCCCGCTTGCACATGCCGGCCTGGCATTGGCGATCAGCCTCGGTGCCTGCCTGAATGCCGGTTTGCTGTTCTGGAAGCTGCGCAGTCagcagc       >  4:38196/1‑149 (MQ=255)
                                                                                                                                            tCGCAGTCTTCACCCTGGTCTGCACCCAGCTGTTCAACCTTGCCCTGGTCGGCCCGCTTGCACATGCCGGCCTGGCATTGGCGATCAGCCTCGGTGCCTGCCTGAATGCCGGTTTGCTGTTCTGGAAGCTGCGCAGTCAGCAGCTGATc  >  3:236128/1‑149 (MQ=255)
                                                                                                                                                aGTCTTCACCCTGGTCTGCACCCAGCTGTTCAACCTTGCCCTGGTCGGCCCGCTTGCACATGCCGGCCTGGCATTGGCGATCAGCCTCGGTGCCTGCCTGAATGCCg                                        >  3:229670/1‑107 (MQ=255)
                                                                                                                                                aGTCTTCACCCTGGTCTGCACCCAGCTGTTCAACCTTGCCCTGGTCGGCCCGCTTGCACATGCCGGCCTGGCATTGGCGATCAGCCTCGGTGCCTGCCTGAATGCCg                                        <  4:229670/107‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                               
TGGCAAGTTCAGCGCGTTTGATGCGGCGATGACCCAGCGGGCGTTGATCGCCTATTCCGTGGGCTTGCTGGCGATCATTCTGATCAAGGTACTGGCACCTGGCTTCTATGCGCAGCAGAATATTCGCACGCCGGTGAAGATCGCAGTCTTCACCCTGGTCTGCACCCAGCTGTTCAACCTTGCCCTGGTCGGCCCGCTTGCACATGCCGGCCTGGCATTGGCGATCAGCCTCGGTGCCTGCCTGAATGCCGGTTTGCTGTTCTGGAAGCTGCGCAGTCAGCAGCTGTTC  >  NC_002947/708615‑708903

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