Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,550,720 T→A 6.7% L305Q (CTG→CAG)  PP_4881 → iron ABC transporter substrate‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,550,7200TA 6.7% 222.1 / 6.8 75L305Q (CTG→CAG) PP_4881iron ABC transporter substrate‑binding protein
Reads supporting (aligned to +/- strand):  ref base T (40/30);  new base A (2/3);  total (42/33)
Fisher's exact test for biased strand distribution p-value = 6.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.34e-01

CAGGCGTGGTCAAGGCCAGCAAACACCCGAAAGAAGCCCAGGCCCTGCTCAACTGGATGGCCAGCGAAGAAGGCCAGCGCGTGATCACCCAGACCACCGCCGAGTACCCACTGCACAAGGGCATGGTTTCTGACCGCGGCCTGAAACCGTTCGACGAGCTGCGCCCGCCGAAGATCTCGCCAGCTGACCTGGGCAATGCCGAGGAAGCCATCGAACTTGAGCGTGAGGTCGGCCTGCTCTGATGACTGCCGCCCTGTCCGAACCGGTGCCGGTACGCTTCGTACCGAGCC  >  NC_002947/5550579‑5550868
                                                                                                                                             |                                                                                                                                                    
cAGGCGTGGTCAAGGCCAGCAAACACCCGAAAGAAGCCCAGGCCCTGCTCAACTGGATGGCCAGCGAAGAAGGCCAGCGCGTGATCACCCAGACCACCGCCGAGTACCCACTGCACAAGGGCATGGTTTCAGGCCGCGGTCAGAAAccat                                                                                                                                              <  1:714104/150‑3 (MQ=255)
           aaGGCCAGCAAACACCCGAAAGAAGCCCAGGCCCTGCTCAACTGGATGGCCAGCGAAGAAGGCCAGCGCGTGATCACCCAGACCACCGCCGAGTACCCACTGCACAAGGGCATGGTTTCTGACCGCGGCCTGAAACCGTTCGACGAGCTg                                                                                                                                   >  1:1370020/1‑150 (MQ=255)
              gCCAGCAAACACCCGAAAGAAGCACAGGCCCTGCTCAACTGGATGGCCAGCGAAGAAGGCCAGCGCGTGATAACCCAGACCACCGCCGAGTACCCACTGCACAAGGGCATGGTTTCTGACCGCGGCCTGAAACCGTTCGACGAGCTgcgc                                                                                                                                >  1:207519/1‑150 (MQ=255)
                cAGCAAACACCCGAAAGAAGCCCAGGCCCTGCTCAACTGGATGGCCAGCGAAGAAGGCCAGCGCGTGATCACCCAGACCACCGCCGAGTACCCACTGCACAAGGGCATGGTTTCTGACCGCGGCCTGAAACCGTTCGACGAGCTGcgccc                                                                                                                              >  2:422458/1‑150 (MQ=255)
                  gCAAACACCCGAAAGAAGCCCAGGCCCTGCTCAACTGGATGGCCAGCGAAGAAGGCCAGCGCGTGATCACCCAGACCACCGCCGAGTACCCACTGCACAAGGGCATGGTTTCTGACCGCGGCCTGAAACCGTTc                                                                                                                                            <  2:321692/134‑1 (MQ=255)
                  gCAAACACCCGAAAGAAGCCCAGGCCCTGCTCAACTGGATGGCCAGCGAAGAAGGCCAGCGCGTGATCACCCAGACCACCGCCGAGTACCCACTGCACAAGGGCATGGTTTCTGACCGCGGCCTGAAACCGTTc                                                                                                                                            >  1:321692/1‑134 (MQ=255)
                     aaCACCCGAAAGAAGCCCAGGCCCTGCTCAACTGGATGGCCAGCGAAGAAGGCCAGCGCGTGATCACCCAGACCACCGCCGAGTACCCACTGCACAAGGGCATGGTTTCTGACCGCGGCCTGAAACCGTTCGACGAGCTGCGCCCGCCGa                                                                                                                         >  1:263165/1‑150 (MQ=255)
                      acacCCGAAAGAAGCCCAGGCCCTGCTCAACTGGATGGCCAGCGAAGAAGGCCAGCGCGTGATCACCCAGACCACCGCCGAGTACCCACTGCACAAGGGCATGGTTTCTGACCGCGGCCTGAAACCGTTCGACGAGCTGCGCCCGCCGaa                                                                                                                        >  1:1276401/1‑150 (MQ=255)
                          ccGAAAGAAGCCCAGGCCCTGCTCAACTGGATGGCCAGCGAAGAAGGCCAGCGCGTGATCACCCAGACCACCGCCGAGTACCCACTGCACAAGGGCATGGTTTCTGACCGCGGCCTGAAACCGTTCGACGAGCTGCGCCCGCCGAAGAtc                                                                                                                    >  2:277619/1‑150 (MQ=255)
                           cGAAAGAAGCCCAGGCCCTGCTCAACTGGATGGCCAGCGAAGAAGGCCAGCGCGTGATCACCCAGACCACCGCCGAGTACCCACTGCACAAGGGCATGGTTTCTGACCGCGGCCTGAAACCGTTCGACGAGCTGCGCCCGCCGAAGAtct                                                                                                                   <  1:1191107/150‑1 (MQ=255)
                           cGAAAGAAGCCCAGGCCCTGCTCAACTGGATGGCCAGCGAAGAAGGCCAGCGCGTGATCACCCAGACCACCGCCGAGTACCCACTGCACAAGGGCATGGTTTCTGACCGCGGCCTGAAACCGTTCGACGAGCTGCGCCCGCCGAAGAtct                                                                                                                   <  1:422458/150‑1 (MQ=255)
                            gAAAGAAGCCCAGGCCCTGCTCAACTGGATGGCCAGCGAAGAAGGCCAGCGCGTGATCACCCAGACCACCGCCGAGTACCCACTGCACAAGGGCATGGTTTCTGACCGCGGCCTGAAACCGTTCGACGAGCTGCGCCCGCCGAAGAtctc                                                                                                                  <  2:1184540/150‑1 (MQ=255)
                            gAAAGAAGCCCAGGCCCTGCTCAACTGGATGGCCAGCGAAGAAGGCCAGCGCGTGATCACCCAGACCACCGCCGAGTACCCACTGCACAAGGGCATGGTTTCTGACCGCGGCCTGAAACCGTTCGACGAGCTGCGCCCGCCGAAGAtctc                                                                                                                  >  2:387656/1‑150 (MQ=255)
                            gAAAGAAGCCCAGGCCCTGCTCAACTGGATGGCCAGCGAAGAAGGCCAGCGCGTGATCACCCAGACCACCGCCGAGTACCCACTGCACAAGGGCATGGTTTCTGACCGCGGCCTGAAACCGTTCGACGAGCTGCGCCCGCCGAAGAtctc                                                                                                                  >  1:285978/1‑150 (MQ=255)
                                gaagCCCAGGCCCTGCTCCGCTGGATGGCCACCGAAGAAGGCCAGCGCGGGATCACCCAGACCACCGCCGAGTACCCACTGCACAAGGGCATGGTTTCTGACCGCGGCCTGAAACCGTTCGACGAGCTGCGCCCGCCGAAGATCTCGCCa                                                                                                              <  2:662810/150‑1 (MQ=255)
                                gaagCCCAGGCCCTGCTCAACTGGATGGCCAGCGAAGAAGGCCAGCGCGTGATCACCCAGACCACCGCCGAGTACCCACTGCACAAGGGCATGGTTTCTGACCGCGGCCTGAAACCGTTCGACGAGCTGCGCCCGCCGAAGATCTCGCCa                                                                                                              >  1:973990/1‑150 (MQ=255)
                                     ccAGGCCCTGCTCAACTGGATGGCCAGCGAAGAAGGCCAGCGCGTGATCACCCAGACCACCGCCGAGTACCCACTGCACAAGGGCATGGTTTCTGACCGCGGCCTGAAACCGTTCGACGAGCTGCTCCCGCCGAAGATCTCGCCAGCTGa                                                                                                         >  1:453948/1‑150 (MQ=255)
                                      cAGGCCCTGCTCAACTGGATGGCCAGCGAAGAAGGCCAGCGCGTGATCACCCAGACCACCGCCGAGTACCCACTGCACAAGGGCATGGTTTCTGACCGCGGCCTGAAACCGTTCGACGAGCTGCGCCCGCCGAAGATCTCGCCAGCTGAc                                                                                                        >  2:885340/1‑150 (MQ=255)
                                        ggCCCTGCTCAACTGGATGGCCAGCGAAGAAGGCCAGCGCGTGATCACCCAGACCACCGCCGAGTACCCACTGCACAAGGGCATGGTTTCTGACCGCGGCCTGAAACCGTTCGACGAGCTGCGCCCGCCGAAGAt                                                                                                                     <  2:1258773/135‑1 (MQ=255)
                                        ggCCCTGCTCAACTGGATGGCCAGCGAAGAAGGCCAGCGCGTGATCACCCAGACCACCGCCGAGTACCCACTGCACAAGGGCATGGTTTCTGACCGCGGCCTGAAACCGTTCGACGAGCTGCGCCCGCCGAAGAt                                                                                                                     >  1:1258773/1‑135 (MQ=255)
                                          cccTGCTCAACTGGATGGCCAGCGAAGAAGGCCAGCGCGTGATCACCCAGACCACCGCCGAGTACCCACTGCACAAGGGCATGGTTTCTGACCGCGGCCTGAAACCGTTCGACGAGCTGCGCCCGCCGAAGATCTCGCCAGCTGACCTgg                                                                                                    <  1:451572/150‑1 (MQ=255)
                                             tGCTCAACTGGATGGCCAGCGAAGAAGGCCAGCGCGTGATCACCCAGACCACCGCCGAGTACCCACTGCACAAGGGCATGGTTTCAGGCCGCGGTCAGAAAccatg                                                                                                                                             >  1:360729/1‑103 (MQ=38)
                                             tGCTCAACTGGATGGCCAGCGAAGAAGGCCAGCGCGTGATCACCCAGACCACCGCCGAGTACCCACTGCACAAGGGCATGGTTTCAGGCCGCGGTCAGAAAccatg                                                                                                                                             <  2:360729/106‑4 (MQ=38)
                                                          ggCCAGCGAAGAAGGCCAGCGCGTGATCACCCAGACCACCGCCGAGTACCCACTGCACAAGGGCATGGTTTCTGACCGCGGCCTGAAACCGTTCGACGAGCTGCGCCCGCCGAAGATCTCGCCAGCTGACCTGGGCAATGCCGAGGAAGc                                                                                    >  1:353545/1‑150 (MQ=255)
                                                             cAGCGAAGAAGGCCAGCGCGTGATCACCCAGACCACCGCCGAGTACCCACTGCACAAGGGCATGGTTTCTGACCGCGGCCTGAAACCGTTCGACGAGCTGCGCCCGCCGAAGATCTCGCCAGCTGACCTGGGCAATGCCGAGGAAGCCAt                                                                                 >  1:272691/1‑150 (MQ=255)
                                                                    gaagGCCAGCGCGTGATCACCCAGACCACCGCCGAGTCCCCACTCCACAAGGGCATGGTTTCTGACCGCGGCCTGAAACCGTTCGACGAGCTGCGCCCGCCGAAGATCTCGCCAGCTGACCTGAGCAATGCCGAGGAAGCCATCGAACtt                                                                          <  2:1028427/150‑1 (MQ=255)
                                                                      agGCCAGCGCGTGATCGCCCAGACCACCGCCGAGTACCCACTGCACAAGGGCATGGTTTCTGACCGCGGCCTGAAACCGTTCGACGAGCTGCGCCCGCCGAAGATCTCGCCAGCTGACCTGGGCAATGCCGAGGAAGCCATCGAACTTGa                                                                        <  2:1370020/150‑1 (MQ=255)
                                                                        gCCAGCGCGTGATCACCCAGACCACCGCCGAGTACCCACTGCACAAGGGCATGGTTTCTGACCGCGGCCTGAAACCGTTCGACGAGCTGCGCCCGCCGAAGATCTCGCCAGCTGACCTGGGCAATGCCGAGGa                                                                                       <  2:191926/133‑1 (MQ=255)
                                                                        gCCAGCGCGTGATCACCCAGACCACCGCCGAGTACCCACTGCACAAGGGCATGGTTTCTGACCGCGGCCTGAAACCGTTCGACGAGCTGCGCCCGCCGAAGATCTCGCCAGCTGACCTGGGCAATGCCGAGGa                                                                                       >  1:191926/1‑133 (MQ=255)
                                                                          cAGCGCGTGATCACCCAGACCACCGCCGAGTACCCACTGCACAAGGGCATGGTTTCTGACCGCGGCCTGAAACCGTTCGACGAGCTGCGCCCGCCGAAGATCTCGCCAGCTGACCTGGGCAATGCCGAGGAAGCCATCGAACTTGAGCGt                                                                    <  2:608865/150‑1 (MQ=255)
                                                                           aGCGCGTGATCACCCAGACCACCGCCGAGTAACCACTGCAAAAGGACATGGTTTATGACCGCGACCTGAAAGCGGTCGACAAGCTTCGCCCTCCGAAGATCTCGCCAGCTGACCTGGGCAATGCCGAGGAAGCCATCGAACTTGAGCGTg                                                                   >  2:1230675/1‑150 (MQ=255)
                                                                             cgcgTGATCACCCAGACCACCGCCGAGTACCCACTGCACAAGGGCATGGTTTCTGACCGCGGCCTGAAACCGTTCGACGAGCTGCGCCCGCCGAAGATCTCGCCAGCTGACCTGGGCAATGCCGAGGAAGCCATCGAACTTGAGCGTGAg                                                                 <  2:1401865/150‑1 (MQ=255)
                                                                               cgTGATCACCCAGACCACCGCCGAGTACCCACTGCACAAGGGCATGGTTTCTGGCAGCGGCTTTAATCCGTTCCACGAGCTCCTTCTGTGGAGGATCTTTCCGGGTTTGCTGGGTAATGCCGAGGAAGCCATCGAACTTGAGCGTGAGGt                                                               <  2:1230167/150‑1 (MQ=255)
                                                                                gTGATCACCCAGACCACCGCCGAGTACCCACTGCACAAGGGCATGGTTTCTGACCGCGGCCTGAAACCGTTCGACGAGCTGCGCCCGCCGAAGATCTCGCCAGCTGACCTGGGCAATGCCGAGGAAGCCATCGAACTTGAGCGTGAGGTc                                                              >  1:1256625/1‑150 (MQ=255)
                                                                                   ataaCCCAGACCACCGCCGAGTACCCACTGCACAAGGGCATGGTTTCTGACCGCGGCCTGAAACCGTTCGACGAGCTGCGCCCGCCGAAGATCTCGCCAGCTGACCTGGGCAATGCCGAGGAAGCCATCGAACTTGAGCGTGAGGTCGGc                                                           <  2:207519/147‑1 (MQ=255)
                                                                                   aTCACCCAGACCACCGCCGAGTACCCACTGCACAAGGGCATGGTTTCTGATCTCGGCCTGAAACCGTTCGACGAGCTGCGCCCGCCGAAGATCTCGCCAGCTAACCTGGGCAATGCCGAGGAAGCCATCGAACTTGAGCGTGAGGTCGGc                                                           <  1:968569/150‑1 (MQ=255)
                                                                                    acACCCAGACCACCGCCGAGTACACACTGGACAACGCCATGGTTTACAAACGCGGCCTGACACCGGTCGAGGTCATGAACCCGCGGAATGGTGCGCCGGCTGACCTTGGCAATGCCGAGGAAGCGATCGACCTTGAGCGTGAGTTCGTcc                                                          >  2:519730/2‑150 (MQ=255)
                                                                                      aCCCAGACCACCGCCGAGTACCCACTGCACAAGGGCATGGTTTCTGACCGCGGCCTGAAACCGTTCGACGAGCTGCGCCCGCCGAAGATCTCGCCAGCTGACCTGGGCAATGCCGAGGAAGCCATCGAACTTGAGCGTGAGGTCGGCCTg                                                        >  1:143479/1‑150 (MQ=255)
                                                                                       cccAAACAAACGCCGAGGACCCACTGCACAAGGGCATGGTTTCTGACCGCGGCCTGAAACCGTTCGACGAGCTGCGCCCGCCGAAGATCTCGCCAGCTGACCTGGGCAATGCCGAGGAAGCCATCGAACTTGAGCGTGAGGTCGGCCTGc                                                       >  2:1348384/1‑150 (MQ=255)
                                                                                        ccAGACCACCGCCGAGTACCCACTGCACAAGGGCATGGTTTCTGACCGCGGCCTGAAACCGTTCGACGAGCTGCGCCCGCCGAAGATCTCGCCAGCTGACCTGGGCAATGCCGAGGAAg                                                                                     <  2:1337630/119‑1 (MQ=255)
                                                                                        ccAGACCACCGCCGAGTACCCACTGCACAAGGGCATGGTTTCTGACCGCGGCCTGAAACCGTTCGACGAGCTGCGCCCGCCGAAGATCTCGCCAGCTGACCTGGGCAATGCCGAGGAAg                                                                                     >  1:1337630/1‑119 (MQ=255)
                                                                                        ccAGACCACCGCCGAGTACCCACTGCACAAGGGCATGGTTTCTGACCGCGGCCTGAAACCGTTCGACGAGCTGCGCCCGCCGAAGATCTCGCCAGCTGACCTGGGCAATGCCGAGGAAGCCATCGAACTTGAGCGTGAGGTCGGCCTGct                                                      >  2:711887/1‑150 (MQ=255)
                                                                                        ccAGACCACCGCCGAGTACCCACTGCACAAGGGCATGGTTTCTGACCGCGGCCTGAAACCGTTCGACGAGCTGCGCCCGCCGAAGATCTCGCCAGCTGACCTGGGCAATGCCGAGGAAGCCATCGAACTTGAGCGTGAGGTCGGCCTGct                                                      >  2:1050398/1‑150 (MQ=255)
                                                                                          aGACCACCGCCGAGTACCCACTGCACAAGGGCATGGTTCCTGACCGCGGCCTGAAACCGTTCGACGAGCTGCGCCCGCCGAAGATCTCGCCAGCTGACCTGGGCAATGCCGAGGAAGCCATCGAACTTGAGCGTGAGGTCGGCCTGctct                                                    >  2:174438/1‑150 (MQ=255)
                                                                                            accaccGCCGAGTACCCACTGCACAAGGGCATGGTTTCTGACCGCGGCCTGAAACCGTTCGACGAGCTGCGCCCGCCGAAGATCTCGCCAGCTGACCTGGGCAATGCCGAGGAAGCCATCGAACTTGAGCGTGAGGTCGGCCTGCTCtga                                                  <  2:285978/150‑1 (MQ=255)
                                                                                             ccaccGCCGAGTACCCACTGCACAAGGGCATGGTTTCTGACCGCGGCCTGAAACCATTCGACGAGCTGCGCCCGCCGAAGATCTCGCCAGCTGACCTGGGCAATGCCGAg                                                                                         <  2:258452/110‑1 (MQ=255)
                                                                                             ccaccGCCGAGTACCCACTGCACAAGGGCATGGTTTCTGACCGCGGCCTGAAACCATTCGACGAGCTGCGCCCGCCGAAGATCTCGCCAGCTGACCTGGGCAATGCCGAg                                                                                         >  1:258452/1‑110 (MQ=255)
                                                                                              caccGCCGAGTACCCACTGCACAAGGGCATGGTTTCTGACCGCGGCCTGAAACCGTTCGACGAGCTGCGCCCGCCGAAGATCTCGCCAGCTGACCTGGGCAATGCCGAGGAAGCCATCGAACTTGAGCGTGAGGTCGGCCTGCTCtgatg                                                >  1:92903/1‑150 (MQ=255)
                                                                                                ccgccgAGTACCCACTGCACAAGGGCATGGTTTCTGACCGCGGCCTGAAACCGTTCGACGAGCTGCGCCCGCCGAAGATCTCGCCAGCTGACCTGGGCAATGCCGAGGAAGCCATCGAACTTGAGCGTGAGGTCGGCCTGCTCTGATGAc                                              <  2:272691/150‑1 (MQ=255)
                                                                                                    cgAGTACCCACTGCACAAGGGCATGGTTTGTGACCGCGGCCTGACTCCGTTCGACGACGCGCTCGCGCCGAAGATCTCGCCAGCTGACCTGGGCAATGCCGAGGAAGCCATCGAACTTGAGCGTGAGGTCGGCCTGCTCTGATGACTgcc                                          <  1:526499/150‑1 (MQ=255)
                                                                                                    cgAGTACCCACTGCACAAGGGCATGGTTTCTGACCGCGGCCTGAAACCGTTCGACGAGCTGCGCCCGCCGAAGATCTCGCCAGCTGACCTGGGCAATGCCGAGGAAGCCATCGAACTTGAGCGTGAGGTCGGCCTGCTCTGATGACTgcc                                          <  1:1348384/150‑1 (MQ=255)
                                                                                                       gTACCCACTGCACAAGGGCATGGTTTCTGACCGCGGCCTGAAACCGTTCGACGAGCTGCGCCCGCCGAAGATCTCGCCAGCTGACCTGGGCAATGCCGAGGAAGCCATCGAACTTGAGCGTGAGGTCGGCCTGCTCTGATGACTgccgcc                                       >  1:551913/1‑150 (MQ=255)
                                                                                                           ccACTGCACAAGGGCATGGTTTCTGACCGCGGCCTGAAACCGTTCGACGAGCTGCGCCCGCCGAAGATCTCGCCAGCTGACCTGGGCAATGCCGAGGAAGCCATCGAACTTGAGCGTGAGGTCGGCCTGCTCTGATGACTGCCGCCCTGt                                   <  2:1256625/150‑1 (MQ=255)
                                                                                                            cACTGCACAAGGGCATGGTTTCTGACCGCGGCCTGAAACCGTTCGACGAGCTGCGCCCGCCGAAGAtctc                                                                                                                  <  1:657256/70‑1 (MQ=255)
                                                                                                            cACTGCACAAGGGCATGGTTTCTGACCGCGGCCTGAAACCGTTCGACGAGCTGCGCCCGCCGAAGAtctc                                                                                                                  >  2:657256/1‑70 (MQ=255)
                                                                                                                 cacaAGGGCATGGTTTCTGACCGCGGCCTGAAACGGTTCGACGAGCTGCGCTCGCCGAAGATCTTGCCAGCTGACGTGAGCGATGCCGAGGAAACGAACGAGCTTGAGCGTGAGGAAGACCTGCTCTGATGGCTGCTGCACTGTCTGAAc                             >  2:203064/1‑150 (MQ=255)
                                                                                                                        gCATGGTTTCTGACCGCGGCCTGAAACCGTTCGACGAGCTGCGCCCGCCGAAGATCTCGCCAGCTGACCTGGGCAATg                                                                                              >  2:758905/1‑78 (MQ=255)
                                                                                                                        gCATGGTTTCTGACCGCGGCCTGAAACCGTTCGACGAGCTGCGCCCGCCGAAGATCTCGCCAGCTGACCTGGGCAATg                                                                                              <  1:758905/78‑1 (MQ=255)
                                                                                                                         cATGGTTTCTGACCGCGGCCTGAAACCGTTCGACGAGCTGCGCCCGCCGAAGATCTCGCCAGCTGACCTGGGCAATGCCGAGGAAGCCATCGAACTTGAGCGTGAGGTCGGCCTGCTCTGATGACTGCCGCCCTGTCCGAACCGGTGCCg                     <  1:277619/150‑1 (MQ=255)
                                                                                                                         aacggTTTCAGGCCGCGGTCAGAAACCGTTCGACGAGCTGCGCCCGCCGAAGATCTCGCCAGCTGACCTGGGCAATGCCGAGGAAGCCATCGAACTTGAGCGTGAGGTCGGCCTGCTCTGATGACTGCCGCCCTGTCCGAACCGGt                         <  1:546065/143‑1 (MQ=255)
                                                                                                                         aacggTTTCAGGCCGCGGTCAGAAACCGTTCGACGAGCTGCGCCCGCCGAAGATCTCGCCAGCTGACCTGGGCAATGCCGAGGAAGCCATCGAACTTGAGCGTGAGGTCGGCCTGCTCTGATGACTGCCGCCCTGTCCGAACCGGt                         >  2:546065/4‑146 (MQ=255)
                                                                                                                          aTGGTTTCTGACCGCGGCCTGAAACCGTTCGACGAGCTGCGCCCGCCGAAGATCTCGCCAGCTGACCTGGGCAATGCCGAGGAAGCCATCGAACTTGAGCGTGAGGTCGGCCTGCTCTGATGACTGCCGCCCTGTCCGAACCGGTGCCgg                    >  1:601585/1‑150 (MQ=255)
                                                                                                                             gTTTCTGACCGCGGCCTGAAACCGTTCGACGAGCTGCGCCCGCCGAAGATCTCGCCAGCTGACCTGGGCAATGCCGAGGAAGCCATCGAACTTGAGCGTGAGGTCGGCCTGCTCTGATGACTGCCGCCCTGTCCGAACCGGTGCCGGTAc                 <  1:711887/150‑1 (MQ=255)
                                                                                                                              tttCTGACCGCGGCCTGAAACCGTTCGACGAGCTGCGCCCGCCGAAGATCTCGCCAGCTGACCTGGGCAATGCCGAg                                                                                         <  1:1245908/77‑1 (MQ=255)
                                                                                                                              tttCTGACCGCGGCCTGAAACCGTTCGACGAGCTGCGCCCGCCGAAGATCTCGCCAGCTGACCTGGGCAATGCCGAg                                                                                         >  2:1245908/1‑77 (MQ=255)
                                                                                                                                  tggccGCGGCCTGAAACCGTTCGACGAGCTGCGCCCGCCGAAGATCTCGCCAGCTGACCTGGGCAATGCCGAGGAAGCCATCGAACTTGAGCGTGAGGTCGGCCTGCTCTGATGACTGCCGCCCTGTCCGAACCGGTGCCGGTACGCTTc            >  1:1422821/4‑150 (MQ=255)
                                                                                                                                  tGACCGCGGCCTGAAACCGTTCGACGAGCTGCGCCCGCCGAAGATCTCGCCAGCTGACCTGGGCAATGCCGAGGAAGCCATCGAACTTGAGCGTGa                                                                  <  2:48183/96‑1 (MQ=255)
                                                                                                                                  tGACCGCGGCCTGAAACCGTTCGACGAGCTGCGCCCGCCGAAGATCTCGCCAGCTGACCTGGGCAATGCCGAGGAAGCCATCGAACTTGAGCGTGa                                                                  >  1:48183/1‑96 (MQ=255)
                                                                                                                                     ccGCGGCCTGAAACCGTTCGACGAGCTGCGCCCGCCGAAGATCTCGCCAGCTGACCTGGGCAATGCCGAGGAAGCCATCGAACTTGAGCGTGAGGTCGGCCTGCTCTGATGACTGCCGCCCTGTCCGAACCGGTGCCGGTACGCTTCGTa         >  1:72012/1‑150 (MQ=255)
                                                                                                                                      cgcgGCCTGAAACCGTTCGACGAGCTGCGCCCGCCGAAGATCTCGCCAGCTGACCTGGGCAATGCCGAGGAAGCCATCGAACTTGAGCGTGAGGTCGGCCTGCTCTGATGACTGCCGCCCTGTCCGAACCGGTGCCGGTACGCTTCGTAc        >  1:1297492/1‑150 (MQ=255)
                                                                                                                                         ggCCTGAAACCGTTCGACGAGCTGCGCCCGCCGAAGATCTCGCCAGCTGACCTGGGCAATGCCGAGGAAGCCATCGAACTTGAGCGTGAGGTCGGCCTGCTCTGATGACTGCCGCCCTGTCCGAACCGGTGCCGGTACGCTTCGTACCGa     <  2:353545/150‑1 (MQ=255)
                                                                                                                                         ggCCTGAAACCGTTCGACGAGCTGCGCCCGCCGAAGATCTCGCCAGCTGACCTGGGCAATGCCGAGGAAGCCATCGAACTTGAGCGTGAGGTCGGCCTGCTCTGATGACTGCCGCCCTGTCCGAACCGGTGCCGGTACGCTTCGTACCGa     <  1:1050398/150‑1 (MQ=255)
                                                                                                                                          gCCTGAAACCGTTCGACGAGCTGCGCCCGCCGAAGATCTCGCCAGCTGACCTGGGCAATGCCGAGGAAGCCATCGAACTTGAGCGTGAGGTCGGCCTGCTCTGATGACTGCCGCCCTGTCCGAACCGGTGCCGGTACGCTTCGTACCGAg    >  2:268020/1‑150 (MQ=255)
                                                                                                                                           ccTGAAACCGTTCGACGAGCTGCGCCCGCCGAAGATCTCGCCAGCTGACCTGGGCAATGCCGAGGAAGCCATCGAACTTGAGCGTGAGGTCGGCCTGCTCTGATGACTGCCGCCCTGGCCGAACCGGTGCCGGTACGCTTCGTACCGAGc   >  1:205123/1‑150 (MQ=255)
                                                                                                                                            cTGAAACCGTTCGACGAGCTGCGCCCGCCGAAGATCTCGCCAGCTGACCTGGGCAATGCCGAGGAAGCCATCGAACTTGAGCGTGAGGTCGGCCTGCTCTGATGACTGCCGCCCTGTCCGAACCGGTGCCGGTACGCTTCGTACCGAgag  >  1:1015239/1‑148 (MQ=255)
                                                                                                                                             |                                                                                                                                                    
CAGGCGTGGTCAAGGCCAGCAAACACCCGAAAGAAGCCCAGGCCCTGCTCAACTGGATGGCCAGCGAAGAAGGCCAGCGCGTGATCACCCAGACCACCGCCGAGTACCCACTGCACAAGGGCATGGTTTCTGACCGCGGCCTGAAACCGTTCGACGAGCTGCGCCCGCCGAAGATCTCGCCAGCTGACCTGGGCAATGCCGAGGAAGCCATCGAACTTGAGCGTGAGGTCGGCCTGCTCTGATGACTGCCGCCCTGTCCGAACCGGTGCCGGTACGCTTCGTACCGAGCC  >  NC_002947/5550579‑5550868

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑