Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,077,797 G→A 6.1% G381E (GGG→GAG)  PP_3588 → Bcr/CflA family multidrug resistance transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,077,7970GA 6.1% 192.9 / 2.1 66G381E (GGG→GAG) PP_3588Bcr/CflA family multidrug resistance transporter
Reads supporting (aligned to +/- strand):  ref base G (27/34);  new base A (2/2);  total (29/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.73e-01

CGCGTCAGGCAGGTGCGGCAGCGTCGCTCAATACGACATTGCTGTTCGCCGGTGCCGGTTTGACGAGCAGTGTGGTCGCAGGTTTGGAAAGCGAGCTACCGGCCGGCTTGGGCGTGCTGTTCATGGTGGCCAGCCTGTGCGGAGGGTGGCTGCTTGCATGCACCAAGGATGATGGAATGCAACGCTGCTGAAACCTTGATCAGGCAAAAAAAATACCCGGGTGGCGACAAGCGACCCGGGCATTCTTGAATGAAACGGGTGAATCAGCTTACAAGGCCGTACACCACG  >  NC_002947/4077653‑4077940
                                                                                                                                                |                                                                                                                                               
cgcgTCAGGCAGGTGCGGCAGCGTCGCTCCATACGCCATGGCTGTTCGCCGGGGCCGGATTGACGAGCAGTGTGGTCGCAGGTTTGGAAAGCGAGCTACCGGCCGGCTTGGGCGTGCTGTTCATGGTGGCCAGCCTGTGCGGAgggggg                                                                                                                                             <  8:292120/149‑1 (MQ=255)
cgcgTCAGGCAGGTGCGGCAGCGTCGCTCAATACGACATTGCTGTTCGCCGGTGCCGGTTTGACGAGCAGTGTGGTCGCAGGTTTGGAAAGCGAGCTACCGGCCGGCTTGGGCGTGCTGTTCATGGTGGCCAGCCTGTGCGGAGGGTgg                                                                                                                                             <  8:15934/149‑1 (MQ=255)
 gcgTCAGGCAGGTGCGGCAGCGTCGCTCAATACGACATTGCTGTTCGCCGGTGCCGGTTTGACGAGCAGTGTGGTCGCAGGTTTGGAAAGCGAGCTACCGGCCGGCTTGGGCGTGCTGTTCATGGTGGCCAGCCTGTGCGGAGGGTg                                                                                                                                              >  8:1114/1‑147 (MQ=255)
 gcgTCAGGCAGGTGCGGCAGCGTCGCTCAATACGACATTGCTGTTCGCCGGTGCCGGTTTGACGAGCAGTGTGGTCGCAGGTTTGGAAAGCGAGCTACCGGCCGGCTTGGGCGTGCTGTTCATGGTGGCCAGCCTGTGCGGAGGGTg                                                                                                                                              <  7:1114/147‑1 (MQ=255)
      aggcaggTGCGGCAGCGTCGCTCAATACGACATTGCTGTTCGCCGGTGCCGGTTTGACGAGCAGTGTGGTCGCAGGTTTGGAAAGCGAGCTACCGGCCGGCTTGGGCGTGCTGTTCATGGTGGCCAGCCTGTGCGGAGGGTGGCTGCtt                                                                                                                                       <  3:73813/149‑1 (MQ=255)
      aggcaggTGCGGCAGCGTCGCTCAATACGACATTGCTGTTCGCCGGTGCCGGTTTGACGAGCAGTGTGGTCGCAGGTTTGGAAAGCGAGCTACCGGCCGGCTTGGGCGTGCTGTTCATGGTGGCCAGCCTGTGCGGAGGGTGGCTGCtt                                                                                                                                       <  4:92060/149‑1 (MQ=255)
      aggcaggTGCGGCAGCGTCGCTCAATACGACATTGCTGTTCGCCGGTGCCGGTTTGACGAGCAGTGTGGTCGCAGGTTTGGAAAGCGAGCTACCGGCCGGCTTGGGCGTGCTGTTCATGGTGGCCAGCCTGTGCGGAGGGTGGCTGCtt                                                                                                                                       <  7:129145/149‑1 (MQ=255)
          aggTGCGGCAGCGTCGCTCAATACGACATGGCTGTTCGCCGGTGCCGGTTTGACGAGCAGTGTGGTCGCAGGTTTGGAAAGCGAGCTACCGGCCGGCTTGGGCGTGCTGTTCATGGTGGCCAGCCTGTGCGGAGGGTGGCTGCTtgcat                                                                                                                                   <  4:256841/149‑1 (MQ=255)
             tGCGGCAGCGTCGCTCAATACGACATTGCTGTTCGCCGGTGCCGGTTTGACGAGCAGTGTGGTCGCAGGTTTGGAAAGCGAGCTACCGGCCGGCTTGGGCGTGCTGTTCATGGTGGCCAGCCTGTGCGGAGGGTGGCTGCTtgcagacg                                                                                                                                >  5:274892/1‑145 (MQ=255)
               cGGCAGCGTCGCTCAATACGACATTGCTGTTCGCCGGTGCCGGTTTGACGAGCAGTGTGGTCGCAGGTTTGGAAAGCGAGCTACCGGCCGGCTTGGGCGTGCTGTTCATGGTGGCCAGCCTGTGCGGAGGGTGGCTGCTTGCATGCAc                                                                                                                               <  5:186919/148‑1 (MQ=255)
                 gCAGCGTCGCTCAATACGACATTGCTGTTCGCCGGTGCCGGTTTGACGAGCAGTGTGGTCGCAGGTTTGGAAAGCGAGCTACCGGCCGGCTTGGGCGTGCTGTTCATGGTGGCCAGCCTGTGCGGAGGGTGGCTGCTTGCATGCAcc                                                                                                                              >  1:53664/1‑147 (MQ=255)
                 gCAGCGTCGCTCAATACGACATTGCTGTTCGCCGGTGCCGGTTTGACGAGCAGTGTGGTCGCAGGTTTGGAAAGCGAGCTACCGGCCGGCTTGGGCGTGCTGTTCATGGTGGCCAGCCTGTGCGGAGGGTGGCTGCTTGCATGCACCaa                                                                                                                            <  7:17252/149‑1 (MQ=255)
                      gTCGCTCAATACGACATTGCTGTTCGCCGGTGCCGGTTTGACGAGCAGTGTGGTCGCAGGTTTGGAAAGCGAGCTACCGGCCGGCTTGGGCGTGCTGTTCATGGTGGCCAGCCTGTGCGGAGGGTGGCTGCTTGCATGCACCAAGgat                                                                                                                        >  6:77114/1‑148 (MQ=255)
                       tCGCTCAATACGACATTGCTGTTCGCCGGTGCCGGTTTGACGAGCAGTGTGGTCGCAGGTTTGGAAAGCGAGCTACCGGCCGGCTTGGGCGTGCTGTTCATGGTGGCCAGCCTGTGCGGAGGGTGGCTGCTTGCATGCACCAAGgatga                                                                                                                      >  1:131179/1‑149 (MQ=255)
                       tCGCTCAATACGACATTGCTGTTCGCCGGTGCCGGTTTGACGAGCAGTGTGGTCGCAGGTTTGGAAAGCGAGCTACCGGCCGGCTTGGGCGTGCTGTTCATGGTGGCCAGCCTGTGCGGAGGGTGGCTGCTTGCATGCACCAAGgatga                                                                                                                      <  5:183671/149‑1 (MQ=255)
                       tCGCTCAATACGACATTGCTGTTCGCCGGTGCCGGTTTGACGAGCAGTGTGGTCGCAGGTTTGGAAAGCGAGCTACCGGCCGGCTTGGGCGTGCTGTTCATGGTGGCCAGCCTGTGCGGAGGGTGGCTGCTTGCATGCACCAAGgatga                                                                                                                      >  2:189589/1‑149 (MQ=255)
                         gCTCAATACGACATTGCTGTTCGCCGGTGCCGGTTTGACGAGCAGTGTGGTCGCAGGTTTGGAAAGCGAGCTACCGGCCGGCTTGGGCGTGCTGTTCATGGTGGCCAGCCTGTGCGGAGGGTGGCTGCTTGCATGCACCAAGGatgatg                                                                                                                    <  2:253038/149‑1 (MQ=255)
                         gCTCAATACGACATTGCTGTTCGCCGGTGCCGGTTTGACGAGCAGTGTGGTCGCAGGTTTGGAAAGCGAGCTACCGGCCGGCTTGGGCGTGCTGTTCATGGTGGCCAGCCTGTGCGGAGGGTGGCTGCTTGCATGCACCAAGGatgatg                                                                                                                    <  2:118427/149‑1 (MQ=255)
                         gCTCAATACGACATTGCTGTTCGCCGGTGCCGGTTTGACGAGCAGTGTGGTCGCAGGTTTGGAAAGCGAGCTACCGGCCGGCTTGGGCGTGCTGTTCATGGTGGCCAGCCTGTGCGGAGGGTGGCTGCTTGCATGCACCAAGGatgatg                                                                                                                    <  1:340165/149‑1 (MQ=255)
                         gCTCAATACGACATTGCTGTTCGCCGGTGCCGGTTTGACGAGCAGTGTGGTCGCAGGTTTGGAAAGCGAGCTACCGGCCGGCTTGGGCGTGCTGTTCATGGTGGCCAGCCTGTGCGGAGGGTGGCTGCTTGCATGCACCAAGGatgatg                                                                                                                    <  7:306699/149‑1 (MQ=255)
                         gCTCAATACGACATTGCTGTTCGCCGGTGCCGGTTTGACGAGCAGTGTGGTCGCAGGTTTGGAAAGCGAGCTACCGGCCGGCTTGGGCGTGCTGTTCATGGTGGCCAGCCTGTGCGGAGCGTGGCTGCTTGCATGCACCAAGgatgat                                                                                                                     <  1:189589/148‑1 (MQ=255)
                           tCAATACGACATTGCTGTTCGCCGGTGCCGGTTTGACGAGCAGTGTGGTCGCAGGTTTGGAAAGCGAGCTACCGGCCGGCTTGGGCGTGCTGTTCATGGTGGCCAGCCTGTGCGGAGGGTGGCTGCTTGCATGCACCAAGGATGATgg                                                                                                                   >  5:81608/1‑148 (MQ=255)
                             aaTACGACATTGCTGTTCGCCGGTGCCGGTTTGACGAGCAGTGTGGTCGCAGGTTTGGAAAGCGAGCTACCGGCCGGCTTGGGCGTGCTGTTCATGGTGGCCAGCCTGTGCGGAGGGTGGCTGCTTGCATGCACCAAGGATGATGGAAt                                                                                                                >  7:33761/1‑149 (MQ=255)
                             aaTACGACATTGCTGTTCGCCGGTGCCGGTTTGACGAGCAGTGTGGTCGCAGGTTTGGAAAGCGAGCTACCGGCCGGCTTGGGCGTGCTGTTCATGGTGGCCAGCCTGTGCGGAGGGTGGCTGCTTGCATGCACCAAGGATGATGGAAt                                                                                                                <  4:124996/149‑1 (MQ=255)
                               tACGACATTGCTGTTCGCCGGTGCCGGTTTGACGAGCAGTGTGGTCGCAGGTTTGGAAAGCGAGCTACCGGCCGGCTTGGGCGTGCTGTTCATGGTGGCCAGCCTGTGCGGAGGGTGGCTGCTTGCATGCACCAAGGATGATGGAATGc                                                                                                              <  5:77114/149‑1 (MQ=255)
                                  gACATTGCTGTTCGCCGGTGCCGGTTTGACGAGCAGTGTGGTCGCAGGTTTGGAAAGCGAGCTACCGGCCGGCTTGGGCGTGCTGTTCATGGTGGCCAGCCTGTGCGGAGGGTGGCTGCTTGCATGCACCAAGGATGATGGAATGCAAc                                                                                                           >  6:120672/1‑149 (MQ=255)
                                     aTTGCTGTTCGCCGGTGCCGGTTTGACGAGCAGTGTGGTCGCAGGTTTGGAAAGCGAGCTACCGGCCGGCTTGGGCGTGCTGTTAATGGTGGCCAGCCTGTGCGGAGGGTGGCTGCTTGCATGCACCAAGGATGATGGAATGCAACgct                                                                                                        >  3:257502/1‑149 (MQ=255)
                                            ttCGCCGGTGCCGGTTTGACGAGCAGTGTGGTCGCAGGTTTGGAAAGCGAGCTACCGGCCGGCTTGGGCGTGCTGTTCATGGTGGCCAGCCTGTGCGGAGGGTGGCTGCTTGCATGCACCAAGGATGATGGAATGCAACGCTGCTGaaa                                                                                                 >  1:332689/1‑149 (MQ=255)
                                              cGCCGGTGCCGGTTTGACGAGCAGTGTGGTCGCAGGTTTGGAAAGCGAGCTACCGGCCGGCTTGGGCGTGCTGTTCATGGTGGCCAGCCTGTGCGGAGGGTGGCTGCTTGCATGCACCAAGGATGATGGAATGCAACGCTGCTGAAAc                                                                                                <  2:172290/148‑1 (MQ=255)
                                                    tgccggtTTGACGAGCAGTGTGGTCGCAGGTTTGGAAAGCGAGCTACCGGCCGGCTTGGGCGTGCTGTTCATGGTGGCCAGCCTGTGCGGAGGGTGGCTGCTTGCATGCACCAAGGATGATGGAATGCAACGCTGCTGAAACCTTGAt                                                                                          >  1:259179/1‑148 (MQ=255)
                                                       cggtTTGACGAGCATTGTGGTCGCAGGTTTGGTTAGCGAGCTACCGGCCGGCTTGGGCGTGCTGTTCATTGTGGCCATCCTGTGCGGAGGGTGGCTGCTTGCATGCTCCAAGGATGATGGAATGCAACGCTGCTGAAACCTTGATCAgg                                                                                      <  3:218484/149‑1 (MQ=255)
                                                       cggtTTGACGAGCAGTGTGGTCGCAGGTTTGGAAAGCGAGCTACCGGCCGGCTTGGGCGTGCTGTTCATGGTGGCCAGCCTGTGCGGAGGGTGGCTGCTTGCATGCACCAAGGATGATGGAATGCAACGCTGCTGAAACCTTGATCAgg                                                                                      >  4:143913/1‑149 (MQ=255)
                                                         gtTTGACGAGCAGTGTGGTCGCAGGTTTGGAAAGCGAGCTACCGGCCGGCTTGGGCGTGCTGTTCATGGTGGCCAGCCTGTGCGGAGGGTGGCTGCTTGCATGCACCAAGGATGATGGAATGCAACGCTGCTGAAACCTTGATCAGGc                                                                                     >  7:35746/1‑148 (MQ=255)
                                                            tGACGAGCAGTGTGGTCGCAGGTTTGGAAAGCGAGCTACCGGCCGGCTTGGGCGTGCTGTTCATGGTGGCCAGCCTGTGCGGAGGGTGGCTGCTTGCATGCACCAAGGATGATGGAATGCAACGCTGCTGAAACCTTGATCAGGCaaa                                                                                  <  8:56529/148‑1 (MQ=255)
                                                              aCGAGCAGTGTGGTCGCAGGTTTGGAAAGCGAGCTACCGGCCGGCTTGGGCGTGCTGTTCATGGTGGCCAGCCTGTGCGGAGGGTGGCTGCTTGCATGCACCAAGGATGATGGAATGCAACGCTGCTGAAACCTTGATCAGGCaaaaaa                                                                               <  1:327966/149‑1 (MQ=255)
                                                              aCGAGCAGTGTGGTCGCAGGTTTGGAAAGCGAGCTACCGGCCGGCTTGGGCGTGCTGTTCATGGTGGCCAGCCTGTGCGGAGGGTGGCTGCTTGCATGCACCAAGGATGATGGAATGCAACGCTGCTGAAACCTTGATCAGGCaaaaa                                                                                <  8:33761/148‑1 (MQ=255)
                                                              aCGAGCAGTGTGGTCGCAGGTTTGGAAAGCGAGCTACCGGCCGGCTTGGGCGTGCTGTTCATGGTGGCCAGCCTGTGCGGAGGGTGGCTGCTTGCATGCACCAAGGATGATGGAATGCAACGCTGCTGAAACCTTGATCAGGCaaaaa                                                                                >  6:108019/1‑148 (MQ=255)
                                                                gAGCAGTGTGGTCACAGGTTTGGAAAGCGAGCTACCGGCCGGCTTGGGCGTGCTGTTCATGGTGGCCAGCCTGTGCGGAGGGTGGCTGCTTGCATGCACCAAGGATGATGGAATGCAACGCTGCTGAAACCTTGATCAGGCaaaaaaaa                                                                             <  3:170039/149‑1 (MQ=255)
                                                                     gtgtGGTCGCAGGTTTGGAAAGCGAGCTACCGGCCGGCTTGGGCGTGCTGTTCATGGTGGCCAGCCTGTGCGGAGGGTGGCTGCTTGCATGCACCAAGGATGATGGAATGCAACGCTGCTGAAACCTTGATCAGGCAAAAAAAATAccc                                                                        <  6:81608/149‑1 (MQ=255)
                                                                     gtgtGGTCGCAGGTTTGGAAAGCGAGCTACCGGCCGGCTTGGGCGTGCTGTTCATGGTGGCCAGCCTGTGCGGAGGGGGGCTGCTTGCATGCACCAAGGATGATGGAATGCAACGCTGCTGAAACCTTGATCAGGCAAAAAAAATAccc                                                                        <  2:288331/149‑1 (MQ=255)
                                                                           tCGCAGGTTTGGAAAGCGAGCTACCGGCCGGCTTGGGCGTGCTGTTCATGGTGGCCAGCCTGTGCGGAGGGTGGCTGCTTGCATGCACCAAGGATGATGGAATGCAACGCTGCTGAAACCTTGATCAGGCAAAAAAAATACCCGGGTgg                                                                  >  5:278693/1‑149 (MQ=255)
                                                                            cGCAGGTTTGGAAAGCGAGCTACCGGCCGGCTTGGGCGTGCTGTTCATGGTGGCCAGCCTGTGCGGAGGGTGGCTGCTTGCATGCACCAAGGATGATGGAATGCAACGCTGCTGAAACCTTGATCAGGCAAAAAAAATACCCGGGTGGc                                                                 >  5:44045/1‑149 (MQ=255)
                                                                                      gAAAGCGAGCTACCGGCCGGCTTGGGCGTGCTGTTCATGGTGGCCAGCCTGTGCGGAGGGTGGCTGCTTGCATGCACCAAGGATGATGGAATGCAACGCTGCTGAAACCTTGATCAGGCAAAAAAAATACCCGGGTGGCGACAAGCGAc                                                       <  5:108019/149‑1 (MQ=255)
                                                                                      gAAAGCGAGCTACCGGCCGGCTTGGGCGTGCTGTTCATGGTGGCCAGCCTGTGCGGAGGGTGGCTGCTTGCATGCACCAAGGATGATGGAATGCAACGCTGCTGAAACCTTGATCAGGCAAAAAAAATACCCGGGTGGCGACAAGCGAc                                                       <  2:259179/149‑1 (MQ=255)
                                                                                      gAAAGCGAGCTACCGGCCGGCTTGGGCGTGCTGTTCATGGTGGCCAGCCTGTGCGGAGGGTGGCTGCTTGCATGCACCAAGGATGATGGAATGCAACGCTGCTGAAACCTTGATCAGGCAAAAAAAATACCCGGGTGGCGACAAGCGAc                                                       >  4:247322/1‑149 (MQ=255)
                                                                                             aGCTACCGGCCGGCTTGGGCGTGCTGTTCATGGTGGCCAGCCTGTGCGGAGGGTGGCTGCTTGCATGCACCAAGGATGATGGAATGCAACGCTGCTGAAACCTTGATCAGGCAAAAAAAATACCCGGGTGGCGACAAGCGACCCGGGCa                                                >  2:63745/1‑149 (MQ=255)
                                                                                                   cggccggcTTGGGCGTGCTGTTCATGGTGGCCAGCCTGTGCGGAGGGTGGCTGCTTGCATGCACCAAGGATGATGGAATGCAACGCTGCTGAAACCTTGATCAGGCAAAAAAAATACCCGGGTGGCGACAAGCGACCCGGGCATTCtt                                           >  5:56699/1‑148 (MQ=255)
                                                                                                    ggccggcTTGGGCGTGCTGTTCATGGTGGCCAGCCTGTGCGGAGGGTGGCTGCTTGCATGCACCAAGGATGATGGAATGCAACGCTGCTGAAACCTTGATCAGGCAAAAAAAATACCCGGGTGGCGACAAGCGACCCGGGCATTCTtga                                         <  7:252993/149‑1 (MQ=255)
                                                                                                     ggcgggTTGGGGGTGGTGTTCATGGTGGGCAGGCTGTGGGGTTGGTGGGTGGTTGCATTCAACAATGATGATGGAATGCAAAGGTGGTGAAAACTTTATCAGGCAAAAAAAATACCCGGGTGGCGACAAGCGACCCGGGCATTCTtgaa                                        <  8:35746/147‑1 (MQ=255)
                                                                                                     gccggcTTGGGCGTGCTGTTCATGGTGGCCAGCCTGTGCGGAGGGTGGCTGCTTGCATGCACCAAGGATGATGGAATGCAACGCTGCTGAAACCTTGATCAGGCAAAAAAAATACCCGGGTGGCGACAAGCGACCCGGGCATTCTtgaa                                        <  2:332689/149‑1 (MQ=255)
                                                                                                     gccggcTTGGGCGTGCTCTTCATGGTGGCCACCCTGTGCGGAGGGTGGCTCCTTGCATGCACCAAGGATGATGGAATGCAACGCTGCTGAAACCTTGATCAGGCAAAAAAAATACCCGGGAGGCGACAACCCACCCGGGCATACTtgaa                                        >  2:286271/1‑149 (MQ=255)
                                                                                                      ccggcTTGGGCGTGCTGTTCATGGTGGCCAGCCTGTGCGGAGAGTGGCTGCTTGCATGCACCAAGGATGATGGAATGCAACGCTGCTGAAACCTTGATCAGGCAAAAAAAATACCCGGGTGGCGACAAGCGACCCGGGCATTCTtgaat                                       >  4:112498/1‑149 (MQ=255)
                                                                                                      ccggcTTGGGCGTGCTGTTCATGGTGGCCAGCCTGTGCGGAGAGTGGCTGCTTGCATGCACCAAGGATGATGGAATGCAACGCTGCTGAAACCTTGATCAGGCAAAAAAAATACCCGGGTGGCGACAAGCGACCCGGGCATTCTtgaat                                       >  6:271734/1‑149 (MQ=255)
                                                                                                       cggcTTGGGCGTGCTGTTCATGGTGGCCAGCCTGTGCGGAGAGTGGCTGCTTGCATGCACCAAGGATGATGGAATGCAACGCTGCTGACACCTTGATCAGGCAAAAAAAATAGCGGGGGGGCGACAAGCGACCCGGGGATTCTtgaatg                                      <  5:271734/149‑1 (MQ=255)
                                                                                                       cggcTTGGGCGTGCTGTTCATGGTGGCCAGCCTGTGCGGAGAGTGGCTGCTTGCATGCACCAAGGATGATGGAATGCAACGCTGCTGAAACCTTGATCAGGCAAAAAAAATACCCGGGTGGCGACAAGCGACCCGGGCATTCTtgaatg                                      <  3:112498/149‑1 (MQ=255)
                                                                                                              ggCGTGCTGTTCATGGTGGCCAGCCTGTGCGGAGGGTGGCTGCTTGCATGCACCAAGGATGATGGAATGCAACGCTGCTGAAACCTTGATCAGGCAAAAAAAATACCCGGGTGGCGACAAGCGACCCGGGCATTCTTGAATGAAACggg                               <  6:56699/149‑1 (MQ=255)
                                                                                                              ggCGTGCTGTTCATGGTGGCCAGCCTGTGCGGAGGGTGGCTGCTTGCATGCACCAAGGATGATGGAATGCAACGCTGCTGAAACCTTGATCAGGCAAAAAAAATACCCGGGTGGCGACAAGCGACCCGGGCATTCTTGAATGAAACgg                                <  2:131179/148‑1 (MQ=255)
                                                                                                              ggCGTGCTGTTCATGGTGGCCAGCCTGTGCGGAGGGTGGCTGCTTGCATGCACCAAGGATGATGGAATGCAACGCTGCTGAAACCTTGATCAGGCAAAAAAAATACCCGGGTGGCGACAAGCGACCCGGGCATTCTTGAATGAAACgg                                <  5:120672/148‑1 (MQ=255)
                                                                                                              ggCGTGCTGTTCATGGTGGCCAGCCTGTGCGGAGGGTGGCTGATTGCATGCACCAAGGATGATGGAATGCAACGCTGCTGAAACCTTGATCAGGCAAAAAAAATACCCGGGTGGCGACAAGCGACCCGGGCATTCTTGAATGAAACggg                               >  1:64769/1‑149 (MQ=255)
                                                                                                               gCGTGCTGTTCATGGTGGCCAGCCTGTGCGGAGGGTGGCTGCTTGCATGCACCAAGGATGATGGAATGCAACGCTGCTGAAACCTTGATCAGGCAAAAAAAATACCCGGGTGGCGACAAGCGACCCGGGCATTCTTGAATGAAACGGGt                              <  6:278693/149‑1 (MQ=255)
                                                                                                                          aTGGTGGCCAGCCTGTGCGGAGGGTGGCTGCTTGCATGCACCAAGGATGATGGAATGCAACGCTGCTGAAACCTTGATCAGGCAAAAAAAATACCCGGGTGGCGACAAGCGACCCGGGCATTCTTGAATGAAACGGGTGAATCAGCTTa                   >  1:228192/1‑149 (MQ=255)
                                                                                                                                gCCAGCCTGTGCGGAGGGTGGCTGCTTGCATGCACCAAGGATGATGGAATGCAACGCTGATGAAACCTTGATCAGGCAAAAAAAATACCCGGGAGGCGACAAGCGACCCGGGCATTCTTGAATGAAACGGGTGAATCAGCTTACAAGGc             >  3:14123/1‑149 (MQ=255)
                                                                                                                                  cAGCCTGTGCGGAGGGTGGCTGCTTGCATGCACCAAGGATGATGGAATGCAACGCTGCTGAAACCTTGATCAGGCAAAAAAAATACCCGGGTGGCGACAAGCGACCCGGGCATTCTTGAATGAAACGGGTGAATCAGCTTACAAGGCCg           >  6:28400/1‑149 (MQ=255)
                                                                                                                                        gtgCGGAGGGTGGCTGCTTGCATGCACCAAGGATGATGGAATGCAACGCTGCTGAAACCTTGATCAGGCAAAAAAAATACCCGGGTGGCGACAAGCGACCCGGGCATTCTTGAATGAAACGGGTGAATCAGCTTACAAGGCCGTacac      >  4:209383/1‑148 (MQ=255)
                                                                                                                                           cGGAGGGTGGCTGCTTGCATGCACCAAGGATGATGGAATGCAACGCTGCTGAAACCTTGATCAGGCAAAAAAAATACCCGGGTGGCGACAAGCGACCCGGGCATTCTTGAATGAAACGGGTGAATCAGCTTACAAGGCCGTACACCACg  <  2:228192/149‑1 (MQ=255)
                                                                                                                                             gAGGGTGGCTGCTTGCATGCACCAAGGATGATGGAATGCAACGCTGCTGAAACCTTGATCAGGCAAAAAAAATACCCGGGTGGCGACaa                                                            <  2:196706/89‑1 (MQ=255)
                                                                                                                                             gAGGGTGGCTGCTTGCATGCACCAAGGATGATGGAATGCAACGCTGCTGAAACCTTGATCAGGCAAAAAAAATACCCGGGTGGCGACaa                                                            >  1:196706/1‑89 (MQ=255)
                                                                                                                                                |                                                                                                                                               
CGCGTCAGGCAGGTGCGGCAGCGTCGCTCAATACGACATTGCTGTTCGCCGGTGCCGGTTTGACGAGCAGTGTGGTCGCAGGTTTGGAAAGCGAGCTACCGGCCGGCTTGGGCGTGCTGTTCATGGTGGCCAGCCTGTGCGGAGGGTGGCTGCTTGCATGCACCAAGGATGATGGAATGCAACGCTGCTGAAACCTTGATCAGGCAAAAAAAATACCCGGGTGGCGACAAGCGACCCGGGCATTCTTGAATGAAACGGGTGAATCAGCTTACAAGGCCGTACACCACG  >  NC_002947/4077653‑4077940

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