Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,173,299 G→A 9.8% intergenic (‑292/+40) fadD‑I ← / ← fadD‑II long‑chain‑fatty acid‑‑CoA ligase/long‑chain‑fatty acid‑‑CoA ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,173,2990GA 9.8% 271.3 / 17.5 102intergenic (‑292/+40)fadD‑I/fadD‑IIlong‑chain‑fatty acid‑‑CoA ligase/long‑chain‑fatty acid‑‑CoA ligase
Reads supporting (aligned to +/- strand):  ref base G (48/44);  new base A (6/4);  total (54/48)
Fisher's exact test for biased strand distribution p-value = 7.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.24e-01

Only 100 of 103 total aligned reads displayed.

AAAACAGGCCTTATGGTCACAAAAATGACTCAATGAACAATACAGGTCATTTTTAGAGTAATTATCCAAAATGTCGCAGATAAGGCTGTATAGAGCCTGTTAGAGCGGTTTAAGGCATTCGCGGGTAAACCCGCTCCTACATTGGGCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCTATATCGCGCAACTCCCGCCGCAAAATCTTGCCCACCGGCGTCATGGGCAACGATTCGCGCAGCACGATGTGCTTGGGTACTTTGTAGCCGGTGAAGTTGG  >  NC_002947/5173151‑5173445
                                                                                                                                                    |                                                                                                                                                  
aaaaCAGGCCTTATGGTCACAAAAATGACTCAATGAACAATACAGGTCATTTTTAGAGTAATTATCCAAAATGTCGCAGATAAGGCTGTATAGAGCCTGTTAGAGCGGTTTAAGGCATTCGCGGGTAAACCCGCTCCTACATTGGGCAgt                                                                                                                                                   >  2:773766/1‑150 (MQ=255)
     aGGCCTTATGGTCACAAAAATGACTCAATGAACAATACAGGTCATTTTTAGAGTAATTATCCAAAATGTCGCAGATAAGGCTGTATAGAGCCTGTTAGAGCGGTTTAAGGCATTCGCGGGTAAACCCGCTCCTACATTGGGCAGTGTAgg                                                                                                                                              >  1:461242/1‑150 (MQ=255)
      ggCCTTATGGTCACAAAAATGACTCAATGAACAATACAGGTCATTTTTAGAGTAATTATCCAAAATGTCGCAGATAAGGCTGTATAGAGCCTGTTAGAGCGGTTTAAGGCATTCGCGGGTAAACCCGCTCCTACATTGGGCAGTGTAGGa                                                                                                                                             <  1:1337401/150‑1 (MQ=255)
       gCCTTATGGTCACAAAAATGACTCAATGAACAATACAGGTCATTTTTAGAGTAATTATCCAAAATGTCGCAGATAAGGCTGTATAGATCCTGTTAGAGCGGTTTAAGGCATTCGCGGGTAAACCCGCTCCTACATTGTGCAGTGTa                                                                                                                                                >  2:401606/1‑146 (MQ=255)
       gCCTTATGGTCACAAAAATGACTCAATGAACAATACAGGTCATTTTTAGAGTAATTATCCAAAATGTCGCAGATAAGGCTGTATAGAGCCTGTTAGAGCGGTTTAAGGCATTCGCGGGTAAACCCGCTCCTACATTGGGCAGTGTa                                                                                                                                                <  1:401606/146‑1 (MQ=255)
       gCCTTATGGTCACAAAAATGACTCAATGAACAATACAGGTCATTTTTAGAGTAATTATCCAAAATGTCGCAGATAAGGCTGTATAGAGCCTGTTAGAGCGGTTTAAGGCATTCGCGGGTAAACCCGCTCCTACACTGCCCAATGTAGGAg                                                                                                                                            >  2:1957354/1‑150 (MQ=39)
          ttATGGTCACAAAAATGACTCAATGAACAATACAGGTCATTTTTAGAGTAATTATCCAAAATGTCGCAGATAAGGCTGTATAGAGCCTGTTAGAGCGGTTTAAGGCATTCGCGGGTAAACCCGCTCCTACATTGGGCAGTGTAGGAGCgg                                                                                                                                         <  1:773766/150‑1 (MQ=255)
          ttATGGTCACAAAAATGACTCAATGAACAATACAGGTCATTTTTAGAGTAATTATCCAAAATGTCGCAGATAAGGCTGTATAGAGCCTGTTAGAGCGGTTTAAGGCATTCGCGGGTAAACCCGCTCCTACATTGGGCAGTGTAGGAGCgg                                                                                                                                         <  2:2148510/150‑1 (MQ=255)
          ttATGGTCACAAAAATGACTCAATGAACAATACAGGTCATTTTTAGAGTAATTATCCAAAATGTCGCAGATAAGGCTGTATAGAGCCTGTTAGAGCGGTTTAAGGCATTCGCGGGTAAACCCGCTCCTACATTGGGCAGTGTAGGAGCgg                                                                                                                                         <  2:348364/150‑1 (MQ=255)
          ttATGGTCACAAAAATGACTCAATGAACAATACAGGTCATTTTTAGAGTAATTATCCAAAATGTCGCAGATAAGGCTGTATAGAGCCTGTTAGAGCGGTTTAAGGCATTCGCGGGTAAACCCGCTCCTACACTGCCCAATGTAGGAGCgg                                                                                                                                         <  2:2081149/150‑1 (MQ=38)
           tATGGTCACAAAAATGACTCAATGAACAATACAGGTCATTTTTAGAGTAATTATCCAAAATGTCGCAGATAAGGCTGTATAGAGCCTGTTAGAGCGGTTTAAGGCATTCGCGGGTAAACCCGCTCCTACACTGCCCAATGTAGGAGCggg                                                                                                                                        <  1:1957354/150‑1 (MQ=38)
             tGGTCACAAAAATGACTCAATGAACAATACAGGTCATTTTTAGAGTAATTATCCAAAATGTCGCAGATAAGGCTGTATAGAGCCTGTTAGAGCGGTTTAAGGCATTCGCGGGTAAACCCGCTCCTACATTGGGCAGTGTAGGAGCGGGtt                                                                                                                                      >  1:1155294/1‑150 (MQ=255)
                tCACAAAAATGACTCAATGAACAATACAGGTCATTTTTAGAGTAATTATCCAAAATGTCGCAGATAAGGCTGTACAGAGCCTGTTAGAGCGGTTTAAGGCATTCGCGGGTAAACCCGCTCCTACATTGGGCAGTGTAGGAGCGGGTTTAc                                                                                                                                   >  2:1190068/1‑150 (MQ=255)
                  acaAAAATGACTCAATGAACAATACAGGTCATTTTTAGAGTAATTATCCAAAATGTCGCAGATAAGGCTGTATAGAGCCTGTTAGAGCGGTTTAAGGCATTCGCGGGTAAACCCGCTCCTACATTGGGCAGTGTAGGAGc                                                                                                                                           >  1:2314422/1‑140 (MQ=255)
                  acaAAAATGACTCAATGAACAATACAGGTCATTTTTAGAGTAATTATCCAAAATGTCGCAGATAAGGCTGTATAGAGCCTGTTAGAGCGGTTTAAGGCATTCGCGGGTAAACCCGCTCCTACATTGGGCAGTGTAGGAGc                                                                                                                                           <  2:2314422/140‑1 (MQ=255)
                   caAAAATGACTCAATGAACAATACAGGTCATTTTTAGAGTAATTATCCAAAATGTCGCAGATAAGGCTGTATAGAGCCTGTTAGAGCGGTTTAAGGCATTCGCGGGTAAACCCGCTCCTACATTGGGCAGTGTAGGAGCGGGTTTACCcg                                                                                                                                <  2:153930/150‑1 (MQ=255)
                           aCTCAATGAACAATACAGGTCATTTTTAGAGTAATTATCCAAAATGTCGCAGATAAGGCTGTATAGAGCCTGTTAGAGCGGTTTAAGGCATTCGCGGGTAAACCCGCTCCTACATTGGGCAGTGTAGGAGCGGGTTTACCCACGaagaag                                                                                                                        <  2:629876/150‑1 (MQ=255)
                               aaTGAACAATACAGGTCATTTTTAGAGTAATTATCCAAAATGTCGCAGATAAGGCTGTATAGAGCCTGTTAGAGCGGTTTAAGGCATTCGCGGGTAAACCCGCTCCTACATTGGGCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCAt                                                                                                                    >  2:1714486/1‑150 (MQ=255)
                                aTGAACAATACAGGTCATTTTTAGAGTAATTATCCAAAATGTCGCAGATAAGGCTGTATAGAGCCTGTTAGAGCGGTTTAAGGCATTCGCGGGTAAACCCGCTCCTACACTGCCCAATGTAGGAGCgg                                                                                                                                         <  1:1904334/128‑1 (MQ=38)
                                   aaCAATACAGGTCATTTTTAGAGTAATTATCCAAAATGTCGCAGATAAGGCTGTATAGAGCCTGTTAGAGCGGTTTAAGGCATTCGCGGGTAAACCCGCTCCTACATTGGGCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCg                                                                                                                >  1:269382/1‑150 (MQ=255)
                                        tACAGGTCATTTTTAGAGTAATTATCCAAAATGTCGCAGATAAGGCTGTATAGAGCCTGTTAGAGCGGTTTAAGGCATTCGCGGGTAAACCCGCTCCTACATTGGGCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGGCATCCCGGTATc                                                                                                           >  1:1156974/1‑150 (MQ=255)
                                            ggTCATTTTTAGAGTAATTATCCAAAATGTCGCAGATAAGGCTGTATAGAGCCTGTTAGAGCGGTTTAAGGCATTCGCGGGTAAACCCGCTCCTACATTGGGCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGc                                                                                                       <  2:269382/150‑1 (MQ=255)
                                              tCATTTTTAGAGTAATTATCCAAAATGTCGCAGATAAGGCTGTATAGAGCCTGTTAGAGCGGTTTAAGGCATTCGGGGGTAAACCCGCTCCTACATTGGGCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCta                                                                                                     <  1:1342022/150‑1 (MQ=255)
                                                       gagTAATTATCCAAAATGTCGCAGATAAGGCTGTATAGAGCCTGTTAGAGCGGTTTAAGGCATTCGCGGGTAAACCCGCTCCTACATTGGGCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCTATATCGCGCa                                                                                            >  2:2032889/1‑150 (MQ=255)
                                                               aTCCAAAATGTCGCAGATAAGGCTGTATAGAGCCTGTTAGAGCGGTTTAAGGCATTCGCGGGTAAACCCGCTCCTACATTGGGCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGGCATCCCGGTATCACGCTATATCGCGCAACTCccgc                                                                                    <  2:1156974/150‑1 (MQ=255)
                                                                   aaaaTGTCGCAGATAAGGCTGTATAGAGCCTGTTAGAGCGGTTTAAGGCATTCGCGGGTAAACCCGCTCCTACATTGGGCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCTATATCGCGCAACTCCCGCCGCa                                                                                <  1:1714486/150‑1 (MQ=255)
                                                                                     cTGTATAGAGCCTGTTAGAGCGGTTTAAGGCATTCGCGGGTAAACCCGCTCCTACATTGGGCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCTATATCgcgc                                                                                             >  1:915037/1‑119 (MQ=255)
                                                                                      tGTATAGAGCCTGTTAGAGCGGTTTAAGGCATTCGCGGGTAAACCCGCTCCTACATTGGGCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCTATATCGCGCAACTCCCGCCACAAAATCTTGCCCACCGGCGt                                                             >  1:1068933/1‑150 (MQ=255)
                                                                                      tGTACAGAGCCTGTTAGAGCGGTTTAAGGCATTCGCGGGTAAACCCGCTCCTACATTGGGCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCTATATCGCGCAACTCCCGCCGCAAAATCTTGCCCACCGGCGt                                                             <  1:1190068/150‑1 (MQ=255)
                                                                                        tataGAGCCTGTTAGAGCGGTTTAAGGCATTCGCGGGTAAACCCGCTCCTACATTGGGCAGTGTAGGAGCGGGTTTACCcg                                                                                                                                >  1:2299421/1‑81 (MQ=255)
                                                                                        tataGAGCCTGTTAGAGCGGTTTAAGGCATTCGCGGGTAAACCCGCTCCTACATTGGGCAGTGTAGGAGCGGGTTTACCcg                                                                                                                                <  2:2299421/81‑1 (MQ=255)
                                                                                         ataGAGCCTGTTAGAGCGGTTTAAGGCATTCGCGGGTAAACCCGCTCCTACATTGGGCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCTATATCGCGCAACTCCCGCCGCAAAATCTTGCCCACCGGCGTCAt                                                          <  1:940314/150‑1 (MQ=255)
                                                                                          taGAGCCTGTTAGAGCGGTTTAAGGCATTCTCGGGTAAACGCGCTCCTACCTTGGTCAGTGTAGTAGCGGGTTTACCCGCAAAAAAGCCATCACGCTATCACGCTATGTCGCGCGACTGCTGCCGCAAGAACTTGCCCAGTGGCGTCATg                                                         >  2:1205726/1‑150 (MQ=255)
                                                                                          taGAGCCTGTTAGAGCGGTTTAAGGCATTCGCGGGTAAACCCGCTCCTACATTGGGCAGTGTAGGAGCggg                                                                                                                                        >  1:2245824/1‑71 (MQ=255)
                                                                                          taGAGCCTGTTAGAGCGGTTTAAGGCATTCGCGGGTAAACCCGCTCCTACATTGGGCAGTGTAGGAGCggg                                                                                                                                        <  2:2245824/71‑1 (MQ=255)
                                                                                            gagCCTGTTAGAGCGGTTTAAGGCATTCGCGGGTAAACCCGCTCCTACATTGGGCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCTATATCGCGCAACTCCCGCCGCAAAATCTTGCCCACCGGCGTCATggg                                                       >  1:541502/1‑150 (MQ=255)
                                                                                            gagCCTGTTAGAGCGGTTTAAGGCATTCGCGGATAAACCCGCTCCTACATTGGGCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCTATATCGCGCAACTCCCGCCGCAAAATCTTGCCCACCGGCGTCATggg                                                       >  1:1489983/1‑150 (MQ=255)
                                                                                               ccTGTTAGAGCGGTTTAAGGCATTCGCGGGTAAACCCGCTCCTACATTGGGCAGTGTAGGAGCGGGTTTACCCGCgaagaa                                                                                                                         >  2:1254898/1‑81 (MQ=255)
                                                                                               ccTGTTAGAGCGGTTTAAGGCATTCGCGGGTAAACCCGCTCCTACATTGGGCAGTGTAGGAGCGGGTTTACCCGCgaagaa                                                                                                                         <  1:1254898/81‑1 (MQ=255)
                                                                                               ccTGTTAGAGCGGTTTAAGGCATTCGCGGGTAAACCCGCTCCTACATTGGGCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCTATATCGCGCAACTCCCGCCGCAAAATCTTGCCCACCGGCGTCATGGGCaa                                                    >  1:1405084/1‑150 (MQ=255)
                                                                                               ccTGTTAGAGCGGTTTAAGGCATTCGCGGGTAAACCCGCTCCTACATTGGGCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCTATATCGCGCAACTCCCGCCGCAAAATCTTGCCCACCGGCGTCATGGGCaa                                                    >  1:1395385/1‑150 (MQ=255)
                                                                                                cTGTTAGAGCGGTTTAAGGCATTCGCGGGTAAACCCGCTCCTACATTGGGCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCTATATCGCGCAACTCCCGCCGCAAAATCTTGCCCACCGGCGTCATGGGCAAc                                                   <  2:1155294/150‑1 (MQ=255)
                                                                                                  gTTAGAGCGGTTTAAGGCATTCGCGGGTAAACCCGCTCCTACATTGGGCAGTGTAGGAGCGGGt                                                                                                                                       >  1:526803/1‑64 (MQ=255)
                                                                                                  gTTAGAGCGGTTTAAGGCATTCGCGGGTAAACCCGCTCCTACATTGGGCAGTGTAGGAGCGGGt                                                                                                                                       <  2:526803/64‑1 (MQ=255)
                                                                                                  gTTAGAGCGGTTTAAGGCATTCGCGGGTAAACCCGCTCCTACATCGGGTCGTCTATGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCTATATCGCGCAACTCCCGCCGCAAAATCTTGCCCACCGGCGTCATGGGCAACGa                                                 <  1:2273039/150‑1 (MQ=39)
                                                                                                  gTTAGAGCGGTTTAAGGCATTCGCGGGTAAACCCACTCCTACATTGGGCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCTATATCGCGCAACTCCCGCCGCAAAATCTTGCCCACCGGCGTCATGGGCAACGa                                                 >  1:1033106/1‑150 (MQ=255)
                                                                                                   ttAGAGTGGTTTAAGGCATTCGCGGGTAACCCCGTTCCTACATTGGGTAGTGCGGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCTATATCGCGCAACTCCCGCCGCAAAATCTTGCCCACCGGCGTCATGGGCAACGAt                                                <  1:68893/150‑1 (MQ=255)
                                                                                                   ttAGAGCGGTTTAAGGCATTCGCGGGTAAACCCGCTCCTACATTGGGCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCTATATCGCGCAACTCCCGCCGCAAAATc                                                                           <  2:347557/123‑1 (MQ=255)
                                                                                                   ttAGAGCGGTTTAAGGCATTCGCGGGTAAACCCGCTCCTACATTGGGCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCTATATCGCGCAACTCCCGCCGCAAAATc                                                                           >  1:347557/1‑123 (MQ=255)
                                                                                                    tAGAGCGGTTTAAGGCATTCGCGGGTAAACCCGCTCCTACATTGGGCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCTATATCGCGCAACTCCCGCCGCAAAATCTTGCCCACCGGCGTCATGGGCAACGAtt                                               <  1:2032889/150‑1 (MQ=255)
                                                                                                       agCGGTTTAAGGCATTCGCGGGTAAACCCGCTCCTACATTGGGCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCTATATCGCGCAACTCCCGCCGCAAAATCTTGCCCACCGGCGTCATGGGCAACGATTcgc                                            >  2:393448/1‑150 (MQ=255)
                                                                                                         cGGTTTAAGGCATTCGCGGATAAACCCGCTCCTACATTGGGCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCTATATCGCGCAACTCCCGCCGCAAAATCTTGCCCACCGGCGTCATGGGCAACGATTCgcgc                                          <  2:1489983/150‑1 (MQ=255)
                                                                                                          ggattAAGGCATTCGCGGGTAAACCCGCTCCTACATTGGGCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCTATATCGCGCAACTCCCGCCGCAAAATCTTGCCCACCGGCGTCATGGGCAACGATTcgcgaa                                         >  1:2334124/4‑148 (MQ=255)
                                                                                                          ggTTTAAGGCATTCGCGGGTAAACCCGCTCCTACATTGGGCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCTATATCGCGCAACTCCCGCCGCAAAATCTTGCCCACCGGCGTCATGGGCAACGATTCGCgca                                         <  2:2124864/150‑1 (MQ=255)
                                                                                                           gTTTAAGGCATTCGCGGGTAAACCCGCTCCTACATTGGGCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCTATATCGCGCAACTCCCGCCGCAAAATCTTGCCCACCGGCGTCATGGGCAACGATTCGCgcag                                        >  2:1922812/1‑150 (MQ=255)
                                                                                                                  gCATTCGCGGGTAAACCCGCTCCTACATTGGGCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCTATATCGCGCAACTCCCGCCGCAAAATCTTGCCCACCGGc                                                               >  2:1342514/1‑120 (MQ=255)
                                                                                                                  gCATTCGCGGGTAAACCCGCTCCTACATTGGGCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCTATATCGCGCAACTCCCGCCGCAAAATCTTGCCCACCGGc                                                               <  1:1342514/120‑1 (MQ=255)
                                                                                                                    aTTCGCGGGTAAACCCGCTTCTACATTGGGCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCTATATCGCGCAACTCCCGCCGCAAAATCTTGCCCACCGGCGTCATGGGCAACGATTCGCGCAGCACGAtgtg                               >  1:671174/1‑150 (MQ=255)
                                                                                                                    aTTCGCGGGTAAACCCGCTCCTACATTGGGCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCTATATCGCGCAACTCCCGCCGCAAAATCTTGCCCACCGGCGTCATGGGCAACGATTCGCGCAGCACGAtgtg                               >  2:588954/1‑150 (MQ=255)
                                                                                                                       cgcgGGTAAACCCGCTCCTACATTGGGCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCTATATCGCGCAACTCCCGCCGCAAAATCTTGCCCACCGGCGTCATGGGCAACGATTCGCGCAGCACGATGTGCtt                            <  2:461242/150‑1 (MQ=255)
                                                                                                                       cgcgGAAAAACCCACTCCTACATTGGGCAGTGGAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCTATATCGCGCAACTCCCGCCGCAAAACCTTGCCCACCGGCGTCATGGGCAACGATTCGCGCAGCACGATGTGCtt                            <  1:1205726/150‑1 (MQ=255)
                                                                                                                        gcgGGTAAACCCGCTCCTACATTGGGCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCTATATCGCGCAACTCCCGCCGCAAAATCTTGCCCACCGGCGTCa                                                           <  1:58744/118‑1 (MQ=255)
                                                                                                                        gcgGGTAAACCCGCTCCTACATTGGGCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCTATATCGCGCAACTCCCGCCGCAAAATCTTGCCCACCGGCGTCa                                                           >  2:58744/1‑118 (MQ=255)
                                                                                                                          gggTAAACCCGCTCCTACATTGGGCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCgg                                                                                                               >  2:1740081/1‑64 (MQ=255)
                                                                                                                          gggTAAACCCGCTCCTACATTGGGCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCgg                                                                                                               <  1:1740081/64‑1 (MQ=255)
                                                                                                                          gggTAAACCCGCTCCTACATTGGGCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCTATATCGCGCAACTCCCGCCGCAAAATCTTGCCCACCGGCGTCATGGGCAACGATTCGCGCAGCACGATGTGCTTggg                         >  2:516812/1‑150 (MQ=255)
                                                                                                                               aaCCCGCTCCTACATTGGGCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCTATATCGCGCAACTCCCGCCGCAAAATCTTGCCCACCGGCGTCATGGGCAACGATTCGCGCAGCACGATGTGCTTGGGTACtt                    <  2:2158674/150‑1 (MQ=255)
                                                                                                                               aaCCCGCTCCTACATTGGGCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCTATATCGCGCAACTCCCGCCGCAAAATCTTGCCCACCGGCGTCATGGGCAACGATTCGCGCAGCACGATGTGCTTGGGTACtt                    <  2:1405084/150‑1 (MQ=255)
                                                                                                                                aCCCGCTCCTACACTGCCCAATGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCTATATCGCGCAACTCCCGCCGCAAAATCTTGc                                                                       >  1:1543389/1‑98 (MQ=25)
                                                                                                                                aCCCGCTCCTACACTGCCCAATGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCTATATCGCGCAACTCCCGCCGCAAAATCTTGc                                                                       <  2:1543389/98‑1 (MQ=25)
                                                                                                                                 cccGCTCCTACATTGGGCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCTATATCGCGCAACTCCCGCCGCAAAATCTTGCCCACCGGCGTCATGGGCAACGATTCGCGCAGCACGATGTGCTTGGGTACTTTg                  >  1:2008878/1‑150 (MQ=255)
                                                                                                                                 cccGCTCCTACACTGCCCAATGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCTATATCGCGCAACTCCCGCCGCAAAATCTTGCCCACCGGCGTCATGGGCAACGATTCGCGCAGCACGATGTGCTTGGGTACTTTg                  >  2:1836033/1‑150 (MQ=38)
                                                                                                                                 cccGCTCCTACACTGCCCAATGTAGGAGCGGGTTTACCCGAGGGGAAGCCATCCCGGTATCACGCTATATCGCGCAACTCCCGCCGCAAAATCTTGCCCACCGGCGTCATGGGCAACGATTCGCGCAGCACGATGTGCTTGGGTACTTTg                  >  2:81951/1‑150 (MQ=38)
                                                                                                                                  ccacTCCTACATTGGGCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCTATATCGCGCAACTCCCGCCGCAAAATc                                                                           <  1:355425/89‑1 (MQ=255)
                                                                                                                                  ccacTCCTACATTGGGCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCTATATCGCGCAACTCCCGCCGCAAAATc                                                                           >  2:355425/4‑92 (MQ=255)
                                                                                                                                  ccGCTCCTACATTGGGCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCTATATCGCGCAACTCCCGCCGCAAAATCTTGCCCACCGGCGTCATGGGCAACGATTCGCGCAGCACGATGTGCtt                            <  2:202630/139‑1 (MQ=255)
                                                                                                                                  ccGCTCCTACATTGGGCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCTATATCGCGCAACTCCCGCCGCAAAATCTTGCCCACCGGCGTCATGGGCAACGATTCGCGCAGCACGATGTGCtt                            >  1:202630/1‑139 (MQ=255)
                                                                                                                                   cGCTCCTACATTGGGCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCTATATCGCGCCACTCCCGCCGCAAAATCTTGCCCACCGGCGTCATGGGCAACGATTCGCGCAGCACGATGTGCTTGGGTACTTTGTa                >  2:1134313/1‑150 (MQ=255)
                                                                                                                                   cGCTCCTACACTGCCCAATGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCTATATCGCGCAACTCCCGCCGCAAAATCTTGCCCACCGGCGTCATGGGCAACGATTCGCGCAGCACGATGTGCTTGGGTACTTTGTa                >  1:617216/1‑150 (MQ=39)
                                                                                                                                     cTCCTACATTGGGCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCTATATCGCGCAACTCCCGCCGCAAAATCTTGCCCACCGGCGTCATGGGCAACGATTCGCGCAGCACGATGTGCTTGGGTACTTTGTAGc              <  2:2008878/150‑1 (MQ=255)
                                                                                                                                      tCCTACATTGGGCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCTATATCGCGCAACTCCCGCCGCAAAATCTTGCCCACCGGCGTCATg                                                         <  1:405593/106‑1 (MQ=255)
                                                                                                                                      tCCTACATTGGGCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCTATATCGCGCAACTCCCGCCGCAAAATCTTGCCCACCGGCGTCATg                                                         >  2:405593/1‑106 (MQ=255)
                                                                                                                                       ccTACATTGGGCAGTGTAGGAGCGGGTTTACCCGCGAGGAAGCCATCCCGGTATCACGCTATATCGCGCAACTCCCGCCGCAAAATCTTGCAGCTTTACAGGAGTGGCAACGATTCGCGCAGCACGATGTGCTTGGGTACTTTGTAGCCg            <  2:2322123/150‑1 (MQ=255)
                                                                                                                                       ccTACATTGGGCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCTATATCGCGCAACTCCCGCCGCAAAATCTTGCCCACCGGCGTCATGGGCAACGATTCGCGCAGCACGATGTGCTTGGGTACTTTGTAGCCg            >  2:861722/1‑150 (MQ=255)
                                                                                                                                       ccTACATTGGGCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCTATATCGCGCAACTCCCGCCGCAAAATCTTGCCCACCGGCGTCATGGGCAACGATTCGCGCAGCACGATGTGCTTGGGTACTTTGTAGCCg            >  2:2099784/1‑150 (MQ=255)
                                                                                                                                        cTACATTGGGCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCTATATCGCGCAACTCCCGCCGCAAAATCTTGc                                                                       >  1:368513/1‑90 (MQ=255)
                                                                                                                                        cTACATTGGGCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCTATATCGCGCAACTCCCGCCGCAAAATCTTGc                                                                       <  2:368513/90‑1 (MQ=255)
                                                                                                                                        cTACATTGGGCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCTATATCGCGCAACTCCCGCCGCAAAATCTTGCCCACCGGCGTCATGGGCAACGATTCGCGCAGCACGATGTGCTTGGGTACTTTGTAGCCgg           >  1:1829178/1‑150 (MQ=255)
                                                                                                                                         tACATTGGGCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCTATATCGCGCAACTCCCGCCGCAAAATCTTGCCCACCGGCGTCATGGGCAACGATTCGCGCAGCACGATGTGCTTGGGTACTTTGTAGCCGGt          >  1:1284875/1‑150 (MQ=255)
                                                                                                                                           cATTGGGCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCTATATCGCGCAACTCCCGCCGCAAAATCTTGCCCACCGGCGTCATGGGCAACGATTCGCGCAGCACGATGTGCtt                            <  1:2047743/130‑1 (MQ=255)
                                                                                                                                           cATTGGGCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCTATATCGCGCAACTCCCGCCGCAAAATCTTGCCCACCGGCGTCATGGGCAACGATTCGCGCAGCACGATGTGCtt                            >  2:2047743/1‑130 (MQ=255)
                                                                                                                                           cATTGGGCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCTATATCGCGCAACTCCCGCCGCAAAATCTTGCCCACCGGCGTCATGGGCAACGATTCGCGCAGCACGATGTGCTTg                           >  2:1327959/1‑131 (MQ=255)
                                                                                                                                           cATTGGGCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCTATATCGCGCAACTCCCGCCGCAAAATCTTGCCCACCGGCGTCATGGGCAACGATTCGCGCAGCACGATGTGCTTg                           <  1:1327959/131‑1 (MQ=255)
                                                                                                                                           cATTGGGCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCTATATCGCGCAACTCCCGCCGCAAAATCTTGCCCACCGGCGTCATGGGCAACGATTCGCGCAGCACGATGTGCTTGGGTACTTTGTAGCCGGTGa        >  2:195950/1‑150 (MQ=255)
                                                                                                                                            aTTGGGCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCTATATCGCGCAACTCCCGCCGCAAAATCTTGCCCACCGGCGTCATGGGCAACGATTCGCGCAGCACGATGTGCTTGGGTACTTTGTAGCCGGTGaa       >  2:1837550/1‑150 (MQ=255)
                                                                                                                                               ggccAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCTATATCGCGCAACTCCCGCCGCAAAATCTTGCCCACCGGCGTCATGGGCAACGATTCGCGCAGCACGATGTGCTTGGGTACTTTGTAGCCGGTGAAGtt    <  1:393448/147‑1 (MQ=255)
                                                                                                                                                ggCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCTATATCGCGCAACTCCCGCCGCAAAATCTTGCCCACCGGCGTCATGGGCAACGATTCGCGCAGCACGATGTGCTTGGGTACTTTGTAGCCGGTGAAGTTg   <  1:195950/150‑1 (MQ=255)
                                                                                                                                                 gCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCTATATCGCGCAACTCCCGCCGCAAAATCTTGCCCACCGGCGTCATGGGCAACGATTCGCGCAGCACGATGTGCTTGGGTACTTTGTAGCCGGTGAAGTTgg  <  2:1829178/150‑1 (MQ=255)
                                                                                                                                                 gCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCTATATCGCGCAACTCCCGCCGCAAAATCTTGCCCACCGGCGTCATGGGCAACGATTCGCGCAGCACGATGTGCTTGGGTACTTTGTAGCCGGTGAAGTTgg  >  2:1336342/1‑150 (MQ=255)
                                                                                                                                                 gCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCTATATCGCGCAACTCCCGCCGCAAAATCTTGCCCACCGGCGTCATGGGCAACGATTCGCGCAGCACGATGTGCTTGGGTACTTTGTAGCCGGTGAAGTTgg  >  1:1250453/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
AAAACAGGCCTTATGGTCACAAAAATGACTCAATGAACAATACAGGTCATTTTTAGAGTAATTATCCAAAATGTCGCAGATAAGGCTGTATAGAGCCTGTTAGAGCGGTTTAAGGCATTCGCGGGTAAACCCGCTCCTACATTGGGCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCTATATCGCGCAACTCCCGCCGCAAAATCTTGCCCACCGGCGTCATGGGCAACGATTCGCGCAGCACGATGTGCTTGGGTACTTTGTAGCCGGTGAAGTTGG  >  NC_002947/5173151‑5173445

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