Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,745,037 A→G 14.9% Q108R (CAA→CGA) ‡ paoC → aromatic aldehyde dehydrogenase molybdopterin‑binding subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,745,0370AG14.9% 158.9 / 18.1 67Q108R (CAA→CGA) ‡paoCaromatic aldehyde dehydrogenase molybdopterin‑binding subunit
Reads supporting (aligned to +/- strand):  ref base A (27/30);  new base G (4/6);  total (31/36)
Fisher's exact test for biased strand distribution p-value = 7.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.20e-01

CGTGCCCGGGGTGATCGCCGTGCTTGATCATGGCCATCGTCCACATATTTCCAGTTATGACGACGACTACAGCGACGCCGACGCCGCGCAGGGGGCGCCGTTTAGGCCCTTGTTCGATGATCGAGTGCTCTACAACGGTCAACCGTTGGCGCTGGTGGTTGCCCAAACCCTCGAACTGGCCCGGTATGCAGGTTCGCTGATCAGGATCGAGTATGCCGAGGAGGCCCACCAGACCGACCTTGGCGCATGCCAACAGTCGCATGCAGCACCCGCAGAGACACCGCCGCCG  >  NC_002947/3744897‑3745185
                                                                                                                                            |                                                                                                                                                     
cgTGCCCGGGGTGATCGCCGTGCTTGATCATGGCCATCGTCCACATATTTCCAGTTATGACGACGACTACAGCGACGCCGACGCCGCGCAGGGGGCGCCGTTTAGGCCCTTGTTCGATGATCGAGTGCTCTACAACGGTCAACCGTTgg                                                                                                                                               <  5:288201/149‑1 (MQ=255)
 gTGCCCGGGGTGATCGCCGTGCTTGATCATGGCCATCGTCCACATATTTCCAGTTATGACGACGACTACAGCGACGCCGACGCCGCGCAGGGGGCGCCGTTTAGGCCCTTGTTCGATGATCGAGTGCTCTACAACGGTTGACCGTTgta                                                                                                                                              <  3:360665/149‑3 (MQ=255)
    cccGGGGTGATCGCGGTGCTTGATCTTGGCCATCGTCCAGATATTTCCAGTTATGACGACGACTACAGCGACGCCGACGCCGCGCAGGGGGCGGCGTTTAGGCCCTTGTTCGATGATCGAGTGCTCTACAACGGTTGACCGTTgtaga                                                                                                                                            <  4:51020/148‑5 (MQ=255)
    cccGGGGTGATCGCCGTGCTTGATCATGGCCATCGTCCACATATTTCCAGTTATGACGACGACTACAGCGACGCCGACGCCGCGCAGGGGGCGCCGTTTAGGCCCTTGTTCGATGATCGAGTGCTCTACAACGGTTGACCGTTgtagag                                                                                                                                           <  1:84475/149‑6 (MQ=255)
    cccGGGGTGATCGCCGTGCTTGATCATGGCCATCGTCCACATATTTCCAGTTATGACGACGACTACAGCGACGCCGACGCCGCGCAGGGGGCGCCGTTTAGGCCCTTGTTCGATGATCGAGTGCTCTACAACGGTTGACCGTTgtagag                                                                                                                                           <  2:179605/149‑6 (MQ=255)
         ggTGATCGCCGTGCTTGATCATGGCCATCGTCCACATATTTCCAGTTATGACGACGACTACAGCGACGCCGACGCCGCGCAGGGGGCGCCGTTTAGGCCCTTGTTCGATGATCGAGTGCTCTACAACGGTCAACCGTTGGCGCtggtgg                                                                                                                                      <  1:92679/149‑1 (MQ=255)
           tGATCGCCGTGCTTGATCATGGCCATCGTCCACATATTTCCAGTTATGACGACGACTACAGCGACGCCGACGCCGCGCAGGGGGCGCCGTTTAGGCCCTTGTTCGATGATCGAGTGCTCTACAACGGTCAACCGTTGGCGCTGGTGGtt                                                                                                                                    <  4:114778/149‑1 (MQ=255)
                   gTGCTTGATCATGGCCATCGTCCACATATTTCCAGTTATGACGACGACTACAGCGACTCAGACGCCTCGCAGGGTGCGCCGTTTAGGCCCTTGTTCGATTACCGAGTGCTCTACAACGGTCAACCGTTGGCGCTGGTGGTTGCCCAAAc                                                                                                                           <  1:88177/149‑1 (MQ=255)
                       ttGATCATGGCCATCGTCCACATATTTCCAGTTATGACGACGACTACAGCGACGCCGACGCCGCGCAGGGGGCGCCGTTTAGGCCCTTGTTCGATGATCGAGTGCTCTACAACGGTCAACCGTTGGCGCtg                                                                                                                                          >  8:238010/1‑131 (MQ=255)
                       ttGATCATGGCCATCGTCCACATATTTCCAGTTATGACGACGACTACAGCGACGCCGACGCCGCGCAGGGGGCGCCGTTTAGGCCCTTGTTCGATGATCGAGTGCTCTACAACGGTCAACCGTTGGCGCtg                                                                                                                                          <  7:238010/131‑1 (MQ=255)
                              tGGCCATCGTCCACATATTTCCAGTTATGACGACGACTACAGCGACGCCGACGCCGCGCAGGGGGCGCCGTTTAGGCCCTTGTTCGATGATCGAGTGCTCTACAACGGTCAACCGTTGGCGCTGGTGGTTGCCCAAACCCTCGAACTgg                                                                                                                <  1:61493/149‑1 (MQ=255)
                              tGGCCATCGTCCACATATTTCCAGTTATGACGACGACTACAGCGACGCCGACGCCGCGCAGGGGGCGCCGTTTAGGCCCTTGTTCGATGATCGAGTGCTCTACAACGGTCAACCGTTGGCGCTGGTGGTTGCCCAAACCCTCGAACTgg                                                                                                                <  3:22031/149‑1 (MQ=255)
                              tGGCCATCGTCCACATATTTCCAGTTATGACGACGACTACAGCGACGCCGACGCCGCGCAGGGGGCGCCGTTTAGGCCCTTGTTCGATGATCGAGTGCTCTACAACGGTCAACCGTTGGCGCTGGTGGTCGCCCAAACCCTCGAACTgg                                                                                                                <  1:261487/149‑1 (MQ=255)
                               ggCCATCGTCCACATATTTCCAGTTATGACGACGACTACAGCGACGCCGACGCCGCGCAGGGGGCGCCGTTTAGGCCCTTGTTCGATGATCGAGTGCTCTACAACGGTCAACCGTTGGCGCTGGTGGTTGCCCAAACCCTCGAACTGGc                                                                                                               >  7:87707/1‑149 (MQ=255)
                                gCCATCGTCCACATATTTCCAGTTATGACGACGACTACAGCGACGCCGACGCCGCGCAGGGGGCGCCGTTTAGGCCCTTGTTCGATGATCGAGTGCTCTACAACGGTCAACCGTTGGCGCTGGTGGTTGCCCAAACCCTCGAACTGGcc                                                                                                              >  1:16003/1‑149 (MQ=255)
                                  cATCGTCCACATATTTCCAGTTATGACGACGACTACAGCGACGCCGACGCCGCGCAGGGGGCGCCGTTTAGGCCCTTGTTCGATGATCGAGTGCTCTACAACGGTCAACCGTTGGCGCTGGTGGTTGCCCAAAccc                                                                                                                         >  4:57875/1‑136 (MQ=255)
                                  cATCGTCCACATATTTCCAGTTATGACGACGACTACAGCGACGCCGACGCCGCGCAGGGGGCGCCGTTTAGGCCCTTGTTCGATGATCGAGTGCTCTACAACGGTCAACCGTTGGCGCTGGTGGTTGCCCAAAccc                                                                                                                         <  3:57875/136‑1 (MQ=255)
                                  cATCGTCCACATATTTCCAGTTATGACGACGACTACAGCGACGCCGACGCCGCGCAGGGGGCGCCGTTTAGGCCCTTGTTCGATGATCGAGTGCTCTACAACGGTCAACCGTTGGCGCTGGTGGTTGCCCAAACCCTCGAACTGGCCCg                                                                                                            >  2:241243/1‑149 (MQ=255)
                                         cacaTATTTCCAGTTATGACGACGACTACAGCGACGCCGACGCCGCGCAGGGGGCGCCGTTTAGGCCCTTGTTCGATGATCGAGTGCTCTACAACGGTCAACCGTTGGCGCTGGTGGTTGCCCAAACCCTCGAACTGGCCCGGTATGc                                                                                                      <  6:181397/148‑1 (MQ=255)
                                          acaTATTTCCAGTTATGACGACGACTACAGCGACGCCGACGCCGCGCAGGGGGCGCCGTTTAGGCCCTTGTTCGATGATCGAGTGCTCTACAACGGTCAACCGTTGGCGCTGGTGGTTGCCCAAACCCTCGAACTGGCCCGGTATGCAg                                                                                                    >  1:340370/1‑149 (MQ=255)
                                            atatTTCCAGTTATGACGACGACTACAGCGACGCCGACGCCGCGCAGGGGGCGCCGTTTAGGCCCTTGTTCGATGATCGAGTGCTCTACAACGGTCAACCGTTGGCGCTGGTGGTTGCCCAAAc                                                                                                                           >  6:263996/1‑124 (MQ=255)
                                            atatTTCCAGTTATGACGACGACTACAGCGACGCCGACGCCGCGCAGGGGGCGCCGTTTAGGCCCGTGTTCGATGATCGAGTGCTCTACAACGGTCAACCGTTGGCGCTGGTGGTTGCCCAAAc                                                                                                                           <  5:263996/124‑1 (MQ=255)
                                                ttCCAGTTATGACGACGACTACAGCGACGCCGACGCCGCGCAGGGGGCGCCGTTTAGGCCCTTGTTCGATGATCGAGTGCTCTACAACGGTCAACCGTTGGCGCTGGTGGTTGCCCAAACCCTCGAACTGGCCCGGTATGCAGGTTCGc                                                                                              >  2:62484/1‑149 (MQ=255)
                                                      ttATGACGACGACTACAGCGACGCCGACGCCGCGCAGGGGGCGCCGTTTAGGCCCTTGTTCGATGATCGAGTGCTCTACAACGGTCAACCGTTGGCGCTGGTGGTTGCAAAAAACCCTCGAAATGGCCCGGTATGCAGGTTCGCTGAt                                                                                          >  1:124937/1‑148 (MQ=255)
                                                      ttATGACGACGACTAAAGCGACGCCGACGCCGCGCAGGGGGCGCCGTTTAGGCACTTGTTCGATGATCGAGTGCTCTACAACGGTCAACAGTTGGCTCTGGTGGTTGCCCAAACACTCGCACTGGCCCGGTATGCAGGTTCGCTGAt                                                                                          >  8:197924/1‑147 (MQ=255)
                                                         tGACGACGACTACAGCGACGCCGACGCCGCGCAGGGGGCGCCGTTTAGGCCCTTGTTCGATGATCGAGTGCTCTACAACGGTCAACCGTTGGCGCTGGTGGTTGCCCAAACCCTCGAACTGGCCCGGTATGCAGGTTCGCTGATCAGGa                                                                                     <  3:339615/149‑1 (MQ=255)
                                                         tGACGACGACTACAGCGACGCCGACGCCGCGCAGGGGGCGCCGTTTAGGCCCTTGTTCGATGATCGAGTGCTCTACAACGGTCAACCGTTGGCGCTGGTGGTTGCCCAAACCCTCGAACTGGCCCGGTATGCAGGTTCGCTGATCAGGa                                                                                     <  6:130119/149‑1 (MQ=255)
                                                           acgacgacTACAGCGACGCCGACGCCGCGCAGGGGGCGCCGTTTAGGCCCTTGTTCGATGATCGAGTGCTCTACAACGGTCAACCGTTGGCGCTGGTGGTTGCCCAAACCCTCGAACTGGCCCGGTATGCAGGTTCGCTGATCAGGATc                                                                                   <  6:86505/149‑1 (MQ=255)
                                                           acgacgacTACAGCGACGCCGACGCCGCGCAGGGGGCGCCGTTTAGGCCCTTGTTCGATGATCGAGTGCTCTACAACGGTCAACCGTTGGCGCTGGTGGTTGCCCAAACCCTCGAACTGGCCCGGCATGCAGGTTCGCTGATCAGGATc                                                                                   >  2:223979/1‑149 (MQ=255)
                                                               cgacTACAGCGACGCCGACGCCGCGCAGGGGGCGCCGTTTAGGCCCTTGTTCGATGATCGAGTGCTCTACAACGGTCAACCGTTGGCGCTGGTGGTTGCCCAAACCCTCGAACTGGCCCGGTATGCAGGTTCGCTGATCAGGATCGAg                                                                                >  3:209499/1‑148 (MQ=255)
                                                                 acTACAGCGACGCCGACGCCGCGCAGGGGGCGCCGTTTAGGCCCTTGTTCGATGATCGAGTGCTCTACAACGGTCAACCGTTGGCGCTGGTGGTTGCCCAAACCCTCGAACTGGCCCGGTATGCAGGTTCGCTGATCAGGATCGAGTAt                                                                             <  2:128082/149‑1 (MQ=255)
                                                                    aCAGCGACGCCGACGCCGCGCAGGGGGCGCCGTTTAGGCCCTTGTTCGATGATCGAGTGCTCTACAACGGTCAACCGTTGGCGCTGGTGGTTGCCCAAACCCTCGAACTGGCCCGGTATGCAGGTTCGCTGATCAGGATCGAGTATGcc                                                                          <  8:158314/149‑1 (MQ=255)
                                                                     cAGCGACGCCGACGCCGCGCAGGGGGCGCCGTTTAGGCCCTTGTTCGATGCTCGAGCGCTCTACAACGGTCAACCGTTGGCGCTGGTGGGTGCCCAAACCCTCGAACTGGCCCGGTATGCAGGTTCGCTGATCAGGATCGAGTATGcc                                                                          <  6:169711/148‑1 (MQ=255)
                                                                     cAGCGACGCCGACGCCGCGCAGGGGGCGCCGTTTAGGCCCTTGTTCGATGATCGAGTGCTCTACAACGGTCAACCGTTGGCGCTGGTGGTTGCCCAACCCCTCGAACTGGCCCGGTATGCAGGTTCGCTGATCAGGATCGAGTATGcc                                                                          <  2:214514/148‑1 (MQ=255)
                                                                     cAGCGACGCCGACGCCGCGCAGGGGGCGCCGTTTAGGCCCTTGTTCGATGATCGAGTGCTCTACAACGGTCAACCGTTGGCGCTGGTGGTTGCCCAAACCCTCGAACTGGCCCGGTATGCAGGTTCGCTGATCAGGATCGa                                                                                 <  3:319235/141‑1 (MQ=255)
                                                                     cAGCGACGCCGACGCCGCGCAGGGGGCGCCGTTTAGGCCCTTGTTCGATGATCGAGTGCTCTACAACGGTCAACCGTTGGCGCTGGTGGTTGCCCAAACCCTCGAACTGGCCCGGTATGCAGGTTCGCTGATCAGGATCGa                                                                                 >  4:319235/1‑141 (MQ=255)
                                                                        cgacgccgacgcCGCGCAGGGGGCGCCGTTTAGGCCCTTGTTCGATGATCGAGTGCTCTACAACGGTCAACCGTTGGCGCTGGTGGTTGCCCAAACCCTCGAACTGGCCCGGTATGCAGGTTCGCTGATCAGGATCGAGTATGCCgag                                                                       >  2:279736/1‑148 (MQ=255)
                                                                        cgacgccgacgcCGCGCAGGGGGCGCCGTTTAGGCCCTTGTTCGATGATCGAGTGCTCTACAACGGTCAACCGTTGGCGCTGGTGGTTGCCCAAACCCTCGAACTGGCCCGGTATGCAGGTTCGCTGATCAGGAGCGAGTATGCCgag                                                                       >  6:365193/1‑148 (MQ=255)
                                                                           cgccgacgccgCGCAGGGGGCGCCGTTTAGGCCCTTGTTCGATGATCGAGTGCTCTACAACGGTCAACCGTTGGCGCTGGTGGTTGCCCAAACCCTCGAACTGGCCCGGTATGCAGGTTCGCTGATCAGGATCGAGTATGCCGaggagg                                                                   <  8:87707/149‑1 (MQ=255)
                                                                                  gccgcGCAGGGGGCGCCGTTTAGGCCCTTGTTCGATGATCGAGTGCTCTACAACGGTCAACCGTTCGCGCTGGTGGTTGCCCAACCCCTCGAACTGGCCCGGTATGCAGGTTCGCTGATCAGGATCGAGTATGCCGAGGAGGCccac                                                              >  5:32491/1‑147 (MQ=255)
                                                                                     gcgcAGGGGGCGCCGTTTAGGCCCTTGTTCGATGATCGAGTGCTCTACAACGGTCAACCGTTGGCGCTGGGGGTTGCCCAAACCCTCGAACTGGCCCGGTATGCAGGTTCGCTGATCAGGATCGAGTATGCCGAGGAGGCCCACCAgac                                                         <  2:274473/149‑1 (MQ=255)
                                                                                              gcgcCGTTTAGGCCCTTGTTCGATGATCGAGTGCTCTACAACGGTCAACCGTTGGCGCTGGTGGTTGCCCAAACCCTCGAACTGGCCCGGTATGCAGGTTCGCTGATCAGGATCGAGTATGCAGAGGAGGCCCACCAGACCGACCTTgg                                                >  6:244749/1‑149 (MQ=255)
                                                                                               cgcCGTTTAGGCCCTTGTTCGATGATCGAGTGCTCTACAACGGTCAACCGTTGGCGCTGGTGGTTGCCCAAACCCTCGAACTGGCCCGGTATGCAGGTTCGATGATCAGGATCGAGTATGCCGAGGAGGCCCACCAGACCGCCCTTGgc                                               >  7:245976/1‑149 (MQ=255)
                                                                                                gcCGTTTAGGCCCTTGTTCGATGATCGAGTGCTCTACAACGGTCAACCGTTGGCGCTGGTGGTTGCCCAAACCCTCGAACTGGCCCGGTATGCAGGTTCGCTGATCAGGATCGAGTATGCCGAGGAGGCCCACCAGACCGACCTTGgcg                                              <  4:209499/149‑1 (MQ=255)
                                                                                                       aGGCCCTTGTTCGATGATCGAGTGCTCTACAACGGTCAACCGTTGGCGCTGGTGGTTGCCCAAACCCTCGAACTGGCCCGGTATGCAGGTTCGCTGATCAGGATCGAGTATGCCGAGGAGGCCCACCAGACCGACCTTGGCGCATGcc                                        >  6:175726/1‑148 (MQ=255)
                                                                                                        ggCCCTTGTTCGATGATCGAGTGCTCTACAACGGTCAACCGTTGGCGCTGGTGGTTGCCCAAACCCTCGAACTGGCCCGGTATGCAGGTTCGCTGATCAGGATCGAGTATGCCGAGGAGGCCCACCAGACCGACCTTGGCGCATGCCaa                                      <  1:62484/149‑1 (MQ=255)
                                                                                                          cccTTGTTCGATGATCGAGTGCTCTACAACGGTCAACCGTTGGCGCTGGTGGTTGCCCAAACCCTCGAACTGGCCCGGTATGCAGGTTCGCTGATCAGGATCGAGTATGCCGAGGAGGCCCACCAGACCGACCTTGGCGCATGCCaa                                      >  3:29065/1‑147 (MQ=255)
                                                                                                                ttCGATGATCGAGTGCTCTACAACGGTCAACCGTTGGCGCTGGTGGTTGCCCAAACCCTCGAACTGGCCCGGTATGCAGGTTCGCTGATCAGGATCGAGTATGCCGAGGAGGCCCACCAGACCGACCTTGGCGCATGCCAACAGTCGCa                              >  4:208523/1‑149 (MQ=255)
                                                                                                                  cGATGATCGAGTGCTCTACAACGGTCAACCGTTGGCGCTGGTGGTTGCCCATACCCTCGAACTGTCCCGGTATGCAGTTTCGCTGATCAGGTTCGAGTATGCCGAGGAGGCCCACCAGACCGACCTTGGGGCATGCCAACAGTCGCAt                             <  7:197924/148‑1 (MQ=255)
                                                                                                                   gatgatCGAGTGCTCTACAACGGTCAACCGTTGGCGCTGGTGGTTGCCCAAACCCTCGAACTGGCCCGGTATGCAGGTTCGCTGATCAGGATCGAGTATGCCGAGGAGGCCCACCAGACCGACCTTGGCGCATGCCAACAGTCGCATg                            >  4:300377/1‑148 (MQ=255)
                                                                                                                    atgatCGAGTGCTCTACAACGGTCAACCGTTGGCGCTGGTGGTTGCCCAAACCCTCGAACTGGCCCGGTATGCAGGTTCGCTGATCAGGATCGAGTATGCCGAGGCGGCCCACCAGACCGACCTTGGCGCATGCCAACAGTCGCATgca                          <  5:365193/149‑1 (MQ=255)
                                                                                                                    atgatCGAGTGCTCTACAACGGTCAACCGTTGGCGCTGGTGGTTGCCCAAACCCTCGAACTGGCCCGGTATGCAGGTTCGCTGATCAGGATCGAGTATGCCGAGGAGGCCCACCAGACCGACCTTGGCGCATGCCAACAGTCGCATgca                          <  3:208523/149‑1 (MQ=255)
                                                                                                                      gatCGAGTGCTCTACAACGGTCAACCGTTGGCGCTGGTGGTTGCCCAAACCCTCGAACTGGCCCGGTATGCAGGTTCGCTGATCAGGATCGAGTATGCCGAGGAGGCCCACCAGACCGACCTTGGCGCATGCCAACAGTCGCATgcag                         >  3:263530/1‑148 (MQ=255)
                                                                                                                        tCGAGTGCTCTACAACGGTCAACCGTTGGCGCTGGTGGTTGCCCAAACCCTCGAACTGGCCCGGTATGCAGGTTCGCTGATCAGGATCGAGTATGCCGAGGAGGCCCACCAGACCGACCTTGGCGCATGCCAACAGTCGCATGCAGCAc                      <  1:241243/149‑1 (MQ=255)
                                                                                                                             tGCTCTACAACGGTCAACCGTTGGCGCTGGTGGTTGCCCAAACCCTCGAACTGGCCCGGTATGCAGGTTCGCTGATCAGGATCGAGTATGCCGAGGAGGCCCACCAGACCGACCTTGGCGCATGcc                                        >  7:203586/1‑126 (MQ=255)
                                                                                                                             tGCTCTACAACGGTCAACCGTTGGCGCTGGTGGTTGCCCAAACCCTCGAACTGGCCCGGTATGCAGGTTCGCTGATCAGGATCGAGTATGCCGAGGAGGCCCACCAGACCGACCTTGGCGCATGcc                                        <  8:203586/126‑1 (MQ=255)
                                                                                                                               agcgccAACGGTTGACCGTTGGCGCTGGTGGTTGCCCAAACCCTCGAACTGGCCCGGTATGCAGGTTCGCTGATCAGGATCGAGTATGCCGAGGAGGCCCACCAGACCGACCTTGGCGCAt                                           >  7:40288/6‑121 (MQ=255)
                                                                                                                               agcgccAACGGTTGACCGTTGGCGCTGGTGGTTGCCCAAACCCTCGAACTGGCCCGGTATGCAGGTTCGCTGATCAGGATCGAGTATGCCGAGGAGGCCCACCAGACCGACCTTGGCGCAt                                           <  8:40288/116‑1 (MQ=255)
                                                                                                                                tctACAACGGTCAACCGTTGGCGCTGGTGGTTGCCCAAACCCTCGAACTGGCCCGGTATGCAGGTTCGCTGATCAGGATCGAGTATGCCGAGGAGGCCCACCAGACCGACCTTGGCGCATGCCAACAGTCGCATGCAGCACCCGCa                 >  6:287472/1‑146 (MQ=255)
                                                                                                                                  tACAACGGTCAACCGTTGGCGCTGGTGGTTGCCCAAACCCTCGAACTGGCCCGGTATGCAGGTTCGCTGATCAGGATCGAGTATGCCGAGGAGGCCCACCAGACCGACCTTGGCGCATGCCAACAGTCGCATGCAGCACCCGCAgaga             >  8:219230/1‑148 (MQ=255)
                                                                                                                                  gccAACGGTTGACCGTTGGCGCTGGTGGTTGCCCAAACCCTCGAACTGGCCCGGTATGCAGGTTCGCTGATCAGGATCGAGTATGCCGAGGAGGCCCACCAGACCGACCTTGGCGCATGc                                         >  3:310698/3‑120 (MQ=255)
                                                                                                                                  gccAACGGTTGACCGTTGGCGCTGGTGGTTGCCCAAACCCTCGAACTGGCCCGGTATGCAGGTTCGCTGATCAGGATCGAGTATGCCGAGGAGGCCCACCAGACCGACCTTGGCGCATGc                                         <  4:310698/118‑1 (MQ=255)
                                                                                                                                  gccAACGGTTGACCGTTGGCGCTGGTGGTTGCCCAAACCCTCGAACTGGCCCGGTATGCAGGTTCGCTGATCAGGATCGAGTATGCCGAGGAGGCCCACCAGACCGACCTTGGCGCATGCCAACAGTCGCATGCAGCACCCGCAgaga             >  3:132514/3‑148 (MQ=255)
                                                                                                                                   ccAACGGTTGACCGTTGGCGCTGGTGGTTGCCCAACCCCTCGAACTGGCCCGGTATGCAGGTTCGCTGATCAGGATCGAGTATGCCGCGGAGGCCCACCAGACCGACCTTGGCGAATGCCAACAGTCGCATGCAGCACCCGCAGAGaca           >  7:72284/2‑149 (MQ=255)
                                                                                                                                    cAACGGTCAACCGTTGGCGCTGGTGGTTGCCCAAACCCTCGAACTGGCCCGGTATGCAGGTTCGCTGATCAGGATCGAGTATGCCGAGGAGGCCCACCAGACCGACCTTGGCGCATGCCAACAGTCGCATGCAGCACCCGCAGAGacac          <  4:263530/149‑1 (MQ=255)
                                                                                                                                        ggTCAACCGTTGGCGCTGGTGGTTGCCCAAACCCTCGAACTGGCCCGGTATGCAGGTTCGCTGATCAGGATCGAGTATGCCGAGGAGGCCCACCAGACCGACCTTGGCGCATGCCAACAGTCGCATGCAGCACCCGCAGAGACAccgcc      <  2:124937/149‑1 (MQ=255)
                                                                                                                                          tCAACCGTTGGCGATGGTGGTTCCCCAAACCCTCGAACTGGCCCGGTATGCAGGTTCGCTGATCAGGATCGAGTATGACGACGAGGCACACCAGACCCACCTTGGCGCATGCCCACAGTCGCATTCAGCACCCGCAGAGACAccgccgc    >  8:301884/1‑149 (MQ=255)
                                                                                                                                            aaCCGTTGGCGCTGGTGGTTGCCCAAACCCTCGAACTGGCCCGGTATGCAGGTTCGCTGATCAGGATCGAGTATGCCGAGGAGGCCCACCAGACCGACCTTGGCGCATGCCAACAGTCGCATGCAGCACCCGCAGAGACAccgccgccg  >  4:5144/1‑149 (MQ=255)
                                                                                                                                            |                                                                                                                                                     
CGTGCCCGGGGTGATCGCCGTGCTTGATCATGGCCATCGTCCACATATTTCCAGTTATGACGACGACTACAGCGACGCCGACGCCGCGCAGGGGGCGCCGTTTAGGCCCTTGTTCGATGATCGAGTGCTCTACAACGGTCAACCGTTGGCGCTGGTGGTTGCCCAAACCCTCGAACTGGCCCGGTATGCAGGTTCGCTGATCAGGATCGAGTATGCCGAGGAGGCCCACCAGACCGACCTTGGCGCATGCCAACAGTCGCATGCAGCACCCGCAGAGACACCGCCGCCG  >  NC_002947/3744897‑3745185

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: