Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,478,890 G→A 8.1% R777C (CGC→TGC)  PP_3090 ← OmpA domain‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,478,8900GA 8.1% 175.0 / 3.9 62R777C (CGC→TGC) PP_3090OmpA domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base G (28/28);  new base A (3/2);  total (31/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

CACCCGGCGGTTGGGCTCCAGGCACTTGATCAACTGCGTTCGCGGTTGTTGCATGTCGCAGTTGACCAATGGCTTGCGGCTGCCTTCGCCCTGCGCGACCACCAGCTCGGCAGGTACGCCCTTGCCAACCAGATAGGTGCGGATGGTTTGCGCCCGTTGTTTCGAAACCTGCAGGTTGCCCTGGCTGTCACCCAGTTGGTCGGCGTGCCCCGAGATGATGATCTTGCCGATATTCGGGGTGTTCAGCAGCTTGCTGGCAATTGCGCTCAGCTGATGCTGGCCTTCGTTT  >  NC_002947/3478750‑3479038
                                                                                                                                            |                                                                                                                                                     
cACCCGGCGGTTGGGCTCCAGGCACTTGATCAACTGCGTTCGCGGTTGTTGCATGTCGCAGTTGACCAATGGCTTGCGGCTGCCTTCGCCCTGCGCGACCACCAGCTCGGCAGGTACGCCCTTGCCAACCAGATAGGTGCGGATGGttt                                                                                                                                               <  4:153660/149‑1 (MQ=255)
cACCCGGCGGTTGGGCTCCAGGCACTTGATCAACTGCGTTCGCGGTTGTTGCATGTCGCAGTTGACCAATGGCTTGCGGCTGCCTTCGCCCTGCGCGACCACCAGCTCGGCAGGTACGCCCTTGCCAACCAGATAGGTGCGGATGGttt                                                                                                                                               >  7:36525/1‑149 (MQ=255)
        ggTTGGGCTCCAGGCACTTGATCAACTGCGTTCGCGGTTGTTGCATGTCGCAGTTGACCAATGGCTTGCGGCTGCCTTCGCCCTGCGCGACCACCTGCTCGGCAGGTAAGCGCTTGCCTAGAAGATAGGTGAGGATGGTTTGCGCCCgt                                                                                                                                       >  5:309381/1‑149 (MQ=255)
        ggTTGGGCTCCAGGCACTTGATCAACTGCGTTCGCGGTTGTTGCATGTCGCAGTTGACCAATGGCTTGCGGCTGCCTTCGCCCTGCGCGACCACCAGCTCGGCAGGTACGCCCTTGCCAACCAGATAGGTGCGGATGGTTTGCGCCCgt                                                                                                                                       >  5:309385/1‑149 (MQ=255)
            gggCTCCAGGCACTTGATCAACTGCGTTCGCGGTTGTTGCATGTCGCAGTTGACCAATGGCTTGCGGCTGCCTTCGCCCTGCGCGACCACCAGCTCGGCAGGTACGCCCTTGCCAACCAGATAGGTGCGGATGGTTTGCGCCCgttgtt                                                                                                                                   >  2:214705/1‑149 (MQ=255)
              gCTCCAGGCACTTGATCAACTGCGTTCGCGGTTGTTGCATGTCGCAGTTGACCAATGGCTTGCGGCTGCCTTCGCCCTGCGCGACCACCAGCTCGGCAGGTACGCCCTTGCCAACCAGATAGGTGCGGATg                                                                                                                                                   >  4:133057/1‑131 (MQ=255)
              gCTCCAGGCACTTGATCAACTGCGTTCGCGGTTGTTGCATGTCGCAGTTGACCAATGGCTTGCGGCTGCCTTCGCCCTGCGCGACCACCAGCTCGGCAGGTACGCCCTTGCCAACCAGATAGGTGCGGATg                                                                                                                                                   <  3:133057/131‑1 (MQ=255)
                  cAGGCACTTGATCAACTGCGTTCGCGGTTGTTGCATGTCGCAGTTGACCAATGGCTTGCGGCTGCCTTCGCCCTGCGCGACCACCAGCTCGGCAGGTACGCCCTTGCCAACCAGATAGGTGCGGATGGTTTGCGCCCGTTGTTTCGaaa                                                                                                                             >  2:294696/1‑149 (MQ=255)
                  cAGGCACTTGATCAACTGCGTTCGCGGTTGTTGCATGTCGCAGTTGACCAATGGCTTGCGGCTGCCTTCGCCCTGCGCGACCACCAGCACGGCAGGTACGCCCTTGCCAACCAGATAGGTGCGGATCGTTTGAGCCCGTTGTTTCGaaa                                                                                                                             >  7:339296/1‑149 (MQ=255)
                   aGGCACTTGATCAACTGCGTTCGCGGTTGTTGCATGTCGCAGTTGACCAATGGCTTGCGGCTGCCTTCGCCCTGCGCGACCACCAGCTCGGCAGGTACGCCCTTGCCAACCAGATAGGTGCGGATTGTTTGCGCCCGTTGTTTCGAAAc                                                                                                                            <  8:228642/149‑1 (MQ=255)
                   aGGCACTTGATCAACTGCGTTCGCGGTTGTTGCATGTCGCAGTTGACCAATGGCTTGCGGCTGCCTTCGCCCTGCGCGACCACCAGCTCGGCAGGTACGCCCTTGCCAACCAGATAGGTGCGGATGGTTTGCGCCCGTTGTTTCGAAAc                                                                                                                            >  8:230208/1‑149 (MQ=255)
                        cTTGATCAACTGCGTTCGCGGTTGTTGCATGTCGCAGTTGACCAATGGCTTGCGGCTGCCTTCGCCCTGCGCGACCACCAGCTCGGCAGGTACGCCCTTGCCAACCAGATAGGTGCGGATGGTTTGCGCCCGTTGTTTCGAAACCTg                                                                                                                         >  4:247926/1‑147 (MQ=255)
                        cTTGATCAACTGCGTTCGCGGTTGTTGCATGTCGCAGTTGACCAATGGCTTGCGGCTGCCTTCGCCCTGCGCGACCACCAGCTCGGCAGGTACGCCCTTGCCAACCAGATAGGTGCGGATGGTTTGCGCCCGTTGTTTCGAAACCTg                                                                                                                         >  7:328087/1‑147 (MQ=255)
                           gATCAACTGCGTTCGCGGTTGTTGCATGTCGCAGTTGACCAATGGCTTGCGGCTGCCTTCGCCCTGCGCGACCACCAGCTCGGCAGGTACGCCCTTGCCAACCAGATAGGTGCGGATGGTTTGCGCCCGTTGTTTCGAAACCTGCAGGt                                                                                                                    >  7:123320/1‑149 (MQ=255)
                               aaCTGCGTTCGCGGTTGTTGCATGTCGCAGTTGACCAATGGCTTGCGGCTGCCTTCGCCCTGCGCGACCACCAGCTCGGCAGGTACGCCCTTGCCAACCAGATAGGTGCGGATGGTTTGCGCCCGTTGTTTCGAAACCTGCAGGTTGcc                                                                                                                <  8:123320/149‑1 (MQ=255)
                               aaCTGCGTTCGCGGTTGTTGCATGTCGCAGTTGACCAATGGCTTGCGGCTGCCTTCGCCCTGCGCGACCACCAGCTCGGCAGGTACGCCCTTGCCAACCAGATAGGTGCGGATGGTTTGCGCCCGTTGTTTCGAAACCTGCAGGTTGcc                                                                                                                >  2:66622/1‑149 (MQ=255)
                               aaCTGCGTTCGCGGTTGTTGCATGTCGCAGTTGACCAATGGCTTGCGGCTGCCTTCGCCCTGCGAGACCACCAGCTCGGCAGGTACTCCCTTGCCAACCAGATAGCTGCAGATGGTTTGCGCCCGTAGATTCGAAACCTGCAGGTTGcc                                                                                                                >  2:285177/1‑149 (MQ=255)
                                 cTGCGTTCGCGGTTGTTGCATGTCGCAGTTGACCAATGGCTTGCGGCTGCCTTCGCCCTGCGCGACCACCAGCTCGGCAGGTACGCCCTTGCCAACCAGATAGGTGCGGATGGTTTGCGCCCGTTGTTTCGAAACCTGCAGGTTGCCCt                                                                                                              >  6:202576/1‑149 (MQ=255)
                                 cTGCGTTCGCGGTTGTTGCATGTCGCAGTTGACCAATGGCTTGCGGCTGCCTTCGCCCTGCGCGACCACCAGCTCGGCAGGTACGCACTTGACAACCAGATAGGTGAGGATGGTTTGCGCCCGTTGTTTAGAAACCTGCAGGTTGCACt                                                                                                              >  1:155673/1‑149 (MQ=255)
                                    cGTTCGCGGTTGTTGCATGTCGCAGTTGACCAATGGCTTGCGGCTGCCTTCGCCCTGCGCGACCACCAGCTCGGCAGGTACGCCCTTGCCAACCAGATAGGTGCGGATGGTTTGCGCCCGTTGTTTCGAAACCTGCAGGTTGCCCTgg                                                                                                            <  5:202576/148‑1 (MQ=255)
                                     gTTCGCGGTTGTTGCATGTCGCAGTTGACCAATGGCTTGCGGCTGCCTTCGCCCTGCGCGACCACCAGCTCGGCAGGTACGCCCTTGCCAACCAGATAGGTGCGGATGGTTTGCGCCCGTTGTTTCGAAACCTGCAGGTTGCCCTGGCt                                                                                                          <  2:92879/149‑1 (MQ=255)
                                          cgGTTGTTGCATGTCGCAGTTGACCAATGGCTTGCGGCTGCCTTCGCCCTGCGCGACCACCAGCTCGGCAGGTACGCCCTTGCCAACCAGATAGGTGCGGATGGTTTGAGCCCGTTGTTTCGAAACCTGCAGGTTGCCCTGGCTGTCAc                                                                                                     >  6:70343/1‑149 (MQ=255)
                                                ttgCATGTCGCAGTTGACCAATGGCTTGCGGCTGCCTTCGCCCTGCGCGACCACCAGCTCGGCAGGTACGCCCTTGCCAACCAGATAGGTGCGGATGGTTTGCGCCCGTTGTTTCGAAACCTGCAGGTTGCCCTGGCTGTCACCCAGt                                                                                                <  3:247926/148‑1 (MQ=255)
                                                  gCATGTCGCAGTTGACCAATGGCTTGCGGCTGCCTTCGCCCTGCGCGACCACCAGCTCGGCAGGTACGCCCTTGCCAACCAGATAGGTGCGGATGGTTTGCGCCCGTTGTTTCGAAACCTGCAGGTTGCCCTGGCTGTCACCCAGTTgg                                                                                             >  1:325640/1‑149 (MQ=255)
                                                         gCAGTTGACCAATGGCTTGCGGCTGCCTTCGCCCTGCGCGACCACCAGCTCGGCAGGTACGCCCTTGCCAACCAGATAGGTGCGGATGGTTTGCGCCCGTTGTTTCGAAACCTGCAGGTTGCCCTGGCTGTCACCCAGTTGGTCGGCGt                                                                                     >  6:175506/1‑149 (MQ=255)
                                                         gCAGTTGACCAATGGCTTGCGGCTGCCTTCGCCCTGCGCGACCACCAGCTCGGCAGGTACGCCCTTGCCAACCAGATAGGTGCGGATGGTTTGCGCCCGTTGTTTCGAAACCTGCAGGTTACCCTGGCTGTCACCCAGTTGGTCGGCGt                                                                                     <  1:214705/149‑1 (MQ=255)
                                                                  cAATGGCTTGCGGCTGCCTTCGCCCTGCGCGACCACCAGCTCGGCAGGTACGCCCTTGCCAACCAGATAGGTGCTGATGGTTTGCGCCCGTTGTTTCGAAACCTGCAGGTTGCCCTGGCTGTCACCCACTTGGTCGGCGTGCCCCgaga                                                                            <  2:1235/149‑1 (MQ=255)
                                                                  cAATGGCTTGCGGCTGCCTTCGCCCTGCGCGACCACCAGCTCGGCAGGTACGCCCTTGCCAACCAGATAGGTGCGGATGGTTTGCGCCCGTTGTTTCGAAACCTGCAGGTTGCCCTGGCTGTCACCCAGTTGGTCGGCGTGCCCCgaga                                                                            <  1:314262/149‑1 (MQ=255)
                                                                  cAATGGCTTGCGGCTGCCTTCGCCCTGCGCGACCACCAGCTCGGCAGGTACGCCCTTGCCAACCAGATAGGTGCGGATGGTTTGCGCCCGTTGTTTCGAAACCTGCAGGTTGCCCTGGCTGTCACCCAGTTGGTCGGCGTGCCCCgaga                                                                            <  8:328087/149‑1 (MQ=255)
                                                                  cAATGGCTTGCGGCTGCCTTCGCCCTGCGCGACCACCAGCTCGGCAGGTACGCCCTTGCCAACCAGATAGGTGCGGATGGTTTGCGCCCGTTGTTTCGAAACCTGCAGGTTGCCCTGGCTGTCACCCAGTTGGTCGGCGTGCCCCgaga                                                                            <  1:66622/149‑1 (MQ=255)
                                                                   aaTGGCTTGCGGCTGCCTTCGCCCTGCGCGACCACCAGCTCGGCAGGTACGCCCTTGCCAACCAGATAGGTGCGGATGGTTTGCGCCCGTTGTTTCGAAACCTGCAGGTTGCCCTgg                                                                                                            <  8:81083/117‑1 (MQ=255)
                                                                    aTGGCTTGCGGCTGCCTTCGCCCTGCGCGACCACCAGCTCGGCAGGTACGCCCTTGCCAACCAGATAGGTGCGGATGGTTTGCGCCCGTTGTTTCGAAACCTGCAGGTTGCCCTGGCTGTCACCCAGTTGGTCGGc                                                                                        >  4:171076/1‑136 (MQ=255)
                                                                    aTGGCTTGCGGCTGCCTTCGCCCTGCGCGACCACCAGCTCGGCAGGTACGCCCTTGCCAACCAGATAGGTGCGGATGGTTTGCGCCCGTTGTTTCGAAACCTGCAGGTTGCCCTGGCTGTCACCCAGTTGGTCGGc                                                                                        <  3:171076/136‑1 (MQ=255)
                                                                    aTGGCTTGCGGCTGCCTTCGCCCTGCGCGACCACCAGCTCGGCAGGTACGCCCTTGCCAACCAGATAGGTGCGGATGGTTTGCGCCCGTTGTTTCGAAACCTGCAGGTTGCCCTGGCTGTCACCCAGTTGGTCGGCGTGCCCCGAgatg                                                                          <  8:224943/149‑1 (MQ=255)
                                                                       gCTTGCGGCTGCCTTCGCCCTGCGCGACCACCAGCTCGGCAGGTACGCCCTTGCCAACCAGATAGGTGCGGATGGTTTGCGCCCGTTGTTTCGAAACCTGCAGGTTGCCCTGGCTGTCACCCAGTTGGTCGGCGTGCCCCGAgatgat                                                                        <  8:36525/148‑1 (MQ=255)
                                                                            cGGCTGCCTTCGCCCTGCGCGACCACCAGCTCGGCAGGTACGCCCTTGCCAACCAGATAGGTGCGGATGGTTTGCGCCCGTTGTTTCGAAACCTGCAGGTTGCCCTGGCTGTCAc                                                                                                     >  4:90823/1‑115 (MQ=255)
                                                                            cGGCTGCCTTCGCCCTGCGCGACCACCAGCTCGGCAGGTACGCCCTTGCCAACCAGATAGGTGCGGATGGTTTGCGCCCGTTGTTTCGAAACCTGCAGGTTGCCCTGGCTGTCAc                                                                                                     <  3:90823/115‑1 (MQ=255)
                                                                                  ccatctcccTGCGCGACCACCAGCTCGGCAGGTACGCCCTTGCCAACCAGATACGTGCGGATGGTTTGCGCCCGATGTTTCGAAACCTGCAGGTTGCCCTGGCTGTCACCCAGTTGGTCGGCGTGCCCCGATATGATGATCTTGCCg                                                              <  6:309381/141‑1 (MQ=255)
                                                                                  ccTTCGCCCTGCGCGACCACCAGCTCGGCAGGTACGCCCTTGCCAACCAGATAGGTGCGGATGGTTTGCGCCCGTTGTTTCGAAACCTGCAGGTTGCCCTGGCTGTCACCCAGTTGGTCGGCGTGCCCCGAGATGATGATCTTGCCg                                                              <  6:309385/147‑1 (MQ=255)
                                                                                   cTTCGCCCTGCGCGTCCACCAGCTCGGCAGGTACGGCCTTGCCAACCAGATAGGTGCGGATGGTTTGCGCCCGTTGTTTCGAAACCTGGAGGTTGCCCTTGCTGTCACCCAGTTGGTCGGCGTGCCCCGAGATGATGATCTTGCCGata                                                           <  2:326998/149‑1 (MQ=255)
                                                                                   cTTCGCCCTGCGCGACCACCAGCTCGGCAGGTACGCCCTTGCCAACCAGATAGGTGCGGATGGTTTGCGCCCGTTGTTTCGAAACCTGCAGGTTGCCCTGGCTGTCACCCAGTTGGTCGGCGTGCCCCGAGATGATGATCTTGCCGata                                                           >  2:232392/1‑149 (MQ=255)
                                                                                    ttCGCCCTGCGCGACCACCAGCTCGGCAGGTACGCCCTTGCCAACCAGATAGGTGCGGATGGTTTGCGCCCGTTGTTTCGAAACCTGCAGGTTGCCCTGGCTGTCACCCAGTTGGTCGGCGTGCCCCg                                                                               >  8:355700/1‑128 (MQ=255)
                                                                                    ttCGCCCTGCGCGACCACCAGCTCGGCAGGTACGCCCTTGCCAACCAGATAGGTGCGGATGGTTTGCGCCCGTTGTTTCGAAACCTGCAGGTTGCCCTGGCTGTCACCCAGTTGGTCGGCGTGCCCCg                                                                               <  7:355700/128‑1 (MQ=255)
                                                                                    ttCGCCCTGCGCGACCACCAGCTCGGCAGGTACGCCCTTGCCAACCAGATAGGTGCGGATGGTTTGCGCCCGTTGTTTCGAAACCTGCAGGTTGCCCTGGCTGTCACCCAGTTGGTCGGCGTGCCCCGAGATGATGATCTTGCCGatat                                                          >  7:79236/1‑149 (MQ=255)
                                                                                    ttCGCCCTGCGCGACCACCAGCTCGGCAGGTACGCCCTTGCCAACCAGATAGGTGCAGATGGTTTGCGCCCGTTGTTTCGAAACCTGCAGGTTGCCCTGGCTGTCACCCAGTTGGTCGGCGTGCCCCGAGATGATGATCTTGCGGatat                                                          >  3:324787/1‑149 (MQ=255)
                                                                                    ttCGCCCTGCGCGACCACCAGCTCGGCAGGTACGCCCTTGCCAACCAGATAGGTGCAGATGGTTTGCGCCCGTTGTTTCGAAACCTGCAGGTTGCCCTGGCTGTCACCCAGTTGGTCGGCGTGCCCCGAGATGATGATCTTGCCGatat                                                          >  3:319326/1‑149 (MQ=255)
                                                                                        cccTGCGCGACCACCAGCTCGGCAGGTACGCCCTTGCCAACCAGATAGGTTCGGATGGTTTGCGCCCGTTGTTTCGAAACCTGCAGGTTGCCCTGGCTGTCACCCAGTTGGTCGGCGTGCCCCGAGATGATGATCTTGCCGATATTCgg                                                      <  7:230208/149‑1 (MQ=255)
                                                                                        cccTGCGCGACCACCAGCTCGGCAGGTACGCCCTTGCCAACCAGATAGGTGCGGATGGGTTGCGCCCGTTGTTTCGAAACCTGCAGGTTGCCCTGGCGGTCACCCAGTTGGTCGGCGTGCCCCGAGATGATGATCTTGCCGATATTCgg                                                      <  5:118144/149‑1 (MQ=255)
                                                                                              gcgACCACCAGCTCGGCAGGTACGCCCTTGCCAACCAGATAGGTGCAGATGGTTTGCGCCCGTTGTTTCGAAACCTGCAGGTTGCCCTGGCTGTCACCCAGTTGGTCGGCGTGCCCCGAGATGATGATCTTGCCGATATTCGGGGTGtt                                                <  4:324787/149‑1 (MQ=255)
                                                                                              gcgACCACCAGCTCGGCAGGTACGCCCTTGCCAACCAGATAGGTGCAGATGGTTTGCGCCCGTTGTTTCGAAACCTGCAGGTTGCCCTGGCTGTCACCCAGTTGGTCGGCGTGCCCCGAGATGATGATCTTGCCGATATTCGGGGTGtt                                                <  4:319326/149‑1 (MQ=255)
                                                                                                   caccaGCTCGGCAGGTACGCCCTTGCCAACCAGATAGGTGCGGATGGTTTGCGCCCGTTGTTTCGAAACCTGCAGGTTGCCCTGGCTGTCACCCAGTTGGTCGGCGTGCCCCGAGATGATGATCTTGCCGATATTCGGGGTGTTcagca                                           <  2:325640/149‑1 (MQ=255)
                                                                                                    accaGCTCGGCAGGTACGCCCTTGCCAACCAGATAGGTGCGGATGGTTTGCGCCCGTTGTTTCGAAACCTGCAGGTTGCCCTGGCTGTCACCCAGTTGGTCGGCGTGCCCCGAGATGATGATCTTGCCGATATTCGGGGTGTTcagcag                                          >  3:346840/1‑149 (MQ=255)
                                                                                                       aGCTCGGCAGGTACGCCCTTGCCAACCAGATAGGTGCGGATGGTTTGCGCCCGTTGTTTCGAAACCTGCAGGTTGCCCTGGCTGTCACCCAGTTGGTCGGCGTGCCCCGAGATGATGATCTTGCCGATATTCGGGGTGTTCAGCAGCtt                                       <  5:175506/149‑1 (MQ=255)
                                                                                                         cTCGGCAGGTACGCCCTTGCCAACCAGATAGGTGCGGATGGTTTGCGCCCGTTGTTTCGAAACCTGCAGGTTGCCCTGGCTGTCACCCAGTTGGTCGGCGTGCCCCGAGATGATGATCTTGCCGATATTCGGGGTGTTCAGCAGCTTGc                                     <  8:79236/149‑1 (MQ=255)
                                                                                                          tCGGCAGGTACGCCCTTGCCAACCAGATAGGTGCGGATGGTTTGCGCCCGTTGTTTCGAAACCTGCAGGTTGCCCTGGCTGTCACCCAGTTGGTCGGCGTGCCCCGAGATGATGATCTTGCCGATATTCGGGGTGTTCAGCAGCTTGCt                                    >  2:38648/1‑149 (MQ=255)
                                                                                                            ggCAGGTAAGCCCTTGCCAACCAGATAGGTGCGGATGGTTTGCGCCCGTTGTTTCGAAACCTGCAGGTTGCCCTGGCTGTCACCCAGTTGGTCGGCGTGCCCCGAGATGATGATCTTGCCGATATTCGGGGTGTTCAGCAGCTTGCTgg                                  <  2:155673/149‑1 (MQ=255)
                                                                                                             gCATGTACGCCCTTGCCTACCAGATAGGTGCGGATGGTTTGGGGCCGTTGTTTCGAAACCTGCAGGTTGCCCTGGCTGTCACCCAGTTGGTCGGCGGTGCCCCGAGATGATGATCTTGCCGATATTCGGGGTGTTCAGCAGCTTGCTgg                                  <  1:285177/149‑1 (MQ=255)
                                                                                                              cAGGTACGCCCTTGCCAACCAGATAGGTGCGGATGGTTTGCGCCCGTTGTTTCGAAACCTGCAGGTTGCCCTGGCTGTCACCCAGTTGGTCGGCGTGCCCCGAGATGATGATCTTGCCGATATTCGGGGTGTTCAGCAGCTTGCTGGc                                 >  1:229967/1‑148 (MQ=255)
                                                                                                              cAGGTACGCCCTTGCCAACCAGATAGGTGCGGATGGTTTGCGCCCGTTGTTTCGAAACCTGCAGGTTGCCCTGGCTGTCACCCAGTTGGTCGGCGTGCCCCGAGATGATGATCTTGCCGATATTCGGGGTGTTCAGCAGCTTGCTGGc                                 >  2:186874/1‑148 (MQ=255)
                                                                                                                                   gATAGGTGCGGATGGTTTGCGCCCGTTGTTTCGAAACCTGCAGGTTGCCCTGGCTGTCACCCAGTTGGTCGGCGTGCCCCGAGATGATGATCTTGCCGATATTCGGGGTGTTCAGCAGCTTGCTGGCAATTGCGCTCAGCTGATGCTgg           >  6:212171/1‑149 (MQ=255)
                                                                                                                                           cGGCTGGTTTCCGGCCGTTGGTTCGAAACCTGCAGGCTGCGCTGGCTGTCACACATTTGGTCGGCGTGCCCCGAGATGATGATCTTGCCGATTTTCGGGTTGTTCAGGAGCTTGCTGGCAATTGCGCTCAGCTGTTGCTGGCCTTCg     <  4:272227/147‑1 (MQ=255)
                                                                                                                                           cGGATGGTTTGCGCCCGTTGTTTCCAAACCTGCAGGTTGCCCTGGCTGTCACCCAGATGGTCGGCGTGCCCCGAGATGATCATCATGCCGATATTCGGGGTGTTAAGCAGCTTGCTGGCAATTGCGCTCAGCTGATGCTGGCCTTCGt  >  3:272227/1‑148 (MQ=255)
                                                                                                                                            ggATGGTTTGCGCCCGTTGTTTCGAAACCTGCAGGTTGCCCTGGCTGTCACCCAGTTGGTCGGCGTGCCCCGAGATGATGATCTTGCCGATATTCGGGGTGTTCAGCAGCTTGCTGGCAATTGCGCTCAGCTGATGCTGGCCTTCGttt  >  1:231917/1‑149 (MQ=255)
                                                                                                                                            ggATGGTTTGCGCCCGTTGTTTCGAAACCTGCAGGATGCCCTGGCTGTCACCCAGTTGGTCGCCGTGCCCCGAGATGAAGATCTTGCCGATATTCGGGGTGTTCAGCAGCTTGCTGGCAATTGCGGTCAGCTGATGCTGGACTGCGttt  >  6:204141/1‑149 (MQ=255)
                                                                                                                                            |                                                                                                                                                     
CACCCGGCGGTTGGGCTCCAGGCACTTGATCAACTGCGTTCGCGGTTGTTGCATGTCGCAGTTGACCAATGGCTTGCGGCTGCCTTCGCCCTGCGCGACCACCAGCTCGGCAGGTACGCCCTTGCCAACCAGATAGGTGCGGATGGTTTGCGCCCGTTGTTTCGAAACCTGCAGGTTGCCCTGGCTGTCACCCAGTTGGTCGGCGTGCCCCGAGATGATGATCTTGCCGATATTCGGGGTGTTCAGCAGCTTGCTGGCAATTGCGCTCAGCTGATGCTGGCCTTCGTTT  >  NC_002947/3478750‑3479038

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: